Line | |
---|
1 | subroutine swap(np,wtp) |
---|
2 | c---swap coordinates of particle np with particle ngood |
---|
3 | c---------------------------------------------------------------------- |
---|
4 | save |
---|
5 | c |
---|
6 | include 'param_sz.h' |
---|
7 | include 'constcom.h' |
---|
8 | include 'coordcom.h' |
---|
9 | include 'ncordscom.h' |
---|
10 | include 'pcordcom.h' |
---|
11 | include 'syscom.h' |
---|
12 | include 'ucom.h' |
---|
13 | c |
---|
14 | common/aswap/acor(7,imaa) |
---|
15 | common/back/bmax,byz0,byzs,pzowmin,pzowmax,iback |
---|
16 | common/jones/ajones,zjones,zcath(imaa) |
---|
17 | common/part/nparticle(imaa) |
---|
18 | common/pnum/ipnum(imaa) |
---|
19 | dimension xyz(7) |
---|
20 | c-------------------------------------------------------------------------- |
---|
21 | c* |
---|
22 | xyz(1)=x |
---|
23 | xyz(2)=bgx |
---|
24 | xyz(3)=y |
---|
25 | xyz(4)=bgy |
---|
26 | xyz(5)=z |
---|
27 | xyz(6)=bgz |
---|
28 | xyz(7)=rne |
---|
29 | c |
---|
30 | c for back bombardment |
---|
31 | if(iback.eq.1) then |
---|
32 | if(xyz(7).le.1.and.xyz(6).le.0.) then |
---|
33 | write(nback1) xyz(1),xyz(2),xyz(3),xyz(4),xyz(5),xyz(6), |
---|
34 | 1 erest*(gamma-1.),phizero(np),wtp |
---|
35 | endif |
---|
36 | endif |
---|
37 | if(np.eq.ngood)return |
---|
38 | c real t |
---|
39 | do 10 i=1,6 |
---|
40 | cord(i,np)=cord(i,ngood) |
---|
41 | cord(i,ngood)=xyz(i) |
---|
42 | 10 continue |
---|
43 | do 20 j=1,7 ! alpha magnet |
---|
44 | atmp=acor(j,np) |
---|
45 | acor(j,np)=acor(j,ngood) |
---|
46 | acor(j,ngood)=atmp |
---|
47 | 20 continue |
---|
48 | cord(7,np)=cord(7,ngood) |
---|
49 | cord(7,ngood)=xyz(7) |
---|
50 | wt=weight(np) |
---|
51 | weight(np)=weight(ngood) |
---|
52 | weight(ngood)=wt |
---|
53 | ph=phizero(np) |
---|
54 | phizero(np)=phizero(ngood) |
---|
55 | phizero(ngood)=ph |
---|
56 | iparticle=nparticle(np) |
---|
57 | nparticle(np)=nparticle(ngood) |
---|
58 | nparticle(ngood)=iparticle |
---|
59 | gam(np)=gam(ngood) |
---|
60 | gam(ngood)=gamma |
---|
61 | it=ipnum(np) |
---|
62 | ipnum(np)=ipnum(ngood) |
---|
63 | ipnum(ngood)=it |
---|
64 | t=zcath(np) |
---|
65 | zcath(np)=zcath(ngood) |
---|
66 | zcath(ngood)=t |
---|
67 | return |
---|
68 | end |
---|
69 | c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++* |
---|
Note: See
TracBrowser
for help on using the repository browser.