====================================================================== File dres_s_06347_351_013.nc netcdf dres_s_06347_351_013 { dimensions: time = UNLIMITED ; // (24 currently) lon = 144 ; lat = 72 ; lev = 57 ; ilev = 57 ; mlon = 81 ; mlat = 97 ; mlev = 60 ; imlev = 60 ; mtimedim = 3 ; latlon = 2 ; dtidedim = 2 ; sdtidedim = 10 ; datelen = 24 ; filelen = 1024 ; variables: double time(time) ; time:long_name = "time" ; time:units = "minutes since 2006-12-13 0:0:0" ; time:initial_year = 2006 ; time:initial_day = 347 ; time:initial_mtime = 347, 0, 0 ; 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" ; int coupled_cmit(time) ; coupled_cmit:long_name = "1 if coupled with CISM/CMIT, 0 otherwise" ; 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 WN(time, lev, lat, lon) ; WN:long_name = "NEUTRAL VERTICAL WIND (plus up)" ; WN:units = "cm/s" ; WN: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 DEN(time, ilev, lat, lon) ; DEN:long_name = "Total Density" ; DEN:units = "g/cm3" ; DEN: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 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 SIGMA_PED(time, lev, lat, lon) ; SIGMA_PED:long_name = "Pedersen Conductivity" ; SIGMA_PED:units = "S/m" ; SIGMA_PED:missing_value = 1.e+36 ; float SIGMA_HAL(time, lev, lat, lon) ; SIGMA_HAL:long_name = "Hall Conductivity" ; SIGMA_HAL:units = "S/m" ; SIGMA_HAL:missing_value = 1.e+36 ; float HMF2(time, lat, lon) ; HMF2:long_name = "Height of the F2 Layer" ; HMF2:units = "km" ; HMF2:missing_value = 1.e+36 ; float NMF2(time, lat, lon) ; NMF2:long_name = "Peak Density of the F2 Layer" ; NMF2:units = "1/cm3" ; NMF2: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 QJOULE_INTEG(time, lat, lon) ; QJOULE_INTEG:long_name = "Height-integrated Joule Heating" ; QJOULE_INTEG:units = "erg/cm2/s" ; QJOULE_INTEG: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 = "Geometric Height ZG" ; ZG:units = "cm" ; ZG:missing_value = 1.e+36 ; float ALFA(time, lat, lon) ; ALFA:long_name = "2D ALFA" ; ALFA:units = "keV" ; ALFA:missing_value = 1.e+36 ; float EFLUX(time, lat, lon) ; EFLUX:long_name = "2D EFLUX" ; EFLUX:units = "mW/m^2" ; EFLUX: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 = "tiegcm_dres" ; :create_date = "10/15/11 18:33:07" ; :logname = "emery" ; :host = "puffin" ; :system = "LINUX" ; :model_name = "tiegcm" ; :model_version = "tiegcm1.94" ; :svn_version = "698" ; :output_file = "dres_s_06347_351_013.nc" ; :history_type = "secondary" ; :run_type = "initial" ; :source_file = "/hao/aim3/emery/tiegcm_dres/tiegcm-linux/dres_p_06328_347_002.nc (initial)" ; :source_mtime = 347, 0, 0 ; :initial_file = "/hao/aim3/emery/tiegcm_dres/tiegcm-linux/dres_p_06328_347_002.nc" ; :initial_mtime = 347, 0, 0 ; :lev_to_hPa_method1 = "p0*exp(-lev(k))" ; :lev_to_hPa_method2 = "p0_model*1.e-3*exp(-lev(k))" ; :nhist = 24 ; :delhist_mins = 15 ; :missing_value = 1.e+36 ; :potential_model = "WEIMER" ; :tuv_lbc_intop = 0 ; :contents = "06350 350 0 0 to 06350 350 5 45 by 15 tiegcm1.94 secondary" ; :contents_desc = "yyddd day hour min to yyddd day hour min by delta_mins" ; data: time = 4320, 4335, 4350, 4365, 4380, 4395, 4410, 4425, 4440, 4455, 4470, 4485, 4500, 4515, 4530, 4545, 4560, 4575, 4590, 4605, 4620, 4635, 4650, 4665 ; lon = -180, -177.5, -175, -172.5, -170, -167.5, -165, -162.5, -160, -157.5, -155, -152.5, -150, -147.5, -145, -142.5, -140, -137.5, -135, -132.5, -130, -127.5, -125, -122.5, -120, -117.5, -115, -112.5, -110, -107.5, -105, -102.5, -100, -97.5, -95, -92.5, -90, -87.5, -85, -82.5, -80, -77.5, -75, -72.5, -70, -67.5, -65, -62.5, -60, -57.5, -55, -52.5, -50, -47.5, -45, -42.5, -40, -37.5, -35, -32.5, -30, -27.5, -25, -22.5, -20, -17.5, -15, -12.5, -10, -7.5, -5, -2.5, 0, 2.5, 5, 7.5, 10, 12.5, 15, 17.5, 20, 22.5, 25, 27.5, 30, 32.5, 35, 37.5, 40, 42.5, 45, 47.5, 50, 52.5, 55, 57.5, 60, 62.5, 65, 67.5, 70, 72.5, 75, 77.5, 80, 82.5, 85, 87.5, 90, 92.5, 95, 97.5, 100, 102.5, 105, 107.5, 110, 112.5, 115, 117.5, 120, 122.5, 125, 127.5, 130, 132.5, 135, 137.5, 140, 142.5, 145, 147.5, 150, 152.5, 155, 157.5, 160, 162.5, 165, 167.5, 170, 172.5, 175, 177.5 ; lat = -88.75, -86.25, -83.75, -81.25, -78.75, -76.25, -73.75, -71.25, -68.75, -66.25, -63.75, -61.25, -58.75, -56.25, -53.75, -51.25, -48.75, -46.25, -43.75, -41.25, -38.75, -36.25, -33.75, -31.25, -28.75, -26.25, -23.75, -21.25, -18.75, -16.25, -13.75, -11.25, -8.75, -6.25, -3.75, -1.25, 1.25, 3.75, 6.25, 8.75, 11.25, 13.75, 16.25, 18.75, 21.25, 23.75, 26.25, 28.75, 31.25, 33.75, 36.25, 38.75, 41.25, 43.75, 46.25, 48.75, 51.25, 53.75, 56.25, 58.75, 61.25, 63.75, 66.25, 68.75, 71.25, 73.75, 76.25, 78.75, 81.25, 83.75, 86.25, 88.75 ; lev = -6.875, -6.625, -6.375, -6.125, -5.875, -5.625, -5.375, -5.125, -4.875, -4.625, -4.375, -4.125, -3.875, -3.625, -3.375, -3.125, -2.875, -2.625, -2.375, -2.125, -1.875, -1.625, -1.375, -1.125, -0.875, -0.625, -0.375, -0.125, 0.125, 0.375, 0.625, 0.875, 1.125, 1.375, 1.625, 1.875, 2.125, 2.375, 2.625, 2.875, 3.125, 3.375, 3.625, 3.875, 4.125, 4.375, 4.625, 4.875, 5.125, 5.375, 5.625, 5.875, 6.125, 6.375, 6.625, 6.875, 7.125 ; ilev = -7, -6.75, -6.5, -6.25, -6, -5.75, -5.5, -5.25, -5, -4.75, -4.5, -4.25, -4, -3.75, -3.5, -3.25, -3, -2.75, -2.5, -2.25, -2, -1.75, -1.5, -1.25, -1, -0.75, -0.5, -0.25, 0, 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 ; 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 = -7.625, -7.375, -7.125, -6.875, -6.625, -6.375, -6.125, -5.875, -5.625, -5.375, -5.125, -4.875, -4.625, -4.375, -4.125, -3.875, -3.625, -3.375, -3.125, -2.875, -2.625, -2.375, -2.125, -1.875, -1.625, -1.375, -1.125, -0.875, -0.625, -0.375, -0.125, 0.125, 0.375, 0.625, 0.875, 1.125, 1.375, 1.625, 1.875, 2.125, 2.375, 2.625, 2.875, 3.125, 3.375, 3.625, 3.875, 4.125, 4.375, 4.625, 4.875, 5.125, 5.375, 5.625, 5.875, 6.125, 6.375, 6.625, 6.875, 7.125 ; imlev = -7.75, -7.5, -7.25, -7, -6.75, -6.5, -6.25, -6, -5.75, -5.5, -5.25, -5, -4.75, -4.5, -4.25, -4, -3.75, -3.5, -3.25, -3, -2.75, -2.5, -2.25, -2, -1.75, -1.5, -1.25, -1, -0.75, -0.5, -0.25, 0, 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 ; mtime = 350, 0, 0, 350, 0, 15, 350, 0, 30, 350, 0, 45, 350, 1, 0, 350, 1, 15, 350, 1, 30, 350, 1, 45, 350, 2, 0, 350, 2, 15, 350, 2, 30, 350, 2, 45, 350, 3, 0, 350, 3, 15, 350, 3, 30, 350, 3, 45, 350, 4, 0, 350, 4, 15, 350, 4, 30, 350, 4, 45, 350, 5, 0, 350, 5, 15, 350, 5, 30, 350, 5, 45 ; 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 ; day = 350, 350, 350, 350, 350, 350, 350, 350, 350, 350, 350, 350, 350, 350, 350, 350, 350, 350, 350, 350, 350, 350, 350, 350 ; 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 ; write_date = "10/15/11 18:33:07 ", "10/15/11 18:33:33 ", "10/15/11 18:34:07 ", "10/15/11 18:34:33 ", "10/15/11 18:34:59 ", "10/15/11 18:35:26 ", "10/15/11 18:35:52 ", "10/15/11 18:36:18 ", "10/15/11 18:36:44 ", "10/15/11 18:37:10 ", "10/15/11 18:37:37 ", "10/15/11 18:38:03 ", "10/15/11 18:38:29 ", "10/15/11 18:38:55 ", "10/15/11 18:39:21 ", "10/15/11 18:39:48 ", "10/15/11 18:40:14 ", "10/15/11 18:40:40 ", "10/15/11 18:41:07 ", "10/15/11 18:41:33 ", "10/15/11 18:42:00 ", "10/15/11 18:42:26 ", "10/15/11 18:42:52 ", "10/15/11 18:43:18 " ; iter = 403200, 403212, 403224, 403236, 403248, 403260, 403272, 403284, 403296, 403308, 403320, 403332, 403344, 403356, 403368, 403380, 403392, 403404, 403416, 403428, 403440, 403452, 403464, 403476 ; 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 ; coupled_cmit = 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ; ut = 0, 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 ; timestep = 75 ; f107d = 81.73125, 81.6814321824062, 81.6318936571663, 81.5826370239258, 81.5336648823302, 81.4849798320252, 81.4365844726562, 81.3884814038689, 81.3406732253086, 81.2931625366211, 81.2459519374518, 81.1990440274462, 81.15244140625, 81.1061466735086, 81.0601624288677, 81.0144912719727, 80.9691358024691, 80.9240986200027, 80.8793823242188, 80.834989514763, 80.7909227912808, 80.747184753418, 80.7037780008198, 80.660705133132 ; f107a = 82.2662808641975, 82.2655455968526, 82.2648112072194, 82.264077702275, 82.2633450889965, 82.2626133743609, 82.2618825653453, 82.2611526689266, 82.260423692082, 82.2596956417884, 82.258968525023, 82.2582423487627, 82.2575171199846, 82.2567928456657, 82.2560695327831, 82.2553471883138, 82.2546258192349, 82.2539054325233, 82.2531860351562, 82.2524676341107, 82.2517502363635, 82.251033848892, 82.2503184786731, 82.2496041326838 ; hpower = 60.2098945412675, 65.4470737653134, 67.5950168390684, 68.0714939944829, 64.5297277932174, 58.3815804648237, 50.9216522651172, 50.3811806061716, 63.9465955959804, 67.2458623125667, 62.079577205206, 48.0715147976467, 57.8455846070054, 48.7139134917093, 48.3002301859591, 43.3901553213673, 42.5493385425949, 42.0677937781801, 42.6417432722056, 44.7193176280679, 45.0940229987272, 44.8835425500543, 44.0529894697593, 44.1458668502527 ; ctpoten = 56.1378103816358, 65.3918564752705, 69.1376992707381, 69.948598681115, 62.483406048654, 50.0891813003721, 45.4164555259802, 40.8789932702744, 60.211636053122, 67.3896768135099, 60.0649114369179, 43.7403250677272, 50.6360133454367, 35.8125912624081, 34.9193384581299, 37.4406499412019, 28.3680952082371, 27.4099813229253, 25.5418248382935, 31.2937041988118, 36.4104492082635, 35.1704383372712, 31.3722726301575, 31.654824984735 ; 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 ; bximf = -0.699333333333333, -0.0993333349586367, 0.571333340328512, 1.07633333333333, 1.22333333774715, -1.11666666849912, -1.55666666666667, -0.648000015010155, -0.182666663639141, 0.896666666666667, 1.73466666222098, -2.13600000683583, -1.73333333333333, -1.36133333403444, 0.100666670204092, -1.422, -2.44266666367101, -1.89999998878222, -0.435333333333333, 1.34266665732914, 2.01533333754, 2.452, 2.63999999759391, 2.88333333532513 ; byimf = -0.827568135095073, -0.144719886067381, -0.109263315775826, 0.178345051906645, 0.179106368308063, -0.482307484285158, -1.20779399559309, -0.804812852028291, -0.325161056405895, 0.831176915621073, 1.27626842234044, -1.28844728148646, 0.810323563568192, -0.0355650328861249, -0.351726978727352, -1.79259042056314, -0.848292763611695, -0.600000204370789, 0.0848024070909798, 0.754686093522636, 0.791397010919975, 0.503453817737722, -0.202414986624121, 0.0568189190864244 ; bzimf = -1.51248235492771, -2.11169473480953, -2.3434397447073, -2.41898787304159, -2.04763244926825, -1.4314014556434, -0.605712791408554, -0.573976807387939, -2.07365437059926, -2.39163211773779, -1.81531003055016, -0.341345243999922, -1.50344427297334, -0.507408054772025, -0.365888106277706, 0.742677187844953, 1.48658037995051, 1.53852638031442, 2.16581614798204, 1.55737073872246, 0.977583340275638, 0.944496381252146, 0.915998061772819, 0.977384344648099 ; swvel = 650.470301198142, 650.668584892827, 651.815622487266, 650.229662012782, 647.464732594287, 640.131431045276, 636.992247522416, 634.412631112735, 640.268777270346, 645.091052651112, 643.564720986726, 631.39907204594, 629.504981419898, 623.79250771284, 631.574536099071, 633.901553213673, 625.493385425949, 620.677937781801, 626.417432722056, 647.193176280679, 650.940229987272, 648.835425500543, 640.529894697593, 641.458668502527 ; swden = 1.29066666666667, 1.59900000118233, 1.36633333267737, 1.33433333333333, 1.43233333124593, 1.66999999956977, 1.682, 2.00666666641172, 1.75333333296684, 1.38333333333333, 1.35599999797634, 1.18199999789667, 1.12266666666667, 1.60433332515901, 1.91599999663785, 1.25933333333333, 1.1833333317399, 1.42766666856285, 1.582, 1.51816666827205, 1.46716666719648, 1.34333333333333, 1.87239999431781, 1.91893333412368 ; 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 ; 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 ; 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 ; 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 ; 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] ", " " ; gpi_ncfile = "/hao/aim3/tgcm/data/gpi_2000001-2011212.nc ", " ", "/hao/aim3/tgcm/data/gpi_2000001-2011212.nc ", " ", "/hao/aim3/tgcm/data/gpi_2000001-2011212.nc ", " ", "/hao/aim3/tgcm/data/gpi_2000001-2011212.nc ", " ", "/hao/aim3/tgcm/data/gpi_2000001-2011212.nc ", " ", "/hao/aim3/tgcm/data/gpi_2000001-2011212.nc ", " ", "/hao/aim3/tgcm/data/gpi_2000001-2011212.nc ", " ", "/hao/aim3/tgcm/data/gpi_2000001-2011212.nc ", " ", "/hao/aim3/tgcm/data/gpi_2000001-2011212.nc ", " ", "/hao/aim3/tgcm/data/gpi_2000001-2011212.nc ", " ", "/hao/aim3/tgcm/data/gpi_2000001-2011212.nc ", " ", "/hao/aim3/tgcm/data/gpi_2000001-2011212.nc ", " ", "/hao/aim3/tgcm/data/gpi_2000001-2011212.nc ", " ", "/hao/aim3/tgcm/data/gpi_2000001-2011212.nc ", " ", "/hao/aim3/tgcm/data/gpi_2000001-2011212.nc ", " ", "/hao/aim3/tgcm/data/gpi_2000001-2011212.nc ", " ", "/hao/aim3/tgcm/data/gpi_2000001-2011212.nc ", " ", "/hao/aim3/tgcm/data/gpi_2000001-2011212.nc ", " ", "/hao/aim3/tgcm/data/gpi_2000001-2011212.nc ", " ", "/hao/aim3/tgcm/data/gpi_2000001-2011212.nc ", " ", "/hao/aim3/tgcm/data/gpi_2000001-2011212.nc ", " ", "/hao/aim3/tgcm/data/gpi_2000001-2011212.nc ", " ", "/hao/aim3/tgcm/data/gpi_2000001-2011212.nc ", " ", "/hao/aim3/tgcm/data/gpi_2000001-2011212.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] ", " " ; 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] ", " " ; 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] ", " " ; imf_ncfile = "/hao/aim3/tgcm/data/imf_OMNI_2006001-2006365.nc ", " ", "/hao/aim3/tgcm/data/imf_OMNI_2006001-2006365.nc ", " ", "/hao/aim3/tgcm/data/imf_OMNI_2006001-2006365.nc ", " ", "/hao/aim3/tgcm/data/imf_OMNI_2006001-2006365.nc ", " ", "/hao/aim3/tgcm/data/imf_OMNI_2006001-2006365.nc ", " ", "/hao/aim3/tgcm/data/imf_OMNI_2006001-2006365.nc ", " ", "/hao/aim3/tgcm/data/imf_OMNI_2006001-2006365.nc ", " ", "/hao/aim3/tgcm/data/imf_OMNI_2006001-2006365.nc ", " ", "/hao/aim3/tgcm/data/imf_OMNI_2006001-2006365.nc ", " ", "/hao/aim3/tgcm/data/imf_OMNI_2006001-2006365.nc ", " ", "/hao/aim3/tgcm/data/imf_OMNI_2006001-2006365.nc ", " ", "/hao/aim3/tgcm/data/imf_OMNI_2006001-2006365.nc ", " ", "/hao/aim3/tgcm/data/imf_OMNI_2006001-2006365.nc ", " ", "/hao/aim3/tgcm/data/imf_OMNI_2006001-2006365.nc ", " ", "/hao/aim3/tgcm/data/imf_OMNI_2006001-2006365.nc ", " ", "/hao/aim3/tgcm/data/imf_OMNI_2006001-2006365.nc ", " ", "/hao/aim3/tgcm/data/imf_OMNI_2006001-2006365.nc ", " ", "/hao/aim3/tgcm/data/imf_OMNI_2006001-2006365.nc ", " ", "/hao/aim3/tgcm/data/imf_OMNI_2006001-2006365.nc ", " ", "/hao/aim3/tgcm/data/imf_OMNI_2006001-2006365.nc ", " ", "/hao/aim3/tgcm/data/imf_OMNI_2006001-2006365.nc ", " ", "/hao/aim3/tgcm/data/imf_OMNI_2006001-2006365.nc ", " ", "/hao/aim3/tgcm/data/imf_OMNI_2006001-2006365.nc ", " ", "/hao/aim3/tgcm/data/imf_OMNI_2006001-2006365.nc ", " " ; gswm_mi_di_ncfile = "/hao/aim3/tgcm/data/gswm_diurn_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_diurn_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_diurn_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_diurn_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_diurn_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_diurn_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_diurn_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_diurn_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_diurn_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_diurn_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_diurn_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_diurn_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_diurn_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_diurn_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_diurn_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_diurn_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_diurn_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_diurn_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_diurn_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_diurn_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_diurn_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_diurn_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_diurn_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_diurn_2.5d_99km.nc ", " " ; gswm_mi_sdi_ncfile = "/hao/aim3/tgcm/data/gswm_semi_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_semi_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_semi_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_semi_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_semi_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_semi_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_semi_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_semi_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_semi_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_semi_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_semi_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_semi_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_semi_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_semi_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_semi_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_semi_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_semi_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_semi_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_semi_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_semi_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_semi_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_semi_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_semi_2.5d_99km.nc ", " ", "/hao/aim3/tgcm/data/gswm_semi_2.5d_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] ", " " ; 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] ", " " ; e1 = 3.15997679580763, 3.26805310097088, 3.30989770050071, 3.31899974982491, 3.24976188847556, 3.12001916637422, 2.94286747731245, 2.92904062171752, 3.23799897574322, 3.30318704347345, 3.19960283330466, 2.86823143472595, 3.1080675997017, 2.88543304175975, 2.87438199175173, 2.73546737646687, 2.71011081222888, 2.69536220459187, 2.71292185988873, 2.77456553756232, 2.78537787419422, 2.77931545049521, 2.75511252858558, 2.75784160148178 ; e2 = 7.62308839953942, 8.19917811418448, 8.43545185229753, 8.48786433939312, 8.09827005725392, 7.4219738511306, 6.60138174916289, 6.54192986667887, 8.03412551555785, 8.39704485438234, 7.82875349257266, 6.28786662774113, 7.3630143067706, 6.35853048408802, 6.3130253204555, 5.7729170853504, 5.68042723968544, 5.62745731559982, 5.69059175994262, 5.91912493908747, 5.96034252985999, 5.93718968050597, 5.84582884167352, 5.8560453535278 ; 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 ; h2 = 4.00524736353169, 4.13617684413284, 4.18987542097671, 4.20178734986207, 4.11324319483044, 3.95953951162059, 3.83421652265117, 3.82881180606172, 4.09866488989951, 4.18114655781417, 4.05198943013015, 3.80571514797647, 3.94613961517514, 3.81213913491709, 3.80800230185959, 3.75890155321367, 3.75049338542595, 3.7456779377818, 3.75141743272206, 3.77219317628068, 3.77594022998727, 3.77383542550054, 3.76552989469759, 3.76645866850253 ; 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 ; ec = 0.128681258685298, 0.13569907884552, 0.138577322564352, 0.139215801952607, 0.134469835242911, 0.126231317822864, 0.116235014035257, 0.11551078201227, 0.133688438098614, 0.138109455498839, 0.131186633454976, 0.112415829828847, 0.125513083373387, 0.11327664407889, 0.112722308449185, 0.106142808130632, 0.105016113647077, 0.104370843662761, 0.105139935984756, 0.107923885621611, 0.108425990818294, 0.108143947017073, 0.107031005889477, 0.107155461579339 ; 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 ; ed = 0.0373259367247605, 0.040468244259188, 0.041757010103441, 0.0420428963966897, 0.0399178366759305, 0.0362289482788942, 0.0317529913590703, 0.0314287083637029, 0.0395679573575883, 0.04154751738754, 0.0384477463231236, 0.030042908878588, 0.0359073507642033, 0.0304283480950256, 0.0301801381115755, 0.0272340931928203, 0.0267296031255569, 0.0264406762669081, 0.0267850459633234, 0.0280315905768407, 0.0282564137992363, 0.0281301255300326, 0.0276317936818556, 0.0276875201101516 ; crit1 = 21.4636046844552, 21.7472780012191, 21.8303836275615, 21.8393574092667, 21.6695965580509, 21.355047139651, 20.9088064867033, 20.8849953335105, 21.6278320995212, 21.7757980579769, 21.5209222377482, 20.661846381156, 21.2091418697115, 20.3312077199706, 20.3727253382502, 20.2898276219217, 19.790976935416, 19.7125514300724, 19.6933613903089, 20.0399319274801, 20.1906478481581, 20.0946377659555, 19.9987311274618, 19.998498011796 ; crit2 = 36.4636046844552, 36.7472780012191, 36.8303836275615, 36.8393574092667, 36.6695965580509, 36.355047139651, 35.9088064867033, 35.8849953335105, 36.6278320995212, 36.7757980579769, 36.5209222377482, 35.661846381156, 36.2091418697115, 35.3312077199706, 35.3727253382502, 35.2898276219217, 34.790976935416, 34.7125514300724, 34.6933613903089, 35.0399319274801, 35.1906478481581, 35.0946377659555, 34.9987311274618, 34.998498011796 ; mag = _, _, _, _ ; p0 = 5e-07 ; p0_model = 0.0005 ; grav = 870 ; LBC = -7 ; }