======================================================================== Begin execution of tiegcm1.94.2 at 08/31/12 16:12:24 Host = puffin System = LINUX Logname = emery ======================================================================== init_timer: level= 1 rtc=F sys=T Current working directory (cwd) = /hao/aim3/emery/tiegcm_cmit/calccloc-2_1_5b/CMITs Process ID (pid) = 11420 &tgcm_input LABEL = 'AGU_CMITs' START_YEAR = 2006 START_DAY = 348 CALENDAR_ADVANCE = 1 STEP = 60 START = 348,22,0 STOP = 348,22,5 HIST = 0,0,5 OUTPUT = 'AGU_CMITs_tie_2006-12-14T22-00-00.nc','AGU_CMITs_old_2006-12-14T22-00-05_pold.nc' MXHIST_PRIM = 1 SECSTART = 348,22,5 SECSTOP = 348,22,5 SECHIST = 0,0,5 SECOUT = 'AGU_CMITs_sech_tie_2006-12-14T22-00-00_testb_old.nc' MXHIST_SECH = 12 SECFLDS = 'TN','UN','VN','O2','O1','NO','NO_COOL','OP','TI','TE','NE','O2P','OMEGA','Z','POTEN','UI_ExB','VI_ExB','WI_ExB','TEC','QJI_TN','gpot','gpotm','geng','gflx','gzigm1','gzigm2','gnsrhs','CUSP','DRIZZLE','ALFA','NFLUX','EFLUX' AURORA = 1 COLFAC = 1.5 POTENTIAL_MODEL = 'HEELIS' GPI_NCFILE = '$TGCMDATA/gpi_2000001-2011059.nc' TIDE = 0.,0.,0.,0.,0.,0.,0.,0.,0.,0. TIDE2 = 0.,0. GSWM_MI_DI_NCFILE = '$TGCMDATA/gswm_diurn_5.0d_99km.nc' GSWM_MI_SDI_NCFILE = '$TGCMDATA/gswm_semi_5.0d_99km.nc' / Reading input data... Completed successful read of namelist inputs. Will use the Heelis potential model Note input: Setting BY to 0 with HEELIS potential model. Input: Expanded secout file AGU_CMITs_sech_tie_2006-12-14T22-00-00_testb_old.nc to AGU_CMITs_sech_tie_2006-12-14T22-00-00_testb_old.nc INPUT NOTE: adding mandatory field ZG to secondary history fields (field 33) INPUT NOTE: adding mandatory field ZMAG to secondary history fields (field 34) ------------------------------------------------------------------------ USER INPUT PARAMETERS: label = AGU_CMITs (optional text label for current run) High-lat electric potential model: potential_model = HEELIS gpi run: gpi_ncfile = $TGCMDATA/gpi_2000001-2011059.nc gswm migrating diurnal file: gswm_mi_di_ncfile = /hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc gswm migrating semi-diurnal file: gswm_mi_sdi_ncfile = /hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc start_year = 2006 (starting calendar day) start_day = 348 (starting calendar year) calendar_advance = 1 (model will be advanced in calendar time starting on this day) step = 60 (model timestep (seconds)) output (primary history output files) = AGU_CMITs_tie_2006-12-14T22-00-00.nc, AGU_CMITs_old_2006-12-14T22-00-05_pold.nc start (model start times) = 348,22, 0 stop (model stop times) = 348,22, 5 hist (primary history disk write frequencies) = 0, 0, 5 Maxmimum number of histories per primary file = 1 secout (secondary history output files)= AGU_CMITs_sech_tie_2006-12-14T22-00-00_testb_old.nc, secstart (secondary history start times) = 348,22, 5 secstop (secondary history stop times) = 348,22, 5 sechist (secondary history disk write frequencies) = 0, 0, 5 secflds (secondary history fields) = TN UN VN O2 O1 NO NO_COOL OP TI TE NE O2P OMEGA Z POTEN UI_ExB VI_ExB WI_ExB TEC QJI_TN gpot gpotm geng gflx gzigm1 gzigm2 gnsrhs CUSP DRIZZLE ALFA NFLUX EFLUX ZG ZMAG Maximum number of histories per secondary file = 12 Number of bytes for values of fields on secondary histories (sech_nbyte) = 4 eddy_dif = 0 (DOY-dependent eddy diffusion flag) tide (amplitudes and phases of semidiurnal tide) = 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 tide2 (amplitude and phase of diurnal tide) = 0.0E+00 0.00 aurora = 1 (0/1 flag for aurora) colfac = 1.500 (collision factor) If any of the following are spval ( 0.1000E+37), they willprogram tiegcm(rank 0) is IC_Waiting remote program ======= IC_Initialize was successful. ======= be calculated during the simulation on a per timestep basis: power = 0.1000E+37 (Hemispheric Power) ctpoten= 0.1000E+37 (Cross-cap potential) kp = 0.1000E+37 (Kp index) bximf = 0.1000E+37 (BX component of IMF) byimf = 0.0000E+00 (BY component of IMF) bzimf = 0.1000E+37 (Bz component of IMF) swvel = 0.1000E+37 (solar wind velocity) swden = 0.1000E+37 (solar wind density) f107 = 0.1000E+37 (F10.7 solar flux) f107a = 0.1000E+37 (81-day ave F10.7 flux) al = 0.1000E+37 (AL, lower auroral mag index) END USER INPUT PARAMETERS ------------------------------------------------------------------------ Initializing TIEGCM InterComm Interface Model version = tiegcm1.94.2 Set constants: nlat= 36 nlon= 72 nlev= 28 dz= 0.50 dlat= 5.00 dlon= 5.00 zbound (cm) = 0.963723E+07 zmbot, zmtop = -6.750 7.250 (bottom,top midpoint levels) zibot, zitop = -7.000 7.000 (bottom,top interface levels) dt = 60.00 secs grav = 870.00 freq_3m3 = 0.3506E-04 freq_semidi= 0.1454E-03 dipmin = 0.170 check_exp = F 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 init: iyear= 2006 iday= 348 gswm_mi_di_ncfile = /hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc gswm_mi_sdi_ncfile = /hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc hist_init: nstep= 5 Allocated 3d data for all fields on root task: nlevp1*nlonp4*nlat*nf4d = 29* 76* 36* 26 = 2062944 words *8 = 16503552 bytes. init_lbc: allocated subdomains tlbc, ulbc, vlbc init_lbc_nm: allocated subdomains tlbc_nm, ulbc_nm, vlbc_nm init_lbc: allocated globals tlbc_glb, ulbc_glb, vlbc_glb: nlonp4= 76 nlat= 36 init_lbc: allocated globals tlbc_nm_glb, ulbc_nm_glb, vlbc_nm_glb: nlonp4= 76 nlat= 36 Model run initialization: nstep = 5 (Number of time steps this run) iter = 502440 (Initial iteration number) iyear = 2006 (Beginning calendar year) iday = 348 (Beginning calendar day) igswm_mi_di = 1 (If > 0, GSWM diurnal tidal database will be used.) igswm_mi_sdi= 1 (If > 0, GSWM semidiurnal tidal database will be used.) igswm_nm_di= 0 (If > 0, GSWM nonmigrating diurnal tidal database will be used.) igswm_nm_sdi= 0 (If > 0, GSWM nonmigrating semidiurnal tidal database will be used.) nlev = 28 (Number of levels (midpoints and interfaces)) zpmid (midpoint levels) = -6.750 -6.250 -5.750 -5.250 -4.750 -4.250 -3.750 -3.250 -2.750 -2.250 -1.750 -1.250 -0.750 -0.250 0.250 0.750 1.250 1.750 2.250 2.750 3.250 3.750 4.250 4.750 5.250 5.750 6.250 6.750 7.250 zpint (interface levels) = -7.000 -6.500 -6.000 -5.500 -5.000 -4.500 -4.000 -3.500 -3.000 -2.500 -2.000 -1.500 -1.000 -0.500 0.000 0.500 1.000 1.500 2.000 2.500 3.000 3.500 4.000 4.500 5.000 5.500 6.000 6.500 7.000 nmlevp1= 32 zpmag_mid (mag midpoint levels) = -8.250 -7.750 -7.250 -6.750 -6.250 -5.750 -5.250 -4.750 -4.250 -3.750 -3.250 -2.750 -2.250 -1.750 -1.250 -0.750 -0.250 0.250 0.750 1.250 1.750 2.250 2.750 3.250 3.750 4.250 4.750 5.250 5.750 6.250 6.750 7.250 nmlevp1= 32 zpmag_int (mag interface levels) = -8.500 -8.000 -7.500 -7.000 -6.500 -6.000 -5.500 -5.000 -4.500 -4.000 -3.500 -3.000 -2.500 -2.000 -1.500 -1.000 -0.500 0.000 0.500 1.000 1.500 2.000 2.500 3.000 3.500 4.000 4.500 5.000 5.500 6.000 6.500 7.000 Primary Histories: nsource = 0 (If > 0, a primary source history was provided) nseries_prim = 1 (Number of primary time series) nhist_total = 1 (Number of primary histories to be written) nfiles_prim = 1 (Number of primary output files to be written) mxhist_prim = 1 (Maximum number of primary histories per file) Secondary Histories: nseries_sech = 1 (Number of secondary time series) nsech_total = 1 (Number of secondary histories to be written) nfiles_sech = 1 (Number of secondary output files to be written) mxhist_sech = 12 (Maximum number of secondary histories per file) nfsech = 34 (Number of requested secondary history fields) secondary history field 1: TN secondary history field 2: UN secondary history field 3: VN secondary history field 4: O2 secondary history field 5: O1 secondary history field 6: NO secondary history field 7: NO_COOL secondary history field 8: OP secondary history field 9: TI secondary history field 10: TE secondary history field 11: NE secondary history field 12: O2P secondary history field 13: OMEGA secondary history field 14: Z secondary history field 15: POTEN secondary history field 16: UI_ExB secondary history field 17: VI_ExB secondary history field 18: WI_ExB secondary history field 19: TEC secondary history field 20: QJI_TN secondary history field 21: gpot secondary history field 22: gpotm secondary history field 23: geng secondary history field 24: gflx secondary history field 25: gzigm1 secondary history field 26: gzigm2 secondary history field 27: gnsrhs secondary history field 28: CUSP secondary history field 29: DRIZZLE secondary history field 30: ALFA secondary history field 31: NFLUX secondary history field 32: EFLUX secondary history field 33: ZG secondary history field 34: ZMAG ------------------------------------------------------------------------------------------------ Table of Available Diagnostic Fields: Shortnames may be added to namelist SECFLDS Field Shortname Units Levels Caller Longname - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 1 CO2_COOL erg/g/s lev newton.F CO2 Cooling 2 NO_COOL erg/g/s lev newton.F NO Cooling 3 DEN g/cm3 ilev dt.F Total Density 4 HEATING erg/g/s lev dt.F Total Heating 5 HMF2 km none elden.F Height of the F2 Layer 6 NMF2 1/cm3 none elden.F Peak Density of the F2 Layer 7 JE13D A/m2 mlev current.F Eastward current density (3d) 8 JE23D A/m2 mlev current.F Downward current density (3d) 9 JQR A/m2 none current.F Upward current density (2d) 10 KQLAM A/m none current.F Height-integrated current density (+north) 11 KQPHI A/m none current.F Height-integrated current density (+east) 12 LAMDA_HAL 1/s lev lamdas.F Hall Ion Drag Coefficient 13 LAMDA_PED 1/s lev lamdas.F Pedersen Ion Drag Coefficient 14 MU_M g/cm/s lev cpktkm.F Molecular Viscosity Coefficient 15 QJOULE erg/g/s lev qjoule.F Joule Heating 16 SCHT km lev addiag.F Pressure Scale Height 17 SIGMA_HAL S/m lev lamdas.F Hall Conductivity 18 SIGMA_PED S/m lev lamdas.F Pedersen Conductivity 19 TEC 1/cm2 none elden.F Total Electron Content 20 UI_ExB cm/s ilev ionvel.F Zonal Ion Drift (ExB) 21 VI_ExB cm/s ilev ionvel.F Meridional Ion Drift (ExB) 22 WI_ExB cm/s ilev ionvel.F Vertical Ion Drift (ExB) 23 WN cm/s ilev swdot.F NEUTRAL VERTICAL WIND (plus up) 24 O_N2 none lev comp.F O/N2 RATIO 25 QJOULE_INTEG erg/cm2/s none qjoule.F Height-integrated Joule Heating ------------------------------------------------------------------------------------------------ ------------------------------------------------------------------------ Getfile: remote=AGU_CMITs_tie_2006-12-14T22-00-00.nc Getfile: Found file AGU_CMITs_tie_2006-12-14T22-00-00.nc ------------------------------------------------------------------------ Acquired source history file AGU_CMITs_tie_2006-12-14T22-00-00.nc (disk file is AGU_CMITs_tie_2006-12-14T22-00-00.nc) Reading source history from diskfile AGU_CMITs_tie_2006-12-14T22-00-00.nc: nc_rdhist: seeking 348 22 0 found 348 22 0 n= 1 Note nc_rdhist: unused variable: calendar_advance Note nc_rdhist: unused variable: write_date Note nc_rdhist: unused variable: crit1 Note nc_rdhist: unused variable: crit2 Note nc_rdhist: unused variable: mag Read field TN 3d subdomain min,max= 0.1548E+03 0.1085E+04 Read field UN 3d subdomain min,max= -0.4869E+05 0.5835E+05 Read field VN 3d subdomain min,max= -0.3627E+05 0.5289E+05 Read field O2 3d subdomain min,max= 0.2359E-04 0.2302E+00 Read field O1 3d subdomain min,max= 0.3861E-02 0.9976E+00 Read field N4S 3d subdomain min,max= 0.1000E-11 0.1755E-01 Read field NO 3d subdomain min,max= 0.1173E-06 0.1141E-02 Read field OP 3d subdomain min,max= 0.0000E+00 0.1185E+07 Read field N2D 3d subdomain min,max= 0.0000E+00 0.3838E-02 Read field TI 3d subdomain min,max= 0.1548E+03 0.3912E+04 Read field TE 3d subdomain min,max= 0.1548E+03 0.3983E+04 Read field NE 3d subdomain min,max= 0.9784E+03 0.1166E+07 Read field O2P 3d subdomain min,max= 0.0000E+00 0.1752E+06 Read field OMEGA 3d subdomain min,max= -0.3766E-03 0.9896E-04 Read field Z 3d subdomain min,max= 0.9578E+07 0.6064E+08 Read field POTEN 3d subdomain min,max= -0.7516E+05 0.1009E+06 Read field TN_NM 3d subdomain min,max= 0.1548E+03 0.1085E+04 Read field UN_NM 3d subdomain min,max= -0.4869E+05 0.5874E+05 Read field VN_NM 3d subdomain min,max= -0.3649E+05 0.5307E+05 Read field O2_NM 3d subdomain min,max= 0.2357E-04 0.2302E+00 Read field O1_NM 3d subdomain min,max= 0.3852E-02 0.9975E+00 Read field N4S_NM 3d subdomain min,max= 0.1000E-11 0.1755E-01 Read field NO_NM 3d subdomain min,max= 0.1169E-06 0.1141E-02 Read field OP_NM 3d subdomain min,max= 0.0000E+00 0.1188E+07 Read LBC from source history: i= 89 sh%lbc= -7.00 ------------------------------------------------------------------------ Read TGCM PRIMARY HISTORY (source history) Diskfile: AGU_CMITs_tie_2006-12-14T22-00-00.nc label = AGU_CMITs model_name = tiegcm model_version = tiegcm1.94.2 create_date= 05/10/12 11:28:38 write_date = logname = ps14 host = unknown system = LINUX run_type = source_file = output_file = source_mtime = 348 20 0 initial_year = 2006 initial_day = 328 initial_mtime= 328 0 0 type = primary ihist = 1 (nth history on history file) delhmins= 0 (delta minutes between histories) calendar year,day = 2006,348 (model IS being advanced in calendar time) modeltime = 348,22, 0, 0 (model time day,hour,minute,seconds) time = ut = 22.00 (ut hours) step = 60 (time step in seconds) iter = 502440 (number of steps from 0,0,0) nlat = 36 (number of latitudes) nlon = 72 (number of longitudes) nlev = 29 (number of levels) zmtop = 7.250 (top midpoint level) zmbot = -6.750 (bottom midpoint level) zitop = 7.000 (top interface level) zibot = -7.000 (bottom interface level) dtcpl-itm_potential recv with tag 2500 succeed cpl-itm_energy recv with tag 2501 succeed cpl-itm_flux recv with tag 2502 succeed ide = 0.0E+00 0.0 (amp/phase of diurnal tide) 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 (amp/phase of semi-diurnal tide) f107d = 0.8826E+02 (daily solar flux) f107a = 0.8234E+02 (average solar flux) hpower = 0.1752E+03 (Gw) ctpoten = 0.1671E+03 (Volts) kp = 0.1000E+37 () byimf = 0.0000E+00 bzimf = 0.1000E+37 swvel = 0.1000E+37 swden = 0.1000E+37 al = 0.1000E+37 e1,e2 = 0.4544E+01 0.2027E+02 (ergs/cm2/s) alfac,alfad = 0.1000E+00 0.5000E+00 (keV) ec,ed = 0.2827E+00 0.1063E+00 (ergs/cm2/s) colfac = 0.15E+01 joulefac = 0.10E+01 p0 = 0.50E-06 p0_model= 0.50E-03 grav = 0.87E+03 nflds = 24 (number of model fields) gswm_mi_di_ncfile = /users/ps14/data/tiegcm1.94_data/gswm_diurn_5.0d_99km.nc gswm_mi_sdi_ncfile = /users/ps14/data/tiegcm1.94_data/gswm_semi_5.0d_99km.nc gswm_nm_di_ncfile = [none] gswm_nm_sdi_ncfile = [none] see_ncfile = [none] gpi_ncfile = /users/ps14/data/tiegcm1.94_data/gpi_2000001-2011059.nc ncep_ncfile = [none] imf_ncfile = [none] saber_ncfile = [none] tidi_ncfile = [none] tuv_lbc_intop= 0 (if 1, then lbc of t,u,v are stored in top k slot (old histories)) LBC = -7.00 (lower boundary interface level) ntask_mpi = 0 (number of MPI tasks) coupled_cmit = 1 (1 if coupled with CISM/CMIT, 0 otherwise) There are 24 fields on this history, as follows: TN UN VN O2 O1 N4S NO OP N2D TI TE NE O2P OMEGA Z POTEN TN_NM UN_NM VN_NM O2_NM O1_NM N4S_NM NO_NM OP_NM ------------------------------------------------------------------------ Enter advance: iter=502440 nstep= 5 ------------------------------------------------------------------------ Getfile: remote=$TGCMDATA/gpi_2000001-2011059.nc Getfile: Found file /hao/aim1/tgcm/data/gpi_2000001-2011059.nc ------------------------------------------------------------------------ ------------------------------------------------------------------------ RDGPI: read GPI data file: Number of days on GPI data file = 4077 First and last year-days = 2000003 2011057 Completed read from GPI data file /hao/aim1/tgcm/data/gpi_2000001-2011059.nc ------------------------------------------------------------------------ ------------------------------------------------------------------------ GETGPI: get geophysical indices from data file /hao/aim1/tgcm/data/gpi_2000001-2011059.nc Initial requested iyear=2006 iday=348 iutsec= 79260 Obtained the following GPI at the requested date and time: f107 = 0.8825E+02 f107a = 0.8234E+02 power = 0.1761E+03 ctpoten = 0.1675E+03 ------------------------------------------------------------------------ ------------------------------------------------------------------------ Getfile: remote=/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc Getfile: Found file /hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ------------------------------------------------------------------------ ------------------------------------------------------------------------ Reading GSWM data file /hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc Completed read from GSWM data file /hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ------------------------------------------------------------------------ ------------------------------------------------------------------------ Getfile: remote=/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc Getfile: Found file /hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ------------------------------------------------------------------------ ------------------------------------------------------------------------ Reading GSWM data file /hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc Completed read from GSWM data file /hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ------------------------------------------------------------------------ inside receive cism_input_analysis byimf,loc,power,ctpoten,kp= itm-cpl_pedersen send with tag 2550 succeed itm-cpl_hall send with tag 2551 succeed itm-cpl_wfac send with tag 2552 succeed cpl-itm_potential recv with tag 2500 succeed cpl-itm_energy recv with tag 2501 succeed cpl-itm_flux recv with tag 2502 succeed itm-cpl_pedersen send with tag 2550 succeed itm-cpl_hall send with tag 2551 succeed itm-cpl_wfac send with tag 2552 succeed cpl-itm_potential recv with tag 2500 succeed cpl-itm_energy recv with tag 2501 succeed cpl-itm_flux recv with tag 2502 succeed itm-cpl_pedersen send with tag 2550 succeed itm-cpl_hall send with tag 2551 succeed itm-cpl_wfac send with tag 2552 succeed cpl-itm_potential recv with tag 2500 succeed cpl-itm_energy recv with tag 2501 succeed cpl-itm_flux recv with tag 2502 succeed 0.00 0.00 0.1761E+03 0.1675E+03 0.7313E+01 MIX HP power (S,N,b), maxS,N mW/m2= 0.2071E+02 0.1638E+02 0.1855E+02 0.3718E+01 0.3093E+01 MIX gpot model_ctpoten (S,N,av) = 0.1628E+03 0.1992E+03 0.1810E+03 MIX ifbad model_ctpoten (S,N,av) = 0 0 0.1618E+03 0.1981E+03 0.1800E+03 CalcCLoc S,N: offc,dskofc,theta0,phid,n= -0.1357E+01 -0.7845E+00 0.8528E+01 0.2049E+01 0.1542E+02 0.1238E+02 -0.1261E+02 -0.1296E+02 -0.9640E+00 -0.1864E+01 After Revised S,N: offc,dskofc,theta0,phid,n= 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.1542E+02 0.1238E+02 -0.1261E+02 -0.1296E+02 -0.9640E+00 -0.1864E+01 crit before colath = 0.1500E+02 0.3000E+02 crit after colath = 0.1890E+02 0.3390E+02 Step 1 of 5 mtime=348 22 1 GPI run: istep= 1 power= 18.546 ctpoten= 179.967 f107= 88.252 f107a= 82.344 kp= 7.313 aurora_cons: cusp: alfac= 0.100 ec= 0.073 fc= 0.2274E+09 drizzle: alfad= 0.500 ed= 0.012 fd= 0.7695E+07 auroral radius = max of rhp,rcp= 20.059 25.225 roth, rote (MLT) = 0.444 -0.074 1/e-widths = h1,h2= 2.350 3.510 energy flux = e1,e2= 1.634 3.040 add_sproton = F theta0,offa,dskofa,phid,phin,rrad,offc,dskofc= 0.1542E+02 0.1238E+02 0.0000E+00 0.0000E+00 -0.2500E+01 -0.2500E+01 -0.6882E+03 -0.6885E+03 -0.6765E+03 -0.6774E+03 0.1742E+02 0.1438E+02 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 Sending at: 348 22 0 0 inside receive cism_input_analysis byimf,loc,power,ctpoten,kp= 0.00 0.00 0.1771E+03 0.1678E+03 0.7327E+01 MIX HP power (S,N,b), maxS,N mW/m2= 0.2271E+02 0.1606E+02 0.1939E+02 0.4076E+01 0.2990E+01 MIX gpot model_ctpoten (S,N,av) = 0.1812E+03 0.2093E+03 0.1953E+03 MIX ifbad model_ctpoten (S,N,av) = 0 0 0.1808E+03 0.2067E+03 0.1937E+03 CalcCLoc S,N: offc,dskofc,theta0,phid,n= -0.1357E+01 -0.7845E+00 0.7561E+01 0.2049E+01 0.1534E+02 0.1247E+02 -0.1261E+02 -0.1325E+02 -0.6479E+00 -0.1848E+01 After Revised S,N: offc,dskofc,theta0,phid,n= 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.1534E+02 0.1247E+02 -0.1261E+02 -0.1325E+02 -0.6479E+00 -0.1848E+01 crit before colath = 0.1890E+02 0.3390E+02 crit after colath = 0.1890E+02 0.3390E+02 Step 2 of 5 mtime=348 22 2 0 secs/step (sys) = 3.15 GPI run: istep= 2 power= 19.385 ctpoten= 193.727 f107= 88.248 f107a= 82.344 kp= 7.327 theta0,offa,dskofa,phid,phin,rrad,offc,dskofc= 0.1534E+02 0.1247E+02 0.0000E+00 0.0000E+00 -0.2500E+01 -0.2500E+01 -0.6882E+03 -0.6888E+03 -0.6762E+03 -0.6774E+03 0.1734E+02 0.1447E+02 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 Sending at: 348 22 1 0 inside receive cism_input_analysis byimf,loc,power,ctpoten,kp= 0.00 0.00 0.1781E+03 0.1682E+03 0.7340E+01 MIX HP power (S,N,b), maxS,N mW/m2= 0.2368E+02 0.1747E+02 0.2058E+02 0.4363E+01 0.3170E+01 MIX gpot model_ctpoten (S,N,av) = 0.1950E+03 0.2283E+03 0.2116E+03 MIX ifbad model_ctpoten (S,N,av) = 0 0 0.1945E+03 0.2260E+03 0.2102E+03 CalcCLoc S,N: offc,dskofc,theta0,phid,n= 0.3653E+00 -0.1482E+01 0.7561E+01 0.2049E+01 0.1544E+02 0.1240E+02 -0.1261E+02 -0.1293E+02 -0.6317E+00 -0.1832E+01 After Revised S,N: offc,dskofc,theta0,phid,n= 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.1544E+02 0.1240E+02 -0.1261E+02 -0.1293E+02 -0.6317E+00 -0.1832E+01 crit before colath = 0.1890E+02 0.3390E+02 crit after colath = 0.1892E+02 0.3392E+02 Step 3 of 5 mtime=348 22 3 0 secs/step (sys) = 43.00 GPI run: istep= 3 power= 20.576 ctpoten= 210.238 f107= 88.244 f107a= 82.344 kp= 7.340 theta0,offa,dskofa,phid,phin,rrad,offc,dskofc= 0.1544E+02 0.1240E+02 0.0000E+00 0.0000E+00 -0.2500E+01 -0.2500E+01 -0.6882E+03 -0.6885E+03 -0.6762E+03 -0.6774E+03 0.1744E+02 0.1440E+02 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 Sending at: 348 22 2 0 inside receive cism_input_analysis byimf,loc,power,ctpoten,kp= 0.00 0.00 0.1791E+03 0.1686E+03 0.7353E+01 MIX HP power (Sitm-cpl_pedersen send with tag 2550 succeed itm-cpl_hall send with tag 2551 succeed itm-cpl_wfac send with tag 2552 succeed cpl-itm_potential recv with tag 2500 succeed cpl-itm_energy recv with tag 2501 succeed cpl-itm_flux recv with tag 2502 succeed ,N,b), maxS,N mW/m2= 0.2431E+02 0.1631E+02 0.2031E+02 0.4510E+01 0.3180E+01 MIX gpot model_ctpoten (S,N,av) = 0.1953E+03 0.2227E+03 0.2090E+03 MIX ifbad model_ctpoten (S,N,av) = 0 0 0.1946E+03 0.2204E+03 0.2075E+03 CalcCLoc S,N: offc,dskofc,theta0,phid,n= -0.1357E+01 -0.7845E+00 0.7561E+01 0.2049E+01 0.1533E+02 0.1247E+02 -0.1261E+02 -0.1292E+02 -0.9156E+00 -0.1816E+01 After Revised S,N: offc,dskofc,theta0,phid,n= 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.1533E+02 0.1247E+02 -0.1261E+02 -0.1292E+02 -0.9156E+00 -0.1816E+01 crit before colath = 0.1892E+02 0.3392E+02 crit after colath = 0.1890E+02 0.3390E+02 Step 4 of 5 mtime=348 22 4 0 secs/step (sys) = 43.47 GPI run: istep= 4 power= 20.308 ctpoten= 207.542 f107= 88.239 f107a= 82.343 kp= 7.353 theta0,offa,dskofa,phid,phin,rrad,offc,dskofc= 0.1533E+02 0.1247E+02 0.0000E+00 0.0000E+00 -0.2500E+01 -0.2500E+01 -0.6882E+03 -0.6885E+03 -0.6765E+03 -0.6774E+03 0.1733E+02 0.1447E+02 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 Sending at: 348 22 3 0 inside receive cism_input_analysis byimf,loc,power,ctpoten,kp= 0.00 0.00 0.1801E+03 0.1689E+03 0.7367E+01 MIX HP power (S,N,b), maxS,N mW/m2= 0.2550E+02 0.1665E+02 0.2107E+02 0.4720E+01 0.3184E+01 MIX gpot model_ctpoten (S,N,av) = 0.1972E+03 0.2210E+03 0.2091E+03 MIX ifbad model_ctpoten (S,N,av) = 0 0 0.1964E+03 0.2187E+03 0.2076E+03 CalcCLoc S,N: offc,dskofc,theta0,phid,n= -0.4793E+00 -0.1482E+01 0.7561E+01 0.2049E+01 0.1546E+02 0.1236E+02 -0.1261E+02 -0.1290E+02 -0.5995E+00 -0.1799E+01 After Revised S,N: offc,dskofc,theta0,phid,n= 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.1546E+02 0.1236E+02 -0.1261E+02 -0.1290E+02 -0.5995E+00 -0.1799E+01 crit before colath = 0.1890E+02 0.3390E+02 crit after colath = 0.1891E+02 0.3391E+02 Step 5 of 5 mtime=348 22 5 0 secs/step (sys) = 43.21 GPI run: istep= 5 power= 21.074 ctpoten= 207.572 f107= 88.235 f107a= 82.343 kp= 7.367 theta0,offa,dskofa,phid,phin,rrad,offc,dskofc= 0.1546E+02 0.1236E+02 0.0000E+00 0.0000E+00 -0.2500E+01 -0.2500E+01 -0.6882E+03 -0.6884E+03 -0.6761E+03 -0.6773E+03 0.1746E+02 0.1436E+02 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 Allocated 3d sech field ZG(lon= 76,lat= 36,ilev= 29) Initialized diagnostic secondary history field ZG (ix= 33): short_name = ZG long_name = Geometric Height ZG units = cm geo = T mag = F dimnames = lon lat ilev dimsizes = 76 36 29 ndims = 3 task0_only = T Allocated 3d sech field UI_ExB(lon= 76,lat= 36,ilev= 29) Initialized diagnostic secondary history field UI_ExB (ix= 16): short_name = UI_ExB long_name = Zonal Ion Drift (ExB) units = cm/s geo = T mag = F dimnames = lon lat ilev dimsizes = 76 36 29 ndims = 3 task0_only = T Allocated 3d sech field VI_ExB(lon= 76,lat= 36,ilev= 29) Initialized diagnostic secondary history field VI_ExB (ix= 17): short_name = VI_ExB long_name = Meridional Ion Drift (ExB) units = cm/s geo = T mag = F dimnames = lon lat ilev dimsizes = 76 36 29 ndims = 3 task0_only = T Allocated 3d sech field WI_ExB(lon= 76,lat= 36,ilev= 29) Initialized diagnostic secondary history field WI_ExB (ix= 18): short_name = WI_ExB long_name = Vertical Ion Drift (ExB) units = cm/s geo = T mag = F dimnames = lon lat ilev dimsizes = 76 36 29 ndims = 3 task0_only = T Allocated 2d sech field TEC(lon= 76,lat= 36) Initialized diagnostic secondary history field TEC (ix= 19): short_name = TEC long_name = Total Electron Content units = 1/cm2 geo = T mag = F dimnames = lon lat dimsizes = 76 36 0 ndims = 2 task0_only = T Allocated 3d sech field NO_COOL(lon= 76,lat= 36,lev= 29) Initialized diagnostic secondary history field NO_COOitm-cpl_pedersen send with tag 2550 succeed itm-cpl_hall send with tag 2551 succeed itm-cpl_wfac send with tag 2552 succeed L (ix= 7): short_name = NO_COOL long_name = NO Cooling units = erg/g/s geo = T mag = F dimnames = lon lat lev dimsizes = 76 36 29 ndims = 3 task0_only = T Allocated 3d sech field QJI_TN(lon= 76,lat= 36,lev= 29) Initialized diagnostic secondary history field QJI_TN (ix= 20): short_name = QJI_TN long_name = QJI_TN units = ergs/s/g geo = T mag = F dimnames = lon lat lev dimsizes = 76 36 29 ndims = 3 task0_only = T Allocated 2d sech field CUSP(lon= 76,lat= 36) Initialized diagnostic secondary history field CUSP (ix= 28): short_name = CUSP long_name = 2D CUSP units = 0-1 geo = T mag = F dimnames = lon lat dimsizes = 76 36 0 ndims = 2 task0_only = T Allocated 2d sech field DRIZZLE(lon= 76,lat= 36) Initialized diagnostic secondary history field DRIZZLE (ix= 29): short_name = DRIZZLE long_name = 2D DRIZZLE units = 0-1 geo = T mag = F dimnames = lon lat dimsizes = 76 36 0 ndims = 2 task0_only = T Allocated 2d sech field ALFA(lon= 76,lat= 36) Initialized diagnostic secondary history field ALFA (ix= 30): short_name = ALFA long_name = 2D ALFA units = keV geo = T mag = F dimnames = lon lat dimsizes = 76 36 0 ndims = 2 task0_only = T Allocated 2d sech field NFLUX(lon= 76,lat= 36) Initialized diagnostic secondary history field NFLUX (ix= 31): short_name = NFLUX long_name = 2D NFLUX units = #/cm2-s geo = T mag = F dimnames = lon lat dimsizes = 76 36 0 ndims = 2 task0_only = T Allocated 2d sech field EFLUX(lon= 76,lat= 36) Initialized diagnostic secondary history field EFLUX (ix= 32): short_name = EFLUX long_name = 2D EFLUX units = mW/m^2 geo = T mag = F dimnames = lon lat dimsizes = 76 36 0 ndims = 2 task0_only = T Allocated 3d sech field ZMAG(mlon= 81,mlat= 97,imlev= 32) Initialized diagnostic secondary history field ZMAG (ix= 34): short_name = ZMAG long_name = GEOPOTENTIAL HEIGHT (MAGNETIC) units = CM geo = F mag = T dimnames = mlon mlat imlev dimsizes = 81 97 32 ndims = 3 task0_only = T Sending at: 348 22 4 0 Allocated 2d sech field gpot(lon= 76,lat= 36) Initialized diagnostic secondary history field gpot (ix= 21): short_name = gpot long_name = Potential from M-I Coupler (geographic) units = geo = T mag = F dimnames = lon lat dimsizes = 76 36 0 ndims = 2 task0_only = T Allocated 2d sech field geng(lon= 76,lat= 36) Initialized diagnostic secondary history field geng (ix= 23): short_name = geng long_name = Energy from M-I Coupler units = geo = T mag = F dimnames = lon lat dimsizes = 76 36 0 ndims = 2 task0_only = T Allocated 2d sech field gflx(lon= 76,lat= 36) Initialized diagnostic secondary history field gflx (ix= 24): short_name = gflx long_name = Number Flux from M-I Coupler units = geo = T mag = F dimnames = lon lat dimsizes = 76 36 0 ndims = 2 task0_only = T Allocated 2d sech field gpotm(lon= 76,lat= 36) Initialized diagnostic secondary history field gpotm (ix= 22): short_name = gpotm long_name = Potential from M-I Coupler (geographic) units = geo = T mag = F dimnames = lon lat dimsizes = 76 36 0 ndims = 2 task0_only = T Allocated 2d sech field gzigm1(lon= 76,lat= 36) Initialized diagnostic secondary history field gzigm1 (ix= 25): short_name = gzigm1 long_name = Pedersen Conductance (geographic) units = geo = T mag = F dimnames = lon lat dimsizes = 76 36 0 ndims = 2 task0_only = T Allocated 2d sech field gzigm2(lon= 76,lat= 36) Initialized diagnostic secondary history field gzigm2 (ix= 26): short_name = gzigm2 long_name = Hall Conductance (geographic) units = geo = T mag = F dimnames = lon lat dimsizes = 76 36 0 ndims = 2 task0_only = T Allocated 2d sech field gnsrhs(lon= 76,lat= 36) Initialized diagnostic secondary history field gnsrhs (ix= 27): short_name = gnsrhs long_name = Height-integrated neutral wind units = geo = T mag = F dimnames = lon lat dimsizes = 76 36 0 ndims = 2 task0_only = T Created netcdf file AGU_CMITs_old_2006-12-14T22-00-05_pold.nc ncid= 6 ------------------------------------------------------------------------ Write TGCM PRIMARY HISTORY Diskfile: AGU_CMITs_old_2006-12-14T22-00-05_pold.nc label = AGU_CMITs model_name = tiegcm model_version = tiegcm1.94.2 create_date= 08/31/12 16:15:48 write_date = 08/31/12 16:15:48 logname = emery host = puffin system = LINUX run_type = continuation source_file = AGU_CMITs_tie_2006-12-14T22-00-00.nc output_file = AGU_CMITs_old_2006-12-14T22-00-05_pold.nc source_mtime = 348 22 0 initial_year = 0 initial_day = 0 initial_mtime= 0 0 0 type = primary ihist = 1 (nth history on history file) delhmins= 5 (delta minutes between histories) calendar year,day = 2006,348 (model IS being advanced in calendar time) modeltime = 348,22, 5, 0 (model time day,hour,minute,seconds) time = 1325.00 (minutes in current day) ut = 22.08 (ut hours) step = 60 (time step in seconds) iter = 502445 (number of steps from 0,0,0) nlat = 36 (number of latitudes) nlon = 72 (number of longitudes) nlev = 29 (number of levels) zmtop = 7.250 (top midpoint level) zmbot = -6.750 (bottom midpoint level) zitop = 7.000 (top interface level) zibot = -7.000 (bottom interface level) dtide = 0.0E+00 0.0 (amp/phase of diurnal tide) 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 (amp/phase of semi-diurnal tide) f107d = 0.8823E+02 (daily solar flux) f107a = 0.8234E+02 (average solar flux) hpower = 0.2107E+02 (Gw) ctpoten = 0.2076E+03 (Volts) kp = 0.7367E+01 () byimf = 0.0000E+00 bzimf = 0.1000E+37 swvel = 0.1000E+37 swden = 0.1000E+37 al = 0.1000E+37 e1,e2 = 0.1800E+01 0.3318E+01 (ergs/cm2/s) alfac,alfad = 0.1000E+00 0.5000E+00 (keV) ec,ed = 0.7624E-01 0.1384E-01 (ergs/cm2/s) colfac = 0.15E+01 joulefac = 0.10E+01 p0 = 0.50E-06 p0_model= 0.50E-03 grav = 0.87E+03 nflds = 24 (number of model fields) gswm_mi_di_ncfile = /hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc gswm_mi_sdi_ncfile = /hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc gswm_nm_di_ncfile = [none] gswm_nm_sdi_ncfile = [none] see_ncfile = [none] gpi_ncfile = /hao/aim1/tgcm/data/gpi_2000001-2011059.nc ncep_ncfile = [none] imf_ncfile = [none] saber_ncfile = [none] tidi_ncfile = [none] tuv_lbc_intop= 0 (if 1, then lbc of t,u,v are stored in top k slot (old histories)) LBC = -7.00 (lower boundary interface level) ntask_mpi = 0 (number of MPI tasks) coupled_cmit = 1 (1 if coupled with CISM/CMIT, 0 otherwise) There are 24 fields on this history, as follows: TN UN VN O2 O1 N4S NO OP N2D TI TE NE O2P OMEGA Z POTEN TN_NM UN_NM VN_NM O2_NM O1_NM N4S_NM NO_NM OP_NM ------------------------------------------------------------------------ Wrote primary history 348,22, 5 to AGU_CMITs_old_2006-12-14T22-00-05_pold.nc ( 1 of======= IC_Finalize was successful. ======= 1) Starting secondary history time series 1 at model time 348 22 5 0 Created netcdf file AGU_CMITs_sech_tie_2006-12-14T22-00-00_testb_old.nc ncid= 6 ------------------------------------------------------------------------ Write TGCM SECONDARY HISTORY Diskfile: AGU_CMITs_sech_tie_2006-12-14T22-00-00_testb_old.nc label = AGU_CMITs model_name = tiegcm model_version = tiegcm1.94.2 create_date= 08/31/12 16:15:49 write_date = 08/31/12 16:15:49 logname = emery host = puffin system = LINUX run_type = continuation source_file = AGU_CMITs_tie_2006-12-14T22-00-00.nc output_file = AGU_CMITs_sech_tie_2006-12-14T22-00-00_testb_old.nc source_mtime = 348 22 0 initial_year = 0 initial_day = 0 initial_mtime= 0 0 0 type = secondary ihist = 1 (nth history on history file) delhmins= 5 (delta minutes between histories) calendar year,day = 2006,348 (model IS being advanced in calendar time) modeltime = 348,22, 5, 0 (model time day,hour,minute,seconds) time = 1325.00 (minutes in current day) ut = 22.08 (ut hours) step = 60 (time step in seconds) iter = 502445 (number of steps from 0,0,0) nlat = 36 (number of latitudes) nlon = 72 (number of longitudes) nlev = 29 (number of levels) zmtop = 7.250 (top midpoint level) zmbot = -6.750 (bottom midpoint level) zitop = 7.000 (top interface level) zibot = -7.000 (bottom interface level) dtide = 0.0E+00 0.0 (amp/phase of diurnal tide) 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 (amp/phase of semi-diurnal tide) f107d = 0.8823E+02 (daily solar flux) f107a = 0.8234E+02 (average solar flux) hpower = 0.2107E+02 (Gw) ctpoten = 0.2076E+03 (Volts) kp = 0.7367E+01 () byimf = 0.0000E+00 bzimf = 0.1000E+37 swvel = 0.1000E+37 swden = 0.1000E+37 al = 0.1000E+37 e1,e2 = 0.1800E+01 0.3318E+01 (ergs/cm2/s) alfac,alfad = 0.1000E+00 0.5000E+00 (keV) ec,ed = 0.7624E-01 0.1384E-01 (ergs/cm2/s) colfac = 0.15E+01 joulefac = 0.10E+01 p0 = 0.50E-06 p0_model= 0.50E-03 grav = 0.87E+03 nflds = 0 (number of model fields) gswm_mi_di_ncfile = /hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc gswm_mi_sdi_ncfile = /hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc gswm_nm_di_ncfile = [none] gswm_nm_sdi_ncfile = [none] see_ncfile = [none] gpi_ncfile = /hao/aim1/tgcm/data/gpi_2000001-2011059.nc ncep_ncfile = [none] imf_ncfile = [none] saber_ncfile = [none] tidi_ncfile = [none] tuv_lbc_intop= 0 (if 1, then lbc of t,u,v are stored in top k slot (old histories)) LBC = -7.00 (lower boundary interface level) ntask_mpi = 0 (number of MPI tasks) coupled_cmit = 1 (1 if coupled with CISM/CMIT, 0 otherwise) ------------------------------------------------------------------------ Wrote secondary history 348,22, 5 to AGU_CMITs_sech_tie_2006-12-14T22-00-00_testb_old.nc ( 1 of 1) nstep= 5 step= 60 Model simulation time = 300 secs (minutes= 5.00, hours= 0.08, days= 0.003472) Cpu time for run = 205.53 ------------------------------------------------------------------------ TIMER (system_clock): Elapsed run time = 205.53 (secs) 3.43 (mins) 0.06 (hrs) Elapsed init time = 28.23 secs, 0.008 hrs, 13.7% (includes source i/o) Elapsed step time = 176.18 secs, 0.049 hrs, 85.7% Elapsed prep time = 348.65 secs, 0.097 hrs, 169.6% Elapsed dynamics = 3.97 secs, 0.001 hrs, 1.9% Elapsed qrj = 0.62 secs, 0.000 hrs, 0.3% Elapsed oplus = 0.38 secs, 0.000 hrs, 0.2% Elapsed cmpminor = 0.33 secs, 0.000 hrs, 0.2% (n4s,n2d,no) Elapsed cmpmajor = 0.25 secs, 0.000 hrs, 0.1% (o2, o) Elapsed dt = 0.12 secs, 0.000 hrs, 0.1% Elapsed duv = 0.18 secs, 0.000 hrs, 0.1% Elapsed dynamo = 0.84 secs, 0.000 hrs, 0.4% Elapsed i/o = 2.22 secs, 0.001 hrs, 1.1% Elapsed Primary io = 1.11 secs, 0.000 hrs, 0.5% Elapsed Secondary = 1.11 secs, 0.000 hrs, 0.5% Average secs/step = 35.24 Average mins/simulated day = 986.53 ------------------------------------------------------------------------ End execution of tiegcm1.94.2 at 08/31/12 16:16:31 NORMAL EXIT