Rev | Line | |
---|
[12] | 1 | subroutine trwimp(wtp,dwtp) |
---|
| 2 | c--apply impulse due to field from many waves |
---|
| 3 | c-- impulse lasts for a phase length of dwtp (degrees) |
---|
| 4 | c----------------------------------------------------------------------- |
---|
| 5 | save |
---|
| 6 | c |
---|
| 7 | include 'param_sz.h' |
---|
| 8 | include 'constcom.h' |
---|
| 9 | include 'coordcom.h' |
---|
| 10 | include 'misccom.h' |
---|
| 11 | include 'pcordcom.h' |
---|
| 12 | include 'syscom.h' |
---|
| 13 | include 'tstepcom.h' |
---|
| 14 | include 'wavescom.h' |
---|
| 15 | include 'ucom.h' |
---|
| 16 | c |
---|
| 17 | c-------------------------------------------------------------------------- |
---|
| 18 | c* |
---|
| 19 | ne=rne |
---|
| 20 | dwtmax=el(7,ne) |
---|
| 21 | if(dwtp.gt.dwtmax) dwtp=dwtmax |
---|
| 22 | zz=z |
---|
| 23 | rr=sqrt(x*x+y*y) |
---|
| 24 | call trwave(zz,rr,wez,wer,wbphi,iz,wtp) |
---|
| 25 | dummy=0. |
---|
| 26 | if(wez.eq.0..and.wer.eq.0.)return |
---|
| 27 | if((ip.eq.999).and.(nupar.eq.1))write(nimp,*) zz,dummy,dummy,wez |
---|
| 28 | if(rr.le.0.)rr=1.e-10 |
---|
| 29 | eomdt=1.e-2/erest*dcon*dwtp |
---|
| 30 | dbgz=wez*eomdt |
---|
| 31 | cliu ezdbgz2=dbgz/dwtp |
---|
| 32 | dgbr=eomdt*(wer-bgz/gamma*wbphi) |
---|
| 33 | bgz=bgz+dbgz |
---|
| 34 | bgx=bgx+dgbr*x/rr |
---|
| 35 | bgy=bgy+dgbr*y/rr |
---|
| 36 | gamma=sqrt(1.+bgx*bgx+bgy*bgy+bgz*bgz) |
---|
| 37 | return |
---|
| 38 | end |
---|
| 39 | c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++* |
---|
Note: See
TracBrowser
for help on using the repository browser.