subroutine maggrd ! use param_module,only: hr,r,hi,ri,mlon,mlat,igeoplot,imagprt use coef_module,only: coefficient,coef_eqcgrd,coef_eqcgrd, | coef_magpot_ex,coef_magpot_tot use nc_module,only: output use map_module,only: prep_mag2geo ! implicit none ! set array for regular geographic grid and prepare for mapping fmag2geo if(igeoplot == 1) call prep_mag2geo ! ! calculate the coefficients for the toroidal and total equivalent current ! c_tot & c_tor call coefficient(hi) ! referenced to hi (110km) ! ptest ! ! calculate EQC at the ground call coef_eqcgrd ! ! calculate magnetic potential external due to ionosphere call coef_magpot_ex ! ! calculate magnetic potential induced & external call coef_magpot_tot ! ! calculate global magnetic perturbation call xyz_glb ! ! calculate magnetic perturbation for a station if(imagprt == 1 .or. imagprt == 2 ) call xyz_stat ! ! calculate b_0*dB call cal_bmagnitude ! ! put plot result in nc-output file if(imagprt == 0 .or. imagprt == 2 ) call output ! return ! end subroutine maggrd