2: 2: ======================================================================== 2: Begin execution of tiegcm1.94 at 08/19/14 12:19:41 2: Host = puffin 2: System = LINUX 2: Logname = emery 2: ======================================================================== 2: 2: init_timer: level= 1 rtc=F sys=T 2: mp_init: ntask= 4 mytid= 2 2: Current working directory (cwd) = /hao/aim3/emery/mackjr/tiegcm-linux 2: Process ID (pid) = 2865 2: 2: Reading input data... 2: 2: Completed successful read of namelist inputs. 2: 2: mkntask: i= 1 j= 4 i*j= 4 ntask= 4 2: mkntask: i= 2 j= 2 i*j= 4 ntask= 4 2: mkntask: i= 4 j= 1 i*j= 4 ntask= 4 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: Expanded secout file s_dres_mackjr_DW1_06075_001.nc to 2: s_dres_mackjr_DW1_06075_001.nc 2: INPUT NOTE: adding mandatory field ZMAG to secondary history fields (field 48) 2: 2: ------------------------------------------------------------------------ 2: USER INPUT PARAMETERS: 2: label = tiegcm res=2.5 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/aim3/emery/mackjr/saber_lbcs_tngph_windiihrdi_unvn_with_ctmttides_DW1Only.nc 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/aim3/emery/mackjr/pso185mar_saber_winhrdi_eddyfix_AllTides_unvntngph_softxray_v3.nc 2: (file or mss path containing source history) 2: source_start = 185, 0, 0 (model time of source history) 2: output (primary history output files) = 2: p_dres_mackjr_DW1_06075_001.nc, 2: start (model start times) = 2: 75, 0, 0 2: stop (model stop times) = 2: 75, 0, 5 2: hist (primary history disk write frequencies) = 2: 0, 0, 5 2: Maxmimum number of histories per primary file = 25 2: secout (secondary history output files)= 2: s_dres_mackjr_DW1_06075_001.nc, 2: secstart (secondary history start times) = 2: 75, 0, 0 2: UI_ExB VI_ExB WI_ExB SIGMA_PED SIGMA_HAL 2: HMF2 NMF2 QJOULE QJOULE_VOL QJOULE_INTEG 2: ZG ALFA EFLUX DEN PHIHM 2: N2 O_N2 HE EX EY 2: EZ ED1 ED2 BMAG OPLUS 2: NOP_diag O2P_diag ZMAG 2: Maximum number of histories per secondary file = 25 2: Number of bytes for values of fields on secondary histories (sech_nbyte) = 4 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: 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.8000E+01 (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.1200E+03 (F10.7 solar flux) 2: f107a = 0.1200E+03 (81-day ave F10.7 flux) 2: al = 0.1000E+37 (AL, lower auroral mag index) 2: END USER INPUT PARAMETERS 2: ------------------------------------------------------------------------ 2: 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= 37 72 (36) lon0,1= 1 74 (74) ncells=2664 2: 2: Task 2: 2: tasks( 2)%mytid = 2 2: tasks( 2)%mytidi= 0 2: tasks( 2)%mytidj= 1 2: tasks( 2)%nlats = 36 2: tasks( 2)%nlons = 74 2: tasks( 2)%lat0 = 37 2: tasks( 2)%lat1 = 72 2: tasks( 2)%lon0 = 1 2: tasks( 2)%lon1 = 74 2: tasks( 2)%mylats= 2: 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 2: 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 2: 67 68 69 70 71 72 2: tasks( 2)%mylons= 2: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 2: 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 2: 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 2: 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 2: 61 62 63 64 65 66 67 68 69 70 71 72 73 74 2: Model version = tiegcm1.94 2: 2: Set constants: 2: nlat= 72 nlon=144 nlev= 56 2: dz= 0.25 2: dlat= 2.50 dlon= 2.50 2: zbound (cm) = 0.963723E+07 2: zmbot, zmtop = -6.875 7.125 (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.240 2: check_exp = F 2: kut (for filtering) = 1 1 2 2 4 4 8 8 10 10 12 12 15 15 18 18 22 22 26 26 30 30 32 32 34 34 34 34 34 34 34 34 34 34 34 34 2: kut (for filtering) = 34 34 34 34 34 34 34 34 34 34 34 34 32 32 30 30 26 26 22 22 18 18 15 15 12 12 10 10 8 8 4 4 2 2 1 1 2: init: iyear= 2006 iday= 75 2: gswm_mi_di_ncfile = /hao/aim3/emery/mackjr/saber_lbcs_tngph_windiihrdi_unvn_with_ctmttides_DW1Only.nc 2: hist_init: nstep= 5 2: allocdata: all tasks allocate fzg: nlevp1= 57 nlonp4= 148 nlat= 72 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=148 nlat= 72 2: init_lbc: allocated globals tlbc_nm_glb, ulbc_nm_glb, vlbc_nm_glb: nlonp4=148 nlat= 72 2: 2: Model run initialization: 2: nstep = 5 (Number of time steps this run) 2: iter = 108000 (Initial iteration number) 2: iyear = 2006 (Beginning calendar year) 2: iday = 75 (Beginning calendar day) 2: igswm_mi_di = 1 (If > 0, GSWM diurnal tidal database will be used.) 2: igswm_mi_sdi= 0 (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: nlev = 56 (Number of levels (midpoints and interfaces)) 2: 2: zpmid (midpoint levels) = 2: -6.875 -6.625 -6.375 -6.125 -5.875 -5.625 -5.375 -5.125 -4.875 -4.625 2: -4.375 -4.125 -3.875 -3.625 -3.375 -3.125 -2.875 -2.625 -2.375 -2.125 2: -1.875 -1.625 -1.375 -1.125 -0.875 -0.625 -0.375 -0.125 0.125 0.375 2: 0.625 0.875 1.125 1.375 1.625 1.875 2.125 2.375 2.625 2.875 2: 3.125 3.375 3.625 3.875 4.125 4.375 4.625 4.875 5.125 5.375 2: 5.625 5.875 6.125 6.375 6.625 6.875 7.125 2: 2: zpint (interface levels) = 2: -7.000 -6.750 -6.500 -6.250 -6.000 -5.750 -5.500 -5.250 -5.000 -4.750 2: -4.500 -4.250 -4.000 -3.750 -3.500 -3.250 -3.000 -2.750 -2.500 -2.250 2: -2.000 -1.750 -1.500 -1.250 -1.000 -0.750 -0.500 -0.250 0.000 0.250 2: 0.500 0.750 1.000 1.250 1.500 1.750 2.000 2.250 2.500 2.750 2: 3.000 3.250 3.500 3.750 4.000 4.250 4.500 4.750 5.000 5.250 2: 5.500 5.750 6.000 6.250 6.500 6.750 7.000 2: 2: nmlevp1= 60 zpmag_mid (mag midpoint levels) = 2: -7.625 -7.375 -7.125 -6.875 -6.625 -6.375 -6.125 -5.875 -5.625 -5.375 2: -5.125 -4.875 -4.625 -4.375 -4.125 -3.875 -3.625 -3.375 -3.125 -2.875 2: -2.625 -2.375 -2.125 -1.875 -1.625 -1.375 -1.125 -0.875 -0.625 -0.375 2: -0.125 0.125 0.375 0.625 0.875 1.125 1.375 1.625 1.875 2.125 2: 2.375 2.625 2.875 3.125 3.375 3.625 3.875 4.125 4.375 4.625 2: 4.875 5.125 5.375 5.625 5.875 6.125 6.375 6.625 6.875 7.125 2: 2: nmlevp1= 60 zpmag_int (mag interface levels) = 2: -7.750 -7.500 -7.250 -7.000 -6.750 -6.500 -6.250 -6.000 -5.750 -5.500 2: -5.250 -5.000 -4.750 -4.500 -4.250 -4.000 -3.750 -3.500 -3.250 -3.000 2: -2.750 -2.500 -2.250 -2.000 -1.750 -1.500 -1.250 -1.000 -0.750 -0.500 2: -0.250 0.000 0.250 0.500 0.750 1.000 1.250 1.500 1.750 2.000 2: 2.250 2.500 2.750 3.000 3.250 3.500 3.750 4.000 4.250 4.500 2: 4.750 5.000 5.250 5.500 5.750 6.000 6.250 6.500 6.750 7.000 2: 2: This is an initial run: 2: start_year = 2006 (Starting year of initial run) 2: start_day = 75 (Starting day of initial run) 2: start_mtime= 75 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 = 25 (Maximum number of primary histories per file) 2: 2: Secondary Histories: 2: nseries_sech = 1 (Number of secondary time series) 2: nsech_total = 2 (Number of secondary histories to be written) 2: nfiles_sech = 1 (Number of secondary output files to be written) 2: mxhist_sech = 25 (Maximum number of secondary histories per file) 2: nfsech = 48 (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: OMEGA 2: secondary history field 5: O1 2: secondary history field 6: O2 2: secondary history field 7: TE 2: secondary history field 8: TI 2: secondary history field 9: NO 2: secondary history field 10: NE 2: secondary history field 11: N2D 2: secondary history field 12: OP 2: secondary history field 13: N4S 2: secondary history field 14: Z 2: secondary history field 15: BARM 2: secondary history field 16: NO_COOL 2: secondary history field 17: CO2_COOL 2: secondary history field 18: POTEN 2: secondary history field 19: TEC 2: secondary history field 20: WN 2: secondary history field 21: UI_ExB 2: secondary history field 22: VI_ExB 2: secondary history field 23: WI_ExB 2: secondary history field 24: SIGMA_PED 2: secondary history field 25: SIGMA_HAL 2: secondary history field 26: HMF2 2: secondary history field 27: NMF2 2: secondary history field 28: QJOULE 2: secondary history field 29: QJOULE_VOL 2: secondary history field 30: QJOULE_INTEG 2: secondary history field 31: ZG 2: secondary history field 32: ALFA 2: secondary history field 33: EFLUX 2: secondary history field 34: DEN 2: secondary history field 35: PHIHM 2: secondary history field 36: N2 2: secondary history field 37: O_N2 2: secondary history field 38: HE 2: secondary history field 39: EX 2: secondary history field 40: EY 2: secondary history field 41: EZ 2: secondary history field 42: ED1 2: secondary history field 43: ED2 2: secondary history field 44: BMAG 2: secondary history field 45: OPLUS 2: secondary history field 46: NOP_diag 2: secondary history field 47: O2P_diag 2: secondary history field 48: 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 Height of the F2 Layer 2: 6 NMF2 1/cm3 none elden.F Peak Density of the F2 Layer 2: 7 FOF2 MHz none elden.F Critical Frequency of the F2 Layer 2: 8 JE13D A/m2 mlev current.F Eastward current density (3d) 2: 9 JE23D A/m2 mlev current.F Downward current density (3d) 2: 10 JQR A/m2 none current.F Upward current density (2d) 2: 11 KQLAM A/m none current.F Height-integrated current density (+north) 2: 12 KQPHI A/m none current.F Height-integrated current density (+east) 2: 13 LAMDA_HAL 1/s lev lamdas.F Hall Ion Drag Coefficient 2: 14 LAMDA_PED 1/s lev lamdas.F Pedersen Ion Drag Coefficient 2: 15 MU_M g/cm/s lev cpktkm.F Molecular Viscosity Coefficient 2: 16 QJOULE erg/g/s lev qjoule.F Joule Heating 2: 17 SCHT km lev addiag.F Pressure Scale Height 2: 18 SIGMA_HAL S/m lev lamdas.F Hall Conductivity 2: 19 SIGMA_PED S/m lev lamdas.F Pedersen Conductivity 2: 20 TEC 1/cm2 none elden.F Total Electron Content 2: 21 UI_ExB cm/s ilev ionvel.F Zonal Ion Drift (ExB) 2: 22 VI_ExB cm/s ilev ionvel.F Meridional Ion Drift (ExB) 2: 23 WI_ExB cm/s ilev ionvel.F Vertical Ion Drift (ExB) 2: 24 WN cm/s ilev swdot.F 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: ------------------------------------------------------------------------------------------------ 2: 2: 2: ------------------------------------------------------------------------ 2: Getfile: remote=/hao/aim3/emery/mackjr/pso185mar_saber_winhrdi_eddyfix_AllTides_unvntngph_softxray_v3.nc 2: Getfile: Found file /hao/aim3/emery/mackjr/pso185mar_saber_winhrdi_eddyfix_AllTides_unvntngph_softxray_v3.nc 2: ------------------------------------------------------------------------ 2: 2: Acquired source history file /hao/aim3/emery/mackjr/pso185mar_saber_winhrdi_eddyfix_AllTides_unvntngph_softxray_v3.nc 2: (disk file is /hao/aim3/emery/mackjr/pso185mar_saber_winhrdi_eddyfix_AllTides_unvntngph_softxray_v3.nc) 2: Reading source history from diskfile /hao/aim3/emery/mackjr/pso185mar_saber_winhrdi_eddyfix_AllTides_unvntngph_softxray_v3.nc: 2: nc_rdhist: seeking 185 0 0 found 184 0 0 n= 1 2: nc_rdhist: seeking 185 0 0 found 185 0 0 n= 2 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.1559E+03 0.1060E+04 2: Read field UN 3d subdomain min,max= -0.1168E+05 0.3544E+05 2: Read field VN 3d subdomain min,max= -0.3489E+05 0.3412E+05 2: Read field O2 3d subdomain min,max= 0.5839E-04 0.2304E+00 2: Read field O1 3d subdomain min,max= 0.3523E-02 0.9947E+00 2: Read field N4S 3d subdomain min,max= 0.1000E-11 0.2289E-01 2: Read field NO 3d subdomain min,max= 0.1430E-06 0.2658E-03 2: Read field OP 3d subdomain min,max= 0.0000E+00 0.2188E+07 2: Read field N2D 3d subdomain min,max= 0.1323E-14 0.2297E-03 2: Read field TI 3d subdomain min,max= 0.1560E+03 0.2017E+04 2: Read field TE 3d subdomain min,max= 0.1559E+03 0.3448E+04 2: Read field NE 3d subdomain min,max= 0.4903E+03 0.2201E+07 2: Read field O2P 3d subdomain min,max= 0.1215E-01 0.7665E+05 2: Read field OMEGA 3d subdomain min,max= -0.4497E-04 0.6910E-04 2: Read field Z 3d subdomain min,max= 0.9466E+07 0.6422E+08 2: Read field POTEN 3d subdomain min,max= -0.1663E+05 0.9459E+04 2: Read field TN_NM 3d subdomain min,max= 0.1559E+03 0.1060E+04 2: Read field UN_NM 3d subdomain min,max= -0.1168E+05 0.3544E+05 2: Read field VN_NM 3d subdomain min,max= -0.3487E+05 0.3413E+05 2: Read field O2_NM 3d subdomain min,max= 0.5839E-04 0.2304E+00 2: Read field O1_NM 3d subdomain min,max= 0.3523E-02 0.9947E+00 2: Read field N4S_NM 3d subdomain min,max= 0.1000E-11 0.2289E-01 2: Read field NO_NM 3d subdomain min,max= 0.1427E-06 0.2659E-03 2: Read field OP_NM 3d subdomain min,max= 0.0000E+00 0.2189E+07 2: Read LBC from source history: i= 89 sh%lbc= -7.00 2: 2: ------------------------------------------------------------------------ 2: Read TGCM PRIMARY HISTORY (source history) 2: Diskfile: /hao/aim3/emery/mackjr/pso185mar_saber_winhrdi_eddyfix_AllTides_unvntngph_softxray_v3.nc 2: label = tiegcm res=2.5 2: model_name = tiegcm 2: model_version = tiegcm1.94 2: create_date= 08/13/14 10:31:11 2: write_date = 2: logname = mcarthur 2: host = ys6326 2: system = LINUX 2: run_type = 2: source_file = 2: output_file = 2: source_mtime = 184 0 0 2: initial_year = 2006 2: initial_day = 75 2: initial_mtime= 75 0 0 2: type = primary 2: ihist = 2 (nth history on history file) 2: delhmins= 0 (delta minutes between histories) 2: calendar year,day = 2006, 74 2: (model IS being advanced in calendar time) 2: modeltime = 185, 0, 0, 0 (model time day,hour,minute,seconds) 2: iter = 266400 (number of steps from 0,0,0) 2: nlat = 72 (number of latitudes) 2: nlon = 144 (number of longitudes) 2: nlev = 57 (number of levels) 2: zmtop = 7.125 (top midpoint level) 2: zmbot = -6.875 (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.1200E+03 (daily solar flux) 2: f107a = 0.1200E+03 (average solar flux) 2: hpower = 0.8000E+01 (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.5445E+00 0.1880E+01 (ergs/cm2/s) 2: alfac,alfad = 0.1000E+00 0.5000E+00 (keV) 2: ec,ed = 0.5872E-01 0.6000E-02 (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 = 24 (number of model fields) 2: gswm_mi_di_ncfile = /glade/p/work/mcarthur/tiegcm_tides_zonalmean/saber_lbcs_tngph_windiihrdi_unvn_with_ctmttides.nc 2: gswm_mi_sdi_ncfile = [none] 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 = 64 (number of MPI tasks) 2: coupled_cmit = 0 (1 if coupled with CISM/CMIT, 0 otherwise) 2: There are 24 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 OP_NM 2: ------------------------------------------------------------------------ 2: 2: 2: Enter advance: iter=108000 nstep= 5 2: 2: ------------------------------------------------------------------------ 2: Getfile: remote=/hao/aim3/emery/mackjr/saber_lbcs_tngph_windiihrdi_unvn_with_ctmttides_DW1Only.nc 2: Getfile: Found file /hao/aim3/emery/mackjr/saber_lbcs_tngph_windiihrdi_unvn_with_ctmttides_DW1Only.nc 2: ------------------------------------------------------------------------ 2: 2: 2: ------------------------------------------------------------------------ 2: Reading GSWM data file /hao/aim3/emery/mackjr/saber_lbcs_tngph_windiihrdi_unvn_with_ctmttides_DW1Only.nc 2: Completed read from GSWM data file /hao/aim3/emery/mackjr/saber_lbcs_tngph_windiihrdi_unvn_with_ctmttides_DW1Only.nc 2: 2: ------------------------------------------------------------------------ 2: Step 1 of 5 mtime= 75 0 1 2: 2: aurora_cons: 2: cusp: alfac= 0.100 ec= 0.059 fc= 0.1833E+09 2: drizzle: alfad= 0.500 ed= 0.006 fd= 0.3745E+07 2: auroral radius = max of rhp,rcp= 18.372 17.905 2: roth, rote (MLT) = 0.549 -0.004 2: 1/e-widths = h1,h2= 2.264 3.405 2: energy flux = e1,e2= 0.545 1.880 2: add_sproton = F 2: 2: Step 2 of 5 mtime= 75 0 2 0 secs/step (sys) = 13.33 2: Step 3 of 5 mtime= 75 0 3 0 secs/step (sys) = 1.52 2: Step 4 of 5 mtime= 75 0 4 0 secs/step (sys) = 1.40 2: Step 5 of 5 mtime= 75 0 5 0 secs/step (sys) = 1.41 2: 2: Allocated 3d sech field ZG(lon=148,lat= 72,ilev= 57) 2: 2: Initialized diagnostic secondary history field ZG (ix= 31): 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 = 148 72 57 2: ndims = 3 2: task0_only = F 2: 2: Allocated 3d sech field WN(lon=148,lat= 72,lev= 57) 2: 2: Initialized diagnostic secondary history field WN (ix= 20): 2: short_name = WN 2: long_name = NEUTRAL VERTICAL WIND (plus up) 2: units = cm/s 2: geo = T 2: mag = F 2: dimnames = lon lat lev 2: dimsizes = 148 72 57 2: ndims = 3 2: task0_only = F 2: 2: Allocated 3d sech field UI_ExB(lon=148,lat= 72,ilev= 57) 2: 2: Initialized diagnostic secondary history field UI_ExB (ix= 21): 2: short_name = UI_ExB 2: long_name = Zonal Ion Drift (ExB) 2: units = cm/s 2: geo = T 2: mag = F 2: dimnames = lon lat ilev 2: dimsizes = 148 72 57 2: ndims = 3 2: task0_only = F 2: 2: Allocated 3d sech field VI_ExB(lon=148,lat= 72,ilev= 57) 2: 2: Initialized diagnostic secondary history field VI_ExB (ix= 22): 2: short_name = VI_ExB 2: long_name = Meridional Ion Drift (ExB) 2: units = cm/s 2: geo = T 2: mag = F 2: dimnames = lon lat ilev 2: dimsizes = 148 72 57 2: ndims = 3 2: task0_only = F 2: 2: Allocated 3d sech field WI_ExB(lon=148,lat= 72,ilev= 57) 2: 2: Initialized diagnostic secondary history field WI_ExB (ix= 23): 2: short_name = WI_ExB 2: long_name = Vertical Ion Drift (ExB) 2: units = cm/s 2: geo = T 2: mag = F 2: dimnames = lon lat ilev 2: dimsizes = 148 72 57 2: ndims = 3 2: task0_only = F 2: 2: Allocated 2d sech field TEC(lon=148,lat= 72) 2: 2: Initialized diagnostic secondary history field TEC (ix= 19): 2: short_name = TEC 2: long_name = Total Electron Content 2: units = 1/cm2 2: geo = T 2: mag = F 2: dimnames = lon lat 2: dimsizes = 148 72 0 2: ndims = 2 2: task0_only = F 2: 2: Allocated 2d sech field HMF2(lon=148,lat= 72) 2: 2: Initialized diagnostic secondary history field HMF2 (ix= 26): 2: short_name = HMF2 2: long_name = Height of the F2 Layer 2: units = km 2: geo = T 2: mag = F 2: dimnames = lon lat 2: dimsizes = 148 72 0 2: ndims = 2 2: task0_only = F 2: 2: Allocated 2d sech field NMF2(lon=148,lat= 72) 2: 2: Initialized diagnostic secondary history field NMF2 (ix= 27): 2: short_name = NMF2 2: long_name = Peak Density of the F2 Layer 2: units = 1/cm3 2: geo = T 2: mag = F 2: dimnames = lon lat 2: dimsizes = 148 72 0 2: ndims = 2 2: task0_only = F 2: 2: Allocated 3d sech field SIGMA_PED(lon=148,lat= 72,lev= 57) 2: 2: Initialized diagnostic secondary history field SIGMA_PED (ix= 24): 2: short_name = SIGMA_PED 2: long_name = Pedersen Conductivity 2: units = S/m 2: geo = T 2: mag = F 2: dimnames = lon lat lev 2: dimsizes = 148 72 57 2: ndims = 3 2: task0_only = F 2: 2: Allocated 3d sech field SIGMA_HAL(lon=148,lat= 72,lev= 57) 2: 2: Initialized diagnostic secondary history field SIGMA_HAL (ix= 25): 2: short_name = SIGMA_HAL 2: long_name = Hall Conductivity 2: units = S/m 2: geo = T 2: mag = F 2: dimnames = lon lat lev 2: dimsizes = 148 72 57 2: ndims = 3 2: task0_only = F 2: 2: Allocated 3d sech field CO2_COOL(lon=148,lat= 72,lev= 57) 2: 2: Initialized diagnostic secondary history field CO2_COOL (ix= 17): 2: short_name = CO2_COOL 2: long_name = CO2 Cooling 2: units = erg/g/s 2: geo = T 2: mag = F 2: dimnames = lon lat lev 2: dimsizes = 148 72 57 2: ndims = 3 2: task0_only = F 2: 2: Allocated 3d sech field NO_COOL(lon=148,lat= 72,lev= 57) 2: 2: Initialized diagnostic secondary history field NO_COOL (ix= 16): 2: short_name = NO_COOL 2: long_name = NO Cooling 2: units = erg/g/s 2: geo = T 2: mag = F 2: dimnames = lon lat lev 2: dimsizes = 148 72 57 2: ndims = 3 2: task0_only = F 2: 2: Allocated 3d sech field QJOULE(lon=148,lat= 72,lev= 57) 2: 2: Initialized diagnostic secondary history field QJOULE (ix= 28): 2: short_name = QJOULE 2: long_name = Joule Heating 2: units = erg/g/s 2: geo = T 2: mag = F 2: dimnames = lon lat lev 2: dimsizes = 148 72 57 2: ndims = 3 2: task0_only = F 2: 2: Allocated 2d sech field QJOULE_INTEG(lon=148,lat= 72) 2: 2: Initialized diagnostic secondary history field QJOULE_INTEG (ix= 30): 2: short_name = QJOULE_INTEG 2: long_name = Height-integrated Joule Heating 2: units = erg/cm2/s 2: geo = T 2: mag = F 2: dimnames = lon lat 2: dimsizes = 148 72 0 2: ndims = 2 2: task0_only = F 2: 2: Allocated 3d sech field DEN(lon=148,lat= 72,ilev= 57) 2: 2: Initialized diagnostic secondary history field DEN (ix= 34): 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 = 148 72 57 2: ndims = 3 2: task0_only = F 2: 2: Allocated 3d sech field O_N2(lon=148,lat= 72,lev= 57) 2: 2: Initialized diagnostic secondary history field O_N2 (ix= 37): 2: short_name = O_N2 2: long_name = O/N2 RATIO 2: units = none 2: geo = T 2: mag = F 2: dimnames = lon lat lev 2: dimsizes = 148 72 57 2: ndims = 3 2: task0_only = F 2: 2: MPI run with ntask = 4 2: nstep= 5 step= 60 2: Model simulation time = 300 secs 2: (minutes= 5.00, hours= 0.08, days= 0.003472) 2: Cpu time for run = 23.51 2: 2: ------------------------------------------------------------------------ 2: TIMER (system_clock): 2: Elapsed run time = 24.10 (secs) 0.40 (mins) 0.01 (hrs) 2: 2: Elapsed init time = 3.44 secs, 0.001 hrs, 14.3% (includes source i/o) 2: Elapsed step time = 19.55 secs, 0.005 hrs, 81.1% 2: Elapsed prep time = 30.79 secs, 0.009 hrs, 127.8% 2: Elapsed dynamics = 6.66 secs, 0.002 hrs, 27.6% 2: Elapsed qrj = 0.69 secs, 0.000 hrs, 2.9% 2: Elapsed oplus = 0.99 secs, 0.000 hrs, 4.1% 2: Elapsed cmpminor = 0.69 secs, 0.000 hrs, 2.9% (n4s,n2d,no) 2: Elapsed cmpmajor = 0.39 secs, 0.000 hrs, 1.6% (o2, o) 2: Elapsed dt = 0.26 secs, 0.000 hrs, 1.1% 2: Elapsed duv = 0.49 secs, 0.000 hrs, 2.1% 2: Elapsed dynamo = 0.95 secs, 0.000 hrs, 3.9% 2: Elapsed i/o = 3.24 secs, 0.001 hrs, 13.4% 2: Elapsed Primary io = 2.13 secs, 0.001 hrs, 8.9% 2: Elapsed Secondary = 1.11 secs, 0.000 hrs, 4.6% 2: 2: Average secs/step = 3.91 2: Average mins/simulated day = 115.66 2: ------------------------------------------------------------------------ 2: 2: End execution of tiegcm1.94 at 08/19/14 12:20:18 2: NORMAL EXIT