add_acc_ptcls.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
add_acc_ptcls.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
add_acc_ptcls.f:       real*8      p_inj, p_cut, distr_array(maxshells, 250)
add_acc_ptcls.f:       real*8      distr_array_hi(maxshells, 250)
add_acc_ptcls.f:       real*8      esc_shell(in, 250), esc_shell_hi(in, 250)                    
add_acc_ptcls.f:       real*8      esc_dens_up(250), esc_dens_down(250),
add_acc_ptcls.f:c       real*8      esc_number_up(250), esc_number_up_hi(250)
add_acc_ptcls.f:       real*8      esc_dens_up_hi(250), esc_dens_down_hi(250)
add_acc_ptcls.f:       real*8      kappa(250, maxshells),kappa_hi(250, maxshells)
add_acc_ptcls.f:       real*8     distr_array_all(250)
add_acc_ptcls.f:       real *8 distr_all(250)
add_acc_ptcls.f:       del_p = (max_p - 2.0)/250.0
add_acc_ptcls.f:       Do i = 1, 250 
add_acc_ptcls.f:       Do i = 1, 250  
add_acc_ptcls.f:c          Do j = 1, 250
add_acc_ptcls.f:       DO j =1,250
addzx1.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
addzx1.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
addzx2.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
addzx2.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
addzx3.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
addzx3.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
alter.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
alter.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
bndyflgs.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
bndyflgs.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
branch.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
branch.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
bsetb1.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
bsetb1.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
bsetb2.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
bsetb2.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
bsetb3.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
bsetb3.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
bsetemf1.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
bsetemf1.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
bsetemf2.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
bsetemf2.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
bsetemf3.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
bsetemf3.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
bvald.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
bvald.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
bvale.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
bvale.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
bvalemf1.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
bvalemf1.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
bvalemf2.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
bvalemf2.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
bvalemf3.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
bvalemf3.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
bvalv1.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
bvalv1.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
bvalv2.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
bvalv2.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
bvalv3.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
bvalv3.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
coolinit.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
coolinit.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
ct.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
ct.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
currents.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
currents.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
data.f:     x     0.5549392504563904,0.563323366325109,0.5716157638237684, 
data.f:     x     0.9933782250569848,0.9936556501704964,0.9939225708887325, 
dataio.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
dataio.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
defaults.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
defaults.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
defaults.f:       do 250 n=1,niov
defaults.f:250    continue
diff_distr_fn_posn.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
diff_distr_fn_posn.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
diff_distr_fn_posn.f:       real*8      p_inj, p_cut, distr_array(maxshells, 250)   
diff_distr_fn_posn.f:       real*8      newdens(in), diffuse_array(maxshells, 250)
diff_distr_fn_posn.f:       real*8      esc_dens_up(250), esc_dens_down(250),
diff_distr_fn_posn.f:       real*8      esc_shell(in, 250)
diffuse_distr.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
diffuse_distr.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
diffuse_distr.f:       parameter (maxshells = 200, max_p = 10.0D0, p_number = 250, 
diffuse_distr.f:       integer escape_shell_num(1:250)
diffuse_distr.f:       integer escape_shell_num_hi(1:250)
diffuse_distr.f:       real*8 distr_all(250)
diffuse_distr.f:       real *8 temp_check(250), temp_check2(250)
diffuse_distr.f:       Do j = 1, 250
diffuse_distr.f:       Do j = 1, 250
diffuse_distr.f:       Do j =1, 250
diffuse_distr.f:c       Do j = 1, 250
diffuse_distr.f:       Do j = 1, 250
diffuse_distr.f:c     do j =1 ,250
diffuse_distr.f:          Do j = 1, 250
diffuse_distr.f:          Do j = 1, 250
diffuse_distr_integ.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
diffuse_distr_integ.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
diffuse_distr_integ.f:       real*8      p_inj, p_cut, distr_array(maxshells, 250)   
diffuse_distr_integ.f:       real*8      new_dens(in), diffuse_array(maxshells, 250)
diffuse_distr_integ.f:       real*8      esc_dens_up(250), esc_dens_down(250),
diffuse_distr_integ.f:       real*8      esc_shell(in, 250), total_dens 
distr_1AU_time.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
distr_1AU_time.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
distr_1AU_time.f:       real*8      kappa(250,maxshells),kappa_hi(250,maxshells)             
distr_1AU_time.f:       real*8      p_inj, p_cut, distr_array(maxshells, 250)   
distr_1AU_time.f:       real*8      esc_dens_up(250), esc_dens_down(250),
distr_1AU_time.f:       real*8      esc_shell(in, 250), arrive_distr(250,10) 
distr_1AU_time.f:       del_p = (max_p - 2.0)/250.0                
distr_1AU_time.f:       Do j = 1, 250
distr_1AU_time.f:        Do i = 1, 250
distr_fn_out.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
distr_fn_out.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
distr_fn_out.f:       real*8      p_inj, p_cut, distr_array(maxshells, 250)   
distr_fn_out.f:       real*8      new_ptcl_num(in), temp, diffuse_array(maxshells, 250)
distr_fn_out.f:       real*8      esc_dens_up(250), esc_dens_down(250),
distr_fn_out.f:       real*8      esc_shell(in, 250), dens_ptcls
distr_fn_out.f:       real*8      kappa(250,maxshells), kappa_hi(250, maxshells)
distr_fn_out.f:      del_p = (max_p - 2.0D0)/250.0D0     
distr_fn_out.f:	  Do j = 1, 250
diverg.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
diverg.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
doloop.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
doloop.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
emfs.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
emfs.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
emfs.f:           do 250 i=ism2,iep3
emfs.f:250        continue
expand_distr.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
expand_distr.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
expand_distr.f:       real*8      p_inj, p_cut, distr_array(maxshells, 250)  
expand_distr.f:       real*8      esc_dens_up(250), esc_dens_down(250), 
expand_distr.f:       real*8   temp_distr_array(maxshells, 250)
expand_distr.f:       real*8   esc_shell(in, 250)         
expand_distr.f:c       real*8  esc_number_up(250), esc_number_up_hi(250)
expand_distr.f:       real*8      temp_distr_array_hi(maxshells, 250)
expand_distr.f:       real*8   distr_array_hi(maxshells,250),esc_dens_up_hi(250),
expand_distr.f:     *      esc_dens_down_hi(250),esc_shell_hi(in,250)
expand_distr.f:       real*8 fo(250,maxshells)
expand_distr.f:       real*8 kappa(250,maxshells),kappa_hi(250,maxshells)
expand_distr.f:       real*8 Iplus(250,maxshells)
expand_distr.f:       real*8  temp_distr_array_all(maxshells, 250)
expand_distr.f:       real*8  distr_array_all(250)
expand_distr.f:c       Do i=1,250 
expand_distr.f:       del_p = (max_p - p0)/250.
expand_distr.f:          Do j = 1, 250
expand_distr.f:             Do j = 1, 250
expand_distr.f:             End Do             ! Do j = 1, 250
expand_distr.f:             Do j = 1, 250
expand_distr.f:             Do j = 1, 250
expand_distr.f:             Do j = 1, 250
expand_distr.f:             Do j = 1, 250
expand_distr.f:             Do j = 1, 250
expand_distr.f:         DO j =1,250
expand_distr.f:          Do j = ipmin(i), 250 ! proton enhenced turbulence.
expand_distr.f:c          Do j = ipmax(i)+1, 250  ! solar wind turbulence again.
expand_distr.f:          Do j = 1, 250
find_forw_shock_3.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
find_forw_shock_3.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
find_forw_shock_3.f:       real*8      esc_shell(in, 250)       
find_forw_shock.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
find_forw_shock.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
find_forw_shock.f:       real*8      esc_shell(in, 250)
floor.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
floor.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
gen_bndy.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
gen_bndy.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
glsave.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
glsave.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
glsave.f:       real*8      p_inj, p_cut, distr_array(maxshells, 250)
glsave.f:       real*8      distr_array_hi(maxshells, 250)
glsave.f:       real*8      esc_shell(in, 250), esc_shell_hi(in, 250)                    
glsave.f:       real*8      esc_dens_up(250), esc_dens_down(250),
glsave.f:c       real*8      esc_number_up(250), esc_number_up_hi(250)
glsave.f:       real*8      esc_dens_up_hi(250), esc_dens_down_hi(250)
glsave.f:       real*8      kappa(250, maxshells),kappa_hi(250, maxshells)
glsave.f:          Do j =1, 250
greeting.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
greeting.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
greeting.f:       if (iolog .gt. 0) write (iolog, 2250)
greeting.f:       if (iotty .gt. 0) write (iotty, 2250)
greeting.f:2250   format('GREETING:       * adiabatic equation of state')
gridx1.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
gridx1.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
gridx2.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
gridx2.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
gridx3.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
gridx3.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
gsetx1.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
gsetx1.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
gsetx2.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
gsetx2.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
gsetx3.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
gsetx3.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
intchk.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
intchk.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
intchk.f:           write (iotty, 2250) msg
intchk.f:2250   format('INTCHK  : ',a3,': Voxel dump switch on at')
lorentz.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
lorentz.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
lorentz.f:           do 250 i=is,iep1
lorentz.f:250        continue
machno.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
machno.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
maxwell_int.f:     x     0.5549392504563904,0.563323366325109,0.5716157638237684, 
maxwell_int.f:     x     0.9933782250569848,0.9936556501704964,0.9939225708887325, 
mget.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
mget.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
mget.f:         if (iolog .gt. 0) write (iolog, 2250) nt03  , iparm(23)
mget.f:         if (iolog .gt. 0) write (iolog, 2500)
mget.f:2250   format('MGET    : **** WARNING **** Data   partly read  : '
mget.f:2500   format('MGET    : **** WARNING **** JSYM   set in current  '
mocemfs.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
mocemfs.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
mocemfs.f:250      continue
momx1.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
momx1.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
momx1.f:           do 250 i=is,iep1
momx1.f:250        continue
momx3.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
momx3.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
momx3.f:           do 250 i=is,ie
momx3.f:250        continue
msave.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
msave.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
mstart.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
mstart.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
newgrid.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
newgrid.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
newvg.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
newvg.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
newx1.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
newx1.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
newx2.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
newx2.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
newx3.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
newx3.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
nmlsts.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
nmlsts.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
nudt.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
nudt.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
observed_ptcls.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
observed_ptcls.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
observed_ptcls.f:       real*8      p_inj, p_cut, distr_array(maxshells, 250)   
observed_ptcls.f:       real*8      esc_dens_up(250), esc_dens_down(250),
observed_ptcls.f:       real*8      esc_shell(in, 250), tot_esc_shell_num
observed_ptcls.f:       real*8      fpnorm,kappa(250,maxshells),kappa_hi(250,maxshells)
observed_ptcls.f:       del_p = (max_p - 2.0D0)/250.0D0   
parker_field.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
parker_field.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
pdvcool.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
pdvcool.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
pdv.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
pdv.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
posn_distr.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
posn_distr.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
posn_distr.f:       real*8      p_inj, p_cut, distr_array(maxshells, 250)   
posn_distr.f:       real*8      esc_dens_up(250), esc_dens_down(250),
posn_distr.f:       real*8      esc_shell(in, 250)  
posn_distr.f:       del_p = (max_p - 2.0)/250.0                
posn_distr.f:           Do j = 1, 250
ppafnc.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
ppafnc.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
pressure.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
pressure.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
resetall.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
resetall.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
resetv.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
resetv.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
setup.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
setup.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
shell_obs.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
shell_obs.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
shell_obs.f:       real*8      p_inj, p_cut, distr_array(maxshells, 250)   
shell_obs.f:       real*8      esc_dens_up(250), esc_dens_down(250),
shell_obs.f:       real*8      esc_shell(in, 250)  
shell_obs.f:       del_p = (max_p - 2.0)/250.0                
shell_obs.f:       Do j = 1, 250
shell_track.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
shell_track.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
shell_track.f:       real*8      esc_shell(in,250)
sndspd.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
sndspd.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
solar.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
solar.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
spenergy.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
spenergy.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
srcstep.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
srcstep.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
stov.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
stov.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
stv1.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
stv1.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
stv2.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
stv2.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
stv3.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
stv3.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
time_distr.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
time_distr.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
time_distr.f:       real*8      p_inj, p_cut, distr_array(maxshells, 250)   
time_distr.f:       real*8      esc_dens_up(250), esc_dens_down(250),
time_distr.f:       real*8      esc_shell(in, 250)  
time_distr.f:       del_p = (max_p - 2.0)/250.0                  
time_distr.f:           Do j = 1, 250
transprt.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
transprt.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
tranx1.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
tranx1.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
tranx3.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
tranx3.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
tslice.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
tslice.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
tslplot.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
tslplot.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
tslprt.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
tslprt.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
userdmp.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
userdmp.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
v1tos1.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
v1tos1.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
v2tos2.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
v2tos2.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
v3tos3.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
v3tos3.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
vecpot1.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
vecpot1.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
vecpot2.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
vecpot2.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
vecpot3.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
vecpot3.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
viscous.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
viscous.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
vshear.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
vshear.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
vtos.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
vtos.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
wave_intens.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
wave_intens.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
wave_intens.f:       real*8 p_inj, p_cut, distr_array(maxshells, 250)
wave_intens.f:       real*8 esc_dens_up(250), esc_dens_down(250)
wave_intens.f:       real*8 Iplus(250,maxshells), Ipluso(250), Ibohm(250)
wave_intens.f:       real *8 tempIplus(250) ! 2009/09/11
wave_intens.f:       real*8 fo(250,maxshells)
wave_intens.f:       real*8 kappa(250,maxshells),kappa_hi(250,maxshells)
wave_intens.f:       real*8 kappa_perp(250,maxshells), kappa_perp_hi(250,maxshells)
wave_intens.f:       real*8 s_ratio, s(250)
wave_intens.f:       real*8 x, dx, xintegral(250,1000), uintegral, u, Uo
wave_intens.f:       real*8   esc_shell(in, 250), R_shell_x(maxshells)
wave_intens.f:       del_p = (max_p - p0)/250.0D0
wave_intens.f:       Do i = 1, 250
wave_intens.f:       Do i = 1, 250
wave_intens.f:       DO i = 1, 250
wave_intens.f:       DO i = 1, 250
wave_intens.f:          DO i = 1, 250
wave_intens.f:             if (index .ge. 250) then ! the corresponding turbulence must be due to solar wind.
wave_intens.f:          Do i = 1, 250
wave_intens.f:          DO i = 1, 250
wave_intens.f:       DO i = 1, 250
wave_intens.f:          Do i = 1, 250
wave_intens.f:       Do j = 1, 250
wiggle.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
wiggle.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
x1fc3d.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
x1fc3d.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
x1zc1d.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
x1zc1d.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
x1zc3d.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
x1zc3d.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
x2fc3d.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
x2fc3d.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
x2zc1d.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
x2zc1d.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
x2zc3d.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
x2zc3d.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
x3fc3d.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
x3fc3d.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
x3zc1d.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
x3zc1d.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
x3zc3d.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
x3zc3d.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
zeus3d.f:     4  ,(insc(250),ntheta  ),(insc(251),nphi    ),(insc(252),icnvlv  )
zeus3d.f:     4  ,(rlsc(250),thetamin),(rlsc(251),thetamax),(rlsc(252),dtheta  )
zeus3d.f:       real*8   esc_shell(in,250)
zeus3d.f:       real*8   p_inj, p_cut, distr_array(maxshells,250)
zeus3d.f:       real*8   esc_dens_up(250),esc_dens_down(250)
zeus3d.f:       real*8   arrive_distr(250,10) 
zeus3d.f:       real*8   distr_array_hi(maxshells,250),esc_dens_up_hi(250),
zeus3d.f:     *      esc_dens_down_hi(250),esc_shell_hi(in,250)
zeus3d.f:       real*8 esc_ptcl_counter(250) ! time-integrated escaped particles.
zeus3d.f:       real*8 esc_ptcl_counter_hi(250) ! time-integrated escaped particles.
zeus3d.f:       real*8   distr_array_all(250)
zeus3d.f:       real*8   distr_all(250)
zeus3d.f:         Do j = 1, 250 
zeus3d.f:         Do j = 1, 250
zeus3d.f:          Do j =1, 250
zeus3d.f:          Do j =1, 250
zeus3d.f:       Do j =1, 250
