3: 3: ======================================================================== 3: Begin execution of tiegcm_trunk at 07/07/15 09:19:42 3: Host = iris 3: System = LINUX 3: Logname = foster 3: ======================================================================== 3: 3: init_timer: level= 1 rtc=F sys=T 3: Current working directory (cwd) = /hao/aim/foster/tiegcm_rholev1/tiegcm_hao 3: Process ID (pid) = 2204 3: 3: Reading input data... 3: 3: Completed successful read of namelist inputs. 3: 3: Input: mkntask chose ntask_lon= 2 ntask_lat= 2 (ntask= 4) 3: Input: mkntask chose ntask_maglon= 2 ntask_maglat= 2 (ntask= 4) 3: Input: Using default joulefac = 1.50 3: Input: enforce_opfloor= 1 3: Will use the Heelis potential model 3: Note input: Setting BY to 0 with HEELIS potential model. 3: INPUT NOTE: adding mandatory field ZMAG to secondary history fields (field 24) 3: 3: ------------------------------------------------------------------------ 3: USER INPUT PARAMETERS: 3: label = tiegcm res=5.0 3: (optional text label for current run) 3: High-lat electric potential model: potential_model = HEELIS 3: gswm migrating diurnal file: gswm_mi_di_ncfile = /hao/aim/tgcm/data/gswm_diurn_5.0d_99km.nc 3: gswm migrating semi-diurnal file: gswm_mi_sdi_ncfile = /hao/aim/tgcm/data/gswm_semi_5.0d_99km.nc 3: calendar_advance = 1 (model will be advanced in calendar time starting on this day) 3: step = 60 (model timestep (seconds)) 3: ntask_lon = 2 (number of mpi tasks in longitude dimension) 3: ntask_lat = 2 (number of mpi tasks in latitude dimension) 3: total tasks = ntask_lon*ntask_lat = 4 3: source = /hao/aim/tgcm/data/tiegcm1.95/TGCM.tiegcm1.95.pcntr_mareqx_smin.nc 3: (file or mss path containing source history) 3: source_start = 80, 0, 0 (model time of source history) 3: output (primary history output files) = 3: tiegcm.p_rholev1_nobugfix.nc, 3: start (model start times) = 3: 80, 0, 0 3: stop (model stop times) = 3: 81, 0, 0 3: hist (primary history disk write frequencies) = 3: 1, 0, 0 3: Maxmimum number of histories per primary file = 10 3: secout (secondary history output files)= 3: tiegcm.s_rholev1_nobugfix.nc, 3: secstart (secondary history start times) = 3: 80, 1, 0 3: secstop (secondary history stop times) = 3: 81, 0, 0 3: sechist (secondary history disk write frequencies) = 3: 0, 1, 0 3: secflds (secondary history fields) = 3: TN UN VN O2 O1 3: N2 NO N4S HE NE 3: TE TI TEC O2P OMEGA 3: POTEN UI_ExB VI_ExB WI_ExB DEN 3: QJOULE Z ZG ZMAG 3: Maximum number of histories per secondary file = 24 3: Number of bytes for values of fields on secondary histories (sech_nbyte) = 4 3: dynamo = 1 (dynamo will be calculated) 3: current_pg = 1 (Add current due to plasma pressure gradient and gravity to rhs) 3: current_kq = 0 (Calculate height-integrated current density) 3: eddy_dif = 0 (DOY-dependent eddy diffusion flag) 3: tide (amplitudes and phases of semidiurnal tide) = 3: 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 3: tide2 (amplitude and phase of diurnal tide) = 3: 0.0E+00 0.00 3: aurora = 1 (0/1 flag for aurora) 3: colfac = 1.500 (collision factor) 3: joulefac = 1.500 (joule heating factor) 3: calc_helium = 1 (0/1 flag for helium) 3: If any of the following are spval ( 0.1000E+37), they will be calculated 3: during the simulation on a per timestep basis: 3: power = 0.1800E+02 (Hemispheric Power) 3: ctpoten= 0.3000E+02 (Cross-cap potential) 3: kp = 0.1000E+37 (Kp index) 3: bximf = 0.1000E+37 (BX component of IMF) 3: byimf = 0.0000E+00 (BY component of IMF) 3: bzimf = 0.1000E+37 (Bz component of IMF) 3: swvel = 0.1000E+37 (solar wind velocity) 3: swden = 0.1000E+37 (solar wind density) 3: f107 = 0.7000E+02 (F10.7 solar flux) 3: f107a = 0.7000E+02 (81-day ave F10.7 flux) 3: al = 0.1000E+37 (AL, lower auroral mag index) 3: END USER INPUT PARAMETERS 3: ------------------------------------------------------------------------ 3: 3: tgcm: dynamo= 1 -- calling apxparm: start_year= 2002 3: Time in apxparm = 0.903 (secs) 3: 3: ntask= 4 ntaski= 2 ntaskj= 2 Task Table: 3: j= -1 itask_table(:,j)= -1 -1 -1 -1 3: j= 0 itask_table(:,j)= -1 0 1 -1 3: j= 1 itask_table(:,j)= -1 2 3 -1 3: j= 2 itask_table(:,j)= -1 -1 -1 -1 3: 3: mytid= 3 mytidi,j= 1 1 lat0,1= 19 36 (18) lon0,1= 39 76 (38) ncells= 684 3: tgcm after distribute_geo 3: 3: ntask= 4 ntask_maglon= 2 ntask_maglat= 2 Task Table: 3: j= -1 itask_table_mag(:,j)= -1 -1 -1 -1 -1 3: j= 0 itask_table_mag(:,j)= 1 0 1 0 -1 3: j= 1 itask_table_mag(:,j)= 3 2 3 2 -1 3: j= 2 itask_table_mag(:,j)= -1 -1 -1 -1 -1 3: 3: mytid= 3 magtidi,j= 1 1 mlat0,1= 50 97 (48) mlon0,1= 42 81 (40) ncells=1920 3: mp_distribute_geo: nmagtaski= 2 tidcol= 1 3: tgcm after distribute_mag 3: 3: Task 3: 3: 3: Subdomain on geographic grid: 3: tasks( 3)%mytid = 3 3: tasks( 3)%mytidi= 1 3: tasks( 3)%mytidj= 1 3: tasks( 3)%nlats = 18 3: tasks( 3)%nlons = 38 3: tasks( 3)%lat0 = 19 3: tasks( 3)%lat1 = 36 3: tasks( 3)%lon0 = 39 3: tasks( 3)%lon1 = 76 3: Number of geo subdomain grid points = 684 3: 3: Subdomain on geomagnetic grid: 3: tasks( 3)%magtidi= 1 3: tasks( 3)%magtidj= 1 3: tasks( 3)%nmaglats = 48 3: tasks( 3)%nmaglons = 40 3: tasks( 3)%mlat0 = 50 3: tasks( 3)%mlat1 = 97 3: tasks( 3)%mlon0 = 42 3: tasks( 3)%mlon1 = 81 3: Number of mag subdomain grid points = 1920 3: mp_exchange_tasks: mxmaglat= 49 mxmaglon= 41 3: tgcm after esmf_init 3: Model version = tiegcm_trunk 3: 3: Set constants: 3: nlat= 36 nlon= 72 nlev= 28 3: dz= 0.50 3: dlat= 5.00 dlon= 5.00 3: zbound (cm) = 0.963723E+07 3: zmbot, zmtop = -6.750 7.250 (bottom,top midpoint levels) 3: zibot, zitop = -7.000 7.000 (bottom,top interface levels) 3: dt = 60.00 secs 3: grav = 870.00 3: freq_3m3 = 0.3506E-04 freq_semidi= 0.1454E-03 3: dipmin = 0.170 3: check_exp = F 3: 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 3: init: iyear= 2002 iday= 80 3: gswm_mi_di_ncfile = /hao/aim/tgcm/data/gswm_diurn_5.0d_99km.nc 3: gswm_mi_sdi_ncfile = /hao/aim/tgcm/data/gswm_semi_5.0d_99km.nc 3: hist_init: nstep= 1440 3: alloc_ar: allocated module data 3: Allocated gswm t,u,v,z (lon0:lon1,lat0:lat1) 3: Allocated private gwm t,u,v,z (lon0:lon1,lat0:lat1,nmonth,nhour) 3: alloc_pdyn: allocate fmeq_in with mlev0,1= -2 29 3: allocdata: all tasks allocate fzg: nlevp1= 29 nlonp4= 76 nlat= 36 3: init_lbc: allocated subdomains tlbc, ulbc, vlbc 3: init_lbc_nm: allocated subdomains tlbc_nm, ulbc_nm, vlbc_nm 3: init_lbc: allocated globals tlbc_glb, ulbc_glb, vlbc_glb: nlonp4= 76 nlat= 36 3: init_lbc: allocated globals tlbc_nm_glb, ulbc_nm_glb, vlbc_nm_glb: nlonp4= 76 nlat= 36 3: 3: Model run initialization: 3: nstep = 1440 (Number of time steps this run) 3: iter = 115200 (Initial iteration number) 3: iyear = 2002 (Beginning calendar year) 3: iday = 80 (Beginning calendar day) 3: igswm_mi_di = 1 (If > 0, GSWM diurnal tidal database will be used.) 3: igswm_mi_sdi= 1 (If > 0, GSWM semidiurnal tidal database will be used.) 3: igswm_nm_di= 0 (If > 0, GSWM nonmigrating diurnal tidal database will be used.) 3: igswm_nm_sdi= 0 (If > 0, GSWM nonmigrating semidiurnal tidal database will be used.) 3: 3: This is an initial run: 3: start_year = 2002 (Starting year of initial run) 3: start_day = 80 (Starting day of initial run) 3: start_mtime= 80 0 0 (Starting mtime of initial run) 3: 3: Primary Histories: 3: nsource = 1 (If > 0, a primary source history was provided) 3: nseries_prim = 1 (Number of primary time series) 3: nhist_total = 2 (Number of primary histories to be written) 3: nfiles_prim = 1 (Number of primary output files to be written) 3: mxhist_prim = 10 (Maximum number of primary histories per file) 3: 3: Secondary Histories: 3: nseries_sech = 1 (Number of secondary time series) 3: nsech_total = 24 (Number of secondary histories to be written) 3: nfiles_sech = 1 (Number of secondary output files to be written) 3: secondary history field 6: N2 3: secondary history field 7: NO 3: secondary history field 8: N4S 3: secondary history field 9: HE 3: secondary history field 10: NE 3: secondary history field 11: TE 3: secondary history field 12: TI 3: secondary history field 13: TEC 3: secondary history field 14: O2P 3: secondary history field 15: OMEGA 3: secondary history field 16: POTEN 3: secondary history field 17: UI_ExB 3: secondary history field 18: VI_ExB 3: secondary history field 19: WI_ExB 3: secondary history field 20: DEN 3: secondary history field 21: QJOULE 3: secondary history field 22: Z 3: secondary history field 23: ZG 3: secondary history field 24: ZMAG 3: 3: ------------------------------------------------------------------------------------------------ 3: Table of Available Diagnostic Fields: 3: Shortnames may be added to namelist SECFLDS 3: 3: Field Shortname Units Levels Caller Longname 3: - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 3: 1 CO2_COOL erg/g/s lev newton.F CO2 Cooling 3: 2 NO_COOL erg/g/s lev newton.F NO Cooling 3: 3 DEN g/cm3 ilev dt.F Total Density 3: 4 HEATING erg/g/s lev dt.F Total Heating 3: 5 HMF2 km none elden.F HMF2: Height of the F2 Layer 3: 6 NMF2 1/cm3 none elden.F NMF2: Peak Density of the F2 Layer 3: 7 FOF2 MHz none elden.F FOF2: Critical Frequency of the F2 Layer 3: 8 JE13D A/m2 mlev current.F JE13D: Eastward current density (3d) 3: 9 JE23D A/m2 mlev current.F JE23D: Downward current density (3d) 3: 10 JQR A/m2 none current.F JQR: Upward current density (2d) 3: 11 KQLAM A/m none current.F KQLAM: Height-integrated current density (+north) 3: 12 KQPHI A/m none current.F KQPHI: Height-integrated current density (+east) 3: 13 LAMDA_HAL 1/s lev lamdas.F LAMDA_HAL: Hall Ion Drag Coefficient 3: 14 LAMDA_PED 1/s lev lamdas.F LAMDA_PED: Pedersen Ion Drag Coefficient 3: 15 MU_M g/cm/s lev cpktkm.F MU_M: Molecular Viscosity Coefficient 3: 16 QJOULE erg/g/s lev qjoule.F QJOULE: Joule Heating 3: 17 SCHT km lev addiag.F SCHT: Pressure Scale Height 3: 18 SIGMA_HAL S/m lev lamdas.F SIGMA_HAL: Hall Conductivity 3: 19 SIGMA_PED S/m lev lamdas.F SIGMA_PED: Pedersen Conductivity 3: 20 TEC 1/cm2 none elden.F TEC: Total Electron Content 3: 21 UI_ExB cm/s ilev ionvel.F UI: Zonal Ion Drift (ExB) 3: 22 VI_ExB cm/s ilev ionvel.F VI: Meridional Ion Drift (ExB) 3: 23 WI_ExB cm/s ilev ionvel.F WI: Vertical Ion Drift (ExB) 3: 24 WN cm/s ilev swdot.F WN: Neutral Vertical Wind (plus up) 3: 25 O_N2 none lev comp.F O/N2 RATIO 3: 26 QJOULE_INTEG erg/cm2/s none qjoule.F Height-integrated Joule Heating 3: 27 BX none oplus.F BX/BMAG: Normalized eastward component of magnetic field 3: 28 BY none oplus.F BY/BMAG: Normalized northward component of magnetic field 3: 29 BZ none oplus.F BZ/BMAG: Normalized upward component of magnetic field 3: 30 BMAG Gauss none oplus.F BMAG: Magnetic field magnitude 3: 31 EX V/m ilev ionvel.F EX: Zonal component of electric field 3: 32 EY V/m ilev ionvel.F EY: Meridional component of electric field 3: 33 EZ V/m ilev ionvel.F EZ: Vertical component of electric field 3: 34 ED1 V/m imlev dynamo.F ED1: Magnetic eastward component of electric field 3: 35 ED2 V/m imlev dynamo.F ED2: Magnetic downward (equatorward) component of electric field 3: 36 PHIM2D V none dynamo.F PHIM2D: 2d Electric Potential on magnetic grid 3: 37 N2 mmr lev comp.F N2: Molecular Nitrogen 3: 38 ZGMID cm lev addiag.F ZGMID: Geometric Height at midpoints 3: 39 CUSP erg/cm2/s none dynamics.F CUSP (cusp2d*ec) 3: 40 DRIZZLE erg/cm2/s none dynamics.F DRIZZLE (drzl2d*ed) 3: 41 ALFA keV none dynamics.F ALFA 3: 42 NFLUX #/cm2/s none dynamics.F NFLUX 3: 43 EFLUX erg/cm2/s none dynamics.F EFLUX 3: ------------------------------------------------------------------------------------------------ 3: 3: 3: ------------------------------------------------------------------------ 3: Getfile: remote=/hao/aim/tgcm/data/tiegcm1.95/TGCM.tiegcm1.95.pcntr_mareqx_smin.nc 3: Getfile: Found file /hao/aim/tgcm/data/tiegcm1.95/TGCM.tiegcm1.95.pcntr_mareqx_smin.nc 3: ------------------------------------------------------------------------ 3: 3: Acquired source history file /hao/aim/tgcm/data/tiegcm1.95/TGCM.tiegcm1.95.pcntr_mareqx_smin.nc 3: (disk file is /hao/aim/tgcm/data/tiegcm1.95/TGCM.tiegcm1.95.pcntr_mareqx_smin.nc) 3: Reading source history from diskfile /hao/aim/tgcm/data/tiegcm1.95/TGCM.tiegcm1.95.pcntr_mareqx_smin.nc: 3: nc_rdhist: seeking 80 0 0 found 80 0 0 n= 1 3: Note nc_rdhist: unused variable: calendar_advance 3: Note nc_rdhist: unused variable: write_date 3: Note nc_rdhist: unused variable: crit1 3: Note nc_rdhist: unused variable: crit2 3: Note nc_rdhist: unused variable: mag 3: Read field TN 3d subdomain min,max= 0.1518E+03 0.7299E+03 3: Read field UN 3d subdomain min,max= -0.2201E+05 0.9620E+04 3: Read field VN 3d subdomain min,max= -0.2566E+05 0.2567E+05 3: Read field O2 3d subdomain min,max= 0.6157E-04 0.2295E+00 3: Read field O1 3d subdomain min,max= 0.4387E-02 0.9953E+00 3: Read field N4S 3d subdomain min,max= 0.1000E-11 0.1311E-01 3: Read field NO 3d subdomain min,max= 0.1370E-06 0.1792E-03 3: Read field OP 3d subdomain min,max= 0.0000E+00 0.8510E+06 3: Read field N2D 3d subdomain min,max= 0.0000E+00 0.2317E-03 3: Read field TI 3d subdomain min,max= 0.1517E+03 0.1486E+04 3: Read field TE 3d subdomain min,max= 0.1517E+03 0.3067E+04 3: Read field NE 3d subdomain min,max= 0.5080E+03 0.8466E+06 3: Read field O2P 3d subdomain min,max= 0.0000E+00 0.8257E+05 3: Read field OMEGA 3d subdomain min,max= -0.4384E-04 0.2841E-04 3: Read field Z 3d subdomain min,max= 0.9583E+07 0.4913E+08 3: Read field POTEN 3d subdomain min,max= -0.1899E+04 0.1249E+05 3: Read field TN_NM 3d subdomain min,max= 0.1517E+03 0.7295E+03 3: Read field UN_NM 3d subdomain min,max= -0.2200E+05 0.9704E+04 3: Read field VN_NM 3d subdomain min,max= -0.2567E+05 0.2564E+05 3: Read field O2_NM 3d subdomain min,max= 0.6155E-04 0.2295E+00 3: Read field O1_NM 3d subdomain min,max= 0.4386E-02 0.9953E+00 3: Read field N4S_NM 3d subdomain min,max= 0.1000E-11 0.1306E-01 3: Read field NO_NM 3d subdomain min,max= 0.1372E-06 0.1793E-03 3: Note nc_rdhist: unused variable: OP_NM 3: Read LBC from source history: i= 89 sh%lbc= -7.00 3: 3: ------------------------------------------------------------------------ 3: Read TGCM PRIMARY HISTORY (source history) 3: Diskfile: /hao/aim/tgcm/data/tiegcm1.95/TGCM.tiegcm1.95.pcntr_mareqx_smin.nc 3: label = tiegcm mareqx smin control run 3: model_name = tiegcm 3: model_version = tiegcm1.95 3: create_date= 07/14/14 13:09:28 3: write_date = 3: logname = foster 3: host = ys6353 3: system = LINUX 3: run_type = 3: source_file = 3: output_file = 3: source_mtime = 80 0 0 3: initial_year = 2002 3: initial_day = 80 3: initial_mtime= 80 0 0 3: type = primary 3: ihist = 1 (nth history on history file) 3: delhmins= 0 (delta minutes between histories) 3: calendar year,day = 2002, 60 3: (model IS being advanced in calendar time) 3: modeltime = 80, 0, 0, 0 (model time day,hour,minute,seconds) 3: time = 3: ut = 0.00 (ut hours) 3: step = 120 (time step in seconds) 3: iter = 57600 (number of steps from 0,0,0) 3: nlat = 36 (number of latitudes) 3: nlon = 72 (number of longitudes) 3: nlev = 29 (number of levels) 3: zmtop = 7.250 (top midpoint level) 3: zmbot = -6.750 (bottom midpoint level) 3: zitop = 7.000 (top interface level) 3: zibot = -7.000 (bottom interface level) 3: dtide = 0.0E+00 0.0 (amp/phase of diurnal tide) 3: 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 3: (amp/phase of semi-diurnal tide) 3: f107d = 0.7000E+02 (daily solar flux) 3: f107a = 0.7000E+02 (average solar flux) 3: hpower = 0.1800E+02 (Gw) 3: ctpoten = 0.3000E+02 (Volts) 3: kp = 0.1000E+37 () 3: byimf = 0.0000E+00 3: bzimf = 0.1000E+37 3: swvel = 0.1000E+37 3: swden = 0.1000E+37 3: al = 0.1000E+37 3: e1,e2 = 0.1595E+01 0.2980E+01 (ergs/cm2/s) 3: alfac,alfad = 0.1000E+00 0.5000E+00 (keV) 3: ec,ed = 0.7212E-01 0.1200E-01 (ergs/cm2/s) 3: colfac = 0.15E+01 3: joulefac = 0.15E+01 3: p0 = 0.50E-06 3: p0_model= 0.50E-03 3: grav = 0.87E+03 3: nflds = 23 (number of model fields) 3: gswm_mi_di_ncfile = /glade/p/hao/tgcm/data/gswm_diurn_5.0d_99km.nc 3: gswm_mi_sdi_ncfile = /glade/p/hao/tgcm/data/gswm_semi_5.0d_99km.nc 3: gswm_nm_di_ncfile = [none] 3: gswm_nm_sdi_ncfile = [none] 3: see_ncfile = [none] 3: gpi_ncfile = [none] 3: ncep_ncfile = [none] 3: imf_ncfile = [none] 3: saber_ncfile = [none] 3: tidi_ncfile = [none] 3: tuv_lbc_intop= 0 (if 1, then lbc of t,u,v are stored in top k slot (old histories)) 3: LBC = -7.00 (lower boundary interface level) 3: ntask_mpi = 16 (number of MPI tasks) 3: coupled_cmit = 0 (1 if coupled with CISM/CMIT, 0 otherwise) 3: There are 23 fields on this history, as follows: 3: TN UN VN O2 O1 N4S NO OP 3: N2D TI TE NE O2P OMEGA Z POTEN 3: TN_NM UN_NM VN_NM O2_NM O1_NM N4S_NM NO_NM 3: ------------------------------------------------------------------------ 3: 3: readsource: Argon apparently not read from source history. Will init AR and AR_NM to ar_glb 3: readsource: Helium apparently not read from source history. Will init HE and HE_NM to pshelb= 0.1154E-05 3: 3: Enter advance: iter= 115200 nstep= 1440 3: 3: ------------------------------------------------------------------------ 3: Getfile: remote=/hao/aim/tgcm/data/gswm_diurn_5.0d_99km.nc 3: Getfile: Found file /hao/aim/tgcm/data/gswm_diurn_5.0d_99km.nc 3: ------------------------------------------------------------------------ 3: 3: 3: ------------------------------------------------------------------------ 3: Reading GSWM data file /hao/aim/tgcm/data/gswm_diurn_5.0d_99km.nc 3: Completed read from GSWM data file /hao/aim/tgcm/data/gswm_diurn_5.0d_99km.nc 3: 3: ------------------------------------------------------------------------ 3: 3: ------------------------------------------------------------------------ 3: Getfile: remote=/hao/aim/tgcm/data/gswm_semi_5.0d_99km.nc 3: Getfile: Found file /hao/aim/tgcm/data/gswm_semi_5.0d_99km.nc 3: ------------------------------------------------------------------------ 3: 3: 3: ------------------------------------------------------------------------ 3: Reading GSWM data file /hao/aim/tgcm/data/gswm_semi_5.0d_99km.nc 3: Completed read from GSWM data file /hao/aim/tgcm/data/gswm_semi_5.0d_99km.nc 3: 3: ------------------------------------------------------------------------ 3: 3: aurora_cons: 3: cusp: alfac= 0.100 ec= 0.072 fc= 0.2251E+09 3: drizzle: alfad= 0.500 ed= 0.012 fd= 0.7491E+07 3: auroral radius = max of rhp,rcp= 19.999 17.905 3: roth, rote (MLT) = 0.448 -0.072 3: 1/e-widths = h1,h2= 2.350 3.505 3: energy flux = e1,e2= 1.595 2.980 3: add_sproton = F 3: 3: mpitime_init: Completed initialization of mpi timing 3: Step 10 of 1440 mtime= 80 0 10 0 secs/step (sys) = 0.39 3: Step 20 of 1440 mtime= 80 0 20 0 secs/step (sys) = 0.40 3: Step 30 of 1440 mtime= 80 0 30 0 secs/step (sys) = 0.40 3: Step 40 of 1440 mtime= 80 0 40 0 secs/step (sys) = 0.44 3: Step 50 of 1440 mtime= 80 0 50 0 secs/step (sys) = 0.44 3: Step 60 of 1440 mtime= 80 1 0 0 secs/step (sys) = 0.41 3: 3: Allocated 3d sech field ZG(lon= 76,lat= 36,ilev= 29) 3: 3: Initialized diagnostic secondary history field ZG (ix= 23): 3: short_name = ZG 3: long_name = Geometric Height ZG 3: units = cm 3: geo = T 3: mag = F 3: dimnames = lon lat ilev 3: dimsizes = 76 36 29 3: ndims = 3 3: task0_only = F 3: 3: Allocated 3d sech field UI_ExB(lon= 76,lat= 36,ilev= 29) 3: 3: Initialized diagnostic secondary history field UI_ExB (ix= 17): 3: short_name = UI_ExB 3: long_name = UI: Zonal Ion Drift (ExB) 3: units = cm/s 3: geo = T 3: mag = F 3: dimnames = lon lat ilev 3: dimsizes = 76 36 29 3: ndims = 3 3: task0_only = F 3: 3: Allocated 3d sech field VI_ExB(lon= 76,lat= 36,ilev= 29) 3: 3: Initialized diagnostic secondary history field VI_ExB (ix= 18): 3: short_name = VI_ExB 3: long_name = VI: Meridional Ion Drift (ExB) 3: units = cm/s 3: geo = T 3: mag = F 3: dimnames = lon lat ilev 3: dimsizes = 76 36 29 3: ndims = 3 3: task0_only = F 3: 3: Allocated 3d sech field WI_ExB(lon= 76,lat= 36,ilev= 29) 3: 3: Initialized diagnostic secondary history field WI_ExB (ix= 19): 3: short_name = WI_ExB 3: long_name = WI: Vertical Ion Drift (ExB) 3: units = cm/s 3: geo = T 3: mag = F 3: dimnames = lon lat ilev 3: dimsizes = 76 36 29 3: ndims = 3 3: task0_only = F 3: 3: Allocated 2d sech field TEC(lon= 76,lat= 36) 3: 3: Initialized diagnostic secondary history field TEC (ix= 13): 3: short_name = TEC 3: long_name = TEC: Total Electron Content 3: units = 1/cm2 3: geo = T 3: mag = F 3: dimnames = lon lat 3: dimsizes = 76 36 0 3: ndims = 2 3: task0_only = F 3: 3: Allocated 3d sech field QJOULE(lon= 76,lat= 36,lev= 29) 3: 3: Initialized diagnostic secondary history field QJOULE (ix= 21): 3: short_name = QJOULE 3: long_name = QJOULE: Joule Heating 3: units = erg/g/s 3: geo = T 3: mag = F 3: dimnames = lon lat lev 3: dimsizes = 76 36 29 3: ndims = 3 3: task0_only = F 3: 3: Allocated 3d sech field DEN(lon= 76,lat= 36,ilev= 29) 3: 3: Initialized diagnostic secondary history field DEN (ix= 20): 3: short_name = DEN 3: long_name = Total Density 3: units = g/cm3 3: geo = T 3: mag = F 3: dimnames = lon lat ilev 3: dimsizes = 76 36 29 3: ndims = 3 3: task0_only = F 3: 3: Allocated 3d sech field N2(lon= 76,lat= 36,lev= 29) 3: 3: Initialized diagnostic secondary history field N2 (ix= 6): 3: short_name = N2 3: long_name = N2: Molecular Nitrogen 3: units = mmr 3: geo = T 3: mag = F 3: dimnames = lon lat lev 3: dimsizes = 76 36 29 3: ndims = 3 3: task0_only = F 3: 3: Allocated 3d sech field ZMAG(mlon= 81,mlat= 97,imlev= 32) 3: 3: Initialized diagnostic secondary history field ZMAG (ix= 24): 3: short_name = ZMAG 3: long_name = ZMAG from pdynamo 3: units = cm 3: geo = F 3: mag = T 3: dimnames = mlon mlat imlev 3: dimsizes = 81 97 32 3: ndims = 3 3: task0_only = F 3: Step 70 of 1440 mtime= 80 1 10 0 secs/step (sys) = 0.40 3: Step 80 of 1440 mtime= 80 1 20 0 secs/step (sys) = 0.46 3: Step 90 of 1440 mtime= 80 1 30 0 secs/step (sys) = 0.51 3: Step 100 of 1440 mtime= 80 1 40 0 secs/step (sys) = 0.39 3: Step 110 of 1440 mtime= 80 1 50 0 secs/step (sys) = 0.40 3: Step 120 of 1440 mtime= 80 2 0 0 secs/step (sys) = 0.40 3: Step 130 of 1440 mtime= 80 2 10 0 secs/step (sys) = 0.42 3: Step 140 of 1440 mtime= 80 2 20 0 secs/step (sys) = 0.39 3: Step 150 of 1440 mtime= 80 2 30 0 secs/step (sys) = 0.39 3: Step 160 of 1440 mtime= 80 2 40 0 secs/step (sys) = 0.40 3: Step 170 of 1440 mtime= 80 2 50 0 secs/step (sys) = 0.39 3: Step 180 of 1440 mtime= 80 3 0 0 secs/step (sys) = 0.39 3: Step 190 of 1440 mtime= 80 3 10 0 secs/step (sys) = 0.41 3: Step 200 of 1440 mtime= 80 3 20 0 secs/step (sys) = 0.40 3: Step 210 of 1440 mtime= 80 3 30 0 secs/step (sys) = 0.35 3: Step 220 of 1440 mtime= 80 3 40 0 secs/step (sys) = 0.35 3: Step 230 of 1440 mtime= 80 3 50 0 secs/step (sys) = 0.52 3: Step 240 of 1440 mtime= 80 4 0 0 secs/step (sys) = 0.41 3: Step 250 of 1440 mtime= 80 4 10 0 secs/step (sys) = 0.40 3: Step 260 of 1440 mtime= 80 4 20 0 secs/step (sys) = 0.43 3: Step 270 of 1440 mtime= 80 4 30 0 secs/step (sys) = 0.39 3: Step 280 of 1440 mtime= 80 4 40 0 secs/step (sys) = 1.31 3: Step 290 of 1440 mtime= 80 4 50 0 secs/step (sys) = 0.45 3: Step 300 of 1440 mtime= 80 5 0 0 secs/step (sys) = 0.55 3: Step 310 of 1440 mtime= 80 5 10 0 secs/step (sys) = 0.41 3: Step 320 of 1440 mtime= 80 5 20 0 secs/step (sys) = 0.41 3: Step 330 of 1440 mtime= 80 5 30 0 secs/step (sys) = 0.43 3: Step 340 of 1440 mtime= 80 5 40 0 secs/step (sys) = 0.43 3: Step 350 of 1440 mtime= 80 5 50 0 secs/step (sys) = 0.40 3: Step 360 of 1440 mtime= 80 6 0 0 secs/step (sys) = 0.56 3: Step 370 of 1440 mtime= 80 6 10 0 secs/step (sys) = 0.41 3: Step 380 of 1440 mtime= 80 6 20 0 secs/step (sys) = 0.39 3: Step 390 of 1440 mtime= 80 6 30 0 secs/step (sys) = 0.40 3: Step 400 of 1440 mtime= 80 6 40 0 secs/step (sys) = 0.39 3: Step 410 of 1440 mtime= 80 6 50 0 secs/step (sys) = 0.39 3: Step 420 of 1440 mtime= 80 7 0 0 secs/step (sys) = 0.45 3: Step 430 of 1440 mtime= 80 7 10 0 secs/step (sys) = 0.42 3: Step 440 of 1440 mtime= 80 7 20 0 secs/step (sys) = 0.42 3: Step 450 of 1440 mtime= 80 7 30 0 secs/step (sys) = 0.41 3: Step 460 of 1440 mtime= 80 7 40 0 secs/step (sys) = 0.44 3: Step 470 of 1440 mtime= 80 7 50 0 secs/step (sys) = 0.40 3: Step 480 of 1440 mtime= 80 8 0 0 secs/step (sys) = 0.41 3: Step 490 of 1440 mtime= 80 8 10 0 secs/step (sys) = 0.46 3: Step 500 of 1440 mtime= 80 8 20 0 secs/step (sys) = 0.40 3: Step 510 of 1440 mtime= 80 8 30 0 secs/step (sys) = 0.45 3: Step 520 of 1440 mtime= 80 8 40 0 secs/step (sys) = 0.40 3: Step 530 of 1440 mtime= 80 8 50 0 secs/step (sys) = 0.39 3: Step 540 of 1440 mtime= 80 9 0 0 secs/step (sys) = 0.44 3: Step 550 of 1440 mtime= 80 9 10 0 secs/step (sys) = 0.44 3: Step 560 of 1440 mtime= 80 9 20 0 secs/step (sys) = 0.46 3: Step 570 of 1440 mtime= 80 9 30 0 secs/step (sys) = 0.42 3: Step 580 of 1440 mtime= 80 9 40 0 secs/step (sys) = 0.46 3: Step 590 of 1440 mtime= 80 9 50 0 secs/step (sys) = 0.41 3: Step 600 of 1440 mtime= 80 10 0 0 secs/step (sys) = 0.43 3: Step 610 of 1440 mtime= 80 10 10 0 secs/step (sys) = 0.39 3: Step 620 of 1440 mtime= 80 10 20 0 secs/step (sys) = 0.41 3: Step 630 of 1440 mtime= 80 10 30 0 secs/step (sys) = 0.39 3: Step 640 of 1440 mtime= 80 10 40 0 secs/step (sys) = 0.45 3: Step 650 of 1440 mtime= 80 10 50 0 secs/step (sys) = 0.40 3: Step 660 of 1440 mtime= 80 11 0 0 secs/step (sys) = 0.44 3: Step 670 of 1440 mtime= 80 11 10 0 secs/step (sys) = 0.41 3: Step 680 of 1440 mtime= 80 11 20 0 secs/step (sys) = 0.40 3: Step 690 of 1440 mtime= 80 11 30 0 secs/step (sys) = 0.43 3: Step 700 of 1440 mtime= 80 11 40 0 secs/step (sys) = 0.40 3: Step 710 of 1440 mtime= 80 11 50 0 secs/step (sys) = 0.40 3: Step 720 of 1440 mtime= 80 12 0 0 secs/step (sys) = 0.39 3: Step 730 of 1440 mtime= 80 12 10 0 secs/step (sys) = 0.41 3: Step 740 of 1440 mtime= 80 12 20 0 secs/step (sys) = 0.50 3: Step 750 of 1440 mtime= 80 12 30 0 secs/step (sys) = 0.40 3: Step 760 of 1440 mtime= 80 12 40 0 secs/step (sys) = 0.39 3: Step 770 of 1440 mtime= 80 12 50 0 secs/step (sys) = 1.01 3: Step 780 of 1440 mtime= 80 13 0 0 secs/step (sys) = 0.43 3: Step 790 of 1440 mtime= 80 13 10 0 secs/step (sys) = 0.44 3: Step 800 of 1440 mtime= 80 13 20 0 secs/step (sys) = 0.40 3: Step 810 of 1440 mtime= 80 13 30 0 secs/step (sys) = 0.44 3: Step 820 of 1440 mtime= 80 13 40 0 secs/step (sys) = 0.48 3: Step 830 of 1440 mtime= 80 13 50 0 secs/step (sys) = 0.98 3: Step 840 of 1440 mtime= 80 14 0 0 secs/step (sys) = 0.51 3: Step 850 of 1440 mtime= 80 14 10 0 secs/step (sys) = 0.44 3: Step 860 of 1440 mtime= 80 14 20 0 secs/step (sys) = 0.43 3: Step 870 of 1440 mtime= 80 14 30 0 secs/step (sys) = 0.40 3: Step 880 of 1440 mtime= 80 14 40 0 secs/step (sys) = 0.41 3: Step 890 of 1440 mtime= 80 14 50 0 secs/step (sys) = 0.87 3: Step 900 of 1440 mtime= 80 15 0 0 secs/step (sys) = 0.40 3: Step 910 of 1440 mtime= 80 15 10 0 secs/step (sys) = 0.39 3: Step 920 of 1440 mtime= 80 15 20 0 secs/step (sys) = 0.40 3: Step 930 of 1440 mtime= 80 15 30 0 secs/step (sys) = 0.46 3: Step 940 of 1440 mtime= 80 15 40 0 secs/step (sys) = 0.41 3: Step 950 of 1440 mtime= 80 15 50 0 secs/step (sys) = 0.78 3: Step 960 of 1440 mtime= 80 16 0 0 secs/step (sys) = 0.50 3: Step 970 of 1440 mtime= 80 16 10 0 secs/step (sys) = 0.45 3: Step 980 of 1440 mtime= 80 16 20 0 secs/step (sys) = 0.43 3: Step 990 of 1440 mtime= 80 16 30 0 secs/step (sys) = 0.42 3: Step 1000 of 1440 mtime= 80 16 40 0 secs/step (sys) = 0.41 3: Step 1010 of 1440 mtime= 80 16 50 0 secs/step (sys) = 0.43 3: Step 1020 of 1440 mtime= 80 17 0 0 secs/step (sys) = 0.45 3: Step 1030 of 1440 mtime= 80 17 10 0 secs/step (sys) = 0.46 3: Step 1040 of 1440 mtime= 80 17 20 0 secs/step (sys) = 0.43 3: Step 1050 of 1440 mtime= 80 17 30 0 secs/step (sys) = 0.39 3: Step 1060 of 1440 mtime= 80 17 40 0 secs/step (sys) = 0.42 3: Step 1070 of 1440 mtime= 80 17 50 0 secs/step (sys) = 0.40 3: Step 1080 of 1440 mtime= 80 18 0 0 secs/step (sys) = 0.41 3: Step 1090 of 1440 mtime= 80 18 10 0 secs/step (sys) = 0.41 3: Step 1100 of 1440 mtime= 80 18 20 0 secs/step (sys) = 0.44 3: Step 1110 of 1440 mtime= 80 18 30 0 secs/step (sys) = 0.45 3: Step 1120 of 1440 mtime= 80 18 40 0 secs/step (sys) = 0.46 3: Step 1130 of 1440 mtime= 80 18 50 0 secs/step (sys) = 0.40 3: Step 1140 of 1440 mtime= 80 19 0 0 secs/step (sys) = 0.46 3: Step 1150 of 1440 mtime= 80 19 10 0 secs/step (sys) = 0.39 3: Step 1160 of 1440 mtime= 80 19 20 0 secs/step (sys) = 0.43 3: Step 1170 of 1440 mtime= 80 19 30 0 secs/step (sys) = 0.46 3: Step 1180 of 1440 mtime= 80 19 40 0 secs/step (sys) = 0.40 3: Step 1190 of 1440 mtime= 80 19 50 0 secs/step (sys) = 0.50 3: Step 1200 of 1440 mtime= 80 20 0 0 secs/step (sys) = 0.52 3: Step 1210 of 1440 mtime= 80 20 10 0 secs/step (sys) = 0.45 3: Step 1220 of 1440 mtime= 80 20 20 0 secs/step (sys) = 0.41 3: Step 1230 of 1440 mtime= 80 20 30 0 secs/step (sys) = 0.47 3: Step 1240 of 1440 mtime= 80 20 40 0 secs/step (sys) = 0.40 3: Step 1250 of 1440 mtime= 80 20 50 0 secs/step (sys) = 0.41 3: Step 1260 of 1440 mtime= 80 21 0 0 secs/step (sys) = 0.36 3: Step 1270 of 1440 mtime= 80 21 10 0 secs/step (sys) = 0.39 3: Step 1280 of 1440 mtime= 80 21 20 0 secs/step (sys) = 0.40 3: Step 1290 of 1440 mtime= 80 21 30 0 secs/step (sys) = 0.40 3: Step 1300 of 1440 mtime= 80 21 40 0 secs/step (sys) = 0.39 3: Step 1310 of 1440 mtime= 80 21 50 0 secs/step (sys) = 0.38 3: Step 1320 of 1440 mtime= 80 22 0 0 secs/step (sys) = 0.44 3: Step 1330 of 1440 mtime= 80 22 10 0 secs/step (sys) = 0.43 3: Step 1340 of 1440 mtime= 80 22 20 0 secs/step (sys) = 0.42 3: Step 1350 of 1440 mtime= 80 22 30 0 secs/step (sys) = 0.41 3: Step 1360 of 1440 mtime= 80 22 40 0 secs/step (sys) = 0.40 3: Step 1370 of 1440 mtime= 80 22 50 0 secs/step (sys) = 0.42 3: Step 1380 of 1440 mtime= 80 23 0 0 secs/step (sys) = 0.52 3: Step 1390 of 1440 mtime= 80 23 10 0 secs/step (sys) = 2.24 3: Step 1400 of 1440 mtime= 80 23 20 0 secs/step (sys) = 3.12 3: Step 1410 of 1440 mtime= 80 23 30 0 secs/step (sys) = 0.42 3: Step 1420 of 1440 mtime= 80 23 40 0 secs/step (sys) = 0.40 3: Step 1430 of 1440 mtime= 80 23 50 0 secs/step (sys) = 0.43 3: Advancing day (previous,present)= 80 2002 81 2002 sfeps= 0.1007E+01 3: Step 1440 of 1440 mtime= 81 0 0 0 secs/step (sys) = 0.41 3: 3: MPI run with ntask = 4 3: nstep= 1440 step= 60 3: Model simulation time = 86400 secs 3: (minutes= 1440.00, hours= 24.00, days= 1.000000) 3: Cpu time for run = 598.93 3: 3: ------------------------------------------------------------------------ 3: Total run time: mins= 0.1126E+02 hours= 0.1876E+00 3: Total MPI timing: mins= 0.2712E+01 hours= 0.4520E-01 %Total runtime= 24.09 3: Subroutine Time (mins) %Total mpi %Total run 3: mp_gather2root_prim 0.1186E-01 0.44 0.11 3: mp_gather2root_sech 0.1075E-01 0.40 0.10 3: mp_gather2root_lbc 0.5541E-04 0.00 0.00 3: mp_bndlats 0.1251E+00 4.61 1.11 3: mp_bndlats_f2d 0.2402E+00 8.86 2.13 3: mp_bndlons 0.6516E-01 2.40 0.58 3: mp_bndlons_f3d 0.2385E+00 8.79 2.12 3: mp_polelats 0.9642E-01 3.56 0.86 3: mp_polelat_f3d 0.8400E-01 3.10 0.75 3: mp_gatherlons_f3d 0.1568E+00 5.78 1.39 3: mp_scatterlons_f3d 0.6937E+00 25.58 6.16 3: mp_periodic_f4d 0.3222E-01 1.19 0.29 3: mp_periodic_f3d 0.1145E+00 4.22 1.02 3: mp_periodic_f2d 0.3281E-02 0.12 0.03 3: mp_bndlats_kmh 0.1106E-01 0.41 0.10 3: mp_bndlons_kmh 0.2435E-02 0.09 0.02 3: mp_mageq 0.1386E-01 0.51 0.12 3: mp_mageq_jpm1 0.4912E-01 1.81 0.44 3: mp_mageq_jpm3 0.9065E-02 0.33 0.08 3: mp_magpole_2d 0.3715E-01 1.37 0.33 3: mp_magpole_3d 0.6868E-04 0.00 0.00 3: mp_magpoles 0.2687E+00 9.91 2.39 3: mp_conjugate_points 0.3537E-06 0.00 0.00 3: mp_foldhem 0.2359E-01 0.87 0.21 3: mp_mag_periodic_f2d 0.2898E-01 1.07 0.26 3: mp_gather_pdyn 0.1406E-01 0.52 0.12 3: mp_mag_halos 0.6683E-01 2.46 0.59 3: mp_geo_halos 0.0000E+00 0.00 0.00 3: mp_geo_halos_f3d 0.1143E+00 4.21 1.02 3: mp_scatter_coeffs 0.0000E+00 0.00 0.00 3: mp_scatter_phim 0.1976E+00 7.29 1.76 3: mp_gather_f2d 0.2861E-03 0.01 0.00 3: mp_scatter_f2d 0.4638E-01 1.71 0.41 3: 3: ------------------------------------------------------------------------ 3: 3: ------------------------------------------------------------------------ 3: Report MPI wall-clock timing (mpi_wtime): mytid= 3 istep= 1440 nstep= 1440 3: Total Time (mins) in 'pefield ' = 0.3922E+00 ( 3.48% of total run time) 3: Total Time (mins) in 'dynamo_inputs ' = 0.2354E+01 ( 20.91% of total run time) 3: Total Time (mins) in 'pdynamo_fieldline_integrals ' = 0.1263E+00 ( 1.12% of total run time) 3: Total Time (mins) in 'pdynamo_complete_integrals ' = 0.1152E+00 ( 1.02% of total run time) 3: Total Time (mins) in 'pdynamo_rhspde ' = 0.6886E-03 ( 0.01% of total run time) 3: Total Time (mins) in 'highlat_poten ' = 0.2831E-03 ( 0.00% of total run time) 3: Total Time (mins) in 'pthreed ' = 0.1310E+01 ( 11.64% of total run time) 3: ------------------------------------------------------------------------ 3: 3: 3: ------------------------------------------------------------------------ 3: TIMER (system_clock): 3: Elapsed run time = 675.43 (secs) 11.26 (mins) 0.19 (hrs) 3: 3: Elapsed init time = 10.13 secs, 0.003 hrs, 1.5% (includes source i/o) 3: Elapsed step time = 643.19 secs, 0.179 hrs, 95.2% 3: Elapsed prep time = 700.03 secs, 0.194 hrs, 103.6% 3: Elapsed apxparm time= 0.90 secs, 0.000 hrs, 100.0% 3: Elapsed dynamics = 323.40 secs, 0.090 hrs, 47.9% 3: Elapsed qrj = 40.40 secs, 0.011 hrs, 6.0% 3: Elapsed oplus = 46.33 secs, 0.013 hrs, 6.9% 3: Elapsed cmpminor = 46.77 secs, 0.013 hrs, 6.9% (n4s,n2d,no) 3: Elapsed cmpmajor = 37.18 secs, 0.010 hrs, 5.5% (o2, o, he) 3: Elapsed dt = 10.63 secs, 0.003 hrs, 1.6% 3: Elapsed duv = 29.66 secs, 0.008 hrs, 4.4% 3: Elapsed pdynamo = 109.61 secs, 0.030 hrs, 16.2% 3: Elapsed magpres_grav= 10.16 secs, 0.003 hrs, 1.5% 3: Elapsed i/o = 1.43 secs, 0.000 hrs, 0.2% 3: Elapsed Primary io = 0.10 secs, 0.000 hrs, 0.0% 3: Elapsed Secondary = 1.33 secs, 0.000 hrs, 0.2% 3: 3: Average secs/step = 0.45 3: Average mins/simulated day = 11.26 3: ------------------------------------------------------------------------ 3: 3: End execution of tiegcm_trunk at 07/07/15 09:30:59 3: NORMAL EXIT