2: 2: ======================================================================== 2: Begin execution of tiegcm_trunk at 07/07/15 09:19:42 2: Host = iris 2: System = LINUX 2: Logname = foster 2: ======================================================================== 2: 2: init_timer: level= 1 rtc=F sys=T 2: Current working directory (cwd) = /hao/aim/foster/tiegcm_rholev1/tiegcm_hao 2: Process ID (pid) = 2204 2: 2: Reading input data... 2: 2: Completed successful read of namelist inputs. 2: 2: Input: mkntask chose ntask_lon= 2 ntask_lat= 2 (ntask= 4) 2: Will use the Heelis potential model 2: Note input: Setting BY to 0 with HEELIS potential model. 2: INPUT NOTE: adding mandatory field ZMAG to secondary history fields (field 24) 2: 2: ------------------------------------------------------------------------ 2: USER INPUT PARAMETERS: 2: label = tiegcm res=5.0 2: (optional text label for current run) 2: High-lat electric potential model: potential_model = HEELIS 2: gswm migrating diurnal file: gswm_mi_di_ncfile = /hao/aim/tgcm/data/gswm_diurn_5.0d_99km.nc 2: gswm migrating semi-diurnal file: gswm_mi_sdi_ncfile = /hao/aim/tgcm/data/gswm_semi_5.0d_99km.nc 2: start_year = 2002 (starting calendar day) 2: start_day = 80 (starting calendar year) 2: calendar_advance = 1 (model will be advanced in calendar time starting on this day) 2: step = 60 (model timestep (seconds)) 2: ntask_lon = 2 (number of mpi tasks in longitude dimension) 2: ntask_lat = 2 (number of mpi tasks in latitude dimension) 2: total tasks = ntask_lon*ntask_lat = 4 2: source = /hao/aim/tgcm/data/tiegcm1.95/TGCM.tiegcm1.95.pcntr_mareqx_smin.nc 2: (file or mss path containing source history) 2: source_start = 80, 0, 0 (model time of source history) 2: output (primary history output files) = 2: tiegcm.p_rholev1_nobugfix.nc, 2: start (model start times) = 2: 80, 0, 0 2: stop (model stop times) = 2: 81, 0, 0 2: hist (primary history disk write frequencies) = 2: 1, 0, 0 2: Maxmimum number of histories per primary file = 10 2: secout (secondary history output files)= 2: tiegcm.s_rholev1_nobugfix.nc, 2: secstart (secondary history start times) = 2: 80, 1, 0 2: secstop (secondary history stop times) = 2: 81, 0, 0 2: sechist (secondary history disk write frequencies) = 2: 0, 1, 0 2: secflds (secondary history fields) = 2: TN UN VN O2 O1 2: N2 NO N4S HE NE 2: TE TI TEC O2P OMEGA 2: POTEN UI_ExB VI_ExB WI_ExB DEN 2: QJOULE Z ZG ZMAG 2: Maximum number of histories per secondary file = 24 2: Number of bytes for values of fields on secondary histories (sech_nbyte) = 4 2: dynamo = 1 (dynamo will be calculated) 2: current_pg = 1 (Add current due to plasma pressure gradient and gravity to rhs) 2: current_kq = 0 (Calculate height-integrated current density) 2: eddy_dif = 0 (DOY-dependent eddy diffusion flag) 2: tide (amplitudes and phases of semidiurnal tide) = 2: 0.0E+00 0.0E+00 0.0E+00 0.0E+00 0.0E+00 0.00 0.00 0.00 0.00 0.00 2: tide2 (amplitude and phase of diurnal tide) = 2: 0.0E+00 0.00 2: aurora = 1 (0/1 flag for aurora) 2: colfac = 1.500 (collision factor) 2: joulefac = 1.500 (joule heating factor) 2: calc_helium = 1 (0/1 flag for helium) 2: If any of the following are spval ( 0.1000E+37), they will be calculated 2: during the simulation on a per timestep basis: 2: power = 0.1800E+02 (Hemispheric Power) 2: ctpoten= 0.3000E+02 (Cross-cap potential) 2: kp = 0.1000E+37 (Kp index) 2: bximf = 0.1000E+37 (BX component of IMF) 2: byimf = 0.0000E+00 (BY component of IMF) 2: bzimf = 0.1000E+37 (Bz component of IMF) 2: swvel = 0.1000E+37 (solar wind velocity) 2: swden = 0.1000E+37 (solar wind density) 2: f107 = 0.7000E+02 (F10.7 solar flux) 2: f107a = 0.7000E+02 (81-day ave F10.7 flux) 2: al = 0.1000E+37 (AL, lower auroral mag index) 2: END USER INPUT PARAMETERS 2: ------------------------------------------------------------------------ 2: 2: tgcm: dynamo= 1 -- calling apxparm: start_year= 2002 2: Time in apxparm = 0.897 (secs) 2: 2: ntask= 4 ntaski= 2 ntaskj= 2 Task Table: 2: j= -1 itask_table(:,j)= -1 -1 -1 -1 2: j= 0 itask_table(:,j)= -1 0 1 -1 2: j= 1 itask_table(:,j)= -1 2 3 -1 2: j= 2 itask_table(:,j)= -1 -1 -1 -1 2: 2: mytid= 2 mytidi,j= 0 1 lat0,1= 19 36 (18) lon0,1= 1 38 (38) ncells= 684 2: tgcm after distribute_geo 2: 2: ntask= 4 ntask_maglon= 2 ntask_maglat= 2 Task Table: 2: j= -1 itask_table_mag(:,j)= -1 -1 -1 -1 -1 2: j= 0 itask_table_mag(:,j)= 1 0 1 0 -1 2: j= 1 itask_table_mag(:,j)= 3 2 3 2 -1 2: j= 2 itask_table_mag(:,j)= -1 -1 -1 -1 -1 2: 2: mytid= 2 magtidi,j= 0 1 mlat0,1= 50 97 (48) mlon0,1= 1 41 (41) ncells=1968 2: mp_distribute_geo: nmagtaski= 2 tidcol= 1 2: tgcm after distribute_mag 2: 2: Task 2: 2: 2: Subdomain on geographic grid: 2: tasks( 2)%mytid = 2 2: tasks( 2)%mytidi= 0 2: tasks( 2)%mytidj= 1 2: tasks( 2)%nlats = 18 2: tasks( 2)%nlons = 38 2: tasks( 2)%lat0 = 19 2: tasks( 2)%lat1 = 36 2: tasks( 2)%lon0 = 1 2: tasks( 2)%lon1 = 38 2: Number of geo subdomain grid points = 684 2: 2: Subdomain on geomagnetic grid: 2: tasks( 2)%magtidi= 0 2: tasks( 2)%magtidj= 1 2: tasks( 2)%nmaglats = 48 2: tasks( 2)%nmaglons = 41 2: tasks( 2)%mlat0 = 50 2: tasks( 2)%mlat1 = 97 2: tasks( 2)%mlon0 = 1 2: tasks( 2)%mlon1 = 41 2: Number of mag subdomain grid points = 1968 2: mp_exchange_tasks: mxmaglat= 49 mxmaglon= 41 2: tgcm after esmf_init 2: Model version = tiegcm_trunk 2: 2: Set constants: 2: nlat= 36 nlon= 72 nlev= 28 2: dz= 0.50 2: dlat= 5.00 dlon= 5.00 2: zbound (cm) = 0.963723E+07 2: zmbot, zmtop = -6.750 7.250 (bottom,top midpoint levels) 2: zibot, zitop = -7.000 7.000 (bottom,top interface levels) 2: dt = 60.00 secs 2: grav = 870.00 2: freq_3m3 = 0.3506E-04 freq_semidi= 0.1454E-03 2: dipmin = 0.170 2: check_exp = F 2: kut (for filtering) = 1 2 3 5 6 7 9 10 11 13 14 15 17 17 17 17 17 17 17 17 17 17 17 17 15 14 13 11 10 9 7 6 5 3 2 1 2: init: iyear= 2002 iday= 80 2: gswm_mi_di_ncfile = /hao/aim/tgcm/data/gswm_diurn_5.0d_99km.nc 2: gswm_mi_sdi_ncfile = /hao/aim/tgcm/data/gswm_semi_5.0d_99km.nc 2: hist_init: nstep= 1440 2: alloc_ar: allocated module data 2: Allocated gswm t,u,v,z (lon0:lon1,lat0:lat1) 2: Allocated private gwm t,u,v,z (lon0:lon1,lat0:lat1,nmonth,nhour) 2: alloc_pdyn: allocate fmeq_in with mlev0,1= -2 29 2: allocdata: all tasks allocate fzg: nlevp1= 29 nlonp4= 76 nlat= 36 2: init_lbc: allocated subdomains tlbc, ulbc, vlbc 2: init_lbc_nm: allocated subdomains tlbc_nm, ulbc_nm, vlbc_nm 2: init_lbc: allocated globals tlbc_glb, ulbc_glb, vlbc_glb: nlonp4= 76 nlat= 36 2: init_lbc: allocated globals tlbc_nm_glb, ulbc_nm_glb, vlbc_nm_glb: nlonp4= 76 nlat= 36 2: 2: Model run initialization: 2: nstep = 1440 (Number of time steps this run) 2: iter = 115200 (Initial iteration number) 2: iyear = 2002 (Beginning calendar year) 2: iday = 80 (Beginning calendar day) 2: igswm_mi_di = 1 (If > 0, GSWM diurnal tidal database will be used.) 2: igswm_mi_sdi= 1 (If > 0, GSWM semidiurnal tidal database will be used.) 2: igswm_nm_di= 0 (If > 0, GSWM nonmigrating diurnal tidal database will be used.) 2: igswm_nm_sdi= 0 (If > 0, GSWM nonmigrating semidiurnal tidal database will be used.) 2: 2: This is an initial run: 2: start_year = 2002 (Starting year of initial run) 2: start_day = 80 (Starting day of initial run) 2: start_mtime= 80 0 0 (Starting mtime of initial run) 2: 2: Primary Histories: 2: nsource = 1 (If > 0, a primary source history was provided) 2: nseries_prim = 1 (Number of primary time series) 2: nhist_total = 2 (Number of primary histories to be written) 2: nfiles_prim = 1 (Number of primary output files to be written) 2: mxhist_prim = 10 (Maximum number of primary histories per file) 2: 2: Secondary Histories: 2: nseries_sech = 1 (Number of secondary time series) 2: nsech_total = 24 (Number of secondary histories to be written) 2: nfiles_sech = 1 (Number of secondary output files to be written) 2: mxhist_sech = 24 (Maximum number of secondary histories per file) 2: nfsech = 24 (Number of requested secondary history fields) 2: secondary history field 1: TN 2: secondary history field 2: UN 2: secondary history field 3: VN 2: secondary history field 4: O2 2: secondary history field 5: O1 2: secondary history field 6: N2 2: secondary history field 7: NO 2: secondary history field 8: N4S 2: secondary history field 9: HE 2: secondary history field 10: NE 2: secondary history field 11: TE 2: secondary history field 12: TI 2: secondary history field 13: TEC 2: secondary history field 14: O2P 2: secondary history field 15: OMEGA 2: secondary history field 16: POTEN 2: secondary history field 17: UI_ExB 2: secondary history field 18: VI_ExB 2: secondary history field 19: WI_ExB 2: secondary history field 20: DEN 2: secondary history field 21: QJOULE 2: secondary history field 22: Z 2: secondary history field 23: ZG 2: secondary history field 24: ZMAG 2: 2: ------------------------------------------------------------------------------------------------ 2: Table of Available Diagnostic Fields: 2: Shortnames may be added to namelist SECFLDS 2: 2: Field Shortname Units Levels Caller Longname 2: - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 2: 1 CO2_COOL erg/g/s lev newton.F CO2 Cooling 2: 2 NO_COOL erg/g/s lev newton.F NO Cooling 2: 3 DEN g/cm3 ilev dt.F Total Density 2: 4 HEATING erg/g/s lev dt.F Total Heating 2: 5 HMF2 km none elden.F HMF2: Height of the F2 Layer 2: 6 NMF2 1/cm3 none elden.F NMF2: Peak Density of the F2 Layer 2: 7 FOF2 MHz none elden.F FOF2: Critical Frequency of the F2 Layer 2: 8 JE13D A/m2 mlev current.F JE13D: Eastward current density (3d) 2: 9 JE23D A/m2 mlev current.F JE23D: Downward current density (3d) 2: 10 JQR A/m2 none current.F JQR: Upward current density (2d) 2: 11 KQLAM A/m none current.F KQLAM: Height-integrated current density (+north) 2: 12 KQPHI A/m none current.F KQPHI: Height-integrated current density (+east) 2: 13 LAMDA_HAL 1/s lev lamdas.F LAMDA_HAL: Hall Ion Drag Coefficient 2: 14 LAMDA_PED 1/s lev lamdas.F LAMDA_PED: Pedersen Ion Drag Coefficient 2: 15 MU_M g/cm/s lev cpktkm.F MU_M: Molecular Viscosity Coefficient 2: 16 QJOULE erg/g/s lev qjoule.F QJOULE: Joule Heating 2: 17 SCHT km lev addiag.F SCHT: Pressure Scale Height 2: 18 SIGMA_HAL S/m lev lamdas.F SIGMA_HAL: Hall Conductivity 2: 19 SIGMA_PED S/m lev lamdas.F SIGMA_PED: Pedersen Conductivity 2: 20 TEC 1/cm2 none elden.F TEC: Total Electron Content 2: 21 UI_ExB cm/s ilev ionvel.F UI: Zonal Ion Drift (ExB) 2: 22 VI_ExB cm/s ilev ionvel.F VI: Meridional Ion Drift (ExB) 2: 23 WI_ExB cm/s ilev ionvel.F WI: Vertical Ion Drift (ExB) 2: 24 WN cm/s ilev swdot.F WN: Neutral Vertical Wind (plus up) 2: 25 O_N2 none lev comp.F O/N2 RATIO 2: 26 QJOULE_INTEG erg/cm2/s none qjoule.F Height-integrated Joule Heating 2: 27 BX none oplus.F BX/BMAG: Normalized eastward component of magnetic field 2: 28 BY none oplus.F BY/BMAG: Normalized northward component of magnetic field 2: 29 BZ none oplus.F BZ/BMAG: Normalized upward component of magnetic field 2: 30 BMAG Gauss none oplus.F BMAG: Magnetic field magnitude 2: 31 EX V/m ilev ionvel.F EX: Zonal component of electric field 2: 32 EY V/m ilev ionvel.F EY: Meridional component of electric field 2: 33 EZ V/m ilev ionvel.F EZ: Vertical component of electric field 2: 34 ED1 V/m imlev dynamo.F ED1: Magnetic eastward component of electric field 2: 35 ED2 V/m imlev dynamo.F ED2: Magnetic downward (equatorward) component of electric field 2: 36 PHIM2D V none dynamo.F PHIM2D: 2d Electric Potential on magnetic grid 2: 37 N2 mmr lev comp.F N2: Molecular Nitrogen 2: 38 ZGMID cm lev addiag.F ZGMID: Geometric Height at midpoints 2: 39 CUSP erg/cm2/s none dynamics.F CUSP (cusp2d*ec) 2: 40 DRIZZLE erg/cm2/s none dynamics.F DRIZZLE (drzl2d*ed) 2: 41 ALFA keV none dynamics.F ALFA 2: 42 NFLUX #/cm2/s none dynamics.F NFLUX 2: 43 EFLUX erg/cm2/s none dynamics.F EFLUX 2: ------------------------------------------------------------------------------------------------ 2: 2: 2: ------------------------------------------------------------------------ 2: Getfile: remote=/hao/aim/tgcm/data/tiegcm1.95/TGCM.tiegcm1.95.pcntr_mareqx_smin.nc 2: Getfile: Found file /hao/aim/tgcm/data/tiegcm1.95/TGCM.tiegcm1.95.pcntr_mareqx_smin.nc 2: ------------------------------------------------------------------------ 2: 2: Acquired source history file /hao/aim/tgcm/data/tiegcm1.95/TGCM.tiegcm1.95.pcntr_mareqx_smin.nc 2: (disk file is /hao/aim/tgcm/data/tiegcm1.95/TGCM.tiegcm1.95.pcntr_mareqx_smin.nc) 2: Reading source history from diskfile /hao/aim/tgcm/data/tiegcm1.95/TGCM.tiegcm1.95.pcntr_mareqx_smin.nc: 2: nc_rdhist: seeking 80 0 0 found 80 0 0 n= 1 2: Note nc_rdhist: unused variable: calendar_advance 2: Note nc_rdhist: unused variable: write_date 2: Note nc_rdhist: unused variable: crit1 2: Note nc_rdhist: unused variable: crit2 2: Note nc_rdhist: unused variable: mag 2: Read field TN 3d subdomain min,max= 0.1506E+03 0.7773E+03 2: Read field UN 3d subdomain min,max= -0.8885E+04 0.2662E+05 2: Read field VN 3d subdomain min,max= -0.2545E+05 0.2588E+05 2: Read field O2 3d subdomain min,max= 0.6742E-04 0.2295E+00 2: Read field O1 3d subdomain min,max= 0.4386E-02 0.9947E+00 2: Read field N4S 3d subdomain min,max= 0.1000E-11 0.1571E-01 2: Read field NO 3d subdomain min,max= 0.1455E-06 0.1794E-03 2: Read field OP 3d subdomain min,max= 0.0000E+00 0.1053E+07 2: Read field N2D 3d subdomain min,max= 0.0000E+00 0.2296E-03 2: Read field TI 3d subdomain min,max= 0.1506E+03 0.1487E+04 2: Read field TE 3d subdomain min,max= 0.1506E+03 0.3095E+04 2: Read field NE 3d subdomain min,max= 0.5164E+03 0.1052E+07 2: Read field O2P 3d subdomain min,max= 0.0000E+00 0.8384E+05 2: Read field OMEGA 3d subdomain min,max= -0.4630E-04 0.4058E-04 2: Read field Z 3d subdomain min,max= 0.9604E+07 0.5080E+08 2: Read field POTEN 3d subdomain min,max= -0.1634E+05 0.9085E+04 2: Read field TN_NM 3d subdomain min,max= 0.1506E+03 0.7772E+03 2: Read field UN_NM 3d subdomain min,max= -0.8887E+04 0.2665E+05 2: Read field VN_NM 3d subdomain min,max= -0.2541E+05 0.2589E+05 2: Read field O2_NM 3d subdomain min,max= 0.6736E-04 0.2295E+00 2: Read field O1_NM 3d subdomain min,max= 0.4385E-02 0.9947E+00 2: Read field N4S_NM 3d subdomain min,max= 0.1000E-11 0.1571E-01 2: Read field NO_NM 3d subdomain min,max= 0.1457E-06 0.1794E-03 2: Note nc_rdhist: unused variable: OP_NM 2: 2: ------------------------------------------------------------------------ 2: Read TGCM PRIMARY HISTORY (source history) 2: Diskfile: /hao/aim/tgcm/data/tiegcm1.95/TGCM.tiegcm1.95.pcntr_mareqx_smin.nc 2: label = tiegcm mareqx smin control run 2: model_name = tiegcm 2: model_version = tiegcm1.95 2: create_date= 07/14/14 13:09:28 2: write_date = 2: logname = foster 2: host = ys6353 2: system = LINUX 2: run_type = 2: source_file = 2: output_file = 2: source_mtime = 80 0 0 2: initial_year = 2002 2: initial_day = 80 2: initial_mtime= 80 0 0 2: type = primary 2: ihist = 1 (nth history on history file) 2: delhmins= 0 (delta minutes between histories) 2: calendar year,day = 2002, 60 2: (model IS being advanced in calendar time) 2: modeltime = 80, 0, 0, 0 (model time day,hour,minute,seconds) 2: time = 2: ut = 0.00 (ut hours) 2: step = 120 (time step in seconds) 2: iter = 57600 (number of steps from 0,0,0) 2: nlat = 36 (number of latitudes) 2: nlon = 72 (number of longitudes) 2: nlev = 29 (number of levels) 2: zmtop = 7.250 (top midpoint level) 2: zmbot = -6.750 (bottom midpoint level) 2: zitop = 7.000 (top interface level) 2: zibot = -7.000 (bottom interface level) 2: dtide = 0.0E+00 0.0 (amp/phase of diurnal tide) 2: sdtide = 0.0E+00 0.0E+00 0.0E+00 0.0E+00 0.0E+00 0.0 0.0 0.0 0.0 0.0 2: (amp/phase of semi-diurnal tide) 2: f107d = 0.7000E+02 (daily solar flux) 2: f107a = 0.7000E+02 (average solar flux) 2: hpower = 0.1800E+02 (Gw) 2: ctpoten = 0.3000E+02 (Volts) 2: kp = 0.1000E+37 () 2: byimf = 0.0000E+00 2: bzimf = 0.1000E+37 2: swvel = 0.1000E+37 2: swden = 0.1000E+37 2: al = 0.1000E+37 2: e1,e2 = 0.1595E+01 0.2980E+01 (ergs/cm2/s) 2: alfac,alfad = 0.1000E+00 0.5000E+00 (keV) 2: ec,ed = 0.7212E-01 0.1200E-01 (ergs/cm2/s) 2: colfac = 0.15E+01 2: joulefac = 0.15E+01 2: p0 = 0.50E-06 2: p0_model= 0.50E-03 2: grav = 0.87E+03 2: nflds = 23 (number of model fields) 2: gswm_mi_di_ncfile = /glade/p/hao/tgcm/data/gswm_diurn_5.0d_99km.nc 2: gswm_mi_sdi_ncfile = /glade/p/hao/tgcm/data/gswm_semi_5.0d_99km.nc 2: gswm_nm_di_ncfile = [none] 2: gswm_nm_sdi_ncfile = [none] 2: see_ncfile = [none] 2: gpi_ncfile = [none] 2: ncep_ncfile = [none] 2: imf_ncfile = [none] 2: saber_ncfile = [none] 2: tidi_ncfile = [none] 2: tuv_lbc_intop= 0 (if 1, then lbc of t,u,v are stored in top k slot (old histories)) 2: LBC = -7.00 (lower boundary interface level) 2: ntask_mpi = 16 (number of MPI tasks) 2: coupled_cmit = 0 (1 if coupled with CISM/CMIT, 0 otherwise) 2: There are 23 fields on this history, as follows: 2: TN UN VN O2 O1 N4S NO OP 2: N2D TI TE NE O2P OMEGA Z POTEN 2: TN_NM UN_NM VN_NM O2_NM O1_NM N4S_NM NO_NM 2: ------------------------------------------------------------------------ 2: 2: readsource: Argon apparently not read from source history. Will init AR and AR_NM to ar_glb 2: readsource: Helium apparently not read from source history. Will init HE and HE_NM to pshelb= 0.1154E-05 2: 2: Enter advance: iter= 115200 nstep= 1440 2: 2: ------------------------------------------------------------------------ 2: Getfile: remote=/hao/aim/tgcm/data/gswm_diurn_5.0d_99km.nc 2: Getfile: Found file /hao/aim/tgcm/data/gswm_diurn_5.0d_99km.nc 2: ------------------------------------------------------------------------ 2: 2: 2: ------------------------------------------------------------------------ 2: Reading GSWM data file /hao/aim/tgcm/data/gswm_diurn_5.0d_99km.nc 2: Completed read from GSWM data file /hao/aim/tgcm/data/gswm_diurn_5.0d_99km.nc 2: 2: ------------------------------------------------------------------------ 2: 2: ------------------------------------------------------------------------ 2: Getfile: remote=/hao/aim/tgcm/data/gswm_semi_5.0d_99km.nc 2: Getfile: Found file /hao/aim/tgcm/data/gswm_semi_5.0d_99km.nc 2: ------------------------------------------------------------------------ 2: 2: 2: ------------------------------------------------------------------------ 2: Reading GSWM data file /hao/aim/tgcm/data/gswm_semi_5.0d_99km.nc 2: Completed read from GSWM data file /hao/aim/tgcm/data/gswm_semi_5.0d_99km.nc 2: 2: ------------------------------------------------------------------------ 2: 2: aurora_cons: 2: cusp: alfac= 0.100 ec= 0.072 fc= 0.2251E+09 2: drizzle: alfad= 0.500 ed= 0.012 fd= 0.7491E+07 2: auroral radius = max of rhp,rcp= 19.999 17.905 2: roth, rote (MLT) = 0.448 -0.072 2: 1/e-widths = h1,h2= 2.350 3.505 2: energy flux = e1,e2= 1.595 2.980 2: add_sproton = F 2: 2: mpitime_init: Completed initialization of mpi timing 2: Step 10 of 1440 mtime= 80 0 10 0 secs/step (sys) = 0.39 2: Step 20 of 1440 mtime= 80 0 20 0 secs/step (sys) = 0.40 2: Step 30 of 1440 mtime= 80 0 30 0 secs/step (sys) = 0.40 2: Step 40 of 1440 mtime= 80 0 40 0 secs/step (sys) = 0.44 2: Step 50 of 1440 mtime= 80 0 50 0 secs/step (sys) = 0.44 2: Step 60 of 1440 mtime= 80 1 0 0 secs/step (sys) = 0.41 2: 2: Allocated 3d sech field ZG(lon= 76,lat= 36,ilev= 29) 2: 2: Initialized diagnostic secondary history field ZG (ix= 23): 2: short_name = ZG 2: long_name = Geometric Height ZG 2: units = cm 2: geo = T 2: mag = F 2: dimnames = lon lat ilev 2: dimsizes = 76 36 29 2: ndims = 3 2: task0_only = F 2: 2: Allocated 3d sech field UI_ExB(lon= 76,lat= 36,ilev= 29) 2: 2: Initialized diagnostic secondary history field UI_ExB (ix= 17): 2: short_name = UI_ExB 2: long_name = UI: Zonal Ion Drift (ExB) 2: units = cm/s 2: geo = T 2: mag = F 2: dimnames = lon lat ilev 2: dimsizes = 76 36 29 2: ndims = 3 2: task0_only = F 2: 2: Allocated 3d sech field VI_ExB(lon= 76,lat= 36,ilev= 29) 2: 2: Initialized diagnostic secondary history field VI_ExB (ix= 18): 2: short_name = VI_ExB 2: long_name = VI: Meridional Ion Drift (ExB) 2: units = cm/s 2: geo = T 2: mag = F 2: dimnames = lon lat ilev 2: dimsizes = 76 36 29 2: ndims = 3 2: task0_only = F 2: 2: Allocated 3d sech field WI_ExB(lon= 76,lat= 36,ilev= 29) 2: 2: Initialized diagnostic secondary history field WI_ExB (ix= 19): 2: short_name = WI_ExB 2: long_name = WI: Vertical Ion Drift (ExB) 2: units = cm/s 2: geo = T 2: mag = F 2: dimnames = lon lat ilev 2: dimsizes = 76 36 29 2: ndims = 3 2: task0_only = F 2: 2: Allocated 2d sech field TEC(lon= 76,lat= 36) 2: 2: Initialized diagnostic secondary history field TEC (ix= 13): 2: short_name = TEC 2: long_name = TEC: Total Electron Content 2: units = 1/cm2 2: geo = T 2: mag = F 2: dimnames = lon lat 2: dimsizes = 76 36 0 2: ndims = 2 2: task0_only = F 2: 2: Allocated 3d sech field QJOULE(lon= 76,lat= 36,lev= 29) 2: 2: Initialized diagnostic secondary history field QJOULE (ix= 21): 2: short_name = QJOULE 2: long_name = QJOULE: Joule Heating 2: units = erg/g/s 2: geo = T 2: mag = F 2: dimnames = lon lat lev 2: dimsizes = 76 36 29 2: ndims = 3 2: task0_only = F 2: 2: Allocated 3d sech field DEN(lon= 76,lat= 36,ilev= 29) 2: 2: Initialized diagnostic secondary history field DEN (ix= 20): 2: short_name = DEN 2: long_name = Total Density 2: units = g/cm3 2: geo = T 2: mag = F 2: dimnames = lon lat ilev 2: dimsizes = 76 36 29 2: ndims = 3 2: task0_only = F 2: 2: Allocated 3d sech field N2(lon= 76,lat= 36,lev= 29) 2: 2: Initialized diagnostic secondary history field N2 (ix= 6): 2: short_name = N2 2: long_name = N2: Molecular Nitrogen 2: units = mmr 2: geo = T 2: mag = F 2: dimnames = lon lat lev 2: dimsizes = 76 36 29 2: ndims = 3 2: task0_only = F 2: 2: Allocated 3d sech field ZMAG(mlon= 81,mlat= 97,imlev= 32) 2: 2: Initialized diagnostic secondary history field ZMAG (ix= 24): 2: short_name = ZMAG 2: long_name = ZMAG from pdynamo 2: units = cm 2: geo = F 2: mag = T 2: dimnames = mlon mlat imlev 2: dimsizes = 81 97 32 2: ndims = 3 2: task0_only = F 2: Step 70 of 1440 mtime= 80 1 10 0 secs/step (sys) = 0.40 2: Step 80 of 1440 mtime= 80 1 20 0 secs/step (sys) = 0.46 2: Step 90 of 1440 mtime= 80 1 30 0 secs/step (sys) = 0.51 2: Step 100 of 1440 mtime= 80 1 40 0 secs/step (sys) = 0.39 2: Step 110 of 1440 mtime= 80 1 50 0 secs/step (sys) = 0.40 2: Step 120 of 1440 mtime= 80 2 0 0 secs/step (sys) = 0.40 2: Step 130 of 1440 mtime= 80 2 10 0 secs/step (sys) = 0.42 2: Step 140 of 1440 mtime= 80 2 20 0 secs/step (sys) = 0.39 2: Step 150 of 1440 mtime= 80 2 30 0 secs/step (sys) = 0.39 2: Step 160 of 1440 mtime= 80 2 40 0 secs/step (sys) = 0.40 2: Step 170 of 1440 mtime= 80 2 50 0 secs/step (sys) = 0.39 2: Step 180 of 1440 mtime= 80 3 0 0 secs/step (sys) = 0.39 2: Step 190 of 1440 mtime= 80 3 10 0 secs/step (sys) = 0.41 2: Step 200 of 1440 mtime= 80 3 20 0 secs/step (sys) = 0.40 2: Step 210 of 1440 mtime= 80 3 30 0 secs/step (sys) = 0.35 2: Step 220 of 1440 mtime= 80 3 40 0 secs/step (sys) = 0.35 2: Step 230 of 1440 mtime= 80 3 50 0 secs/step (sys) = 0.52 2: Step 240 of 1440 mtime= 80 4 0 0 secs/step (sys) = 0.41 2: Step 250 of 1440 mtime= 80 4 10 0 secs/step (sys) = 0.40 2: Step 260 of 1440 mtime= 80 4 20 0 secs/step (sys) = 0.43 2: Step 270 of 1440 mtime= 80 4 30 0 secs/step (sys) = 0.39 2: Step 280 of 1440 mtime= 80 4 40 0 secs/step (sys) = 1.31 2: Step 290 of 1440 mtime= 80 4 50 0 secs/step (sys) = 0.45 2: Step 300 of 1440 mtime= 80 5 0 0 secs/step (sys) = 0.55 2: Step 310 of 1440 mtime= 80 5 10 0 secs/step (sys) = 0.41 2: Step 320 of 1440 mtime= 80 5 20 0 secs/step (sys) = 0.41 2: Step 330 of 1440 mtime= 80 5 30 0 secs/step (sys) = 0.43 2: Step 340 of 1440 mtime= 80 5 40 0 secs/step (sys) = 0.43 2: Step 350 of 1440 mtime= 80 5 50 0 secs/step (sys) = 0.40 2: Step 360 of 1440 mtime= 80 6 0 0 secs/step (sys) = 0.56 2: Step 370 of 1440 mtime= 80 6 10 0 secs/step (sys) = 0.41 2: Step 380 of 1440 mtime= 80 6 20 0 secs/step (sys) = 0.39 2: Step 390 of 1440 mtime= 80 6 30 0 secs/step (sys) = 0.40 2: Step 400 of 1440 mtime= 80 6 40 0 secs/step (sys) = 0.39 2: Step 410 of 1440 mtime= 80 6 50 0 secs/step (sys) = 0.39 2: Step 420 of 1440 mtime= 80 7 0 0 secs/step (sys) = 0.45 2: Step 430 of 1440 mtime= 80 7 10 0 secs/step (sys) = 0.42 2: Step 440 of 1440 mtime= 80 7 20 0 secs/step (sys) = 0.42 2: Step 450 of 1440 mtime= 80 7 30 0 secs/step (sys) = 0.41 2: Step 460 of 1440 mtime= 80 7 40 0 secs/step (sys) = 0.44 2: Step 470 of 1440 mtime= 80 7 50 0 secs/step (sys) = 0.40 2: Step 480 of 1440 mtime= 80 8 0 0 secs/step (sys) = 0.41 2: Step 490 of 1440 mtime= 80 8 10 0 secs/step (sys) = 0.46 2: Step 500 of 1440 mtime= 80 8 20 0 secs/step (sys) = 0.40 2: Step 510 of 1440 mtime= 80 8 30 0 secs/step (sys) = 0.45 2: Step 520 of 1440 mtime= 80 8 40 0 secs/step (sys) = 0.40 2: Step 530 of 1440 mtime= 80 8 50 0 secs/step (sys) = 0.39 2: Step 540 of 1440 mtime= 80 9 0 0 secs/step (sys) = 0.45 2: Step 550 of 1440 mtime= 80 9 10 0 secs/step (sys) = 0.44 2: Step 560 of 1440 mtime= 80 9 20 0 secs/step (sys) = 0.46 2: Step 570 of 1440 mtime= 80 9 30 0 secs/step (sys) = 0.42 2: Step 580 of 1440 mtime= 80 9 40 0 secs/step (sys) = 0.46 2: Step 590 of 1440 mtime= 80 9 50 0 secs/step (sys) = 0.41 2: Step 600 of 1440 mtime= 80 10 0 0 secs/step (sys) = 0.43 2: Step 610 of 1440 mtime= 80 10 10 0 secs/step (sys) = 0.39 2: Step 620 of 1440 mtime= 80 10 20 0 secs/step (sys) = 0.41 2: Step 630 of 1440 mtime= 80 10 30 0 secs/step (sys) = 0.39 2: Step 640 of 1440 mtime= 80 10 40 0 secs/step (sys) = 0.45 2: Step 650 of 1440 mtime= 80 10 50 0 secs/step (sys) = 0.40 2: Step 660 of 1440 mtime= 80 11 0 0 secs/step (sys) = 0.44 2: Step 670 of 1440 mtime= 80 11 10 0 secs/step (sys) = 0.41 2: Step 680 of 1440 mtime= 80 11 20 0 secs/step (sys) = 0.40 2: Step 690 of 1440 mtime= 80 11 30 0 secs/step (sys) = 0.43 2: Step 700 of 1440 mtime= 80 11 40 0 secs/step (sys) = 0.40 2: Step 710 of 1440 mtime= 80 11 50 0 secs/step (sys) = 0.40 2: Step 720 of 1440 mtime= 80 12 0 0 secs/step (sys) = 0.39 2: Step 730 of 1440 mtime= 80 12 10 0 secs/step (sys) = 0.41 2: Step 740 of 1440 mtime= 80 12 20 0 secs/step (sys) = 0.50 2: Step 750 of 1440 mtime= 80 12 30 0 secs/step (sys) = 0.40 2: Step 760 of 1440 mtime= 80 12 40 0 secs/step (sys) = 0.39 2: Step 770 of 1440 mtime= 80 12 50 0 secs/step (sys) = 1.01 2: Step 780 of 1440 mtime= 80 13 0 0 secs/step (sys) = 0.43 2: Step 790 of 1440 mtime= 80 13 10 0 secs/step (sys) = 0.44 2: Step 800 of 1440 mtime= 80 13 20 0 secs/step (sys) = 0.40 2: Step 810 of 1440 mtime= 80 13 30 0 secs/step (sys) = 0.44 2: Step 820 of 1440 mtime= 80 13 40 0 secs/step (sys) = 0.48 2: Step 830 of 1440 mtime= 80 13 50 0 secs/step (sys) = 0.98 2: Step 840 of 1440 mtime= 80 14 0 0 secs/step (sys) = 0.49 2: Step 850 of 1440 mtime= 80 14 10 0 secs/step (sys) = 0.44 2: Step 860 of 1440 mtime= 80 14 20 0 secs/step (sys) = 0.43 2: Step 870 of 1440 mtime= 80 14 30 0 secs/step (sys) = 0.40 2: Step 880 of 1440 mtime= 80 14 40 0 secs/step (sys) = 0.41 2: Step 890 of 1440 mtime= 80 14 50 0 secs/step (sys) = 0.87 2: Step 900 of 1440 mtime= 80 15 0 0 secs/step (sys) = 0.40 2: Step 910 of 1440 mtime= 80 15 10 0 secs/step (sys) = 0.39 2: Step 920 of 1440 mtime= 80 15 20 0 secs/step (sys) = 0.40 2: Step 930 of 1440 mtime= 80 15 30 0 secs/step (sys) = 0.46 2: Step 940 of 1440 mtime= 80 15 40 0 secs/step (sys) = 0.41 2: Step 950 of 1440 mtime= 80 15 50 0 secs/step (sys) = 0.78 2: Step 960 of 1440 mtime= 80 16 0 0 secs/step (sys) = 0.50 2: Step 970 of 1440 mtime= 80 16 10 0 secs/step (sys) = 0.45 2: Step 980 of 1440 mtime= 80 16 20 0 secs/step (sys) = 0.43 2: Step 990 of 1440 mtime= 80 16 30 0 secs/step (sys) = 0.42 2: Step 1000 of 1440 mtime= 80 16 40 0 secs/step (sys) = 0.41 2: Step 1010 of 1440 mtime= 80 16 50 0 secs/step (sys) = 0.43 2: Step 1020 of 1440 mtime= 80 17 0 0 secs/step (sys) = 0.45 2: Step 1030 of 1440 mtime= 80 17 10 0 secs/step (sys) = 0.46 2: Step 1040 of 1440 mtime= 80 17 20 0 secs/step (sys) = 0.43 2: Step 1050 of 1440 mtime= 80 17 30 0 secs/step (sys) = 0.39 2: Step 1060 of 1440 mtime= 80 17 40 0 secs/step (sys) = 0.42 2: Step 1070 of 1440 mtime= 80 17 50 0 secs/step (sys) = 0.40 2: Step 1080 of 1440 mtime= 80 18 0 0 secs/step (sys) = 0.41 2: Step 1090 of 1440 mtime= 80 18 10 0 secs/step (sys) = 0.41 2: Step 1100 of 1440 mtime= 80 18 20 0 secs/step (sys) = 0.44 2: Step 1110 of 1440 mtime= 80 18 30 0 secs/step (sys) = 0.45 2: Step 1120 of 1440 mtime= 80 18 40 0 secs/step (sys) = 0.46 2: Step 1130 of 1440 mtime= 80 18 50 0 secs/step (sys) = 0.40 2: Step 1140 of 1440 mtime= 80 19 0 0 secs/step (sys) = 0.46 2: Step 1150 of 1440 mtime= 80 19 10 0 secs/step (sys) = 0.39 2: Step 1160 of 1440 mtime= 80 19 20 0 secs/step (sys) = 0.43 2: Step 1170 of 1440 mtime= 80 19 30 0 secs/step (sys) = 0.46 2: Step 1180 of 1440 mtime= 80 19 40 0 secs/step (sys) = 0.40 2: Step 1190 of 1440 mtime= 80 19 50 0 secs/step (sys) = 0.50 2: Step 1200 of 1440 mtime= 80 20 0 0 secs/step (sys) = 0.52 2: Step 1210 of 1440 mtime= 80 20 10 0 secs/step (sys) = 0.45 2: Step 1220 of 1440 mtime= 80 20 20 0 secs/step (sys) = 0.41 2: Step 1230 of 1440 mtime= 80 20 30 0 secs/step (sys) = 0.47 2: Step 1240 of 1440 mtime= 80 20 40 0 secs/step (sys) = 0.40 2: Step 1250 of 1440 mtime= 80 20 50 0 secs/step (sys) = 0.41 2: Step 1260 of 1440 mtime= 80 21 0 0 secs/step (sys) = 0.36 2: Step 1270 of 1440 mtime= 80 21 10 0 secs/step (sys) = 0.39 2: Step 1280 of 1440 mtime= 80 21 20 0 secs/step (sys) = 0.40 2: Step 1290 of 1440 mtime= 80 21 30 0 secs/step (sys) = 0.40 2: Step 1300 of 1440 mtime= 80 21 40 0 secs/step (sys) = 0.39 2: Step 1310 of 1440 mtime= 80 21 50 0 secs/step (sys) = 0.38 2: Step 1320 of 1440 mtime= 80 22 0 0 secs/step (sys) = 0.44 2: Step 1330 of 1440 mtime= 80 22 10 0 secs/step (sys) = 0.43 2: Step 1340 of 1440 mtime= 80 22 20 0 secs/step (sys) = 0.42 2: Step 1350 of 1440 mtime= 80 22 30 0 secs/step (sys) = 0.42 2: Step 1360 of 1440 mtime= 80 22 40 0 secs/step (sys) = 0.40 2: Step 1370 of 1440 mtime= 80 22 50 0 secs/step (sys) = 0.42 2: Step 1380 of 1440 mtime= 80 23 0 0 secs/step (sys) = 0.52 2: Step 1390 of 1440 mtime= 80 23 10 0 secs/step (sys) = 2.24 2: Step 1400 of 1440 mtime= 80 23 20 0 secs/step (sys) = 3.12 2: Step 1410 of 1440 mtime= 80 23 30 0 secs/step (sys) = 0.42 2: Step 1420 of 1440 mtime= 80 23 40 0 secs/step (sys) = 0.40 2: Step 1430 of 1440 mtime= 80 23 50 0 secs/step (sys) = 0.43 2: Advancing day (previous,present)= 80 2002 81 2002 sfeps= 0.1007E+01 2: Step 1440 of 1440 mtime= 81 0 0 0 secs/step (sys) = 0.41 2: 2: MPI run with ntask = 4 2: nstep= 1440 step= 60 2: Model simulation time = 86400 secs 2: (minutes= 1440.00, hours= 24.00, days= 1.000000) 2: Cpu time for run = 625.10 2: 2: ------------------------------------------------------------------------ 2: Total run time: mins= 0.1126E+02 hours= 0.1876E+00 2: Total MPI timing: mins= 0.2253E+01 hours= 0.3755E-01 %Total runtime= 20.01 2: Subroutine Time (mins) %Total mpi %Total run 2: mp_gather2root_prim 0.9546E-02 0.42 0.08 2: mp_gather2root_sech 0.1247E-01 0.55 0.11 2: mp_gather2root_lbc 0.3515E-04 0.00 0.00 2: mp_bndlats 0.1184E+00 5.26 1.05 2: mp_bndlats_f2d 0.2386E+00 10.59 2.12 2: mp_bndlons 0.7328E-01 3.25 0.65 2: mp_bndlons_f3d 0.2222E+00 9.86 1.97 2: mp_polelats 0.8874E-01 3.94 0.79 2: mp_polelat_f3d 0.8432E-01 3.74 0.75 2: mp_gatherlons_f3d 0.2707E+00 12.02 2.40 2: mp_scatterlons_f3d 0.1463E+00 6.49 1.30 2: mp_periodic_f4d 0.3523E-01 1.56 0.31 2: mp_periodic_f3d 0.1180E+00 5.24 1.05 2: mp_periodic_f2d 0.2692E-02 0.12 0.02 2: mp_bndlats_kmh 0.1057E-01 0.47 0.09 2: mp_bndlons_kmh 0.3183E-02 0.14 0.03 2: mp_mageq 0.1425E-01 0.63 0.13 2: mp_mageq_jpm1 0.2112E-01 0.94 0.19 2: mp_mageq_jpm3 0.7643E-02 0.34 0.07 2: mp_magpole_2d 0.6305E-01 2.80 0.56 2: mp_magpole_3d 0.6601E-04 0.00 0.00 2: mp_magpoles 0.2503E+00 11.11 2.22 2: mp_conjugate_points 0.3020E-06 0.00 0.00 2: mp_foldhem 0.1729E-01 0.77 0.15 2: mp_mag_periodic_f2d 0.1695E-01 0.75 0.15 2: mp_gather_pdyn 0.1724E-01 0.77 0.15 2: mp_mag_halos 0.8844E-01 3.93 0.79 2: mp_geo_halos 0.0000E+00 0.00 0.00 2: mp_geo_halos_f3d 0.1276E+00 5.66 1.13 2: mp_scatter_coeffs 0.0000E+00 0.00 0.00 2: mp_scatter_phim 0.1874E+00 8.32 1.66 2: mp_gather_f2d 0.4641E-03 0.02 0.00 2: mp_scatter_f2d 0.4136E-01 1.84 0.37 2: 2: ------------------------------------------------------------------------ 2: 2: ------------------------------------------------------------------------ 2: Report MPI wall-clock timing (mpi_wtime): mytid= 2 istep= 1440 nstep= 1440 2: Total Time (mins) in 'pefield ' = 0.3968E+00 ( 3.53% of total run time) 2: Total Time (mins) in 'dynamo_inputs ' = 0.2356E+01 ( 20.93% of total run time) 2: Total Time (mins) in 'pdynamo_fieldline_integrals ' = 0.1347E+00 ( 1.20% of total run time) 2: Total Time (mins) in 'pdynamo_complete_integrals ' = 0.9480E-01 ( 0.84% of total run time) 2: Total Time (mins) in 'pdynamo_rhspde ' = 0.6928E-03 ( 0.01% of total run time) 2: Total Time (mins) in 'highlat_poten ' = 0.2634E-03 ( 0.00% of total run time) 2: Total Time (mins) in 'pthreed ' = 0.1311E+01 ( 11.65% of total run time) 2: ------------------------------------------------------------------------ 2: 2: ------------------------------------------------------------------------ 2: TIMER (system_clock): 2: Elapsed run time = 675.43 (secs) 11.26 (mins) 0.19 (hrs) 2: 2: Elapsed init time = 10.12 secs, 0.003 hrs, 1.5% (includes source i/o) 2: Elapsed step time = 643.26 secs, 0.179 hrs, 95.2% 2: Elapsed prep time = 700.14 secs, 0.194 hrs, 103.7% 2: Elapsed apxparm time= 0.90 secs, 0.000 hrs, 100.0% 2: Elapsed dynamics = 323.25 secs, 0.090 hrs, 47.9% 2: Elapsed qrj = 39.32 secs, 0.011 hrs, 5.8% 2: Elapsed oplus = 47.13 secs, 0.013 hrs, 7.0% 2: Elapsed cmpminor = 47.04 secs, 0.013 hrs, 7.0% (n4s,n2d,no) 2: Elapsed cmpmajor = 37.23 secs, 0.010 hrs, 5.5% (o2, o, he) 2: Elapsed dt = 13.81 secs, 0.004 hrs, 2.0% 2: Elapsed duv = 29.80 secs, 0.008 hrs, 4.4% 2: Elapsed pdynamo = 109.89 secs, 0.031 hrs, 16.3% 2: Elapsed magpres_grav= 9.75 secs, 0.003 hrs, 1.4% 2: Elapsed i/o = 1.39 secs, 0.000 hrs, 0.2% 2: Elapsed Primary io = 0.09 secs, 0.000 hrs, 0.0% 2: Elapsed Secondary = 1.30 secs, 0.000 hrs, 0.2% 2: 2: Average secs/step = 0.45 2: Average mins/simulated day = 11.26 2: ------------------------------------------------------------------------ 2: 2: End execution of tiegcm_trunk at 07/07/15 09:30:59 2: NORMAL EXIT