Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/callphys.def =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/callphys.def (revision 1491) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/callphys.def (revision 1491) @@ -0,0 +1,193 @@ +##General options +##~~~~~~~~~~~~~~~ +#Run with or without tracer transport ? +tracer=.false. + +##Directory where external input files are: +#datadir=/planeto/emlmd/run_mcd5.v0_new/datadir +##datadir = /san/home/spiga/DATAFILE +##datadir = /d5/aslmd/datafile +##datadir=/users/forget/WWW/datagcm/datafile +datadir = ./datadir/ + +#Diurnal cycle ? if diurnal=False, diurnal averaged solar heating +diurnal=.true. + +#Seasonal cycle ? if season=False, Ls stays constant, to value set in "start" +season = .true. + +#write some more output on the screen ? +lwrite = .false. + +#Save statistics in file "stats.nc" ? +callstats =.true. +callstats =.false. + +#Save EOF profiles in file "profiles" for Climate Database? +calleofdump = .false. + +## Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +## ~~~~~~~~~~~~~ +# =1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario, +# =4 Mars Year 24 from TES assimilation (old version of MY24; dust_tes.nc file) +# =6 "cold" (low dust) scenario ; =7 "warm" (high dust) scenario +# =24 Mars Year 24 ; =25 Mars Year 25 (year with a global dust storm) ; ... +# =30 Mars Year 30 +iaervar = 24 +# Dust opacity at 610 Pa (when constant, i.e. for the iaervar=1 case) +tauvis=0.2 +# Dust vertical distribution: +# (=0: old distrib. (Pollack90), =1: top set by "topdustref", +# =2: Viking scenario; =3 MGS scenario) +iddist = 3 +# Dust top altitude (km). (Matters only if iddist=1) +topdustref = 55. + +## Physical Parameterizations : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call radiative transfer ? +callrad = .true. +# call NLTE radiative schemes ? matters only if callrad=T +callnlte = .true. +callnlte = .false. +# NLTE 15um scheme to use. +# 0-> Old scheme, static oxygen +# 1-> Old scheme, dynamic oxygen +# 2-> New scheme +nltemodel = 2 +nltemodel = 0 +# call CO2 NIR absorption ? matters only if callrad=T +callnirco2 = .true. +# NIR NLTE correction ? matters only if callnirco2=T +nircorr=1 +nircorr=0 +# call turbulent vertical diffusion ? +calldifv = .true. +# call convective adjustment ? +calladj = .true. +calladj = .false. +# Thermals +calltherm = .true. +calltherm = .false. +callrichsl = .true. +# call CO2 condensation ? +callcond =.true. +# call thermal conduction in the soil ? +callsoil = .true. +# call Lott's gravity wave/subgrid topography scheme ? +calllott = .true. +calllott = .false. +# Impose polar cap surface albedos as observed by TES? +TESicealbedo = .true. +TESicealbedo = .false. +## Coefficient for Northern cap albedoes +TESice_Ncoef=1.6 +## Coefficient for Southern cap albedoes +TESice_Scoef=1.6 + + +## Radiative transfer options : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# the rad.transfer is computed every "iradia" physical timestep +iradia = 1 +# Output of the exchange coefficient mattrix ? for diagnostic only +callg2d = .false. +# Rayleigh scattering : (should be .false. for now) +rayleigh = .false. + +## Tracer (dust water, ice and/or chemical species) options (used if tracer=T): +## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# DUST: Transported dust ? (if >0, use 'dustbin' dust bins) +dustbin = 2 +dustbin = 0 +# DUST: Radiatively active dust ? (matters if dustbin>0) +active = .true. +active = .false. +# DUST: use mass and number mixing ratios to predict dust size ? +doubleq = .true. +doubleq = .false. +# DUST: lifted by GCM surface winds ? +lifting = .true. +lifting = .false. +# DUST: lifted by dust devils ? +callddevil = .false. +# DUST: Scavenging by H2O/CO2 snowfall ? +scavenging = .true. +scavenging = .false. +# DUST/WATERICE: Gravitationnal sedimentation ? +sedimentation = .true. +sedimentation = .false. +# WATERICE: Radiatively active transported atmospheric water ice ? +activice = .true. +activice = .false. +# WATER: Compute water cycle +water = .true. +water = .false. +# WATER: Microphysical scheme for water-ice clouds? +microphys = .true. +microphys = .false. +# WATER: microphysics sub-timestep +# imicro = 2 5 10 +imicro = 1 +# WATER: parameter contact +mteta = 0.95 +# WATER: Effective variance for sedimentation for the log-normal +# distribution of ice particles ? +nuice_sed=0.1 +# WATER: current permanent caps at both poles. True IS RECOMMENDED +# (with .true., North cap is a source of water and South pole +# is a cold trap) +caps = .true. +caps = .false. +# WATER: Water ice albedo ? +albedo_h2o_ice = 0.35 +# WATER: Water ice thermal inertia +inert_h2o_ice = 800 +# WATER: Frost thickness threshold for albedo +frost_albedo_threshold = 0.005 +# PHOTOCHEMISTRY: include chemical species +photochem = .true. +photochem = .false. + +# SCATTERERS: set number of scatterers. must be compliant with preceding options. +naerkind = 1 + + +## Thermospheric options (relevant if tracer=T) : +##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call thermosphere ? +callthermos = .true. +callthermos = .false. +# WATER: included without cycle (only if water=.false.) +thermoswater = .false. +# call thermal conduction ? (only if callthermos=.true.) +callconduct = .true. +callconduct = .false. +# call molecular viscosity ? (only if callthermos=.true.) +callmolvis = .true. +callmolvis = .false. +# call molecular diffusion ? (only if callthermos=.true.) +callmoldiff = .true. +callmoldiff = .false. +# call thermospheric photochemistry ? (only if callthermos=.true.) +thermochem = .true. +thermochem = .false. +# call EUV heating ? (only if callthermos=.true.) +calleuv=.true. +calleuv=.false. +#Method to include solar variability? +#0-> Old method (fixed EUV input) 1-> Variability with E10.7 as observed +solvarmod=0 +# date for solar flux calculation: (1985 < date < 2002) +## (Solar min=1996.4 ave=1993.4 max=1990.6) ; Only used if solvarmod=0 +solarcondate = 1993.4 +#Solar variability as observed for MY? (must bebetween MY23 and MY30) +# (only matters if solvarmod=1) +solvaryear=24 +# value for the UV heating efficiency +##(experimental values between 0.19 and 0.23, lower values may +## be used to compensate for low 15 um cooling) +euveff = 0.21 + + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/input_coord =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/input_coord (revision 1491) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/input_coord (revision 1491) @@ -0,0 +1,49 @@ +136.0 lon +4.4 lat +231.2 ls +7. lct +20 ztop +100. nz + +326.0 lon +-26.1 lat +348. ls +08. lct +20. ztop +100. nz + + +-195.5 lon +73. lat +90.0 ls +08. lct +60. ztop +100. nz + +204.6 lon +13.7 lat +140.0 ls +08. lct +20. ztop +100. nz + +04.6 lon +13.7 lat +140.0 ls +08. lct +20. ztop +100. nz + +-195.5 lon +73. lat +90.0 ls +08. lct +20. ztop +100. nz + +353.87 lon +-1.88 lat +244. ls +06. lct +20 ztop +100. nz Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/input_dust =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/input_dust (revision 1491) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/input_dust (revision 1491) @@ -0,0 +1,109 @@ + 2.484180e-05 1.091621e+09 + 2.484180e-05 1.091621e+09 + 2.484180e-05 1.091621e+09 + 2.484180e-05 1.091621e+09 + 2.479834e-05 1.090773e+09 + 2.426126e-05 1.080098e+09 + 2.318709e-05 1.057635e+09 + 2.089747e-05 1.006888e+09 + 1.800252e-05 9.381032e+08 + 1.397773e-05 8.315937e+08 + 1.172798e-05 7.742898e+08 + 1.107991e-05 7.604942e+08 + 1.057517e-05 7.501051e+08 + 1.060382e-05 7.540547e+08 + 1.063248e-05 7.580213e+08 + 1.066774e-05 7.620605e+08 + 1.085635e-05 7.673091e+08 + 1.104496e-05 7.723900e+08 + 1.123356e-05 7.773068e+08 + 1.142217e-05 7.820637e+08 + 1.158144e-05 7.856579e+08 + 1.171280e-05 7.881999e+08 + 1.184415e-05 7.906679e+08 + 1.197550e-05 7.930630e+08 + 1.210686e-05 7.953867e+08 + 1.223821e-05 7.976401e+08 + 1.236956e-05 7.998245e+08 + 1.250092e-05 8.019414e+08 + 1.260200e-05 8.045168e+08 + 1.268613e-05 8.073556e+08 + 1.277025e-05 8.101715e+08 + 1.285437e-05 8.129645e+08 + 1.293849e-05 8.157349e+08 + 1.302261e-05 8.184828e+08 + 1.310673e-05 8.212083e+08 + 1.319086e-05 8.239118e+08 + 1.327498e-05 8.265931e+08 + 1.336635e-05 8.284246e+08 + 1.346006e-05 8.299587e+08 + 1.355377e-05 8.314628e+08 + 1.364748e-05 8.329368e+08 + 1.374119e-05 8.343816e+08 + 1.383490e-05 8.357973e+08 + 1.392861e-05 8.371840e+08 + 1.402232e-05 8.385426e+08 + 1.411604e-05 8.398728e+08 + 1.420975e-05 8.411755e+08 + 1.428466e-05 8.422249e+08 + 1.429673e-05 8.425085e+08 + 1.430880e-05 8.427918e+08 + 1.432087e-05 8.430745e+08 + 1.433294e-05 8.433571e+08 + 1.434501e-05 8.436390e+08 + 1.435708e-05 8.439205e+08 + 1.436915e-05 8.442015e+08 + 1.438122e-05 8.444822e+08 + 1.439329e-05 8.447624e+08 + 1.440535e-05 8.450421e+08 + 1.441742e-05 8.453215e+08 + 1.442949e-05 8.456004e+08 + 1.442034e-05 8.459682e+08 + 1.438056e-05 8.464628e+08 + 1.434078e-05 8.469544e+08 + 1.430100e-05 8.474427e+08 + 1.426122e-05 8.479281e+08 + 1.422145e-05 8.484104e+08 + 1.418167e-05 8.488896e+08 + 1.414189e-05 8.493653e+08 + 1.410211e-05 8.498380e+08 + 1.406234e-05 8.503075e+08 + 1.402256e-05 8.507734e+08 + 1.398278e-05 8.512362e+08 + 1.394300e-05 8.516955e+08 + 1.390322e-05 8.521513e+08 + 1.386345e-05 8.526038e+08 + 1.381027e-05 8.530063e+08 + 1.373205e-05 8.533095e+08 + 1.365382e-05 8.535956e+08 + 1.357559e-05 8.538641e+08 + 1.349737e-05 8.541148e+08 + 1.341914e-05 8.543475e+08 + 1.334091e-05 8.545617e+08 + 1.326268e-05 8.547572e+08 + 1.318446e-05 8.549335e+08 + 1.310623e-05 8.550906e+08 + 1.302800e-05 8.552280e+08 + 1.294977e-05 8.553454e+08 + 1.287155e-05 8.554424e+08 + 1.279332e-05 8.555188e+08 + 1.271509e-05 8.555740e+08 + 1.263687e-05 8.556080e+08 + 1.255864e-05 8.556200e+08 + 1.248085e-05 8.556503e+08 + 1.240883e-05 8.561926e+08 + 1.233682e-05 8.567182e+08 + 1.226480e-05 8.572269e+08 + 1.219278e-05 8.577181e+08 + 1.212076e-05 8.581918e+08 + 1.204875e-05 8.586473e+08 + 1.197673e-05 8.590846e+08 + 1.190471e-05 8.595031e+08 + 1.183270e-05 8.599026e+08 + 1.176068e-05 8.602826e+08 + 1.168866e-05 8.606429e+08 + 1.161665e-05 8.609829e+08 + 1.154463e-05 8.613024e+08 + 1.147261e-05 8.616011e+08 + 1.140059e-05 8.618784e+08 + 1.132858e-05 8.621342e+08 Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/input_more =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/input_more (revision 1491) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/input_more (revision 1491) @@ -0,0 +1,1 @@ + -2652.65 201.81 Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/input_sounding =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/input_sounding (revision 1491) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/input_sounding (revision 1491) @@ -0,0 +1,110 @@ + 7.89 187.87 1.50 + -2652.65 187.87 1.50 0.99 1.38 + -2652.25 187.60 1.50 1.40 1.94 + -2651.75 187.49 1.50 1.57 2.18 + -2650.75 187.38 1.50 1.74 2.41 + -2647.75 187.23 1.50 1.96 2.72 + -2642.75 187.23 1.50 2.11 2.88 + -2632.75 187.23 1.50 2.40 3.21 + -2602.75 187.59 1.50 2.88 3.91 + -2552.75 188.68 1.50 3.79 4.77 + -2452.75 191.66 1.50 5.88 5.85 + -2252.75 194.66 1.50 5.56 4.01 + -2052.75 196.96 1.50 1.86 1.65 + -1852.75 199.23 1.50 -1.99 -0.69 + -1652.75 201.48 1.50 -3.53 -1.44 + -1452.75 203.75 1.50 -5.08 -2.19 + -1252.75 205.99 1.50 -6.55 -2.89 + -1052.75 207.43 1.50 -6.36 -2.57 + -852.75 208.87 1.50 -6.16 -2.25 + -652.75 210.32 1.50 -5.97 -1.92 + -452.75 211.78 1.50 -5.78 -1.60 + -252.75 213.07 1.50 -5.73 -1.36 + -52.75 214.22 1.50 -5.81 -1.20 + 147.25 215.36 1.50 -5.90 -1.04 + 347.25 216.52 1.50 -5.98 -0.88 + 547.25 217.68 1.50 -6.06 -0.72 + 747.25 218.84 1.50 -6.15 -0.57 + 947.25 220.02 1.50 -6.23 -0.41 + 1147.25 221.20 1.50 -6.32 -0.25 + 1347.25 222.44 1.50 -6.04 -0.40 + 1547.25 223.73 1.50 -5.55 -0.72 + 1747.25 225.02 1.50 -5.07 -1.04 + 1947.25 226.32 1.50 -4.58 -1.37 + 2147.25 227.63 1.50 -4.10 -1.69 + 2347.25 228.94 1.50 -3.61 -2.01 + 2547.25 230.26 1.50 -3.13 -2.34 + 2747.25 231.59 1.50 -2.64 -2.66 + 2947.25 232.93 1.50 -2.16 -2.98 + 3147.25 234.20 1.50 -1.55 -2.91 + 3347.25 235.46 1.50 -0.92 -2.72 + 3547.25 236.72 1.50 -0.28 -2.52 + 3747.25 237.99 1.50 0.36 -2.32 + 3947.25 239.27 1.50 1.00 -2.13 + 4147.25 240.56 1.50 1.64 -1.93 + 4347.25 241.85 1.50 2.28 -1.74 + 4547.25 243.15 1.50 2.92 -1.54 + 4747.25 244.45 1.50 3.56 -1.34 + 4947.25 245.76 1.50 4.20 -1.15 + 5147.25 247.03 1.50 4.64 -0.84 + 5347.25 248.14 1.50 4.41 -0.14 + 5547.25 249.26 1.50 4.19 0.56 + 5747.25 250.38 1.50 3.96 1.26 + 5947.25 251.51 1.50 3.73 1.96 + 6147.25 252.64 1.50 3.51 2.66 + 6347.25 253.78 1.50 3.28 3.36 + 6547.25 254.92 1.50 3.05 4.06 + 6747.25 256.07 1.50 2.83 4.76 + 6947.25 257.22 1.50 2.60 5.46 + 7147.25 258.38 1.50 2.38 6.16 + 7347.25 259.54 1.50 2.15 6.86 + 7547.25 260.71 1.50 1.92 7.56 + 7747.25 261.86 1.50 1.76 8.00 + 7947.25 262.97 1.50 1.68 8.06 + 8147.25 264.09 1.50 1.61 8.13 + 8347.25 265.21 1.50 1.54 8.19 + 8547.25 266.34 1.50 1.46 8.25 + 8747.25 267.47 1.50 1.39 8.31 + 8947.25 268.61 1.50 1.31 8.37 + 9147.25 269.75 1.50 1.24 8.44 + 9347.25 270.90 1.50 1.16 8.50 + 9547.25 272.06 1.50 1.09 8.56 + 9747.25 273.21 1.50 1.01 8.62 + 9947.25 274.37 1.50 0.94 8.68 + 10147.25 275.54 1.50 0.87 8.75 + 10347.25 276.71 1.50 0.79 8.81 + 10547.25 277.89 1.50 0.72 8.87 + 10747.25 279.07 1.50 0.62 8.86 + 10947.25 280.25 1.50 0.48 8.69 + 11147.25 281.43 1.50 0.35 8.53 + 11347.25 282.62 1.50 0.21 8.37 + 11547.25 283.81 1.50 0.08 8.21 + 11747.25 285.01 1.50 -0.06 8.05 + 11947.25 286.21 1.50 -0.20 7.89 + 12147.25 287.42 1.50 -0.33 7.73 + 12347.25 288.64 1.50 -0.47 7.57 + 12547.25 289.85 1.50 -0.60 7.41 + 12747.25 291.08 1.50 -0.74 7.25 + 12947.25 292.31 1.50 -0.88 7.09 + 13147.25 293.54 1.50 -1.01 6.93 + 13347.25 294.78 1.50 -1.15 6.76 + 13547.25 296.02 1.50 -1.29 6.60 + 13747.25 297.27 1.50 -1.42 6.44 + 13947.25 298.53 1.50 -1.56 6.28 + 14147.25 299.78 1.50 -1.70 6.12 + 14347.25 300.93 1.50 -1.91 5.89 + 14547.25 302.08 1.50 -2.12 5.66 + 14747.25 303.24 1.50 -2.32 5.43 + 14947.25 304.40 1.50 -2.53 5.20 + 15147.25 305.56 1.50 -2.74 4.97 + 15347.25 306.73 1.50 -2.95 4.74 + 15547.25 307.90 1.50 -3.16 4.51 + 15747.25 309.08 1.50 -3.36 4.28 + 15947.25 310.26 1.50 -3.57 4.05 + 16147.25 311.45 1.50 -3.78 3.82 + 16347.25 312.64 1.50 -3.99 3.59 + 16547.25 313.83 1.50 -4.20 3.36 + 16747.25 315.03 1.50 -4.41 3.13 + 16947.25 316.23 1.50 -4.61 2.90 + 17147.25 317.44 1.50 -4.82 2.67 + 17347.25 318.66 1.50 -5.03 2.44 Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/input_therm =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/input_therm (revision 1491) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/input_therm (revision 1491) @@ -0,0 +1,109 @@ + 190.51 747.35 7.893504e+02 2.070555e-02 200.71 + 190.51 746.99 7.893195e+02 2.070474e-02 200.42 + 190.51 746.83 7.892809e+02 2.070372e-02 200.30 + 190.51 746.67 7.892036e+02 2.070170e-02 200.17 + 190.51 746.45 7.889719e+02 2.069583e-02 200.00 + 190.51 746.42 7.885858e+02 2.068834e-02 199.97 + 190.51 746.36 7.878140e+02 2.067338e-02 199.92 + 190.51 746.65 7.855054e+02 2.058946e-02 200.15 + 190.51 747.81 7.816891e+02 2.039966e-02 201.07 + 190.51 751.17 7.741555e+02 1.994401e-02 203.73 + 190.51 753.90 7.594910e+02 1.936387e-02 205.92 + 190.50 755.68 7.452648e+02 1.887299e-02 207.34 + 190.50 757.38 7.313261e+02 1.839771e-02 208.71 + 190.50 759.05 7.178578e+02 1.794510e-02 210.07 + 190.50 760.71 7.046375e+02 1.750082e-02 211.42 + 190.50 762.33 6.916686e+02 1.706815e-02 212.74 + 190.50 762.91 6.791204e+02 1.672212e-02 213.21 + 190.50 763.48 6.667998e+02 1.638236e-02 213.69 + 190.50 764.06 6.547028e+02 1.604877e-02 214.17 + 190.50 764.64 6.428251e+02 1.572124e-02 214.64 + 190.50 765.01 6.312129e+02 1.541434e-02 214.95 + 190.50 765.19 6.198571e+02 1.512684e-02 215.10 + 190.50 765.37 6.087055e+02 1.484450e-02 215.25 + 190.50 765.55 5.977546e+02 1.456725e-02 215.40 + 190.50 765.73 5.870007e+02 1.429499e-02 215.54 + 190.50 765.91 5.764402e+02 1.402762e-02 215.69 + 190.50 766.09 5.660698e+02 1.376506e-02 215.84 + 190.50 766.27 5.558859e+02 1.350723e-02 215.99 + 190.50 766.53 5.459407e+02 1.325292e-02 216.20 + 190.50 766.83 5.362040e+02 1.300253e-02 216.45 + 190.50 767.13 5.266410e+02 1.275661e-02 216.70 + 190.50 767.43 5.172485e+02 1.251507e-02 216.95 + 190.50 767.73 5.080236e+02 1.227784e-02 217.20 + 190.50 768.03 4.989632e+02 1.204485e-02 217.44 + 190.49 768.33 4.900644e+02 1.181600e-02 217.69 + 190.49 768.63 4.813242e+02 1.159124e-02 217.94 + 190.49 768.93 4.727400e+02 1.137049e-02 218.19 + 190.49 769.16 4.643743e+02 1.115970e-02 218.38 + 190.49 769.36 4.561774e+02 1.095457e-02 218.55 + 190.49 769.57 4.481252e+02 1.075306e-02 218.73 + 190.49 769.78 4.402152e+02 1.055511e-02 218.90 + 190.49 769.99 4.324447e+02 1.036065e-02 219.07 + 190.49 770.19 4.248115e+02 1.016962e-02 219.24 + 190.49 770.40 4.173130e+02 9.981969e-03 219.42 + 190.49 770.61 4.099467e+02 9.797628e-03 219.59 + 190.49 770.81 4.027106e+02 9.616540e-03 219.76 + 190.49 771.02 3.956022e+02 9.438649e-03 219.93 + 190.49 771.17 3.886273e+02 9.266076e-03 220.07 + 190.49 771.16 3.818018e+02 9.103757e-03 220.05 + 190.49 771.15 3.750963e+02 8.944290e-03 220.04 + 190.49 771.13 3.685084e+02 8.787623e-03 220.03 + 190.49 771.12 3.620363e+02 8.633709e-03 220.02 + 190.49 771.11 3.556779e+02 8.482497e-03 220.01 + 190.49 771.09 3.494312e+02 8.333941e-03 219.99 + 190.49 771.08 3.432941e+02 8.187993e-03 219.98 + 190.49 771.06 3.372648e+02 8.044610e-03 219.97 + 190.49 771.05 3.313414e+02 7.903744e-03 219.96 + 190.49 771.04 3.255221e+02 7.765353e-03 219.95 + 190.49 771.02 3.198050e+02 7.629392e-03 219.94 + 190.49 771.01 3.141882e+02 7.495819e-03 219.92 + 190.49 770.97 3.086673e+02 7.365176e-03 219.89 + 190.49 770.89 3.032396e+02 7.237664e-03 219.82 + 190.48 770.81 2.979072e+02 7.112395e-03 219.76 + 190.48 770.73 2.926686e+02 6.989328e-03 219.69 + 190.48 770.66 2.875222e+02 6.868426e-03 219.62 + 190.48 770.58 2.824662e+02 6.749649e-03 219.56 + 190.48 770.50 2.774992e+02 6.632961e-03 219.49 + 190.48 770.42 2.726195e+02 6.518325e-03 219.42 + 190.48 770.34 2.678256e+02 6.405705e-03 219.36 + 190.48 770.26 2.631160e+02 6.295065e-03 219.29 + 190.47 770.18 2.584892e+02 6.186371e-03 219.22 + 190.47 770.11 2.539438e+02 6.079588e-03 219.16 + 190.47 770.03 2.494783e+02 5.974683e-03 219.09 + 190.47 769.95 2.450913e+02 5.871622e-03 219.02 + 190.47 769.87 2.407815e+02 5.770374e-03 218.96 + 190.47 769.79 2.365421e+02 5.670684e-03 218.89 + 190.47 769.69 2.323676e+02 5.572339e-03 218.81 + 190.47 769.60 2.282667e+02 5.475729e-03 218.73 + 190.46 769.51 2.242383e+02 5.380824e-03 218.66 + 190.46 769.42 2.202809e+02 5.287594e-03 218.58 + 190.46 769.33 2.163933e+02 5.196009e-03 218.50 + 190.46 769.24 2.125744e+02 5.106041e-03 218.43 + 190.46 769.15 2.088229e+02 5.017660e-03 218.35 + 190.46 769.05 2.051375e+02 4.930839e-03 218.27 + 190.46 768.96 2.015172e+02 4.845551e-03 218.19 + 190.46 768.87 1.979608e+02 4.761768e-03 218.12 + 190.46 768.78 1.944672e+02 4.679463e-03 218.04 + 190.45 768.69 1.910352e+02 4.598611e-03 217.96 + 190.45 768.60 1.876638e+02 4.519185e-03 217.89 + 190.45 768.50 1.843519e+02 4.441162e-03 217.81 + 190.45 768.41 1.810984e+02 4.364515e-03 217.73 + 190.45 768.32 1.779024e+02 4.289221e-03 217.65 + 190.45 768.22 1.747608e+02 4.215284e-03 217.57 + 190.45 768.01 1.716499e+02 4.143034e-03 217.40 + 190.45 767.81 1.685944e+02 4.072070e-03 217.23 + 190.45 767.60 1.655932e+02 4.002369e-03 217.06 + 190.45 767.40 1.626455e+02 3.933908e-03 216.89 + 190.45 767.19 1.597503e+02 3.866667e-03 216.72 + 190.44 766.99 1.569066e+02 3.800622e-03 216.55 + 190.44 766.78 1.541135e+02 3.735753e-03 216.38 + 190.44 766.58 1.513701e+02 3.672039e-03 216.21 + 190.44 766.37 1.486756e+02 3.609459e-03 216.04 + 190.44 766.17 1.460290e+02 3.547993e-03 215.86 + 190.44 765.96 1.434296e+02 3.487621e-03 215.69 + 190.44 765.75 1.408764e+02 3.428324e-03 215.52 + 190.44 765.55 1.383686e+02 3.370082e-03 215.35 + 190.44 765.34 1.359056e+02 3.312878e-03 215.18 + 190.44 765.13 1.334863e+02 3.256691e-03 215.01 + 190.44 764.93 1.311101e+02 3.201504e-03 214.84 Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/mpi.ll =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/mpi.ll (revision 1491) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/mpi.ll (revision 1491) @@ -0,0 +1,1 @@ +link newphys10m/mpi.ll Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/namelist.input (revision 1491) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/namelist.input (revision 1491) @@ -0,0 +1,1 @@ +link newphys10m/namelist.input Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/newphys10m/mpi.ll =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/newphys10m/mpi.ll (revision 1491) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/newphys10m/mpi.ll (revision 1491) @@ -0,0 +1,22 @@ +# Nom du travail LoadLeveler +# @ job_name= LES +# Fichier de sortie standard du travail +# @ output = $(job_name).$(jobid) +# Fichier de sortie d'erreur du travail +# @ error = $(job_name).$(jobid) +# Type de travail +# @ job_type= parallel +# Nombre de processus demandes +# @ total_tasks = 1024 +# Temps ELAPSED max. pour l'ensemble du job en hh:mm:ss +# @ wall_clock_limit = 20:00:00 +# @ queue +cd $LOADL_STEP_INITDIR +\rm -rf core* +\rm -rf *rsl.* +\rm -rf wrfi* +\rm -rf wrfo* +poe ./ideal.exe +mv rsl.error.0000 ideal.rsl.error.0000 +mv rsl.out.0000 ideal.rsl.out.0000 +poe ./wrf.exe Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/newphys10m/namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/newphys10m/namelist.input (revision 1491) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/newphys10m/namelist.input (revision 1491) @@ -0,0 +1,83 @@ + &time_control + start_month = 01, !! (p1) Start Martian Month + start_day = 01, !! (p1) Start Martian Day + start_hour = 00, !! (p1) Start Martian Hour (at longitude 0) + start_minute = 00, !! two digit minute of starting time + start_second = 00, !! two digit second of starting time + run_seconds = 37000 !! Duration of simulation [LES: !=0] + history_interval_s = 1 !! Frequency of outputs in seconds [LES: !=0] + frames_per_outfile = 3700 !! Size of time dimension in files + start_year = 9999 !! Idealized simulation [LES: 9999] + end_year = 9999 !! Idealized simulation [LES: 9999] + restart = .false. !! (*) restart run ? + restart_interval = 8880 !! (*) Frequency of output restart files ? + io_form_history = 2 !! (*) Choice of NETCDF for ouputs + io_form_restart = 2 !! (*) Choice of NETCDF for ouputs + io_form_input = 2 !! (*) Choice of NETCDF for ouputs + io_form_boundary = 2 !! (*) Choice of NETCDF for ouputs + debug_level = 0 !! (*) Verbose level + / + + &domains + ztop = 6000. !! Altitude above ground of model top [LES: !=0] + time_step = 0 !! Dynamical timestep + time_step_fract_num = 1 !! Additional fraction to time_step: numerator + time_step_fract_den = 10 !! Additional fraction to time_step: denominator + dx = 10 !! (p2) Horizontal resolution + dy = 10 !! (p2) Horizontal resolution (should be equal to dx) + e_we = 417 !! (r)(p2) Number of longitude grid points + e_sn = 417 !! (r)(p2) Number of latitude grid points + e_vert = 301 !! (r)(p2) Number of vertical levels + / + + &physics + isfflx = 1 !! 0 : UST/HFX from tke_... namelist values + !! 1 : UST/HFX from LMD physics + !! 2 : UST from LMD / HFX from namelist.input + init_TI = 260. !! Define constant thermal inertia value [LES: !=0] + init_AL = 0.26 !! Define constant albedo value [LES: !=0] + num_soil_layers = 18 !!! NOUVELLE PHYSIQUE + !! + !! OPTIONAL + !! + radt = 64 !! Ratio between physical and dynamical time step + mars = 0 !! (r)(p2) Configuration of tracers: + !! 0: no tracers, 1: water vapor + ice, 2: dust + init_MU = 0.01 !! (p3) Multiply ini & bdy zonal wind by init_U + init_MV = 0.01 !! (p3) Multiply ini & bdy meridional wind by init_V + init_U = 10. + / + + &dynamics + diff_opt = 2 !! (*d) Diffusion option [LES: 2] + km_opt = 2 !! (*d) Eddy coefficient option [LES: 2] + diff_6th_opt = 0 !! (*d) Knievel numerical diffusion [LES: 0] + mix_full_fields = T !! (*d) Set true if mixing ref + pert profiles [LES: T] + !! + !! OPTIONAL + !! +! c_k = 0.15 !! Mixing coefficient constant for km_opt=2 [default: 0.15] +! c_s = 0.18 !! Mixing coefficient constant for km_opt=3 [default: 0.25] +! pd_tke = T !! Positive definite advection of tke +! tke_heat_flux = 1. !! Constant surface thermal flux (H/(rho*cp), K m/s) +! tke_drag_coefficient = 0.0025, !! Constant surface drag coefficient (Cd, dimensionless) +! mix_upper_bound = 100. !! Non-dimensional upper limit for diffusion coeffs +! mix_isotropic = 1 !! 0=anistropic vertical/horizontal diffusion coeffs, 1=isotropic +! pert_coriolis = T !! Coriolis only acts on wind perturbation (idealized) + / + + &bdy_control + periodic_x = T !! (p3) Periodic boundary conditions over x [LES: T] + periodic_y = T !! (p3) Periodic boundary conditions over y [LES: T] + / + + &grib2 + / + + &fdda + / + + &namelist_quilt !! (*) + nio_tasks_per_group = 0, !! (*) + nio_groups = 1, !! (*) + / !! (*) Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/newphys50m/mpi.ll =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/newphys50m/mpi.ll (revision 1491) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/newphys50m/mpi.ll (revision 1491) @@ -0,0 +1,22 @@ +# Nom du travail LoadLeveler +# @ job_name= LES +# Fichier de sortie standard du travail +# @ output = $(job_name).$(jobid) +# Fichier de sortie d'erreur du travail +# @ error = $(job_name).$(jobid) +# Type de travail +# @ job_type= parallel +# Nombre de processus demandes +# @ total_tasks = 512 +# Temps ELAPSED max. pour l'ensemble du job en hh:mm:ss +# @ wall_clock_limit = 20:00:00 +# @ queue +cd $LOADL_STEP_INITDIR +\rm -rf core* +\rm -rf *rsl.* +\rm -rf wrfi* +\rm -rf wrfo* +poe ./ideal.exe +mv rsl.error.0000 ideal.rsl.error.0000 +mv rsl.out.0000 ideal.rsl.out.0000 +poe ./wrf.exe Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/newphys50m/namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/newphys50m/namelist.input (revision 1491) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/newphys50m/namelist.input (revision 1491) @@ -0,0 +1,83 @@ + &time_control + start_month = 01, !! (p1) Start Martian Month + start_day = 01, !! (p1) Start Martian Day + start_hour = 00, !! (p1) Start Martian Hour (at longitude 0) + start_minute = 00, !! two digit minute of starting time + start_second = 00, !! two digit second of starting time + run_seconds = 88800 !! Duration of simulation [LES: !=0] + history_interval_s = 6 !! Frequency of outputs in seconds [LES: !=0] + frames_per_outfile = 100 !! Size of time dimension in files + start_year = 9999 !! Idealized simulation [LES: 9999] + end_year = 9999 !! Idealized simulation [LES: 9999] + restart = .false. !! (*) restart run ? + restart_interval = 8880 !! (*) Frequency of output restart files ? + io_form_history = 2 !! (*) Choice of NETCDF for ouputs + io_form_restart = 2 !! (*) Choice of NETCDF for ouputs + io_form_input = 2 !! (*) Choice of NETCDF for ouputs + io_form_boundary = 2 !! (*) Choice of NETCDF for ouputs + debug_level = 0 !! (*) Verbose level + / + + &domains + ztop = 8000. !! Altitude above ground of model top [LES: !=0] + time_step = 0 !! Dynamical timestep + time_step_fract_num = 3 !! Additional fraction to time_step: numerator + time_step_fract_den = 4 !! Additional fraction to time_step: denominator + dx = 50 !! (p2) Horizontal resolution + dy = 50 !! (p2) Horizontal resolution (should be equal to dx) + e_we = 289 !! (r)(p2) Number of longitude grid points + e_sn = 289 !! (r)(p2) Number of latitude grid points + e_vert = 151 !! (r)(p2) Number of vertical levels + / + + &physics + isfflx = 1 !! 0 : UST/HFX from tke_... namelist values + !! 1 : UST/HFX from LMD physics + !! 2 : UST from LMD / HFX from namelist.input + init_TI = 260. !! Define constant thermal inertia value [LES: !=0] + init_AL = 0.26 !! Define constant albedo value [LES: !=0] + num_soil_layers = 18 !!! NOUVELLE PHYSIQUE + !! + !! OPTIONAL + !! + radt = 8 !! Ratio between physical and dynamical time step + mars = 0 !! (r)(p2) Configuration of tracers: + !! 0: no tracers, 1: water vapor + ice, 2: dust + init_MU = 0.01 !! (p3) Multiply ini & bdy zonal wind by init_U + init_MV = 0.01 !! (p3) Multiply ini & bdy meridional wind by init_V + init_U = 10. + / + + &dynamics + diff_opt = 2 !! (*d) Diffusion option [LES: 2] + km_opt = 2 !! (*d) Eddy coefficient option [LES: 2] + diff_6th_opt = 0 !! (*d) Knievel numerical diffusion [LES: 0] + mix_full_fields = T !! (*d) Set true if mixing ref + pert profiles [LES: T] + !! + !! OPTIONAL + !! +! c_k = 0.15 !! Mixing coefficient constant for km_opt=2 [default: 0.15] +! c_s = 0.18 !! Mixing coefficient constant for km_opt=3 [default: 0.25] +! pd_tke = T !! Positive definite advection of tke +! tke_heat_flux = 1. !! Constant surface thermal flux (H/(rho*cp), K m/s) +! tke_drag_coefficient = 0.0025, !! Constant surface drag coefficient (Cd, dimensionless) +! mix_upper_bound = 100. !! Non-dimensional upper limit for diffusion coeffs +! mix_isotropic = 1 !! 0=anistropic vertical/horizontal diffusion coeffs, 1=isotropic +! pert_coriolis = T !! Coriolis only acts on wind perturbation (idealized) + / + + &bdy_control + periodic_x = T !! (p3) Periodic boundary conditions over x [LES: T] + periodic_y = T !! (p3) Periodic boundary conditions over y [LES: T] + / + + &grib2 + / + + &fdda + / + + &namelist_quilt !! (*) + nio_tasks_per_group = 0, !! (*) + nio_groups = 1, !! (*) + / !! (*) Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/newphys50m/namelist.input.lesspoints =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/newphys50m/namelist.input.lesspoints (revision 1491) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/newphys50m/namelist.input.lesspoints (revision 1491) @@ -0,0 +1,83 @@ + &time_control + start_month = 01, !! (p1) Start Martian Month + start_day = 01, !! (p1) Start Martian Day + start_hour = 00, !! (p1) Start Martian Hour (at longitude 0) + start_minute = 00, !! two digit minute of starting time + start_second = 00, !! two digit second of starting time + run_seconds = 88800 !! Duration of simulation [LES: !=0] + history_interval_s = 6 !! Frequency of outputs in seconds [LES: !=0] + frames_per_outfile = 100 !! Size of time dimension in files + start_year = 9999 !! Idealized simulation [LES: 9999] + end_year = 9999 !! Idealized simulation [LES: 9999] + restart = .false. !! (*) restart run ? + restart_interval = 8880 !! (*) Frequency of output restart files ? + io_form_history = 2 !! (*) Choice of NETCDF for ouputs + io_form_restart = 2 !! (*) Choice of NETCDF for ouputs + io_form_input = 2 !! (*) Choice of NETCDF for ouputs + io_form_boundary = 2 !! (*) Choice of NETCDF for ouputs + debug_level = 0 !! (*) Verbose level + / + + &domains + ztop = 12000. !! Altitude above ground of model top [LES: !=0] + time_step = 0 !! Dynamical timestep + time_step_fract_num = 3 !! Additional fraction to time_step: numerator + time_step_fract_den = 4 !! Additional fraction to time_step: denominator + dx = 50 !! (p2) Horizontal resolution + dy = 50 !! (p2) Horizontal resolution (should be equal to dx) + e_we = 145 !! (r)(p2) Number of longitude grid points + e_sn = 145 !! (r)(p2) Number of latitude grid points + e_vert = 201 !! (r)(p2) Number of vertical levels + / + + &physics + isfflx = 1 !! 0 : UST/HFX from tke_... namelist values + !! 1 : UST/HFX from LMD physics + !! 2 : UST from LMD / HFX from namelist.input + init_TI = 260. !! Define constant thermal inertia value [LES: !=0] + init_AL = 0.26 !! Define constant albedo value [LES: !=0] + num_soil_layers = 18 !!! NOUVELLE PHYSIQUE + !! + !! OPTIONAL + !! + radt = 8 !! Ratio between physical and dynamical time step + mars = 0 !! (r)(p2) Configuration of tracers: + !! 0: no tracers, 1: water vapor + ice, 2: dust + init_MU = 0.01 !! (p3) Multiply ini & bdy zonal wind by init_U + init_MV = 0.01 !! (p3) Multiply ini & bdy meridional wind by init_V + init_U = 10. + / + + &dynamics + diff_opt = 2 !! (*d) Diffusion option [LES: 2] + km_opt = 2 !! (*d) Eddy coefficient option [LES: 2] + diff_6th_opt = 0 !! (*d) Knievel numerical diffusion [LES: 0] + mix_full_fields = T !! (*d) Set true if mixing ref + pert profiles [LES: T] + !! + !! OPTIONAL + !! +! c_k = 0.15 !! Mixing coefficient constant for km_opt=2 [default: 0.15] +! c_s = 0.18 !! Mixing coefficient constant for km_opt=3 [default: 0.25] +! pd_tke = T !! Positive definite advection of tke +! tke_heat_flux = 1. !! Constant surface thermal flux (H/(rho*cp), K m/s) +! tke_drag_coefficient = 0.0025, !! Constant surface drag coefficient (Cd, dimensionless) +! mix_upper_bound = 100. !! Non-dimensional upper limit for diffusion coeffs +! mix_isotropic = 1 !! 0=anistropic vertical/horizontal diffusion coeffs, 1=isotropic +! pert_coriolis = T !! Coriolis only acts on wind perturbation (idealized) + / + + &bdy_control + periodic_x = T !! (p3) Periodic boundary conditions over x [LES: T] + periodic_y = T !! (p3) Periodic boundary conditions over y [LES: T] + / + + &grib2 + / + + &fdda + / + + &namelist_quilt !! (*) + nio_tasks_per_group = 0, !! (*) + nio_groups = 1, !! (*) + / !! (*) Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/rename.py =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/rename.py (revision 1491) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/rename.py (revision 1491) @@ -0,0 +1,7 @@ +#! /usr/bin/env python +import os +path = os.getcwd() +filenames = os.listdir(path) +for filename in filenames: + print filename + os.rename(filename, filename.replace(":", "-")) Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/run.def =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/run.def (revision 1491) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/run.def (revision 1491) @@ -0,0 +1,2 @@ +# some definitions for the physics, in file 'callphys.def' +INCLUDEDEF=callphys.def Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/transfer =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/transfer (revision 1491) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/2015_spiga_LES4InSight/transfer (revision 1491) @@ -0,0 +1,1 @@ +find wrfout_d01_9999-01-01_* -exec mfput -v {} \;