c
c
      PARAmeter (
     Q    NIP1=NI+1, NIP2=NI+2, NIP3=NI+3, NIP4=NI+4,  
     Q    NIM1=NI-1, NIM2=NI-2, NIM3=NI-3, NIM4=NI-4)  
      PARAmeter (                        
     Q     NJP1=NJ+1, NJP2=NJ+2, NJP3=NJ+3, NJP4=NJ+4,  
     Q     NJM1=NJ-1, NJM2=NJ-2, NJM3=NJ-3, NJM4=NJ-4,  
     Q     NJQ2=2*NJ, NJQ2P1=NJQ2+1, NJQ2P2=NJQ2+2, NJ2=NJ/2)     
      PARAmeter (
     Q     NKP1=NK+1, NKP2=NK+2, NKP3=NK+3, NKP4=NK+4,  
     Q     NKM1=NK-1, NKM2=NK-2, NKM3=NK-3, NKM4=NK-4,  
     Q     NK2=NK/2, NK2P1=NK2+1, nk2p2=nk2+2,
     Q     nk2p3=nk2+3,  nk4=nk/4, nk8=nk/8)            
      Parameter (nion_i = njp1,  nion_j=nk2p3)            
      PARAmeter (NOUT=ni,  NOUTP1=NOUT+1)                   
      parameter (nx=njp3, ny=nk2p3, kk=9, nxm2=nx-2, nym2=ny-2)
      real rearth,rion,rionm
      parameter( rearth=6.38e8,rion=6.5e8, rionm=6.5e6 )
      real erg2ev,effmass,boltzman
      parameter( erg2ev=1.6022e-12,effmass=2.24e-24,boltzman=4.8e-9)