WARNING: Unable to read mpd.hosts or list of hosts isn't provided. MPI job will be run on the current machine only. ======================================================================== Begin execution of tiegcm1.94 at 01/05/12 11:16:45 Host = puffin System = LINUX Logname = emery ======================================================================== init_timer: level= 1 rtc=F sys=T mp_init: ntask= 1 mytid= 0 Current working directory (cwd) = /hao/aim3/emery/tiegcm_1.94.2/tiegcm-linux Process ID (pid) = 25562 &tgcm_input LABEL = 'tiegcm_sres_w01' START_YEAR = 2006 START_DAY = 348 CALENDAR_ADVANCE = 1 START = 348,6,0 STOP = 348,7,0 STEP = 120 HIST = 0,1,0 OUTPUT = 'p_06348_testw01.nc' MXHIST_PRIM = 13 SECSTART = 348,6,0 SECSTOP = 348,7,0 SECHIST = 0,1,0 SECOUT = 's_06348_testw01_002.nc' MXHIST_SECH = 24 SECFLDS = 'TN','UN','VN','WN','O1','NO','DEN','NE','TE','TI','TEC', 'O2','POTEN','UI_ExB','VI_ExB','WI_ExB','SIGMA_PED','SIGMA_HAL', 'HMF2','NMF2','QJOULE','QJOULE_INTEG','Z','ZG','ALFA','EFLUX' 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' POTENTIAL_MODEL = 'WEIMER01' GPI_NCFILE = '$TGCMDATA/gpi_2000001-2011212.nc' IMF_NCFILE = '$TGCMDATA/imf_OMNI_2006001-2006365.nc' AURORA = 1 COLFAC = 1.5 / Reading input data... Completed successful read of namelist inputs. mkntask: i= 1 j= 1 i*j= 1 ntask= 1 Input: mkntask chose ntask_lon= 1 ntask_lat= 1 (ntask= 1) Will use the Weimer 2001 potential model weimer_ncfile=/hao/aim1/tgcm/data/weimer2001_coeffs.nc Note INPUT: Both gpi_ncfile and imf_ncfile are specified. In this case, only f10.7 flux will be used from the gpi file. Note INPUT: gpi_ncfile has been specified with Weimer convection model. In this case, only f10.7 flux will be used from the gpi file. Input: Expanded secout file s_06348_testw01_002.nc to s_06348_testw01_002.nc INPUT NOTE: adding mandatory field ZMAG to secondary history fields (field 27) ------------------------------------------------------------------------ USER INPUT PARAMETERS: label = tiegcm_sres_w01 (optional text label for current run) High-lat electric potential model: potential_model = WEIMER01 weimer coefs: weimer_ncfile = /hao/aim1/tgcm/data/weimer2001_coeffs.nc gpi run: gpi_ncfile = $TGCMDATA/gpi_2000001-2011212.nc imf run: imf_ncfile = $TGCMDATA/imf_OMNI_2006001-2006365.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 = 120 (model timestep (seconds)) ntask_lon = 1 (number of mpi tasks in longitude dimension) ntask_lat = 1 (number of mpi tasks in latitude dimension) total tasks = ntask_lon*ntask_lat = 1 output (primary history output files) = p_06348_testw01.nc, start (model start times) = 348, 6, 0 stop (model stop times) = 348, 7, 0 hist (primary history disk write frequencies) = 0, 1, 0 Maxmimum number of histories per primary file = 13 secout (secondary history output files)= s_06348_testw01_002.nc, secstart (secondary history start times) = 348, 6, 0 secstop (secondary history stop times) = 348, 7, 0 sechist (secondary history disk write frequencies) = 0, 1, 0 secflds (secondary history fields) = TN UN VN WN O1 NO DEN NE TE TI TEC O2 POTEN UI_ExB VI_ExB WI_ExB SIGMA_PED SIGMA_HAL HMF2 NMF2 QJOULE QJOULE_INTEGZ ZG ALFA EFLUX ZMAG Maximum number of histories per secondary file = 24 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 will 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.1000E+37 (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 ------------------------------------------------------------------------ ntask= 1 ntaski= 1 ntaskj= 1 Task Table: j= -1 itask_table(:,j)= -1 -1 -1 j= 0 itask_table(:,j)= -1 0 -1 j= 1 itask_table(:,j)= -1 -1 -1 mytid= 0 mytidi,j= 0 0 lat0,1= 1 36 (36) lon0,1= 1 76 (76) ncells=2736 Task 0: tasks( 0)%mytid = 0 tasks( 0)%mytidi= 0 tasks( 0)%mytidj= 0 tasks( 0)%nlats = 36 tasks( 0)%nlons = 76 tasks( 0)%lat0 = 1 tasks( 0)%lat1 = 36 tasks( 0)%lon0 = 1 tasks( 0)%lon1 = 76 tasks( 0)%mylats= 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 tasks( 0)%mylons= 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 Model version = tiegcm1.94 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 = 120.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= 30 Allocated 3d data for all fields on root task: nlevp1*nlonp4*nlat*nf4d = 29* 76* 36* 26 = 2062944 words *8 = 16503552 bytes. allocdata: all tasks allocate fzg: nlevp1= 29 nlonp4= 76 nlat= 36 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 = 30 (Number of time steps this run) iter = 250740 (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 = 13 (Maximum number of primary histories per file) Secondary Histories: nseries_sech = 1 (Number of secondary time series) nsech_total = 2 (Number of secondary histories to be written) nfiles_sech = 1 (Number of secondary output files to be written) mxhist_sech = 24 (Maximum number of secondary histories per file) nfsech = 27 (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: WN secondary history field 5: O1 secondary history field 6: NO secondary history field 7: DEN secondary history field 8: NE secondary history field 9: TE secondary history field 10: TI secondary history field 11: TEC secondary history field 12: O2 secondary history field 13: POTEN secondary history field 14: UI_ExB secondary history field 15: VI_ExB secondary history field 16: WI_ExB secondary history field 17: SIGMA_PED secondary history field 18: SIGMA_HAL secondary history field 19: HMF2 secondary history field 20: NMF2 secondary history field 21: QJOULE secondary history field 22: QJOULE_INTEG secondary history field 23: Z secondary history field 24: ZG secondary history field 25: ALFA secondary history field 26: EFLUX secondary history field 27: 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=p_06348_testw01.nc Getfile: Found file p_06348_testw01.nc ------------------------------------------------------------------------ Acquired source history file p_06348_testw01.nc (disk file is p_06348_testw01.nc) Reading source history from diskfile p_06348_testw01.nc: nc_rdhist: seeking 348 6 0 found 348 0 0 n= 1 nc_rdhist: seeking 348 6 0 found 348 1 0 n= 2 nc_rdhist: seeking 348 6 0 found 348 2 0 n= 3 nc_rdhist: seeking 348 6 0 found 348 3 0 n= 4 nc_rdhist: seeking 348 6 0 found 348 4 0 n= 5 nc_rdhist: seeking 348 6 0 found 348 5 0 n= 6 nc_rdhist: seeking 348 6 0 found 348 6 0 n= 7 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.1544E+03 0.1034E+04 Read field UN 3d subdomain min,max= -0.3644E+05 0.2816E+05 Read field VN 3d subdomain min,max= -0.3399E+05 0.3208E+05 Read field O2 3d subdomain min,max= 0.3172E-04 0.2308E+00 Read field O1 3d subdomain min,max= 0.3138E-02 0.9978E+00 Read field N4S 3d subdomain min,max= 0.1000E-11 0.1977E-01 Read field NO 3d subdomain min,max= 0.1248E-06 0.1232E-02 Read field OP 3d subdomain min,max= 0.0000E+00 0.1642E+07 Read field N2D 3d subdomain min,max= 0.0000E+00 0.1821E-02 Read field TI 3d subdomain min,max= 0.1544E+03 0.1915E+04 Read field TE 3d subdomain min,max= 0.1544E+03 0.4002E+04 Read field NE 3d subdomain min,max= 0.2839E+03 0.1633E+07 Read field O2P 3d subdomain min,max= 0.0000E+00 0.1184E+06 Read field OMEGA 3d subdomain min,max= -0.1213E-03 0.3884E-03 Read field Z 3d subdomain min,max= 0.9578E+07 0.5857E+08 Read field POTEN 3d subdomain min,max= -0.2064E+05 0.1771E+05 Read field TN_NM 3d subdomain min,max= 0.1544E+03 0.1033E+04 Read field UN_NM 3d subdomain min,max= -0.3649E+05 0.2796E+05 Read field VN_NM 3d subdomain min,max= -0.3432E+05 0.3191E+05 Read field O2_NM 3d subdomain min,max= 0.3162E-04 0.2308E+00 Read field O1_NM 3d subdomain min,max= 0.3135E-02 0.9973E+00 Read field N4S_NM 3d subdomain min,max= 0.1000E-11 0.1977E-01 Read field NO_NM 3d subdomain min,max= 0.1244E-06 0.1232E-02 Read field OP_NM 3d subdomain min,max= 0.0000E+00 0.1644E+07 Read LBC from source history: i= 89 sh%lbc= -7.00 ------------------------------------------------------------------------ Read TGCM PRIMARY HISTORY (source history) Diskfile: p_06348_testw01.nc label = tiegcm1.93 model_name = tiegcm model_version = tiegcm1.94 create_date= 01/05/12 10:43:45 write_date = logname = foster host = be0901en.ucar.ed system = AIX run_type = source_file = output_file = source_mtime = 348 6 0 initial_year = 2006 initial_day = 348 initial_mtime= 348 0 0 type = primary ihist = 7 (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, 6, 0, 0 (model time day,hour,minute,seconds) time = 360.00 (minutes in current day) ut = 6.00 (ut hours) step = 120 (time step in seconds) iter = 250740 (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.9056E+02 (daily solar flux) f107a = 0.8240E+02 (average solar flux) hpower = 0.4073E+02 (Gw) ctpoten = 0.2572E+02 (Volts) kp = 0.1000E+37 () byimf = -0.1017E+01 bzimf = 0.3450E+01 swvel = 0.6073E+03 swden = 0.1217E+01 al = 0.1000E+37 e1,e2 = 0.2653E+01 0.5480E+01 (ergs/cm2/s) alfac,alfad = 0.1000E+00 0.5000E+00 (keV) ec,ed = 0.1026E+00 0.2564E-01 (ergs/cm2/s) colfac = 0.15E+01 joulefac = 0.15E+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-2011212.nc ncep_ncfile = [none] imf_ncfile = /hao/aim1/tgcm/data/imf_OMNI_2006001-2006365.nc 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 = 4 (number of MPI tasks) coupled_cmit = 0 (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=250740 nstep= 30 Starting secondary history time series 1 at model time 348 6 0 0 Created netcdf file s_06348_testw01_002.nc ncid= 7 ------------------------------------------------------------------------ Write TGCM SECONDARY HISTORY Diskfile: s_06348_testw01_002.nc label = tiegcm_sres_w01 model_name = tiegcm model_version = tiegcm1.94 create_date= 01/05/12 11:16:46 write_date = 01/05/12 11:16:46 logname = emery host = puffin system = LINUX run_type = continuation source_file = p_06348_testw01.nc output_file = s_06348_testw01_002.nc source_mtime = 348 6 0 initial_year = 0 initial_day = 0 initial_mtime= 0 0 0 type = secondary ihist = 1 (nth history on history file) delhmins= 60 (delta minutes between histories) calendar year,day = 2006,348 (model IS being advanced in calendar time) modeltime = 348, 6, 0, 0 (model time day,hour,minute,seconds) time = 360.00 (minutes in current day) ut = 6.00 (ut hours) step = 120 (time step in seconds) iter = 250740 (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.9056E+02 (daily solar flux) f107a = 0.8240E+02 (average solar flux) hpower = 0.4073E+02 (Gw) ctpoten = 0.2572E+02 (Volts) kp = 0.1000E+37 () byimf = -0.1017E+01 bzimf = 0.3450E+01 swvel = 0.6073E+03 swden = 0.1217E+01 al = 0.1000E+37 e1,e2 = 0.0000E+00 0.0000E+00 (ergs/cm2/s) alfac,alfad = 0.1000E+00 0.5000E+00 (keV) ec,ed = 0.0000E+00 0.0000E+00 (ergs/cm2/s) colfac = 0.15E+01 joulefac = 0.15E+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 = $TGCMDATA/gpi_2000001-2011212.nc ncep_ncfile = [none] imf_ncfile = $TGCMDATA/imf_OMNI_2006001-2006365.nc 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 = 1 (number of MPI tasks) coupled_cmit = 0 (1 if coupled with CISM/CMIT, 0 otherwise) ------------------------------------------------------------------------ Wrote secondary history 348, 6, 0 to s_06348_testw01_002.nc ( 1 of 2) ------------------------------------------------------------------------ Getfile: remote=$TGCMDATA/gpi_2000001-2011212.nc Getfile: Found file /hao/aim1/tgcm/data/gpi_2000001-2011212.nc ------------------------------------------------------------------------ ------------------------------------------------------------------------ RDGPI: read GPI data file: Number of days on GPI data file = 4230 First and last year-days = 2000003 2011210 Completed read from GPI data file /hao/aim1/tgcm/data/gpi_2000001-2011212.nc ------------------------------------------------------------------------ ------------------------------------------------------------------------ GETGPI: get geophysical indices from data file /hao/aim1/tgcm/data/gpi_2000001-2011212.nc Initial requested iyear=2006 iday=348 iutsec= 21720 Obtained the following GPI at the requested date and time: f107 = 0.9056E+02 f107a = 0.8240E+02 power = 0.2857E+02 ctpoten = 0.5089E+02 ------------------------------------------------------------------------ ------------------------------------------------------------------------ Getfile: remote=$TGCMDATA/imf_OMNI_2006001-2006365.nc Getfile: Found file /hao/aim1/tgcm/data/imf_OMNI_2006001-2006365.nc ------------------------------------------------------------------------ ------------------------------------------------------------------------ RDIMF: read IMF data file: Opened netcdf imf data file /hao/aim1/tgcm/data/imf_OMNI_2006001-2006365.nc rdimf: date min,max= 2006001.0000 2007001.0035 rdimf: bx min,max= -18.37 15.56 rdimf: by min,max= -17.90 17.29 rdimf: bz min,max= -17.74 17.41 rdimf: swvel min,max= 264.37 916.35 rdimf: swden min,max= 0.23 69.58 ------------------------------------------------------------------------ ------------------------------------------------------------------------ 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 ------------------------------------------------------------------------ Step 1 of 30 mtime=348 6 2 GPI+IMF run: istep= 1 GPI f107= 90.560 f107a= 82.403 power from bzimf,swvel= 40.784 aurora_cons: cusp: alfac= 0.100 ec= 0.103 fc= 0.3204E+09 drizzle: alfad= 0.500 ed= 0.026 fd= 0.1602E+08 auroral radius = max of rhp,rcp= 21.640 19.815 roth, rote (MLT) = 0.345 -0.140 1/e-widths = h1,h2= 2.350 3.733 energy flux = e1,e2= 2.655 5.486 add_sproton = F Copyright 2001 Dan Weimer/MHC potential model ------------------------------------------------------------------------ Getfile: remote=/hao/aim1/tgcm/data/weimer2001_coeffs.nc Getfile: Found file /hao/aim1/tgcm/data/weimer2001_coeffs.nc ------------------------------------------------------------------------ Opened weimer coefficients file /hao/aim1/tgcm/data/weimer2001_coeffs.nc Read Weimer coefs, MaxJ = 3 calccloc: ih,bz,y,of18,24,crad= 1 3.21 -1.05 -1.07 -0.78 14.76 nogood=1: ih,bz,y,mlt06,18,of18,24,crad,0 = 2 3.21 -1.05 14.96 19.76 461.33 0.00 93.69 -371.33 calccloc: ih,bz,y,of18,24,crad= 2 3.21 -1.05 0.00 4.20 10.00 weimer01: CPS,N ctpoten= 25.6 25.8 25.7 Step 2 of 30 mtime=348 6 4 0 secs/step (sys) = 2.72 GPI+IMF run: istep= 2 GPI f107= 90.561 f107a= 82.403 power from bzimf,swvel= 40.840 calccloc: ih,bz,y,of18,24,crad= 1 2.97 -1.09 -2.11 -0.78 13.72 nogood=1: ih,bz,y,mlt06,18,of18,24,crad,0 = 2 2.97 -1.09 14.99 19.79 461.33 1.07 93.69 -371.33 calccloc: ih,bz,y,of18,24,crad= 2 2.97 -1.09 0.00 4.20 10.00 weimer01: CPS,N ctpoten= 26.0 25.2 25.6 Step 3 of 30 mtime=348 6 6 0 secs/step (sys) = 0.49 GPI+IMF run: istep= 3 GPI f107= 90.561 f107a= 82.403 power from bzimf,swvel= 40.744 calccloc: ih,bz,y,of18,24,crad= 1 2.73 -1.05 -2.11 -0.78 13.73 nogood=1: ih,bz,y,mlt06,18,of18,24,crad,0 = 2 2.73 -1.05 15.03 19.83 461.33 1.07 93.69 -371.33 calccloc: ih,bz,y,of18,24,crad= 2 2.73 -1.05 0.00 4.20 10.00 weimer01: CPS,N ctpoten= 26.6 24.3 25.4 Step 4 of 30 mtime=348 6 8 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 4 GPI f107= 90.561 f107a= 82.403 power from bzimf,swvel= 40.495 calccloc: ih,bz,y,of18,24,crad= 1 2.48 -0.93 -2.11 -0.78 13.73 nogood=1: ih,bz,y,mlt06,18,of18,24,crad,0 = 2 2.48 -0.93 15.06 19.86 461.33 1.07 93.69 -371.33 calccloc: ih,bz,y,of18,24,crad= 2 2.48 -0.93 0.00 4.20 10.00 weimer01: CPS,N ctpoten= 27.2 22.8 25.0 Step 5 of 30 mtime=348 6 10 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 5 GPI f107= 90.562 f107a= 82.403 power from bzimf,swvel= 40.246 calccloc: ih,bz,y,of18,24,crad= 1 2.23 -0.82 -2.11 -0.78 13.73 nogood=1: ih,bz,y,mlt06,18,of18,24,crad,0 = 2 2.23 -0.82 15.10 19.90 461.33 1.07 93.69 -371.33 calccloc: ih,bz,y,of18,24,crad= 2 2.23 -0.82 0.00 4.20 10.00 weimer01: CPS,N ctpoten= 27.9 21.3 24.6 Step 6 of 30 mtime=348 6 12 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 6 GPI f107= 90.562 f107a= 82.402 power from bzimf,swvel= 40.061 calccloc: ih,bz,y,of18,24,crad= 1 2.04 -0.39 -1.07 -0.78 14.76 calccloc: ih,bz,y,of18,24,crad= 2 2.04 -0.39 0.00 -1.96 13.81 weimer01: CPS,N ctpoten= 29.2 19.7 24.5 Step 7 of 30 mtime=348 6 14 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 7 GPI f107= 90.562 f107a= 82.402 power from bzimf,swvel= 39.877 calccloc: ih,bz,y,of18,24,crad= 1 1.84 0.03 -1.07 -0.94 14.77 calccloc: ih,bz,y,of18,24,crad= 2 1.84 0.03 0.00 -3.03 13.99 weimer01: CPS,N ctpoten= 31.0 21.6 26.3 Step 8 of 30 mtime=348 6 16 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 8 GPI f107= 90.562 f107a= 82.402 power from bzimf,swvel= 39.718 calccloc: ih,bz,y,of18,24,crad= 1 1.85 0.36 -1.07 -0.94 14.77 calccloc: ih,bz,y,of18,24,crad= 2 1.85 0.36 0.00 3.33 14.06 weimer01: CPS,N ctpoten= 31.8 22.9 27.3 Step 9 of 30 mtime=348 6 18 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 9 GPI f107= 90.563 f107a= 82.402 power from bzimf,swvel= 39.585 calccloc: ih,bz,y,of18,24,crad= 1 2.07 0.58 -1.07 -0.94 14.77 calccloc: ih,bz,y,of18,24,crad= 2 2.07 0.58 0.00 1.11 13.71 weimer01: CPS,N ctpoten= 31.5 23.0 27.3 Step 10 of 30 mtime=348 6 20 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 10 GPI f107= 90.563 f107a= 82.402 power from bzimf,swvel= 39.453 calccloc: ih,bz,y,of18,24,crad= 1 2.30 0.80 0.00 -2.05 13.82 calccloc: ih,bz,y,of18,24,crad= 2 2.30 0.80 -1.07 2.18 14.90 weimer01: CPS,N ctpoten= 31.3 23.2 27.2 Step 11 of 30 mtime=348 6 22 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 11 GPI f107= 90.563 f107a= 82.402 power from bzimf,swvel= 39.573 calccloc: ih,bz,y,of18,24,crad= 1 2.55 0.86 -1.07 -0.94 14.77 calccloc: ih,bz,y,of18,24,crad= 2 2.55 0.86 -1.07 1.11 14.78 weimer01: CPS,N ctpoten= 30.6 22.5 26.5 Step 12 of 30 mtime=348 6 24 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 12 GPI f107= 90.563 f107a= 82.402 power from bzimf,swvel= 39.694 calccloc: ih,bz,y,of18,24,crad= 1 2.81 0.92 -1.07 -2.01 14.88 calccloc: ih,bz,y,of18,24,crad= 2 2.81 0.92 -1.07 0.04 14.74 weimer01: CPS,N ctpoten= 30.0 21.8 25.9 Step 13 of 30 mtime=348 6 26 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 13 GPI f107= 90.563 f107a= 82.401 power from bzimf,swvel= 39.791 calccloc: ih,bz,y,of18,24,crad= 1 3.03 0.71 -1.07 -2.01 14.88 calccloc: ih,bz,y,of18,24,crad= 2 3.03 0.71 0.00 1.19 15.86 weimer01: CPS,N ctpoten= 28.9 20.4 24.7 Step 14 of 30 mtime=348 6 28 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 14 GPI f107= 90.564 f107a= 82.401 power from bzimf,swvel= 39.863 calccloc: ih,bz,y,of18,24,crad= 1 3.23 0.22 -1.07 -2.01 14.88 nogood=1: ih,bz,y,mlt06,18,of18,24,crad,0 = 2 3.23 0.22 15.43 19.63 461.33 -1.07 93.68 -371.33 calccloc: ih,bz,y,of18,24,crad= 2 3.23 0.22 0.00 4.20 10.00 weimer01: CPS,N ctpoten= 27.6 20.2 23.9 Step 15 of 30 mtime=348 6 30 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 15 GPI f107= 90.564 f107a= 82.401 power from bzimf,swvel= 39.935 calccloc: ih,bz,y,of18,24,crad= 1 3.43 -0.27 -1.07 -0.86 14.76 nogood=1: ih,bz,y,mlt06,18,of18,24,crad,0 = 2 3.43 -0.27 15.16 19.66 461.33 0.00 93.69 -371.33 calccloc: ih,bz,y,of18,24,crad= 2 3.43 -0.27 0.00 4.20 10.00 weimer01: CPS,N ctpoten= 26.2 22.7 24.4 Step 16 of 30 mtime=348 6 32 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 16 GPI f107= 90.564 f107a= 82.401 power from bzimf,swvel= 39.868 calccloc: ih,bz,y,of18,24,crad= 1 3.30 -0.43 -1.07 -0.86 14.76 nogood=1: ih,bz,y,mlt06,18,of18,24,crad,0 = 2 3.30 -0.43 15.20 19.70 461.33 0.00 93.69 -371.33 calccloc: ih,bz,y,of18,24,crad= 2 3.30 -0.43 0.00 4.20 10.00 weimer01: CPS,N ctpoten= 26.2 22.9 24.5 Step 17 of 30 mtime=348 6 34 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 17 GPI f107= 90.564 f107a= 82.401 power from bzimf,swvel= 39.802 calccloc: ih,bz,y,of18,24,crad= 1 3.18 -0.59 -1.07 -0.78 14.76 nogood=1: ih,bz,y,mlt06,18,of18,24,crad,0 = 2 3.18 -0.59 15.23 19.73 461.33 0.00 93.69 -371.33 calccloc: ih,bz,y,of18,24,crad= 2 3.18 -0.59 0.00 4.20 10.00 weimer01: CPS,N ctpoten= 26.1 23.2 24.6 Step 18 of 30 mtime=348 6 36 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 18 GPI f107= 90.565 f107a= 82.401 power from bzimf,swvel= 39.740 calccloc: ih,bz,y,of18,24,crad= 1 2.87 -0.87 -1.07 -0.78 14.76 nogood=1: ih,bz,y,mlt06,18,of18,24,crad,0 = 2 2.87 -0.87 14.97 19.77 461.33 0.00 93.69 -371.33 calccloc: ih,bz,y,of18,24,crad= 2 2.87 -0.87 0.00 4.20 10.00 weimer01: CPS,N ctpoten= 26.5 23.4 24.9 Step 19 of 30 mtime=348 6 38 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 19 GPI f107= 90.565 f107a= 82.401 power from bzimf,swvel= 39.683 calccloc: ih,bz,y,of18,24,crad= 1 2.37 -1.26 -2.11 0.26 13.70 nogood=1: ih,bz,y,mlt06,18,of18,24,crad,0 = 2 2.37 -1.26 14.71 19.81 461.33 1.07 93.69 -371.33 calccloc: ih,bz,y,of18,24,crad= 2 2.37 -1.26 0.00 4.20 10.00 weimer01: CPS,N ctpoten= 27.1 24.0 25.5 Step 20 of 30 mtime=348 6 40 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 20 GPI f107= 90.565 f107a= 82.400 power from bzimf,swvel= 39.626 nogood=1: ih,bz,y,mlt06,18,of18,24,crad,0 = 1 1.88 -1.65 2.14 18.94 -3.12 5.35 14.02 12.69 calccloc: ih,bz,y,of18,24,crad= 1 1.88 -1.65 0.00 4.20 10.00 nogood=1: ih,bz,y,mlt06,18,of18,24,crad,0 = 2 1.88 -1.65 14.14 19.84 6.83 2.02 13.87 6.83 calccloc: ih,bz,y,of18,24,crad= 2 1.88 -1.65 0.00 4.20 10.00 weimer01: CPS,N ctpoten= 29.1 26.0 27.6 Step 21 of 30 mtime=348 6 42 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 21 GPI f107= 90.565 f107a= 82.400 power from bzimf,swvel= 39.322 calccloc: ih,bz,y,of18,24,crad= 1 1.46 -1.80 -3.12 4.40 13.63 nogood=1: ih,bz,y,mlt06,18,of18,24,crad,0 = 2 1.46 -1.80 13.28 19.88 6.83 2.02 13.88 6.83 calccloc: ih,bz,y,of18,24,crad= 2 1.46 -1.80 0.00 4.20 10.00 weimer01: CPS,N ctpoten= 31.6 27.7 29.7 Step 22 of 30 mtime=348 6 44 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 22 GPI f107= 90.565 f107a= 82.400 power from bzimf,swvel= 39.019 calccloc: ih,bz,y,of18,24,crad= 1 1.04 -1.96 -7.91 3.45 16.26 nogood=1: ih,bz,y,mlt06,18,of18,24,crad,0 = 2 1.04 -1.96 12.71 19.91 5.80 2.99 12.11 5.80 calccloc: ih,bz,y,of18,24,crad= 2 1.04 -1.96 0.00 4.20 10.00 weimer01: CPS,N ctpoten= 35.2 31.0 33.1 Step 23 of 30 mtime=348 6 46 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 23 GPI f107= 90.565 f107a= 82.400 power from bzimf,swvel= 38.765 calccloc: ih,bz,y,of18,24,crad= 1 0.65 -2.19 -6.97 3.45 14.57 nogood=1: ih,bz,y,mlt06,18,of18,24,crad,0 = 2 0.65 -2.19 12.15 19.65 5.80 4.03 12.50 5.80 calccloc: ih,bz,y,of18,24,crad= 2 0.65 -2.19 0.00 4.20 10.00 weimer01: CPS,N ctpoten= 40.7 35.9 38.3 Step 24 of 30 mtime=348 6 48 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 24 GPI f107= 90.566 f107a= 82.400 power from bzimf,swvel= 38.562 calccloc: ih,bz,y,of18,24,crad= 1 0.29 -2.50 -6.97 2.26 14.24 calccloc: ih,bz,y,of18,24,crad= 2 0.29 -2.50 5.80 4.03 12.41 weimer01: CPS,N ctpoten= 47.3 41.5 44.4 Step 25 of 30 mtime=348 6 50 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 25 GPI f107= 90.566 f107a= 82.400 power from bzimf,swvel= 40.926 calccloc: ih,bz,y,of18,24,crad= 1 -0.07 -2.81 -6.97 2.26 14.24 calccloc: ih,bz,y,of18,24,crad= 2 -0.07 -2.81 2.05 2.99 12.05 weimer01: CPS,N ctpoten= 53.6 47.2 50.4 Step 26 of 30 mtime=348 6 52 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 26 GPI f107= 90.566 f107a= 82.399 power from bzimf,swvel= 41.081 calccloc: ih,bz,y,of18,24,crad= 1 -0.10 -2.83 -6.97 2.26 14.23 calccloc: ih,bz,y,of18,24,crad= 2 -0.10 -2.83 2.05 2.99 12.05 weimer01: CPS,N ctpoten= 54.1 47.6 50.8 Step 27 of 30 mtime=348 6 54 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 27 GPI f107= 90.566 f107a= 82.399 power from bzimf,swvel= 41.235 calccloc: ih,bz,y,of18,24,crad= 1 -0.14 -2.85 -6.97 2.26 14.23 calccloc: ih,bz,y,of18,24,crad= 2 -0.14 -2.85 2.05 2.99 12.04 weimer01: CPS,N ctpoten= 54.5 48.0 51.3 Step 28 of 30 mtime=348 6 56 0 secs/step (sys) = 0.45 GPI+IMF run: istep= 28 GPI f107= 90.566 f107a= 82.399 power from bzimf,swvel= 41.723 calccloc: ih,bz,y,of18,24,crad= 1 -0.20 -2.86 -6.97 2.26 14.23 calccloc: ih,bz,y,of18,24,crad= 2 -0.20 -2.86 2.05 2.99 12.04 weimer01: CPS,N ctpoten= 55.3 48.7 52.0 Step 29 of 30 mtime=348 6 58 0 secs/step (sys) = 0.45 GPI+IMF run: istep= 29 GPI f107= 90.566 f107a= 82.399 power from bzimf,swvel= 42.545 calccloc: ih,bz,y,of18,24,crad= 1 -0.30 -2.87 -6.97 2.26 14.22 calccloc: ih,bz,y,of18,24,crad= 2 -0.30 -2.87 2.05 2.99 12.03 weimer01: CPS,N ctpoten= 56.5 49.7 53.1 Step 30 of 30 mtime=348 7 0 0 secs/step (sys) = 0.44 GPI+IMF run: istep= 30 GPI f107= 90.567 f107a= 82.399 power from bzimf,swvel= 43.367 Allocated 3d sech field ZG(lon= 76,lat= 36,ilev= 29) Initialized diagnostic secondary history field ZG (ix= 24): 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 calccloc: ih,bz,y,of18,24,crad= 1 -0.39 -2.89 -6.97 2.26 14.22 calccloc: ih,bz,y,of18,24,crad= 2 -0.39 -2.89 2.05 2.99 12.03 weimer01: CPS,N ctpoten= 57.6 50.7 54.2 Allocated 3d sech field WN(lon= 76,lat= 36,lev= 29) Initialized diagnostic secondary history field WN (ix= 4): short_name = WN long_name = NEUTRAL VERTICAL WIND (plus up) units = cm/s geo = T mag = F dimnames = lon lat lev 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= 14): 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= 15): 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= 16): 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= 11): 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 2d sech field HMF2(lon= 76,lat= 36) Initialized diagnostic secondary history field HMF2 (ix= 19): short_name = HMF2 long_name = Height of the F2 Layer units = km geo = T mag = F dimnames = lon lat dimsizes = 76 36 0 ndims = 2 task0_only = T Allocated 2d sech field NMF2(lon= 76,lat= 36) Initialized diagnostic secondary history field NMF2 (ix= 20): short_name = NMF2 long_name = Peak Density of the F2 Layer units = 1/cm3 geo = T mag = F dimnames = lon lat dimsizes = 76 36 0 ndims = 2 task0_only = T Allocated 3d sech field SIGMA_PED(lon= 76,lat= 36,lev= 29) Initialized diagnostic secondary history field SIGMA_PED (ix= 17): short_name = SIGMA_PED long_name = Pedersen Conductivity units = S/m geo = T mag = F dimnames = lon lat lev dimsizes = 76 36 29 ndims = 3 task0_only = T Allocated 3d sech field SIGMA_HAL(lon= 76,lat= 36,lev= 29) Initialized diagnostic secondary history field SIGMA_HAL (ix= 18): short_name = SIGMA_HAL long_name = Hall Conductivity units = S/m geo = T mag = F dimnames = lon lat lev dimsizes = 76 36 29 ndims = 3 task0_only = T Allocated 3d sech field QJOULE(lon= 76,lat= 36,lev= 29) Initialized diagnostic secondary history field QJOULE (ix= 21): short_name = QJOULE long_name = Joule Heating units = erg/g/s geo = T mag = F dimnames = lon lat lev dimsizes = 76 36 29 ndims = 3 task0_only = T Allocated 2d sech field QJOULE_INTEG(lon= 76,lat= 36) Initialized diagnostic secondary history field QJOULE_INTEG (ix= 22): short_name = QJOULE_INTEG long_name = Height-integrated Joule Heating units = erg/cm2/s geo = T mag = F dimnames = lon lat dimsizes = 76 36 0 ndims = 2 task0_only = T Allocated 3d sech field DEN(lon= 76,lat= 36,ilev= 29) Initialized diagnostic secondary history field DEN (ix= 7): short_name = DEN long_name = Total Density units = g/cm3 geo = T mag = F dimnames = lon lat ilev dimsizes = 76 36 29 ndims = 3 task0_only = T Allocated 3d sech field ZMAG(mlon= 81,mlat= 97,imlev= 32) Initialized diagnostic secondary history field ZMAG (ix= 27): 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 ------------------------------------------------------------------------ Write TGCM PRIMARY HISTORY Diskfile: p_06348_testw01.nc label = tiegcm_sres_w01 model_name = tiegcm model_version = tiegcm1.94 create_date= write_date = 01/05/12 11:17:02 logname = emery host = puffin system = LINUX run_type = continuation source_file = p_06348_testw01.nc output_file = p_06348_testw01.nc source_mtime = 348 6 0 initial_year = 0 initial_day = 0 initial_mtime= 0 0 0 type = primary ihist = 8 (nth history on history file) delhmins= 60 (delta minutes between histories) calendar year,day = 2006,348 (model IS being advanced in calendar time) modeltime = 348, 7, 0, 0 (model time day,hour,minute,seconds) time = 420.00 (minutes in current day) ut = 7.00 (ut hours) step = 120 (time step in seconds) iter = 250770 (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.9057E+02 (daily solar flux) f107a = 0.8240E+02 (average solar flux) hpower = 0.4337E+02 (Gw) ctpoten = 0.5416E+02 (Volts) kp = 0.1000E+37 () byimf = -0.2888E+01 bzimf = -0.3879E+00 swvel = 0.5815E+03 swden = 0.1454E+01 al = 0.1000E+37 e1,e2 = 0.2735E+01 0.5770E+01 (ergs/cm2/s) alfac,alfad = 0.1000E+00 0.5000E+00 (keV) ec,ed = 0.1061E+00 0.2722E-01 (ergs/cm2/s) colfac = 0.15E+01 joulefac = 0.15E+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-2011212.nc ncep_ncfile = [none] imf_ncfile = /hao/aim1/tgcm/data/imf_OMNI_2006001-2006365.nc 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 = 1 (number of MPI tasks) coupled_cmit = 0 (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, 7, 0 to p_06348_testw01.nc ( 1 of 1) >>> WARNING: Field ALFA will not be defined on secondary histories because data was not allocated. (addfld was probably not called for this field) >>> WARNING: Field EFLUX will not be defined on secondary histories because data was not allocated. (addfld was probably not called for this field) Wrote secondary history 348, 7, 0 to s_06348_testw01_002.nc ( 2 of 2) MPI run with ntask = 1 nstep= 30 step= 120 Model simulation time = 3600 secs (minutes= 60.00, hours= 1.00, days= 0.041667) Cpu time for run = 16.57 ------------------------------------------------------------------------ TIMER (system_clock): Elapsed run time = 17.54 (secs) 0.29 (mins) 0.00 (hrs) Elapsed init time = 1.06 secs, 0.000 hrs, 6.1% (includes source i/o) Elapsed step time = 15.55 secs, 0.004 hrs, 88.7% Elapsed prep time = 2.32 secs, 0.001 hrs, 13.2% Elapsed dynamics = 11.15 secs, 0.003 hrs, 63.5% Elapsed qrj = 1.96 secs, 0.001 hrs, 11.2% Elapsed oplus = 1.00 secs, 0.000 hrs, 5.7% Elapsed cmpminor = 1.02 secs, 0.000 hrs, 5.8% (n4s,n2d,no) Elapsed cmpmajor = 0.67 secs, 0.000 hrs, 3.8% (o2, o) Elapsed dt = 0.59 secs, 0.000 hrs, 3.3% Elapsed duv = 0.48 secs, 0.000 hrs, 2.7% Elapsed dynamo = 2.08 secs, 0.001 hrs, 11.9% Elapsed i/o = 2.35 secs, 0.001 hrs, 13.4% Elapsed Primary io = 1.42 secs, 0.000 hrs, 8.1% Elapsed Secondary = 0.92 secs, 0.000 hrs, 5.3% Average secs/step = 0.52 Average mins/simulated day = 7.02 ------------------------------------------------------------------------ End execution of tiegcm1.94 at 01/05/12 11:17:03 NORMAL EXIT