Line | |
---|
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.