Line | |
---|
1 | subroutine prch(nc) |
---|
2 | c |
---|
3 | c write of the field ez,er,bt for cell |
---|
4 | c ip=3 ez(r,z),er(r,z),bt(r,z) rebuilt of the field |
---|
5 | c |
---|
6 | c--------------------------------------------------------------------- |
---|
7 | c |
---|
8 | include 'param_sz.h' |
---|
9 | include 'cfldscom.h' |
---|
10 | include 'constcom.h' |
---|
11 | include 'misccom.h' |
---|
12 | include 'syscom.h' |
---|
13 | include 'ucom.h' |
---|
14 | c |
---|
15 | common/pr/nz,nr,zbm(numz),rbm(numr) |
---|
16 | c-------------------------------------------------------------------------- |
---|
17 | c* |
---|
18 | ne=nel |
---|
19 | if(el(8,ne).eq.0.)then |
---|
20 | zmin=-hcll(nc) |
---|
21 | zmax=hcll(nc) |
---|
22 | endif |
---|
23 | if((el(8,ne).eq.-1.) .or. (el(8,ne).eq.-2.))then |
---|
24 | zmin=0. |
---|
25 | zmax=hcll(nc) |
---|
26 | endif |
---|
27 | if(el(8,ne).eq.1.) then |
---|
28 | zmin=-hcll(nc) |
---|
29 | zmax=0. |
---|
30 | endif |
---|
31 | zpas=(zmax-zmin)/99. |
---|
32 | r=0. |
---|
33 | rpas=el(2,ne)/3. |
---|
34 | do 9401 j=1,4 ! loop on r |
---|
35 | z=zmin |
---|
36 | do 9400 i=1,100 ! loop on z |
---|
37 | call cfield(r,z,nc,ez,er,bt) |
---|
38 | write(nuch,*)nc,el(8,ne),r,z,ez,er,bt |
---|
39 | z=z+zpas |
---|
40 | 9400 continue |
---|
41 | r=r+rpas |
---|
42 | 9401 continue |
---|
43 | return |
---|
44 | end |
---|
45 | c+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
46 | |
---|
Note: See
TracBrowser
for help on using the repository browser.