source: PSPA/parmelaPSPA/trunk/prch.f @ 445

Last change on this file since 445 was 12, checked in by lemeur, 12 years ago

parmela pspa initial

File size: 1.3 KB
Line 
1      subroutine prch(nc)
2c
3c     write of the field ez,er,bt for cell
4c     ip=3 ez(r,z),er(r,z),bt(r,z) rebuilt of the field
5c
6c---------------------------------------------------------------------
7c
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'
14c
15      common/pr/nz,nr,zbm(numz),rbm(numr)
16c--------------------------------------------------------------------------
17c*
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
45c+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
46 
Note: See TracBrowser for help on using the repository browser.