====================================================================== File TGCM.tiegcm1.94.s_dec2006_heelis_gpi_001.nc netcdf TGCM.tiegcm1.94.s_dec2006_heelis_gpi_001 { dimensions: time = UNLIMITED ; // (48 currently) lon = 72 ; lat = 36 ; lev = 29 ; ilev = 29 ; mlon = 81 ; mlat = 97 ; mlev = 32 ; imlev = 32 ; mtimedim = 3 ; latlon = 2 ; dtidedim = 2 ; sdtidedim = 10 ; datelen = 24 ; filelen = 80 ; variables: double time(time) ; time:long_name = "Total Electron Content" ; time:units = "1/cm2" ; time:initial_year = 2006 ; time:initial_day = 330 ; time:initial_mtime = 330, 0, 0 ; time:missing_value = 1.e+36 ; double lon(lon) ; lon:long_name = "geographic longitude (-west, +east)" ; lon:units = "degrees_east" ; double lat(lat) ; lat:long_name = "geographic latitude (-south, +north)" ; lat:units = "degrees_north" ; double lev(lev) ; lev:long_name = "midpoint levels" ; lev:short_name = "ln(p0/p)" ; lev:units = "" ; lev:positive = "up" ; lev:standard_name = "atmosphere_ln_pressure_coordinate" ; lev:formula_terms = "p0: p0 lev: lev" ; lev:formula = "p(k) = p0 * exp(-lev(k))" ; double ilev(ilev) ; ilev:long_name = "interface levels" ; ilev:short_name = "ln(p0/p)" ; ilev:units = "" ; ilev:positive = "up" ; ilev:standard_name = "atmosphere_ln_pressure_coordinate" ; ilev:formula_terms = "p0: p0 lev: ilev" ; ilev:formula = "p(k) = p0 * exp(-ilev(k))" ; double mlon(mlon) ; mlon:long_name = "magnetic longitude (-west, +east)" ; mlon:units = "degrees_east" ; double mlat(mlat) ; mlat:long_name = "magnetic latitude (-south, +north)" ; mlat:units = "degrees_north" ; double mlev(mlev) ; mlev:long_name = "magnetic midpoint levels" ; mlev:short_name = "ln(p0/p)" ; mlev:units = "" ; mlev:positive = "up" ; mlev:standard_name = "atmosphere_ln_pressure_coordinate" ; mlev:formula_terms = "p0: p0 lev: mlev" ; mlev:formula = "p(k) = p0 * exp(-mlev(k))" ; double imlev(imlev) ; imlev:long_name = "magnetic interface levels" ; imlev:short_name = "ln(p0/p)" ; imlev:units = "" ; imlev:positive = "up" ; imlev:standard_name = "atmosphere_ln_pressure_coordinate" ; imlev:formula_terms = "p0: p0 lev: imlev" ; imlev:formula = "p(k) = p0 * exp(-imlev(k))" ; int mtime(time, mtimedim) ; mtime:long_name = "model times (day, hour, minute)" ; mtime:units = "day, hour, minute" ; int year(time) ; year:long_name = "calendar year" ; int day(time) ; day:long_name = "calendar day" ; int calendar_advance(time) ; calendar_advance:long_name = "calendar advance flag (1 if advancing calendar time)" ; char write_date(time, datelen) ; write_date:long_name = "Date and time each history was written" ; int iter(time) ; iter:long_name = "number of time steps from model time 0,0,0" ; int ntask_mpi(time) ; ntask_mpi:long_name = "Number of MPI tasks" ; double ut(time) ; ut:long_name = "universal time (from model time hour and minute)" ; ut:units = "hours" ; int timestep ; timestep:long_name = "timestep" ; timestep:units = "seconds" ; double f107d(time) ; f107d:long_name = "Daily 10.7 cm solar flux" ; f107d:units = "1.e-22 W/m2/Hz" ; double f107a(time) ; f107a:long_name = "81-day average 10.7 cm solar flux" ; f107a:units = "1.e-22 W/m2/Hz" ; double hpower(time) ; hpower:long_name = "Hemispheric power" ; hpower:units = "GW" ; double ctpoten(time) ; ctpoten:long_name = "Cross-tail potential" ; ctpoten:units = "kV" ; double Kp(time) ; Kp:long_name = "Kp index" ; double bximf(time) ; bximf:long_name = "BX component of IMF" ; bximf:units = "nT" ; double byimf(time) ; byimf:long_name = "BY component of IMF" ; byimf:units = "nT" ; double bzimf(time) ; bzimf:long_name = "BZ component of IMF" ; bzimf:units = "nT" ; double swvel(time) ; swvel:long_name = "Solar wind velocity" ; swvel:units = "km/s" ; double swden(time) ; swden:long_name = "Solar wind density" ; swden:units = "cm3" ; double al(time) ; al:long_name = "Lower magnetic auroral activity index" ; al:units = "nT" ; double colfac(time) ; colfac:long_name = "ion/neutral collision factor" ; double dtide(time, dtidedim) ; dtide:long_name = "amplitude and phase of diurnal tide mode (1,1)" ; double sdtide(time, sdtidedim) ; sdtide:long_name = "amplitudes and phases of semi-diurnal tide" ; char ncep_ncfile(time, filelen) ; ncep_ncfile:long_name = "ncep data file" ; ncep_ncfile:data = "10 mb lower boundary for TN and Z" ; char gpi_ncfile(time, filelen) ; gpi_ncfile:long_name = "GeoPhysical Indices data file" ; gpi_ncfile:data = "Kp, f107, f107a, ctpoten, power" ; char saber_ncfile(time, filelen) ; saber_ncfile:long_name = "SABER lbc data file" ; saber_ncfile:data = "T,Z" ; char tidi_ncfile(time, filelen) ; tidi_ncfile:long_name = "TIDI lbc data file" ; tidi_ncfile:data = "U,V" ; char see_ncfile(time, filelen) ; see_ncfile:long_name = "TIMED SEE data file" ; see_ncfile:data = "Solar flux data" ; char imf_ncfile(time, filelen) ; imf_ncfile:long_name = "IMF data file" ; imf_ncfile:data = "Bx,By,Bz,Swvel,Swden,Kp,f107,f107a" ; char gswm_mi_di_ncfile(time, filelen) ; gswm_mi_di_ncfile:long_name = "GSWM migrating diurnal tides data file" ; gswm_mi_di_ncfile:data = "lower boundary perturbations for TN, UN, VN, and Z" ; char gswm_mi_sdi_ncfile(time, filelen) ; gswm_mi_sdi_ncfile:long_name = "GSWM migrating semi-diurnal tides data file" ; gswm_mi_sdi_ncfile:data = "lower boundary perturbations for TN, UN, VN, and Z" ; char gswm_nm_di_ncfile(time, filelen) ; gswm_nm_di_ncfile:long_name = "GSWM non-migrating diurnal tides data file" ; gswm_nm_di_ncfile:data = "lower boundary perturbations for TN, UN, VN, and Z" ; char gswm_nm_sdi_ncfile(time, filelen) ; gswm_nm_sdi_ncfile:long_name = "GSWM non-migrating semi-diurnal tides data file" ; gswm_nm_sdi_ncfile:data = "lower boundary perturbations for TN, UN, VN, and Z" ; double e1(time) ; e1:long_name = "Peak energy flux in noon sector of the aurora" ; e1:units = "ergs/cm2/s" ; double e2(time) ; e2:long_name = "Peak energy flux in midnight sector of the aurora" ; e2:units = "ergs/cm2/s" ; double h1(time) ; h1:long_name = "Gaussian half-width of the noon auroral oval" ; h1:units = "degrees" ; double h2(time) ; h2:long_name = "Gaussian half-width of the midnight auroral oval" ; h2:units = "degrees" ; double alfac(time) ; alfac:long_name = "Characteristic Maxwellian energy of polar cusp electrons" ; alfac:units = "keV" ; double ec(time) ; ec:long_name = "Column energy input of polar cusp electrons" ; ec:units = "ergs/cm**2/s" ; double alfad(time) ; alfad:long_name = "Characteristic Maxwellian energy of drizzle electrons" ; alfad:units = "keV" ; double ed(time) ; ed:long_name = "Column energy input of drizzle electrons" ; ed:units = "ergs/cm**2/s" ; double crit1(time) ; crit1:long_name = "critical cross-over latitude 1" ; crit1:units = "deg" ; double crit2(time) ; crit2:long_name = "critical cross-over latitude 2" ; crit2:units = "deg" ; double mag(latlon, latlon) ; mag:long_name = "lat,lon coordinates of S,N magnetic poles" ; mag:units = "degrees" ; double p0 ; p0:long_name = "Reference pressure" ; p0:units = "millibars" ; double p0_model ; p0_model:long_name = "Reference pressure (as used by the model)" ; p0_model:units = "microbars" ; double grav ; grav:long_name = "gravitational acceleration" ; grav:units = "cm/s" ; grav:info = "constant used in the model, independent of height" ; float TN(time, lev, lat, lon) ; TN:long_name = "NEUTRAL TEMPERATURE" ; TN:units = "K" ; TN:missing_value = 1.e+36 ; float UN(time, lev, lat, lon) ; UN:long_name = "NEUTRAL ZONAL WIND (+EAST)" ; UN:units = "cm/s" ; UN:missing_value = 1.e+36 ; float VN(time, lev, lat, lon) ; VN:long_name = "NEUTRAL MERIDIONAL WIND (+NORTH)" ; VN:units = "cm/s" ; VN:missing_value = 1.e+36 ; float O1(time, lev, lat, lon) ; O1:long_name = "ATOMIC OXYGEN" ; O1:units = "mmr" ; O1:missing_value = 1.e+36 ; float NO(time, lev, lat, lon) ; NO:long_name = "NITRIC OXIDE" ; NO:units = "mmr" ; NO:missing_value = 1.e+36 ; float N4S(time, lev, lat, lon) ; N4S:long_name = "N4S" ; N4S:units = "mmr" ; N4S:missing_value = 1.e+36 ; float NE(time, ilev, lat, lon) ; NE:long_name = "ELECTRON DENSITY" ; NE:units = "cm-3" ; NE:missing_value = 1.e+36 ; float TE(time, lev, lat, lon) ; TE:long_name = "ELECTRON TEMPERATURE" ; TE:units = "K" ; TE:missing_value = 1.e+36 ; float TI(time, lev, lat, lon) ; TI:long_name = "ION TEMPERATURE" ; TI:units = "K" ; TI:missing_value = 1.e+36 ; float TEC(time, lat, lon) ; TEC:long_name = "Total Electron Content" ; TEC:units = "1/cm2" ; TEC:missing_value = 1.e+36 ; float O2(time, lev, lat, lon) ; O2:long_name = "MOLECULAR OXYGEN" ; O2:units = "mmr" ; O2:missing_value = 1.e+36 ; float O2P(time, lev, lat, lon) ; O2P:long_name = "O2+ ION" ; O2P:units = "cm-3" ; O2P:missing_value = 1.e+36 ; float OMEGA(time, ilev, lat, lon) ; OMEGA:long_name = "VERTICAL MOTION" ; OMEGA:units = "s-1" ; OMEGA:former_name = "W" ; OMEGA:missing_value = 1.e+36 ; float POTEN(time, ilev, lat, lon) ; POTEN:long_name = "ELECTRIC POTENTIAL" ; POTEN:units = "volts" ; POTEN:missing_value = 1.e+36 ; float UI_ExB(time, ilev, lat, lon) ; UI_ExB:long_name = "Zonal Ion Drift (ExB)" ; UI_ExB:units = "cm/s" ; UI_ExB:missing_value = 1.e+36 ; float VI_ExB(time, ilev, lat, lon) ; VI_ExB:long_name = "Meridional Ion Drift (ExB)" ; VI_ExB:units = "cm/s" ; VI_ExB:missing_value = 1.e+36 ; float WI_ExB(time, ilev, lat, lon) ; WI_ExB:long_name = "Vertical Ion Drift (ExB)" ; WI_ExB:units = "cm/s" ; WI_ExB:missing_value = 1.e+36 ; float DEN(time, ilev, lat, lon) ; DEN:long_name = "Total Density" ; DEN:units = "g/cm3" ; DEN:missing_value = 1.e+36 ; float QJOULE(time, lev, lat, lon) ; QJOULE:long_name = "Joule Heating" ; QJOULE:units = "erg/g/s" ; QJOULE:missing_value = 1.e+36 ; float Z(time, ilev, lat, lon) ; Z:long_name = "GEOPOTENTIAL HEIGHT" ; Z:units = "cm" ; Z:missing_value = 1.e+36 ; float ZG(time, ilev, lat, lon) ; ZG:long_name = "Geopotential calculated with varying gravity" ; ZG:units = "cm" ; ZG:missing_value = 1.e+36 ; float ZMAG(time, imlev, mlat, mlon) ; ZMAG:long_name = "GEOPOTENTIAL HEIGHT (MAGNETIC)" ; ZMAG:units = "CM" ; ZMAG:missing_value = 1.e+36 ; float TLBC(time, lat, lon) ; TLBC:long_name = "Lower boundary condition of TN" ; TLBC:units = "K" ; TLBC:coordinates = "LBC" ; float ULBC(time, lat, lon) ; ULBC:long_name = "Lower boundary condition of UN" ; ULBC:units = "cm/s" ; ULBC:coordinates = "LBC" ; float VLBC(time, lat, lon) ; VLBC:long_name = "Lower boundary condition of VN" ; VLBC:units = "cm/s" ; VLBC:coordinates = "LBC" ; float TLBC_NM(time, lat, lon) ; TLBC_NM:long_name = "Lower boundary condition of TN (TIME N-1)" ; TLBC_NM:units = "K" ; TLBC_NM:coordinates = "LBC" ; float ULBC_NM(time, lat, lon) ; ULBC_NM:long_name = "Lower boundary condition of UN (TIME N-1)" ; ULBC_NM:units = "cm/s" ; ULBC_NM:coordinates = "LBC" ; float VLBC_NM(time, lat, lon) ; VLBC_NM:long_name = "Lower boundary condition of VN (TIME N-1)" ; VLBC_NM:units = "cm/s" ; VLBC_NM:coordinates = "LBC" ; double LBC ; LBC:long_name = "Interface level of t,u,v lower boundary condition" ; // global attributes: :Conventions = "CF-1.0" ; :label = "tiegcm1.94 dec2006 heelis gpi" ; :create_date = "05/20/11 14:33:00" ; :logname = "foster" ; :host = "arc" ; :system = "LINUX" ; :model_name = "tiegcm" ; :model_version = "tiegcm1.94" ; :svn_version = "600" ; :output_file = "TGCM.tiegcm1.94.s_dec2006_heelis_gpi_001.nc" ; :history_type = "secondary" ; :run_type = "initial" ; :source_file = "/hao/aim1/tgcm/data/tiegcm1.94/dec2006/tiegcm1.94.p_dec2006_heelis_gpi_equil.nc (initial)" ; :source_mtime = 330, 0, 0 ; :initial_file = "/hao/aim1/tgcm/data/tiegcm1.94/dec2006/tiegcm1.94.p_dec2006_heelis_gpi_equil.nc" ; :initial_mtime = 330, 0, 0 ; :lev_to_hPa_method1 = "p0*exp(-lev(k))" ; :lev_to_hPa_method2 = "p0_model*1.e-3*exp(-lev(k))" ; :nhist = 48 ; :delhist_mins = 60 ; :missing_value = 1.e+36 ; :potential_model = "HEELIS" ; :tuv_lbc_intop = 0 ; :contents = "06330 330 1 0 to 06332 332 0 0 by 60 tiegcm1.94 secondary" ; :contents_desc = "yyddd day hour min to yyddd day hour min by delta_mins" ; data: time = 60, 120, 180, 240, 300, 360, 420, 480, 540, 600, 660, 720, 780, 840, 900, 960, 1020, 1080, 1140, 1200, 1260, 1320, 1380, 1440, 1500, 1560, 1620, 1680, 1740, 1800, 1860, 1920, 1980, 2040, 2100, 2160, 2220, 2280, 2340, 2400, 2460, 2520, 2580, 2640, 2700, 2760, 2820, 2880 ; lon = -180, -175, -170, -165, -160, -155, -150, -145, -140, -135, -130, -125, -120, -115, -110, -105, -100, -95, -90, -85, -80, -75, -70, -65, -60, -55, -50, -45, -40, -35, -30, -25, -20, -15, -10, -5, 0, 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85, 90, 95, 100, 105, 110, 115, 120, 125, 130, 135, 140, 145, 150, 155, 160, 165, 170, 175 ; lat = -87.5, -82.5, -77.5, -72.5, -67.5, -62.5, -57.5, -52.5, -47.5, -42.5, -37.5, -32.5, -27.5, -22.5, -17.5, -12.5, -7.5, -2.5, 2.5, 7.5, 12.5, 17.5, 22.5, 27.5, 32.5, 37.5, 42.5, 47.5, 52.5, 57.5, 62.5, 67.5, 72.5, 77.5, 82.5, 87.5 ; lev = -6.75, -6.25, -5.75, -5.25, -4.75, -4.25, -3.75, -3.25, -2.75, -2.25, -1.75, -1.25, -0.75, -0.25, 0.25, 0.75, 1.25, 1.75, 2.25, 2.75, 3.25, 3.75, 4.25, 4.75, 5.25, 5.75, 6.25, 6.75, 7.25 ; ilev = -7, -6.5, -6, -5.5, -5, -4.5, -4, -3.5, -3, -2.5, -2, -1.5, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5, 5, 5.5, 6, 6.5, 7 ; mlon = -180, -175.5, -171, -166.5, -162, -157.5, -153, -148.5, -144, -139.5, -135, -130.5, -126, -121.5, -117, -112.5, -108, -103.5, -99, -94.5, -90, -85.5, -81, -76.5, -72, -67.5, -63, -58.5, -54, -49.5, -45, -40.5, -36, -31.5, -27, -22.5, -18, -13.5, -9, -4.5, 0, 4.5, 9, 13.5, 18, 22.5, 27, 31.5, 36, 40.5, 45, 49.5, 54, 58.5, 63, 67.5, 72, 76.5, 81, 85.5, 90, 94.5, 99, 103.5, 108, 112.5, 117, 121.5, 126, 130.5, 135, 139.5, 144, 148.5, 153, 157.5, 162, 166.5, 171, 175.5, 180 ; mlat = -90, -88.1238292398491, -86.2386359278657, -84.3344382773342, -82.4013318763435, -80.4295344892688, -78.4094552099168, -76.331796630125, -74.1876988925388, -71.9689341802758, -69.6681589022773, -67.2792279882741, -64.7975706790533, -62.2206194320588, -59.5482728298363, -56.7833601290164, -53.9320608459732, -51.0042204168578, -48.0134966005524, -44.9772754602266, -41.916313892128, -38.8540980954293, -35.8159497801506, -32.8279553674349, -29.9158266703621, -27.1038148776609, -24.4137889090065, -21.8645574169981, -19.4714697638694, -17.2462861630082, -15.1972697734841, -13.3294282264571, -11.6448185129562, -10.142824406667, -8.82031765103987, -7.67162666281269, -6.68827297583048, -5.85851734698832, -5.16689314460211, -4.5940469432968, -4.11722526306697, -3.71151170575937, -3.35148255039153, -3.01257883277328, -2.67136426606314, -2.3036287214954, -1.87754943767857, -1.32687203939232, -7.72840966450717e-08, 1.32687203939232, 1.87754943767857, 2.3036287214954, 2.67136426606314, 3.01257883277328, 3.35148255039153, 3.71151170575936, 4.11722526306697, 4.59404694329679, 5.16689314460211, 5.85851734698832, 6.68827297583048, 7.67162666281268, 8.82031765103987, 10.142824406667, 11.6448185129562, 13.3294282264571, 15.1972697734841, 17.2462861630082, 19.4714697638694, 21.8645574169981, 24.4137889090064, 27.1038148776609, 29.9158266703621, 32.8279553674348, 35.8159497801506, 38.8540980954293, 41.916313892128, 44.9772754602266, 48.0134966005524, 51.0042204168578, 53.9320608459731, 56.7833601290163, 59.5482728298363, 62.2206194320588, 64.7975706790533, 67.2792279882741, 69.6681589022773, 71.9689341802758, 74.1876988925387, 76.331796630125, 78.4094552099168, 80.4295344892687, 82.4013318763434, 84.3344382773342, 86.2386359278657, 88.123829239849, 90 ; mlev = -8.25, -7.75, -7.25, -6.75, -6.25, -5.75, -5.25, -4.75, -4.25, -3.75, -3.25, -2.75, -2.25, -1.75, -1.25, -0.75, -0.25, 0.25, 0.75, 1.25, 1.75, 2.25, 2.75, 3.25, 3.75, 4.25, 4.75, 5.25, 5.75, 6.25, 6.75, 7.25 ; imlev = -8.5, -8, -7.5, -7, -6.5, -6, -5.5, -5, -4.5, -4, -3.5, -3, -2.5, -2, -1.5, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5, 5, 5.5, 6, 6.5, 7 ; mtime = 330, 1, 0, 330, 2, 0, 330, 3, 0, 330, 4, 0, 330, 5, 0, 330, 6, 0, 330, 7, 0, 330, 8, 0, 330, 9, 0, 330, 10, 0, 330, 11, 0, 330, 12, 0, 330, 13, 0, 330, 14, 0, 330, 15, 0, 330, 16, 0, 330, 17, 0, 330, 18, 0, 330, 19, 0, 330, 20, 0, 330, 21, 0, 330, 22, 0, 330, 23, 0, 331, 0, 0, 331, 1, 0, 331, 2, 0, 331, 3, 0, 331, 4, 0, 331, 5, 0, 331, 6, 0, 331, 7, 0, 331, 8, 0, 331, 9, 0, 331, 10, 0, 331, 11, 0, 331, 12, 0, 331, 13, 0, 331, 14, 0, 331, 15, 0, 331, 16, 0, 331, 17, 0, 331, 18, 0, 331, 19, 0, 331, 20, 0, 331, 21, 0, 331, 22, 0, 331, 23, 0, 332, 0, 0 ; year = 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006, 2006 ; day = 330, 330, 330, 330, 330, 330, 330, 330, 330, 330, 330, 330, 330, 330, 330, 330, 330, 330, 330, 330, 330, 330, 330, 331, 331, 331, 331, 331, 331, 331, 331, 331, 331, 331, 331, 331, 331, 331, 331, 331, 331, 331, 331, 331, 331, 331, 331, 332 ; calendar_advance = 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ; write_date = "05/20/11 14:33:00 ", "05/20/11 14:33:08 ", "05/20/11 14:33:16 ", "05/20/11 14:33:24 ", "05/20/11 14:33:32 ", "05/20/11 14:33:39 ", "05/20/11 14:33:47 ", "05/20/11 14:33:55 ", "05/20/11 14:34:03 ", "05/20/11 14:34:13 ", "05/20/11 14:34:22 ", "05/20/11 14:34:30 ", "05/20/11 14:34:38 ", "05/20/11 14:34:46 ", "05/20/11 14:34:55 ", "05/20/11 14:35:03 ", "05/20/11 14:35:11 ", "05/20/11 14:35:19 ", "05/20/11 14:35:27 ", "05/20/11 14:35:36 ", "05/20/11 14:35:44 ", "05/20/11 14:35:52 ", "05/20/11 14:36:00 ", "05/20/11 14:36:09 ", "05/20/11 14:36:17 ", "05/20/11 14:36:25 ", "05/20/11 14:36:34 ", "05/20/11 14:36:42 ", "05/20/11 14:36:50 ", "05/20/11 14:36:58 ", "05/20/11 14:37:06 ", "05/20/11 14:37:14 ", "05/20/11 14:37:22 ", "05/20/11 14:37:30 ", "05/20/11 14:37:38 ", "05/20/11 14:37:46 ", "05/20/11 14:37:53 ", "05/20/11 14:38:01 ", "05/20/11 14:38:08 ", "05/20/11 14:38:16 ", "05/20/11 14:38:24 ", "05/20/11 14:38:33 ", "05/20/11 14:38:41 ", "05/20/11 14:38:49 ", "05/20/11 14:38:56 ", "05/20/11 14:39:04 ", "05/20/11 14:39:12 ", "05/20/11 14:39:21 " ; iter = 237630, 237660, 237690, 237720, 237750, 237780, 237810, 237840, 237870, 237900, 237930, 237960, 237990, 238020, 238050, 238080, 238110, 238140, 238170, 238200, 238230, 238260, 238290, 238320, 238350, 238380, 238410, 238440, 238470, 238500, 238530, 238560, 238590, 238620, 238650, 238680, 238710, 238740, 238770, 238800, 238830, 238860, 238890, 238920, 238950, 238980, 239010, 239040 ; ntask_mpi = 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4 ; ut = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 0 ; timestep = 120 ; f107d = 76.1416473765432, 76.1238040123457, 76.109375, 76.0987654320988, 76.0923804012346, 76.090625, 76.0939043209877, 76.1026234567901, 76.1171875, 76.1380015432099, 76.1654706790123, 76.2, 76.3175576292438, 76.4421971450617, 76.57353515625, 76.7111882716049, 76.8547730999229, 77.00390625, 77.1582043306327, 77.3172839506173, 77.48076171875, 77.6482542438272, 77.8193781346451, 77.99375, 78.1709864486883, 78.3507040895062, 78.53251953125, 78.7160493827161, 78.9009102527006, 79.08671875, 79.2730914834105, 79.4596450617284, 79.64599609375, 79.8317611882716, 80.0165569540895, 80.2, 80.362296248071, 80.5226755401235, 80.68095703125, 80.8369598765432, 80.9905032310957, 81.14140625, 81.2894880883488, 81.4345679012346, 81.57646484375, 81.7149980709877, 81.8499867380401, 81.98125 ; f107a = 81.2834120715687, 81.2924337467611, 81.3014443479938, 81.3104404816339, 81.3194187540485, 81.3283757716049, 81.3373081406702, 81.3462124676116, 81.3550853587963, 81.3639234205914, 81.372723259364, 81.3814814814815, 81.3897239755992, 81.397922977633, 81.406080005787, 81.4141965782655, 81.4222742132725, 81.4303144290123, 81.438318743689, 81.4462886755068, 81.4542257426697, 81.4621314633821, 81.470007355848, 81.4778549382716, 81.4856757288571, 81.4934712458086, 81.5012430073302, 81.5089925316263, 81.5167213369008, 81.524430941358, 81.5321228632021, 81.5397986206371, 81.5474597318673, 81.5551077150968, 81.5627440885297, 81.5703703703704, 81.5779966373266, 81.5856157597927, 81.5932291666666, 81.6008382868465, 81.6084445492303, 81.616049382716, 81.6236542162018, 81.6312604785856, 81.6388695987654, 81.6464830056394, 81.6541021281054, 81.6617283950617 ; hpower = 39.0687346829768, 39.7772474249664, 41.2287393549136, 42.7274301813932, 43.4949528720734, 43.4949528720734, 43.4949528720734, 40.9835870655604, 36.3453727556912, 32.1764287050081, 30.2529024657998, 30.2529024657998, 30.2529024657998, 29.147077154032, 27.0388083899093, 25.0612417488077, 26.7003626864626, 32.5735984172433, 39.5398158191161, 42.7274301813932, 41.2287393549136, 39.7772474249664, 37.4589645972178, 34.4142329825722, 31.5885613258783, 29.6955666030302, 28.6072937025108, 27.5532942352594, 27.5532942352594, 28.6072937025108, 29.6955666030302, 31.5885613258783, 34.4142329825722, 37.4589645972178, 35.2610842490009, 28.6072937025108, 23.0569860121119, 21.6070877095079, 23.6589464621026, 25.8698753922253, 28.2522064020339, 30.8192274233656, 33.5852565338579, 31.5885613258783, 25.5438320439339, 20.5015772290971, 18.6768823973199, 19.4422430955406 ; ctpoten = 67.2, 68.192, 70.188, 72.2, 73.212, 73.212, 73.212, 69.8542222222222, 63.272, 56.8675555555555, 53.732, 53.732, 53.732, 51.872, 48.2, 44.592, 47.5942222222222, 57.5, 67.8608888888889, 72.2, 70.188, 68.192, 64.9008888888889, 60.368, 55.9222222222222, 52.8, 50.948, 49.112, 49.112, 50.948, 52.8, 55.9222222222222, 60.368, 64.9008888888889, 61.6542222222222, 50.948, 40.7555555555556, 37.8555555555556, 41.928, 46.0875555555556, 50.3342222222222, 54.668, 59.0888888888889, 55.9222222222222, 45.488, 35.5675555555556, 31.632, 33.308 ; Kp = 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36 ; bximf = 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36 ; byimf = 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ; bzimf = 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36 ; swvel = 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36 ; swden = 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36 ; al = 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36, 1e+36 ; colfac = 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5 ; dtide = 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ; sdtide = 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ; ncep_ncfile = "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] " ; gpi_ncfile = "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc ", "/hao/aim1/tgcm/data/gpi_2000001-2009031.nc " ; saber_ncfile = "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] " ; tidi_ncfile = "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] " ; see_ncfile = "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] " ; imf_ncfile = "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] " ; gswm_mi_di_ncfile = "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_99km.nc " ; gswm_mi_sdi_ncfile = "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_99km.nc " ; gswm_nm_di_ncfile = "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] " ; gswm_nm_sdi_ncfile = "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] ", "[none] " ; e1 = 2.59952492484319, 2.62281375958925, 2.66925587203836, 2.71552310644461, 2.73859326846107, 2.73859326846107, 2.73859326846107, 2.66152785138876, 2.50589608762278, 2.34802522024202, 2.26814948096119, 2.26814948096119, 2.26814948096119, 2.21989714071186, 2.12260658895468, 2.02418969455426, 2.10628464977344, 2.36392201567739, 2.61505591326344, 2.71552310644461, 2.66925587203836, 2.62281375958925, 2.54500223989631, 2.43514972183375, 2.32413185391248, 2.24405490005291, 2.19567484365717, 2.14703101901786, 2.14703101901786, 2.19567484365717, 2.24405490005291, 2.32413185391248, 2.43514972183375, 2.54500223989631, 2.46665022890178, 2.19567484365717, 1.91617997662324, 1.83202112437331, 1.94957601454563, 2.06534000079113, 2.17949010390976, 2.2921822080566, 2.40355416400309, 2.32413185391248, 2.04890498159407, 1.76396625989738, 1.64317787735049, 1.69521934898627 ; e2 = 5.29756081512745, 5.3754972167463, 5.5351613290405, 5.70001731995326, 5.78444481592808, 5.78444481592808, 5.78444481592808, 5.50819457721165, 4.99799100312603, 4.53940715755089, 4.32781927123798, 4.32781927123798, 4.32781927123798, 4.20617848694352, 3.97426892289002, 3.75673659236885, 3.93703989551088, 4.58309582589676, 5.34937974010277, 5.70001731995326, 5.5351613290405, 5.3754972167463, 5.12048610569396, 4.78556562808294, 4.47474174584661, 4.26651232633332, 4.14680230727619, 4.03086236587854, 4.03086236587854, 4.14680230727619, 4.26651232633332, 4.47474174584661, 4.78556562808294, 5.12048610569396, 4.8787192673901, 4.14680230727619, 3.53626846133231, 3.37677964804586, 3.60248411083129, 3.84568629314479, 4.10774270422373, 4.39011501657022, 4.69437821872437, 4.47474174584661, 3.80982152483273, 3.25517349520068, 3.05445706370518, 3.13864674050947 ; h1 = 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35, 2.35 ; h2 = 3.71568734682977, 3.72277247424966, 3.73728739354914, 3.75227430181393, 3.75994952872073, 3.75994952872073, 3.75994952872073, 3.7348358706556, 3.68845372755691, 3.64676428705008, 3.627529024658, 3.627529024658, 3.627529024658, 3.61647077154032, 3.59538808389909, 3.57561241748808, 3.59200362686463, 3.65073598417243, 3.72039815819116, 3.75227430181393, 3.73728739354914, 3.72277247424966, 3.69958964597218, 3.66914232982572, 3.64088561325878, 3.6219556660303, 3.61107293702511, 3.60053294235259, 3.60053294235259, 3.61107293702511, 3.6219556660303, 3.64088561325878, 3.66914232982572, 3.69958964597218, 3.67761084249001, 3.61107293702511, 3.55556986012112, 3.54107087709508, 3.56158946462103, 3.58369875392225, 3.60752206402034, 3.63319227423366, 3.66085256533858, 3.64088561325878, 3.58043832043934, 3.53001577229097, 3.5117688239732, 3.51942243095541 ; alfac = 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1 ; ec = 0.100352104475189, 0.101301511549455, 0.103246510735584, 0.105254756443067, 0.106283236848578, 0.106283236848578, 0.106283236848578, 0.102918006667851, 0.0967027994926262, 0.0911164144647108, 0.0885388893041718, 0.0885388893041718, 0.0885388893041718, 0.0870570833864029, 0.0842320032424784, 0.0815820639434023, 0.0837784859998598, 0.091648621879106, 0.100983353197616, 0.105254756443067, 0.103246510735584, 0.101301511549455, 0.0981950125602719, 0.0941150721966468, 0.0903286721766769, 0.0877920592480605, 0.0863337735613645, 0.0849214142752476, 0.0849214142752476, 0.0863337735613645, 0.0877920592480605, 0.0903286721766769, 0.0941150721966468, 0.0981950125602719, 0.0952498528936612, 0.0863337735613645, 0.07889636125623, 0.0769534975307405, 0.0797029882592175, 0.0826656330255819, 0.0858579565787254, 0.0892977647473099, 0.0930042437553696, 0.0903286721766769, 0.0822287349388715, 0.0754721134869901, 0.0730270224124086, 0.0740526057480245 ; alfad = 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5 ; ed = 0.0246412408097861, 0.0250663484549798, 0.0259372436129482, 0.0268364581088359, 0.0272969717232441, 0.0272969717232441, 0.0272969717232441, 0.0257901522393363, 0.0230072236534147, 0.0205058572230049, 0.0193517414794799, 0.0193517414794799, 0.0193517414794799, 0.0186882462924192, 0.0174232850339456, 0.0162367450492846, 0.0172202176118775, 0.020744159050346, 0.0249238894914696, 0.0268364581088359, 0.0259372436129482, 0.0250663484549798, 0.0236753787583307, 0.0218485397895433, 0.0201531367955269, 0.0190173399618181, 0.0183643762215065, 0.0177319765411557, 0.0177319765411557, 0.0183643762215065, 0.0190173399618181, 0.0201531367955269, 0.0218485397895433, 0.0236753787583307, 0.0223566505494005, 0.0183643762215065, 0.0150341916072672, 0.0141642526257047, 0.0153953678772616, 0.0167219252353352, 0.0181513238412203, 0.0196915364540194, 0.0213511539203147, 0.0201531367955269, 0.0165262992263604, 0.0135009463374583, 0.0124061294383919, 0.0128653458573244 ; crit1 = 19.8649826667366, 19.9163375888544, 20.0178559584696, 20.1178412525175, 20.1672788059236, 20.1672788059236, 20.1672788059236, 20.0010444537101, 19.6553812067271, 19.2897663448114, 19.0984136524971, 19.0984136524971, 19.0984136524971, 18.9805744723556, 18.7374793202015, 18.4834934997635, 18.6959396157452, 19.327318821135, 19.8992641571994, 20.1178412525175, 20.0178559584696, 19.9163375888544, 19.7435486572887, 19.4935134898649, 19.2329978615085, 19.0397889702913, 18.9207537674767, 18.7992244912701, 18.7992244912701, 18.9207537674767, 19.0397889702913, 19.2329978615085, 19.4935134898649, 19.7435486572887, 19.565970468105, 18.9207537674767, 18.194401739401, 17.9608159364553, 18.2850156031872, 18.5907291453687, 18.8805280935904, 19.1564642241002, 19.4202022795288, 19.2329978615085, 18.5480837541961, 17.7660314201362, 17.4057673357028, 17.5634060189633 ; crit2 = 34.8649826667366, 34.9163375888544, 35.0178559584696, 35.1178412525175, 35.1672788059236, 35.1672788059236, 35.1672788059236, 35.0010444537101, 34.6553812067271, 34.2897663448114, 34.0984136524971, 34.0984136524971, 34.0984136524971, 33.9805744723556, 33.7374793202015, 33.4834934997635, 33.6959396157452, 34.327318821135, 34.8992641571994, 35.1178412525175, 35.0178559584696, 34.9163375888544, 34.7435486572887, 34.4935134898649, 34.2329978615085, 34.0397889702913, 33.9207537674767, 33.7992244912701, 33.7992244912701, 33.9207537674767, 34.0397889702913, 34.2329978615085, 34.4935134898649, 34.7435486572887, 34.565970468105, 33.9207537674767, 33.194401739401, 32.9608159364553, 33.2850156031872, 33.5907291453687, 33.8805280935904, 34.1564642241002, 34.4202022795288, 34.2329978615085, 33.5480837541961, 32.7660314201362, 32.4057673357028, 32.5634060189633 ; mag = _, _, _, _ ; p0 = 5e-07 ; p0_model = 0.0005 ; grav = 870 ; LBC = -7 ; }