!
      module eqvcur_module
!      
      use param_module,only: mlat,mlon,rmlat,rmlon
!      
      implicit none
!       
      real :: vallat(mlat),     ! TIEGCM mag. quasi-dipole latitude from south to north pole
     |        vallon(mlon),     ! TIEGCM mag. longitude
     |        reglat(rmlat),    ! reg. grid mag. quasi-dipole latitude from south to north pole
     |        reglon(rmlon),    ! reg grid mag. longitude
     |        regmlt(rmlon),    ! reg grid mag. local time
     |        reg_glat(rmlat),  ! reg. grid geographic latitude from south to north pole
     |        reg_glon(rmlon)   ! reg grid geographic longitude
!      
! kqlam = northward height-integrated current density, A/m
! kqphi = eastward height-integrated current density, A/m      
! jmr   = current along the field line, A/m^2
      real :: kqlam_org(mlon,mlat),     ! TIEGCM: K_qlam irreg. grid
     |        kqphi_org(mlon,mlat),     ! TIEGCM: K_qphi irreg. grid
     |        kqlam(rmlon,rmlat),       ! K_qlam reg. grid
     |        kqphi(rmlon,rmlat),       ! K_qphi
     |        dB_phi(rmlon,rmlat),   	! magnetic perturbation deltaB_phi
     |        dB_lam(rmlon,rmlat),   	! magnetic perturbation deltaB_lam
     |        dB_dwn(rmlon,rmlat),    	! magnetic perturbation deltaB_downward
     |        dB_phi_geo(rmlon,rmlat),  ! magnetic perturbation deltaB_phi (geo.direction)
     |        dB_lam_geo(rmlon,rmlat),  ! magnetic perturbation deltaB_lam (geo.direction)
     |        dB_dwn_geo(rmlon,rmlat),  ! magnetic perturbation deltaB_downward (geo.direction)
     |        bhat(rmlon,rmlat,3),	! unit vector along geomagnetic field direction
     |        babs_mag(rmlon,rmlat),  	! total dB (mag.coord.)
     |        babs_geo(rmlon,rmlat)  	! total dB (geo.coord.)
!----------------------------------------------------------------------------

      end module eqvcur_module