====================================================================== File timegcm1.42.scntr_eqnx_smin_s001.nc netcdf timegcm1.42.scntr_eqnx_smin_s001 { dimensions: time = UNLIMITED ; // (96 currently) lon = 72 ; lat = 36 ; lev = 49 ; ilev = 49 ; mlon = 81 ; mlat = 97 ; mlev = 49 ; imlev = 49 ; mtimedim = 3 ; latlon = 2 ; dtidedim = 2 ; sdtidedim = 10 ; datelen = 24 ; filelen = 80 ; variables: double time(time) ; time:long_name = "2D ELECTRIC POTENTIAL" ; time:units = "VOLTS" ; time:initial_year = 2002 ; time:initial_mtime = 80, 0, 0 ; time:initial_day = 80 ; 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 ilev: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 mlev: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 imlev: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" ; 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 byimf(time) ; byimf:long_name = "BY component of IMF" ; double colfac(time) ; colfac:long_name = "ion/neutral collision factor" ; double dtide(time, dtidedim) ; dtide:long_name = "amplitude and phase of diurnal tide" ; double sdtide(time, sdtidedim) ; sdtide:long_name = "amplitude and phase 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 ncep_reanalysis(time, filelen) ; ncep_reanalysis:long_name = "NCEP reanalysis data file" ; ncep_reanalysis:data = "10 mb lower boundary for TN and Z" ; char ecmwf_ncfile(time, filelen) ; ecmwf_ncfile:long_name = "ECMWF ERA40 data file" ; ecmwf_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 = "f107, f107a, ctpoten, power" ; 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 alfa30(time) ; alfa30:long_name = "Characteristic energy of high-energy auroral electrons" ; alfa30:units = "KeV" ; double e30(time) ; e30:long_name = "Energy flux of high-energy auroral electrons" ; e30:units = "ergs/cm2/s" ; double alfad2(time) ; alfad2:long_name = "Characteristic energy of solar protons in the polar cap" ; alfad2:units = "KeV" ; double ed2(time) ; ed2:long_name = "Energy flux of solar protons in the polar cap" ; ed2:units = "ergs/cm2/s" ; 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" ; UN:units = "CM/S" ; UN:missing_value = 1.e+36 ; float VN(time, lev, lat, lon) ; VN:long_name = "NEUTRAL MERIDIONAL WIND" ; 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)" ; 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 H(time, lev, lat, lon) ; H:long_name = "HYDROGEN (H)" ; H:units = "MMR" ; H: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 NE(time, ilev, lat, lon) ; NE:long_name = "ELECTRON DENSITY" ; NE:units = "CM^3" ; NE: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_VEL(time, lev, lat, lon) ; UI_VEL:long_name = "UI" ; UI_VEL:units = "cm/s" ; UI_VEL:missing_value = 1.e+36 ; float VI_VEL(time, lev, lat, lon) ; VI_VEL:long_name = "VI" ; VI_VEL:units = "cm/s" ; VI_VEL:missing_value = 1.e+36 ; float WI_VEL(time, lev, lat, lon) ; WI_VEL:long_name = "WI" ; WI_VEL:units = "cm/s" ; WI_VEL:missing_value = 1.e+36 ; float SIGMAPED(time, lev, lat, lon) ; SIGMAPED:long_name = "Pedersen Conductivity" ; SIGMAPED:units = "S/m" ; SIGMAPED:missing_value = 1.e+36 ; float SIGMAHAL(time, lev, lat, lon) ; SIGMAHAL:long_name = "Hall Conductivity" ; SIGMAHAL:units = "S/m" ; SIGMAHAL:missing_value = 1.e+36 ; float SIGMA_PEDM(time, mlev, mlat, mlon) ; SIGMA_PEDM:missing_value = 1.e+36 ; float SIGMA_HALM(time, mlev, mlat, mlon) ; SIGMA_HALM:missing_value = 1.e+36 ; float ZN_ExB(time, mlev, mlat, mlon) ; ZN_ExB:long_name = "Zonal ExB drift" ; ZN_ExB:units = "m/s" ; ZN_ExB:missing_value = 1.e+36 ; float UP_ExB(time, mlev, mlat, mlon) ; UP_ExB:long_name = "Upward ExB drift" ; UP_ExB:units = "m/s" ; UP_ExB:missing_value = 1.e+36 ; float PHIM2D(time, mlat, mlon) ; PHIM2D:long_name = "2D ELECTRIC POTENTIAL" ; PHIM2D:units = "VOLTS" ; PHIM2D:missing_value = 1.e+36 ; float PHIM3D(time, mlev, mlat, mlon) ; PHIM3D:long_name = "ELECTRIC POTENTIAL (MAG)" ; PHIM3D:units = "VOLTS" ; PHIM3D: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 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: :label = "" ; :create_date = "11/03/11 11:40:56" ; :logname = "foster" ; :host = "iris" ; :system = "LINUX" ; :model_name = "time-gcm" ; :model_version = "timegcm_trunk" ; :output_file = "timegcm1.42.scntr_eqnx_smin_s001.nc" ; :history_type = "secondary" ; :run_type = "initial" ; :source_file = "/hao/aim1/tgcm/data/timegcm1.42/TGCM.timegcm1.42.pcntr_eqnx_smin.nc (initial)" ; :source_mtime = 80, 0, 0 ; :initial_file = "" ; :initial_mtime = 80, 0, 0 ; :lev_to_hPa_method1 = "p0*exp(-lev(k))" ; :lev_to_hPa_method2 = "p0_model*1.e-3*exp(-lev(k))" ; :nhist = 96 ; :delhist_mins = 15 ; :missing_value = 1.e+36 ; :potential_model = "HEELIS" ; :tuv_lbc_intop = 0 ; :contents = "02080 80 0 15 to 02081 81 0 0 by 15 timegcm_trunk secondary" ; :contents_desc = "yyddd day hour min to yyddd day hour min by delta_mins" ; data: time = 80.0104166666667, 80.0208333333333, 80.03125, 80.0416666666667, 80.0520833333333, 80.0625, 80.0729166666667, 80.0833333333333, 80.09375, 80.1041666666667, 80.1145833333333, 80.125, 80.1354166666667, 80.1458333333333, 80.15625, 80.1666666666667, 80.1770833333333, 80.1875, 80.1979166666667, 80.2083333333333, 80.21875, 80.2291666666667, 80.2395833333333, 80.25, 80.2604166666667, 80.2708333333333, 80.28125, 80.2916666666667, 80.3020833333333, 80.3125, 80.3229166666667, 80.3333333333333, 80.34375, 80.3541666666667, 80.3645833333333, 80.375, 80.3854166666667, 80.3958333333333, 80.40625, 80.4166666666667, 80.4270833333333, 80.4375, 80.4479166666667, 80.4583333333333, 80.46875, 80.4791666666667, 80.4895833333333, 80.5, 80.5104166666667, 80.5208333333333, 80.53125, 80.5416666666667, 80.5520833333333, 80.5625, 80.5729166666667, 80.5833333333333, 80.59375, 80.6041666666667, 80.6145833333333, 80.625, 80.6354166666667, 80.6458333333333, 80.65625, 80.6666666666667, 80.6770833333333, 80.6875, 80.6979166666667, 80.7083333333333, 80.71875, 80.7291666666667, 80.7395833333333, 80.75, 80.7604166666667, 80.7708333333333, 80.78125, 80.7916666666667, 80.8020833333333, 80.8125, 80.8229166666667, 80.8333333333333, 80.84375, 80.8541666666667, 80.8645833333333, 80.875, 80.8854166666667, 80.8958333333333, 80.90625, 80.9166666666667, 80.9270833333333, 80.9375, 80.9479166666667, 80.9583333333333, 80.96875, 80.9791666666667, 80.9895833333333, 81 ; 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 = -16.75, -16.25, -15.75, -15.25, -14.75, -14.25, -13.75, -13.25, -12.75, -12.25, -11.75, -11.25, -10.75, -10.25, -9.75, -9.25, -8.75, -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 ; ilev = -17, -16.5, -16, -15.5, -15, -14.5, -14, -13.5, -13, -12.5, -12, -11.5, -11, -10.5, -10, -9.5, -9, -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 ; 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.1238292398346, -86.2386359278075, -84.3344382772017, -82.4013318761043, -80.4295344888877, -78.409455209355, -76.331796629339, -74.1876988914796, -71.9689341788881, -69.6681589004986, -67.2792279860339, -64.7975706762728, -62.2206194286502, -59.5482728257033, -56.7833601240544, -53.9320608400711, -51.004220409899, -48.013496592417, -44.977275450794, -41.916313881279, -38.8540980830483, -35.8159497661267, -32.8279553516623, -29.9158266527393, -27.1038148580883, -24.4137888873841, -21.8645573932226, -19.4714697378313, -17.2462861345908, -15.1972697425655, -13.3294281929159, -11.6448184766842, -10.1428243675907, -8.82031760915328, -7.67162661822114, -6.68827292880132, -5.85851729799549, -5.16689309434438, -4.5940468926696, -4.117225213084, -3.71151165744767, -3.35148250470276, -3.01257879054767, -2.67136422806136, -2.30362868850279, -1.87754941072085, -1.32687202032625, -7.72840955341401e-08, 1.32687202032624, 1.87754941072084, 2.30362868850279, 2.67136422806136, 3.01257879054767, 3.35148250470276, 3.71151165744767, 4.117225213084, 4.5940468926696, 5.16689309434438, 5.85851729799548, 6.68827292880132, 7.67162661822113, 8.82031760915327, 10.1428243675907, 11.6448184766842, 13.3294281929159, 15.1972697425655, 17.2462861345908, 19.4714697378312, 21.8645573932226, 24.4137888873841, 27.1038148580883, 29.9158266527392, 32.8279553516623, 35.8159497661267, 38.8540980830483, 41.916313881279, 44.977275450794, 48.013496592417, 51.004220409899, 53.9320608400711, 56.7833601240544, 59.5482728257032, 62.2206194286502, 64.7975706762727, 67.2792279860339, 69.6681589004986, 71.9689341788881, 74.1876988914795, 76.331796629339, 78.409455209355, 80.4295344888877, 82.4013318761043, 84.3344382772017, 86.2386359278074, 88.1238292398346, 90 ; mlev = -16.75, -16.25, -15.75, -15.25, -14.75, -14.25, -13.75, -13.25, -12.75, -12.25, -11.75, -11.25, -10.75, -10.25, -9.75, -9.25, -8.75, -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 = -17, -16.5, -16, -15.5, -15, -14.5, -14, -13.5, -13, -12.5, -12, -11.5, -11, -10.5, -10, -9.5, -9, -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 = 80, 0, 15, 80, 0, 30, 80, 0, 45, 80, 1, 0, 80, 1, 15, 80, 1, 30, 80, 1, 45, 80, 2, 0, 80, 2, 15, 80, 2, 30, 80, 2, 45, 80, 3, 0, 80, 3, 15, 80, 3, 30, 80, 3, 45, 80, 4, 0, 80, 4, 15, 80, 4, 30, 80, 4, 45, 80, 5, 0, 80, 5, 15, 80, 5, 30, 80, 5, 45, 80, 6, 0, 80, 6, 15, 80, 6, 30, 80, 6, 45, 80, 7, 0, 80, 7, 15, 80, 7, 30, 80, 7, 45, 80, 8, 0, 80, 8, 15, 80, 8, 30, 80, 8, 45, 80, 9, 0, 80, 9, 15, 80, 9, 30, 80, 9, 45, 80, 10, 0, 80, 10, 15, 80, 10, 30, 80, 10, 45, 80, 11, 0, 80, 11, 15, 80, 11, 30, 80, 11, 45, 80, 12, 0, 80, 12, 15, 80, 12, 30, 80, 12, 45, 80, 13, 0, 80, 13, 15, 80, 13, 30, 80, 13, 45, 80, 14, 0, 80, 14, 15, 80, 14, 30, 80, 14, 45, 80, 15, 0, 80, 15, 15, 80, 15, 30, 80, 15, 45, 80, 16, 0, 80, 16, 15, 80, 16, 30, 80, 16, 45, 80, 17, 0, 80, 17, 15, 80, 17, 30, 80, 17, 45, 80, 18, 0, 80, 18, 15, 80, 18, 30, 80, 18, 45, 80, 19, 0, 80, 19, 15, 80, 19, 30, 80, 19, 45, 80, 20, 0, 80, 20, 15, 80, 20, 30, 80, 20, 45, 80, 21, 0, 80, 21, 15, 80, 21, 30, 80, 21, 45, 80, 22, 0, 80, 22, 15, 80, 22, 30, 80, 22, 45, 80, 23, 0, 80, 23, 15, 80, 23, 30, 80, 23, 45, 81, 0, 0 ; year = 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002, 2002 ; day = 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 81 ; 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, 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 = "11/03/11 11:40:56 ", "11/03/11 11:41:03 ", "11/03/11 11:41:09 ", "11/03/11 11:41:16 ", "11/03/11 11:41:22 ", "11/03/11 11:41:30 ", "11/03/11 11:41:36 ", "11/03/11 11:41:42 ", "11/03/11 11:41:50 ", "11/03/11 11:41:57 ", "11/03/11 11:42:03 ", "11/03/11 11:42:09 ", "11/03/11 11:42:15 ", "11/03/11 11:42:20 ", "11/03/11 11:42:27 ", "11/03/11 11:42:33 ", "11/03/11 11:42:39 ", "11/03/11 11:42:44 ", "11/03/11 11:42:50 ", "11/03/11 11:42:57 ", "11/03/11 11:43:03 ", "11/03/11 11:43:09 ", "11/03/11 11:43:15 ", "11/03/11 11:43:20 ", "11/03/11 11:43:26 ", "11/03/11 11:43:33 ", "11/03/11 11:43:39 ", "11/03/11 11:43:45 ", "11/03/11 11:43:51 ", "11/03/11 11:43:57 ", "11/03/11 11:44:03 ", "11/03/11 11:44:09 ", "11/03/11 11:44:15 ", "11/03/11 11:44:21 ", "11/03/11 11:44:27 ", "11/03/11 11:44:33 ", "11/03/11 11:44:39 ", "11/03/11 11:44:45 ", "11/03/11 11:44:51 ", "11/03/11 11:44:57 ", "11/03/11 11:45:03 ", "11/03/11 11:45:09 ", "11/03/11 11:45:15 ", "11/03/11 11:45:20 ", "11/03/11 11:45:26 ", "11/03/11 11:45:32 ", "11/03/11 11:45:38 ", "11/03/11 11:45:44 ", "11/03/11 11:45:49 ", "11/03/11 11:45:55 ", "11/03/11 11:46:01 ", "11/03/11 11:46:06 ", "11/03/11 11:46:12 ", "11/03/11 11:46:18 ", "11/03/11 11:46:23 ", "11/03/11 11:46:29 ", "11/03/11 11:46:35 ", "11/03/11 11:46:40 ", "11/03/11 11:46:46 ", "11/03/11 11:46:52 ", "11/03/11 11:46:58 ", "11/03/11 11:47:03 ", "11/03/11 11:47:09 ", "11/03/11 11:47:16 ", "11/03/11 11:47:21 ", "11/03/11 11:47:27 ", "11/03/11 11:47:33 ", "11/03/11 11:47:39 ", "11/03/11 11:47:45 ", "11/03/11 11:47:52 ", "11/03/11 11:47:58 ", "11/03/11 11:48:04 ", "11/03/11 11:48:10 ", "11/03/11 11:48:16 ", "11/03/11 11:48:21 ", "11/03/11 11:48:27 ", "11/03/11 11:48:33 ", "11/03/11 11:48:39 ", "11/03/11 11:48:45 ", "11/03/11 11:48:50 ", "11/03/11 11:48:56 ", "11/03/11 11:49:02 ", "11/03/11 11:49:08 ", "11/03/11 11:49:14 ", "11/03/11 11:49:19 ", "11/03/11 11:49:25 ", "11/03/11 11:49:31 ", "11/03/11 11:49:37 ", "11/03/11 11:49:44 ", "11/03/11 11:49:50 ", "11/03/11 11:49:56 ", "11/03/11 11:50:02 ", "11/03/11 11:50:07 ", "11/03/11 11:50:14 ", "11/03/11 11:50:20 ", "11/03/11 11:50:29 " ; iter = 38405, 38410, 38415, 38420, 38425, 38430, 38435, 38440, 38445, 38450, 38455, 38460, 38465, 38470, 38475, 38480, 38485, 38490, 38495, 38500, 38505, 38510, 38515, 38520, 38525, 38530, 38535, 38540, 38545, 38550, 38555, 38560, 38565, 38570, 38575, 38580, 38585, 38590, 38595, 38600, 38605, 38610, 38615, 38620, 38625, 38630, 38635, 38640, 38645, 38650, 38655, 38660, 38665, 38670, 38675, 38680, 38685, 38690, 38695, 38700, 38705, 38710, 38715, 38720, 38725, 38730, 38735, 38740, 38745, 38750, 38755, 38760, 38765, 38770, 38775, 38780, 38785, 38790, 38795, 38800, 38805, 38810, 38815, 38820, 38825, 38830, 38835, 38840, 38845, 38850, 38855, 38860, 38865, 38870, 38875, 38880 ; ut = 0.25, 0.5, 0.75, 1, 1.25, 1.5, 1.75, 2, 2.25, 2.5, 2.75, 3, 3.25, 3.5, 3.75, 4, 4.25, 4.5, 4.75, 5, 5.25, 5.5, 5.75, 6, 6.25, 6.5, 6.75, 7, 7.25, 7.5, 7.75, 8, 8.25, 8.5, 8.75, 9, 9.25, 9.5, 9.75, 10, 10.25, 10.5, 10.75, 11, 11.25, 11.5, 11.75, 12, 12.25, 12.5, 12.75, 13, 13.25, 13.5, 13.75, 14, 14.25, 14.5, 14.75, 15, 15.25, 15.5, 15.75, 16, 16.25, 16.5, 16.75, 17, 17.25, 17.5, 17.75, 18, 18.25, 18.5, 18.75, 19, 19.25, 19.5, 19.75, 20, 20.25, 20.5, 20.75, 21, 21.25, 21.5, 21.75, 22, 22.25, 22.5, 22.75, 23, 23.25, 23.5, 23.75, 0 ; timestep = 180 ; f107d = 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70 ; f107a = 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70 ; hpower = 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8 ; ctpoten = 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, 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 ; 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, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ; 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, 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, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 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, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 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] ", "[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] " ; ncep_reanalysis = "[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] ", "[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] " ; ecmwf_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] ", "[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 = "[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] ", "[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_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_diurn_5.0d_32km.nc " ; gswm_mi_sdi_ncfile = "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.nc ", "/hao/aim1/tgcm/data/gswm_semi_5.0d_32km.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] ", "[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] ", "[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] " ; alfa30 = 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40 ; e30 = 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05 ; alfad2 = 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10 ; ed2 = 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20, 1e-20 ; mag = 1e+36, 1e+36, 1e+36, 1e+36 ; p0 = 0.0005 ; p0_model = 1e+36 ; grav = 1e+36 ; LBC = -17 ; }