Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/LES_DDphoenix/callphys.def =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/LES_DDphoenix/callphys.def (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/LES_DDphoenix/callphys.def (revision 1291) @@ -0,0 +1,94 @@ +General options +~~~~~~~~~~~~~~~ +tracer (Run with or without tracer transport ?) +F +diurnal (Diurnal cycle ? if diurnal=F, diurnal averaged solar heating) +T +season (Seasonal cycle ? if season=F, Ls stays constant like in "start") +T +lwrite (want some more output on the screen ?) +F +stats (Saving statistics in file "cumul" ?) -- noWRF +F +calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF +F +Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +~~~~~~~~~~~~~ +iaervar (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario +1 (=4 Mars Year 24 from TES assimilation) +iddist (Dust vertical distribution: =0: old distrib. (Pollack90) +3 (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario ) +topdustref (Dust top altitude (km). Matter only if iddist=1) +55. +Physical Parameterizations : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +callrad (call radiative transfer ?) +T +callnlte (call NLTE radiative schemes ? matter only if callrad=T) +F +callnirco2 (call CO2 NIR absorption ? matter only if callrad=T) +F +calldifv (call turbulent vertical diffusion ?) +T +calladj (call convective adjustment ?) +F +callcond (call CO2 condensation ?) +F +callsoil (call thermal conduction in the soil ?) +T +calllott (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF +F +Radiative transfer options : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +iradia (the rad.transfer is computed every "iradia" physical timestep) +1 +callg2d (Output of the exchange coefficient mattrix ? for diagnostic only) +F +rayleigh (Rayleigh scattering : should be =F for now) +F +Tracer (dust water, ice and/or chemical species) options (use if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +dustbin (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin)) +0 +active (DUST: Radiatively active dust ? (uses q(1) to q(dustbin)) +F +doubleq (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?) +F +lifting (DUST: lifted by GCM surface winds ?) +F +dustdevil (DUST: lifted by dust devils ?) +F +scavenging (DUST: Scavenging by CO2 snowfall ?) +F +sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?) +F +iceparty (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1)) +F +activice (WATERICE: Radiatively active transported atmospheric water ice ?) +F +water (WATER: Compute water cycle using q(nqmx) ) +F +caps (WATER: put the current permanent caps at both poles) +F +photochem (PHOTOCHEMISTRY: chemical species included) +F +Thermospheric options (relevant if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +callthermos (call thermosphere ?) +F +thermoswater (WATER: included without cycle only if water=F) +F +callconduct (call thermal conduction ? matter only if callthermos=T) +F +calleuv (call EUV heating ? matter only if callthermos=T) +F +callmolvis (call molecular viscosity ? matter only if callthermos=T) +F +callmoldiff (call molecular diffusion ? matter only if callthermos=T) +F +thermochem (call thermospheric photochemistry ? matter only if callthermos=T) +F +solarcondate (date for solar flux calculation: 1985 < date < 2002)) +1993.4 (Solar min=1996.4 ave=1993.4 max=1990.6) + + Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/LES_DDphoenix/dustopacity.def =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/LES_DDphoenix/dustopacity.def (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/LES_DDphoenix/dustopacity.def (revision 1291) @@ -0,0 +1,1 @@ +0.3 Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/LES_DDphoenix/input_coord =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/LES_DDphoenix/input_coord (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/LES_DDphoenix/input_coord (revision 1291) @@ -0,0 +1,101 @@ +233.0 lon +68.0 lat +90.0 ls +09. lct +20. ztop +100. nz + + + +0.0 lon +85.0 lat +0.0 ls +00. lct +20. ztop +100. nz + +16. + +353.87 lon +-1.88 lat +68.5 ls +06. lct +20. ztop +100. nz + + + +205. lon +21.8 lat +47.1 ls +10. lct +20. ztop +100. nz + + + +236.6 lon +-10.2 lat +67.0 ls +13. lct +20. ztop +100. nz + + +205. lon +21.8 lat +47.1 ls +12. lct +60. ztop +200. nz + + +60.0 lon +-65.0 lat +33.0 ls +15. lct +60. ztop +100. nz + + +236.6 lon +-10.2 lat +67.0 ls +08. lct +20. ztop +100. nz + +74. lon +20.6 lat +47.8 ls +08. lct +60. ztop +200. nz + +237.2 lon +12.3 lat +52.1 ls +17. lct +60. ztop +200. nz + +204.6 lon +13.7 lat +51.2 ls +17. lct +60. ztop +200. nz + +205. lon +21.8 lat +47.1 ls +08. lct +60. ztop +200. nz + +237.2 lon +12.3 lat +52.1 ls +12. lct +20. ztop +100. nz Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/LES_DDphoenix/input_more =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/LES_DDphoenix/input_more (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/LES_DDphoenix/input_more (revision 1291) @@ -0,0 +1,1 @@ + -4121.25 235.39 Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/LES_DDphoenix/input_sounding =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/LES_DDphoenix/input_sounding (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/LES_DDphoenix/input_sounding (revision 1291) @@ -0,0 +1,108 @@ + 8.19 218.28 0.00 + -4121.25 218.28 0.00 0.00 0.00 + -4120.25 215.90 0.00 -3.69 -1.01 + -4116.25 206.37 0.00 -4.98 -1.36 + -4111.25 203.46 0.00 -5.43 -1.49 + -4101.25 203.37 0.00 -6.16 -1.69 + -4071.25 203.38 0.00 -7.15 -1.87 + -4021.25 203.39 0.00 -7.94 -1.90 + -3921.25 203.46 0.00 -9.10 -1.72 + -3721.25 203.97 0.00 -11.19 -0.89 + -3521.25 204.73 0.00 -13.09 0.11 + -3321.25 205.74 0.00 -14.04 0.86 + -3121.25 206.75 0.00 -14.99 1.62 + -2921.25 207.68 0.00 -15.50 2.04 + -2721.25 208.52 0.00 -15.48 2.06 + -2521.25 209.36 0.00 -15.45 2.08 + -2321.25 210.20 0.00 -15.43 2.11 + -2121.25 211.05 0.00 -15.40 2.13 + -1921.25 211.70 0.00 -15.25 2.12 + -1721.25 212.22 0.00 -15.01 2.08 + -1521.25 212.73 0.00 -14.76 2.04 + -1321.25 213.25 0.00 -14.52 2.00 + -1121.25 213.77 0.00 -14.28 1.96 + -921.25 214.29 0.00 -14.04 1.92 + -721.25 214.81 0.00 -13.80 1.88 + -521.25 215.34 0.00 -13.58 1.85 + -321.25 215.92 0.00 -13.47 1.89 + -121.25 216.50 0.00 -13.36 1.94 + 78.75 217.08 0.00 -13.25 1.98 + 278.75 217.66 0.00 -13.14 2.02 + 478.75 218.24 0.00 -13.03 2.06 + 678.75 218.82 0.00 -12.92 2.11 + 878.75 219.40 0.00 -12.81 2.15 + 1078.75 219.98 0.00 -12.70 2.19 + 1278.75 220.56 0.00 -12.60 2.24 + 1478.75 221.15 0.00 -12.49 2.28 + 1678.75 221.73 0.00 -12.38 2.32 + 1878.75 222.32 0.00 -12.28 2.36 + 2078.75 223.04 0.00 -12.32 2.36 + 2278.75 223.75 0.00 -12.37 2.37 + 2478.75 224.47 0.00 -12.41 2.38 + 2678.75 225.18 0.00 -12.45 2.38 + 2878.75 225.90 0.00 -12.49 2.39 + 3078.75 226.62 0.00 -12.54 2.39 + 3278.75 227.34 0.00 -12.58 2.40 + 3478.75 228.07 0.00 -12.62 2.40 + 3678.75 228.79 0.00 -12.66 2.41 + 3878.75 229.52 0.00 -12.71 2.41 + 4078.75 230.25 0.00 -12.75 2.42 + 4278.75 230.98 0.00 -12.79 2.42 + 4478.75 231.71 0.00 -12.83 2.43 + 4678.75 232.44 0.00 -12.88 2.43 + 4878.75 233.18 0.00 -12.92 2.44 + 5078.75 233.91 0.00 -12.96 2.44 + 5278.75 234.65 0.00 -13.00 2.45 + 5478.75 235.45 0.00 -13.12 2.43 + 5678.75 236.28 0.00 -13.26 2.40 + 5878.75 237.11 0.00 -13.40 2.38 + 6078.75 237.94 0.00 -13.54 2.35 + 6278.75 238.77 0.00 -13.68 2.32 + 6478.75 239.61 0.00 -13.81 2.30 + 6678.75 240.45 0.00 -13.95 2.27 + 6878.75 241.29 0.00 -14.09 2.24 + 7078.75 242.13 0.00 -14.23 2.22 + 7278.75 242.98 0.00 -14.37 2.19 + 7478.75 243.82 0.00 -14.51 2.16 + 7678.75 244.68 0.00 -14.65 2.14 + 7878.75 245.53 0.00 -14.79 2.11 + 8078.75 246.38 0.00 -14.93 2.08 + 8278.75 247.24 0.00 -15.07 2.06 + 8478.75 248.10 0.00 -15.21 2.03 + 8678.75 248.96 0.00 -15.35 2.00 + 8878.75 249.83 0.00 -15.49 1.98 + 9078.75 250.69 0.00 -15.62 1.95 + 9278.75 251.56 0.00 -15.76 1.92 + 9478.75 252.43 0.00 -15.90 1.90 + 9678.75 253.31 0.00 -16.04 1.87 + 9878.75 254.20 0.00 -16.26 1.84 + 10078.75 255.11 0.00 -16.53 1.81 + 10278.75 256.02 0.00 -16.80 1.78 + 10478.75 256.94 0.00 -17.07 1.75 + 10678.75 257.86 0.00 -17.34 1.72 + 10878.75 258.78 0.00 -17.61 1.68 + 11078.75 259.70 0.00 -17.88 1.65 + 11278.75 260.62 0.00 -18.15 1.62 + 11478.75 261.55 0.00 -18.42 1.59 + 11678.75 262.48 0.00 -18.69 1.56 + 11878.75 263.41 0.00 -18.96 1.53 + 12078.75 264.35 0.00 -19.23 1.50 + 12278.75 265.29 0.00 -19.50 1.47 + 12478.75 266.23 0.00 -19.77 1.43 + 12678.75 267.17 0.00 -20.04 1.40 + 12878.75 268.11 0.00 -20.31 1.37 + 13078.75 269.06 0.00 -20.58 1.34 + 13278.75 270.01 0.00 -20.85 1.31 + 13478.75 270.97 0.00 -21.12 1.28 + 13678.75 271.92 0.00 -21.39 1.25 + 13878.75 272.88 0.00 -21.66 1.22 + 14078.75 273.84 0.00 -21.93 1.18 + 14278.75 274.81 0.00 -22.20 1.15 + 14478.75 275.77 0.00 -22.47 1.12 + 14678.75 276.74 0.00 -22.74 1.09 + 14878.75 277.71 0.00 -23.01 1.06 + 15078.75 278.69 0.00 -23.31 1.03 + 15278.75 279.66 0.00 -23.66 0.99 + 15478.75 280.64 0.00 -24.00 0.96 + 15678.75 281.62 0.00 -24.35 0.93 + 15878.75 282.60 0.00 -24.69 0.90 Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/LES_DDphoenix/input_therm =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/LES_DDphoenix/input_therm (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/LES_DDphoenix/input_therm (revision 1291) @@ -0,0 +1,107 @@ + 191.14 839.14 8.187230e+02 1.962268e-02 235.39 + 191.14 839.14 8.186497e+02 1.962093e-02 232.81 + 191.14 839.14 8.183566e+02 1.961390e-02 222.52 + 191.14 839.14 8.179904e+02 1.960760e-02 219.35 + 191.14 839.14 8.172585e+02 1.959660e-02 219.21 + 191.14 839.14 8.150653e+02 1.955333e-02 219.07 + 191.14 839.14 8.114193e+02 1.948710e-02 218.82 + 191.14 839.14 8.041653e+02 1.934815e-02 218.40 + 191.14 839.14 7.898306e+02 1.903689e-02 217.94 + 191.14 839.14 7.757439e+02 1.870710e-02 217.75 + 191.14 839.14 7.619154e+02 1.836278e-02 217.81 + 191.14 839.14 7.483335e+02 1.802459e-02 217.87 + 191.14 839.14 7.349937e+02 1.770081e-02 217.85 + 191.14 839.14 7.218920e+02 1.739271e-02 217.72 + 191.14 839.14 7.090237e+02 1.709011e-02 217.59 + 191.14 839.14 6.963849e+02 1.679290e-02 217.46 + 191.14 839.14 6.839713e+02 1.650098e-02 217.34 + 191.14 839.14 6.717127e+02 1.622870e-02 216.99 + 191.14 839.14 6.596332e+02 1.597024e-02 216.51 + 191.14 839.14 6.477710e+02 1.571642e-02 216.04 + 191.14 839.14 6.361221e+02 1.546716e-02 215.56 + 191.14 839.14 6.246826e+02 1.522239e-02 215.08 + 191.14 839.14 6.134489e+02 1.498202e-02 214.60 + 191.14 839.14 6.024171e+02 1.474598e-02 214.13 + 191.14 839.14 5.915521e+02 1.451246e-02 213.65 + 191.14 839.14 5.807031e+02 1.427334e-02 213.21 + 191.14 839.14 5.700531e+02 1.403862e-02 212.77 + 191.14 839.14 5.595984e+02 1.380820e-02 212.33 + 191.14 839.14 5.493354e+02 1.358200e-02 211.89 + 191.14 839.14 5.392608e+02 1.335995e-02 211.45 + 191.14 839.14 5.293708e+02 1.314198e-02 211.01 + 191.14 839.14 5.196622e+02 1.292800e-02 210.57 + 191.14 839.14 5.101317e+02 1.271795e-02 210.12 + 191.14 839.14 5.007760e+02 1.251175e-02 209.68 + 191.14 839.14 4.915919e+02 1.230932e-02 209.24 + 191.14 839.14 4.825761e+02 1.211062e-02 208.80 + 191.14 839.14 4.737077e+02 1.191429e-02 208.37 + 191.14 839.14 4.647972e+02 1.170716e-02 208.02 + 191.14 839.14 4.560543e+02 1.150393e-02 207.67 + 191.14 839.14 4.474759e+02 1.130452e-02 207.32 + 191.14 839.14 4.390588e+02 1.110886e-02 206.97 + 191.14 839.14 4.308001e+02 1.091688e-02 206.63 + 191.14 839.14 4.226967e+02 1.072851e-02 206.28 + 191.14 839.14 4.147457e+02 1.054369e-02 205.93 + 191.14 839.14 4.069443e+02 1.036234e-02 205.58 + 191.14 839.14 3.992897e+02 1.018440e-02 205.23 + 191.14 839.14 3.917790e+02 1.000981e-02 204.89 + 191.14 839.14 3.844096e+02 9.838507e-03 204.54 + 191.14 839.14 3.771788e+02 9.670424e-03 204.19 + 191.14 839.14 3.700840e+02 9.505502e-03 203.84 + 191.14 839.14 3.631227e+02 9.343682e-03 203.49 + 191.14 839.14 3.562924e+02 9.184906e-03 203.15 + 191.14 839.14 3.495905e+02 9.029117e-03 202.80 + 191.14 839.14 3.430146e+02 8.876259e-03 202.45 + 191.14 839.14 3.364218e+02 8.717966e-03 202.13 + 191.14 839.14 3.299084e+02 8.559917e-03 201.83 + 191.14 839.14 3.235212e+02 8.404926e-03 201.52 + 191.14 839.14 3.172577e+02 8.252937e-03 201.22 + 191.14 839.14 3.111154e+02 8.103890e-03 200.91 + 191.14 839.14 3.050920e+02 7.957730e-03 200.61 + 191.14 839.14 2.991852e+02 7.814398e-03 200.30 + 191.14 839.14 2.933928e+02 7.673842e-03 200.00 + 191.14 839.14 2.877126e+02 7.536007e-03 199.69 + 191.14 839.14 2.821423e+02 7.400841e-03 199.39 + 191.14 839.14 2.766798e+02 7.268291e-03 199.09 + 191.14 839.14 2.713232e+02 7.138308e-03 198.78 + 191.14 839.14 2.660702e+02 7.010841e-03 198.48 + 191.14 839.14 2.609189e+02 6.885842e-03 198.17 + 191.14 839.14 2.558674e+02 6.763263e-03 197.87 + 191.14 839.14 2.509136e+02 6.643058e-03 197.56 + 191.14 839.14 2.460558e+02 6.525179e-03 197.26 + 191.14 839.14 2.412920e+02 6.409583e-03 196.95 + 191.14 839.14 2.366204e+02 6.296224e-03 196.65 + 191.14 839.14 2.320393e+02 6.185061e-03 196.34 + 191.14 839.14 2.275469e+02 6.076050e-03 196.04 + 191.14 839.14 2.231415e+02 5.969149e-03 195.73 + 191.14 839.14 2.187314e+02 5.859995e-03 195.42 + 191.14 839.14 2.143450e+02 5.749932e-03 195.10 + 191.14 839.14 2.100465e+02 5.642075e-03 194.79 + 191.14 839.14 2.058342e+02 5.536382e-03 194.47 + 191.14 839.14 2.017064e+02 5.432808e-03 194.15 + 191.14 839.14 1.976614e+02 5.331311e-03 193.84 + 191.14 839.14 1.936975e+02 5.231849e-03 193.52 + 191.14 839.14 1.898131e+02 5.134382e-03 193.20 + 191.14 839.14 1.860066e+02 5.038870e-03 192.88 + 191.14 839.14 1.822764e+02 4.945273e-03 192.57 + 191.14 839.14 1.786210e+02 4.853553e-03 192.25 + 191.14 839.14 1.750390e+02 4.763673e-03 191.93 + 191.14 839.14 1.715287e+02 4.675594e-03 191.62 + 191.14 839.14 1.680889e+02 4.589282e-03 191.30 + 191.14 839.14 1.647180e+02 4.504702e-03 190.98 + 191.14 839.14 1.614148e+02 4.421817e-03 190.67 + 191.14 839.14 1.581778e+02 4.340594e-03 190.35 + 191.14 839.14 1.550057e+02 4.261001e-03 190.03 + 191.14 839.14 1.518972e+02 4.183003e-03 189.71 + 191.14 839.14 1.488511e+02 4.106570e-03 189.40 + 191.14 839.14 1.458660e+02 4.031669e-03 189.08 + 191.14 839.14 1.429408e+02 3.958271e-03 188.76 + 191.14 839.14 1.400743e+02 3.886344e-03 188.45 + 191.14 839.14 1.372652e+02 3.815860e-03 188.13 + 191.14 839.14 1.345125e+02 3.746789e-03 187.81 + 191.14 839.14 1.318150e+02 3.679104e-03 187.50 + 191.14 839.14 1.291146e+02 3.610238e-03 187.16 + 191.14 839.14 1.264022e+02 3.539765e-03 186.79 + 191.14 839.14 1.237467e+02 3.470772e-03 186.43 + 191.14 839.14 1.211470e+02 3.403229e-03 186.06 + 191.14 839.14 1.186019e+02 3.337105e-03 185.70 Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/LES_DDphoenix/launch.ciclad =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/LES_DDphoenix/launch.ciclad (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/LES_DDphoenix/launch.ciclad (revision 1291) @@ -0,0 +1,20 @@ +#PBS -S /bin/bash +#PBS -j oe +#PBS -q infini +#PBS -N BIGLES_DD +#PBS -l "nodes=1:ppn=16" +#PBS -l "mem=64gb" + +####### prep run +cd $PBS_O_WORKDIR +\rm *rsl.* 2> /dev/null +\rm wrfo* 2> /dev/null + +## LE LANCER EN INTERACTIF +#\rm wrfi* 2> /dev/null +#./ideal.exe +#mv rsl.out.0000 ideal_rsl.out.0000 + +####### run +ulimit -s unlimited +mpirun wrf.exe Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/LES_DDphoenix/namelist.input =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/LES_DDphoenix/namelist.input (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/LES_DDphoenix/namelist.input (revision 1291) @@ -0,0 +1,76 @@ + &time_control + run_seconds = 37000 !! Duration of simulation [LES: !=0] + history_interval_s = 50 !! Frequency of outputs in seconds [LES: !=0] + frames_per_outfile = 74 !! Size of time dimension in files + start_year = 9999 !! Idealized simulation [LES: 9999] + end_year = 9999 !! Idealized simulation [LES: 9999] + restart = .false. !! (*) Output restart files ? + 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. !!12000. !! Altitude above ground of model top [LES: !=0] + time_step = 0 !! Dynamical timestep + time_step_fract_num = 1 !1 !!1 !! 3 !! Additional fraction to time_step: numerator + time_step_fract_den = 8 !10 !!4 !! 4 !! Additional fraction to time_step: denominator + dx = 10 !!25 !!50 !! (p2) Horizontal resolution + dy = 10 !!25 !!50 !! (p2) Horizontal resolution (should be equal to dx) + e_we = 369 !! (r)(p2) Number of longitude grid points + e_sn = 369 !! (r)(p2) Number of latitude grid points + e_vert = 101 !!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 = 200. !! Define constant thermal inertia value [LES: !=0] + init_AL = 0.2 !! Define constant albedo value [LES: !=0] + !! + !! OPTIONAL + !! + radt = 320 !!160 !!40 !! 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 + / + + &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.10 !! 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/tests_highresles/highresLES/callphys.def =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/highresLES/callphys.def (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/highresLES/callphys.def (revision 1291) @@ -0,0 +1,122 @@ +##General options +##~~~~~~~~~~~~~~~ + +callrichsl=.true. + +# datafile is in: +datadir = /san/work/colaitis/datafile/ + +#Run with or without tracer transport ? +tracer= .false. + +#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" ? +stats = .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 (same as =24 for now) +# =24 Mars Year 24 from TES assimilation (ie: MCD reference case) +# =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm) +# =26 Mars Year 26 from TES assimilation +iaervar = 1 +# Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case) +tauvis=0.3 +# 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 = .false. +# call CO2 NIR absorption ? matters only if callrad=T +callnirco2 = .false. +# call turbulent vertical diffusion ? +calldifv = .true. +# call convective adjustment ? +calladj = .false. +# call CO2 condensation ? +callcond =.false. +# call thermal conduction in the soil ? +callsoil = .true. +# call Lott's gravity wave/subgrid topography scheme ? +calllott = .false. + +## 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 = 0 +# DUST: Radiatively active dust ? (matters if dustbin>0) +active = .false. +# DUST: use mass and number mixing ratios to predict dust size ? +# (must also have dustbin=1) +doubleq = .false. +# DUST: lifted by GCM surface winds ? +lifting = .false. +# DUST: lifted by dust devils ? +callddevil = .false. +# DUST: Scavenging by CO2 snowfall ? +scavenging = .false. +# DUST/WATERICE: Gravitationnal sedimentation ? +sedimentation = .false. +# WATERICE: Radiatively active transported atmospheric water ice ? +activice = .false. +# WATER: Compute water cycle +water = .false. +# 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 = .false. +#caps = .true. +# PHOTOCHEMISTRY: include chemical species +photochem = .false. + +## Thermospheric options (relevant if tracer=T) : +##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call thermosphere ? +callthermos = .false. +# WATER: included without cycle (only if water=.false.) +thermoswater = .false. +# call thermal conduction ? (only if callthermos=.true.) +callconduct = .false. +# call EUV heating ? (only if callthermos=.true.) +calleuv=.false. +# call molecular viscosity ? (only if callthermos=.true.) +callmolvis = .false. +# call molecular diffusion ? (only if callthermos=.true.) +callmoldiff = .false. +# call thermospheric photochemistry ? (only if callthermos=.true.) +thermochem = .false. +# date for solar flux calculation: (1985 < date < 2002) +## (Solar min=1996.4 ave=1993.4 max=1990.6) +solarcondate = 1993.4 + + + Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/highresLES/input_coord =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/highresLES/input_coord (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/highresLES/input_coord (revision 1291) @@ -0,0 +1,6 @@ +353.87 lon +-1.88 lat +244. ls +10. lct +20. ztop +100. nz Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/highresLES/input_more =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/highresLES/input_more (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/highresLES/input_more (revision 1291) @@ -0,0 +1,1 @@ + -1412.19 264.84 Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/highresLES/input_sounding =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/highresLES/input_sounding (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/highresLES/input_sounding (revision 1291) @@ -0,0 +1,108 @@ + 7.16 254.21 0.00 + -1412.19 254.21 0.00 0.00 0.00 + -1411.19 248.72 0.00 -3.51 -0.13 + -1407.19 226.77 0.00 -4.74 -0.17 + -1402.19 218.53 0.00 -5.14 -0.19 + -1392.19 218.13 0.00 -5.70 -0.21 + -1362.19 217.99 0.00 -6.47 -0.28 + -1312.19 217.99 0.00 -7.11 -0.40 + -1212.19 217.99 0.00 -7.89 -0.66 + -1012.19 217.99 0.00 -9.77 -1.55 + -812.19 218.03 0.00 -12.23 -2.78 + -612.19 218.92 0.00 -13.90 -3.70 + -412.19 219.82 0.00 -15.56 -4.62 + -212.19 220.79 0.00 -16.47 -5.23 + -12.19 221.93 0.00 -15.30 -4.98 + 187.81 223.09 0.00 -14.14 -4.73 + 387.81 224.24 0.00 -12.98 -4.48 + 587.81 225.41 0.00 -11.81 -4.23 + 787.81 226.53 0.00 -10.73 -3.88 + 987.81 227.53 0.00 -9.86 -3.23 + 1187.81 228.54 0.00 -9.00 -2.58 + 1387.81 229.55 0.00 -8.14 -1.93 + 1587.81 230.56 0.00 -7.28 -1.28 + 1787.81 231.58 0.00 -6.42 -0.64 + 1987.81 232.60 0.00 -5.56 0.01 + 2187.81 233.63 0.00 -4.70 0.66 + 2387.81 234.66 0.00 -3.84 1.31 + 2587.81 235.70 0.00 -3.28 1.87 + 2787.81 236.72 0.00 -3.65 2.17 + 2987.81 237.75 0.00 -4.01 2.47 + 3187.81 238.78 0.00 -4.37 2.76 + 3387.81 239.82 0.00 -4.73 3.06 + 3587.81 240.86 0.00 -5.09 3.36 + 3787.81 241.90 0.00 -5.45 3.65 + 3987.81 242.95 0.00 -5.81 3.95 + 4187.81 244.01 0.00 -6.17 4.24 + 4387.81 245.07 0.00 -6.53 4.54 + 4587.81 246.13 0.00 -6.89 4.84 + 4787.81 247.20 0.00 -7.25 5.13 + 4987.81 248.27 0.00 -7.61 5.43 + 5187.81 249.35 0.00 -7.97 5.72 + 5387.81 250.36 0.00 -7.91 5.48 + 5587.81 251.37 0.00 -7.78 5.17 + 5787.81 252.37 0.00 -7.66 4.87 + 5987.81 253.39 0.00 -7.54 4.56 + 6187.81 254.40 0.00 -7.42 4.25 + 6387.81 255.42 0.00 -7.30 3.94 + 6587.81 256.45 0.00 -7.18 3.63 + 6787.81 257.47 0.00 -7.06 3.32 + 6987.81 258.51 0.00 -6.94 3.01 + 7187.81 259.54 0.00 -6.81 2.70 + 7387.81 260.58 0.00 -6.69 2.39 + 7587.81 261.63 0.00 -6.57 2.08 + 7787.81 262.67 0.00 -6.45 1.77 + 7987.81 263.73 0.00 -6.33 1.46 + 8187.81 264.78 0.00 -6.21 1.15 + 8387.81 265.84 0.00 -6.09 0.85 + 8587.81 266.91 0.00 -5.96 0.54 + 8787.81 267.98 0.00 -5.84 0.23 + 8987.81 269.04 0.00 -5.71 -0.07 + 9187.81 269.94 0.00 -5.30 -0.23 + 9387.81 270.84 0.00 -4.89 -0.38 + 9587.81 271.75 0.00 -4.48 -0.53 + 9787.81 272.65 0.00 -4.07 -0.68 + 9987.81 273.56 0.00 -3.66 -0.83 + 10187.81 274.47 0.00 -3.25 -0.99 + 10387.81 275.39 0.00 -2.84 -1.14 + 10587.81 276.30 0.00 -2.43 -1.29 + 10787.81 277.22 0.00 -2.02 -1.44 + 10987.81 278.14 0.00 -1.61 -1.59 + 11187.81 279.07 0.00 -1.20 -1.75 + 11387.81 280.00 0.00 -0.80 -1.90 + 11587.81 280.93 0.00 -0.39 -2.05 + 11787.81 281.86 0.00 0.02 -2.20 + 11987.81 282.79 0.00 0.43 -2.35 + 12187.81 283.73 0.00 0.84 -2.51 + 12387.81 284.67 0.00 1.25 -2.66 + 12587.81 285.61 0.00 1.66 -2.81 + 12787.81 286.56 0.00 2.07 -2.96 + 12987.81 287.51 0.00 2.48 -3.11 + 13187.81 288.46 0.00 2.89 -3.27 + 13387.81 289.41 0.00 3.30 -3.42 + 13587.81 290.37 0.00 3.71 -3.57 + 13787.81 291.33 0.00 4.12 -3.72 + 13987.81 292.26 0.00 4.06 -3.68 + 14187.81 293.18 0.00 3.81 -3.57 + 14387.81 294.10 0.00 3.56 -3.45 + 14587.81 295.03 0.00 3.32 -3.34 + 14787.81 295.95 0.00 3.07 -3.23 + 14987.81 296.88 0.00 2.82 -3.11 + 15187.81 297.81 0.00 2.58 -3.00 + 15387.81 298.74 0.00 2.33 -2.88 + 15587.81 299.68 0.00 2.08 -2.77 + 15787.81 300.62 0.00 1.84 -2.65 + 15987.81 301.56 0.00 1.59 -2.54 + 16187.81 302.50 0.00 1.34 -2.43 + 16387.81 303.44 0.00 1.10 -2.31 + 16587.81 304.39 0.00 0.85 -2.20 + 16787.81 305.34 0.00 0.60 -2.08 + 16987.81 306.29 0.00 0.36 -1.97 + 17187.81 307.24 0.00 0.11 -1.86 + 17387.81 308.19 0.00 -0.14 -1.74 + 17587.81 309.15 0.00 -0.38 -1.63 + 17787.81 310.11 0.00 -0.63 -1.51 + 17987.81 311.07 0.00 -0.88 -1.40 + 18187.81 312.04 0.00 -1.12 -1.29 + 18387.81 313.00 0.00 -1.37 -1.17 + 18587.81 313.97 0.00 -1.62 -1.06 Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/highresLES/input_therm =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/highresLES/input_therm (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/highresLES/input_therm (revision 1291) @@ -0,0 +1,107 @@ + 191.14 839.14 7.156943e+02 1.644131e-02 264.84 + 191.14 839.14 7.156329e+02 1.643990e-02 259.12 + 191.14 839.14 7.153877e+02 1.643427e-02 236.23 + 191.14 839.14 7.150811e+02 1.643886e-02 227.62 + 191.14 839.14 7.144677e+02 1.645781e-02 227.16 + 191.14 839.14 7.126275e+02 1.643700e-02 226.86 + 191.14 839.14 7.095670e+02 1.638459e-02 226.61 + 191.14 839.14 7.034742e+02 1.627987e-02 226.10 + 191.14 839.14 6.914027e+02 1.607142e-02 225.11 + 191.14 839.14 6.795093e+02 1.586280e-02 224.14 + 191.14 839.14 6.677845e+02 1.559250e-02 224.06 + 191.14 839.14 6.562619e+02 1.532687e-02 223.98 + 191.14 839.14 6.449444e+02 1.506191e-02 223.96 + 191.14 839.14 6.338391e+02 1.479088e-02 224.13 + 191.14 839.14 6.229249e+02 1.452450e-02 224.29 + 191.14 839.14 6.121988e+02 1.426272e-02 224.45 + 191.14 839.14 6.016573e+02 1.400544e-02 224.61 + 191.14 839.14 5.913035e+02 1.375578e-02 224.73 + 191.14 839.14 5.811459e+02 1.351970e-02 224.72 + 191.14 839.14 5.711629e+02 1.328768e-02 224.72 + 191.14 839.14 5.613513e+02 1.305965e-02 224.71 + 191.14 839.14 5.517083e+02 1.283554e-02 224.70 + 191.14 839.14 5.422309e+02 1.261527e-02 224.69 + 191.14 839.14 5.329163e+02 1.239879e-02 224.68 + 191.14 839.14 5.237617e+02 1.218602e-02 224.68 + 191.14 839.14 5.147644e+02 1.197692e-02 224.67 + 191.14 839.14 5.059217e+02 1.177166e-02 224.66 + 191.14 839.14 4.972308e+02 1.157073e-02 224.63 + 191.14 839.14 4.886892e+02 1.137325e-02 224.61 + 191.14 839.14 4.802943e+02 1.117916e-02 224.58 + 191.14 839.14 4.720437e+02 1.098841e-02 224.56 + 191.14 839.14 4.639348e+02 1.080093e-02 224.53 + 191.14 839.14 4.559651e+02 1.061667e-02 224.51 + 191.14 839.14 4.481324e+02 1.043558e-02 224.48 + 191.14 839.14 4.404343e+02 1.025760e-02 224.46 + 191.14 839.14 4.328684e+02 1.008267e-02 224.43 + 191.14 839.14 4.254324e+02 9.910755e-03 224.41 + 191.14 839.14 4.181242e+02 9.741791e-03 224.38 + 191.14 839.14 4.109416e+02 9.575727e-03 224.36 + 191.14 839.14 4.038823e+02 9.412517e-03 224.33 + 191.14 839.14 3.969221e+02 9.253722e-03 224.24 + 191.14 839.14 3.900792e+02 9.097873e-03 224.14 + 191.14 839.14 3.833541e+02 8.944712e-03 224.04 + 191.14 839.14 3.767451e+02 8.794190e-03 223.93 + 191.14 839.14 3.702499e+02 8.646264e-03 223.83 + 191.14 839.14 3.638667e+02 8.500888e-03 223.73 + 191.14 839.14 3.575936e+02 8.358018e-03 223.63 + 191.14 839.14 3.514287e+02 8.217611e-03 223.52 + 191.14 839.14 3.453700e+02 8.079626e-03 223.42 + 191.14 839.14 3.394158e+02 7.944019e-03 223.32 + 191.14 839.14 3.335642e+02 7.810749e-03 223.22 + 191.14 839.14 3.278135e+02 7.679777e-03 223.11 + 191.14 839.14 3.221619e+02 7.551063e-03 223.01 + 191.14 839.14 3.166078e+02 7.424568e-03 222.91 + 191.14 839.14 3.111494e+02 7.300254e-03 222.81 + 191.14 839.14 3.057852e+02 7.178084e-03 222.70 + 191.14 839.14 3.005134e+02 7.058019e-03 222.60 + 191.14 839.14 2.953325e+02 6.940025e-03 222.50 + 191.14 839.14 2.902368e+02 6.824114e-03 222.39 + 191.14 839.14 2.851480e+02 6.711194e-03 222.12 + 191.14 839.14 2.801484e+02 6.600253e-03 221.85 + 191.14 839.14 2.752365e+02 6.491258e-03 221.59 + 191.14 839.14 2.704107e+02 6.384173e-03 221.32 + 191.14 839.14 2.656696e+02 6.278966e-03 221.05 + 191.14 839.14 2.610115e+02 6.175603e-03 220.78 + 191.14 839.14 2.564352e+02 6.074053e-03 220.52 + 191.14 839.14 2.519390e+02 5.974284e-03 220.25 + 191.14 839.14 2.475217e+02 5.876264e-03 219.98 + 191.14 839.14 2.431819e+02 5.779962e-03 219.71 + 191.14 839.14 2.389181e+02 5.685349e-03 219.45 + 191.14 839.14 2.347291e+02 5.592394e-03 219.18 + 191.14 839.14 2.306136e+02 5.501069e-03 218.91 + 191.14 839.14 2.265702e+02 5.411346e-03 218.65 + 191.14 839.14 2.225977e+02 5.323195e-03 218.38 + 191.14 839.14 2.186948e+02 5.236591e-03 218.11 + 191.14 839.14 2.148604e+02 5.151505e-03 217.84 + 191.14 839.14 2.110932e+02 5.067910e-03 217.58 + 191.14 839.14 2.073921e+02 4.985781e-03 217.31 + 191.14 839.14 2.037558e+02 4.905093e-03 217.04 + 191.14 839.14 2.001833e+02 4.825818e-03 216.77 + 191.14 839.14 1.966735e+02 4.747934e-03 216.51 + 191.14 839.14 1.932252e+02 4.671416e-03 216.24 + 191.14 839.14 1.898373e+02 4.596239e-03 215.97 + 191.14 839.14 1.864248e+02 4.519545e-03 215.66 + 191.14 839.14 1.830415e+02 4.443143e-03 215.32 + 191.14 839.14 1.797197e+02 4.368127e-03 214.99 + 191.14 839.14 1.764581e+02 4.294473e-03 214.65 + 191.14 839.14 1.732557e+02 4.222156e-03 214.32 + 191.14 839.14 1.701115e+02 4.151151e-03 213.98 + 191.14 839.14 1.670243e+02 4.081434e-03 213.65 + 191.14 839.14 1.639931e+02 4.012983e-03 213.31 + 191.14 839.14 1.610169e+02 3.945774e-03 212.98 + 191.14 839.14 1.580948e+02 3.879785e-03 212.64 + 191.14 839.14 1.552257e+02 3.814993e-03 212.31 + 191.14 839.14 1.524086e+02 3.751378e-03 211.97 + 191.14 839.14 1.496427e+02 3.688917e-03 211.64 + 191.14 839.14 1.469270e+02 3.627589e-03 211.30 + 191.14 839.14 1.442605e+02 3.567374e-03 210.97 + 191.14 839.14 1.416425e+02 3.508252e-03 210.63 + 191.14 839.14 1.390719e+02 3.450203e-03 210.30 + 191.14 839.14 1.365480e+02 3.393208e-03 209.97 + 191.14 839.14 1.340699e+02 3.337247e-03 209.63 + 191.14 839.14 1.316368e+02 3.282301e-03 209.30 + 191.14 839.14 1.292478e+02 3.228353e-03 208.96 + 191.14 839.14 1.269023e+02 3.175383e-03 208.63 + 191.14 839.14 1.245992e+02 3.123375e-03 208.29 + 191.14 839.14 1.223380e+02 3.072311e-03 207.96 Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/highresLES/namelist.input =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/highresLES/namelist.input (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/highresLES/namelist.input (revision 1291) @@ -0,0 +1,85 @@ + &time_control + run_seconds = 37000 !! Duration of simulation [LES: !=0] + history_interval_s = 50 !! Frequency of outputs in seconds [LES: !=0] + frames_per_outfile = 74 !! Size of time dimension in files + start_year = 9999 !! Idealized simulation [LES: 9999] + start_month = 01 + start_day = 01 + start_hour = 01 + start_minute = 40 + start_second = 00 + end_year = 9999 !! Idealized simulation [LES: 9999] + restart = .true. !! (*) restart run ? + restart_interval = 10 !! (*) 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 = 5000. !! 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 = 12 !15 !20 !15 !! Additional fraction to time_step: denominator + dx = 12 !! (p2) Horizontal resolution + dy = 12 !! (p2) Horizontal resolution (should be equal to dx) + e_we = 249 !! (r)(p2) Number of longitude grid points + e_sn = 249 !! (r)(p2) Number of latitude grid points + e_vert = 401 !! (r)(p2) Number of vertical levels + !!!!! OPTIONAL + force_sfc_in_vinterp = 1 !! (p3) Number of levels hardwired in the PBL + / + + &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 = 40 !! 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 = 5. + / + + &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/tests_highresles/highresLES/run.def =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/highresLES/run.def (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_highresles/highresLES/run.def (revision 1291) @@ -0,0 +1,2 @@ +# some definitions for the physics, in file 'callphys.def' +INCLUDEDEF=callphys.def Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/LMD_MMM_d1_10km_PSFC_UV_Ls244.8.sh =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/LMD_MMM_d1_10km_PSFC_UV_Ls244.8.sh (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/LMD_MMM_d1_10km_PSFC_UV_Ls244.8.sh (revision 1291) @@ -0,0 +1,1 @@ +/donnees/aslmd/MODELES//MESOSCALE//LMD_MM_MARS//SRC/PYTHON/pp.py -f wrfout_d01_9999-09-09_09:36:00 -v PSFC -W -s 2 -i 4 -l 1. --div 25 -z 33 -m 300. -M 500. -S png Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/callphys.def =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/callphys.def (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/callphys.def (revision 1291) @@ -0,0 +1,94 @@ +General options +~~~~~~~~~~~~~~~ +tracer (Run with or without tracer transport ?) +F +diurnal (Diurnal cycle ? if diurnal=F, diurnal averaged solar heating) +T +season (Seasonal cycle ? if season=F, Ls stays constant like in "start") +F T +lwrite (want some more output on the screen ?) +F +stats (Saving statistics in file "cumul" ?) -- noWRF +F +calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF +F +Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +~~~~~~~~~~~~~ +iaervar (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario +1 (=4 Mars Year 24 from TES assimilation) +iddist (Dust vertical distribution: =0: old distrib. (Pollack90) +3 (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario ) +topdustref (Dust top altitude (km). Matter only if iddist=1) +55. +Physical Parameterizations : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +callrad (call radiative transfer ?) +T +callnlte (call NLTE radiative schemes ? matter only if callrad=T) +T +callnirco2 (call CO2 NIR absorption ? matter only if callrad=T) +T +calldifv (call turbulent vertical diffusion ?) +T +calladj (call convective adjustment ?) +T +callcond (call CO2 condensation ?) +F T +callsoil (call thermal conduction in the soil ?) +T +calllott (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF +F +Radiative transfer options : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +iradia (the rad.transfer is computed every "iradia" physical timestep) +1 37 75 +callg2d (Output of the exchange coefficient mattrix ? for diagnostic only) +F +rayleigh (Rayleigh scattering : should be =F for now) +F +Tracer (dust water, ice and/or chemical species) options (use if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +dustbin (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin)) +0 +active (DUST: Radiatively active dust ? (uses q(1) to q(dustbin)) +F +doubleq (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?) +F +lifting (DUST: lifted by GCM surface winds ?) +F +dustdevil (DUST: lifted by dust devils ?) +F +scavenging (DUST: Scavenging by CO2 snowfall ?) +F +sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?) +F +iceparty (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1)) +F +activice (WATERICE: Radiatively active transported atmospheric water ice ?) +F +water (WATER: Compute water cycle using q(nqmx) ) +F +caps (WATER: put the current permanent caps at both poles) +F +photochem (PHOTOCHEMISTRY: chemical species included) +F +Thermospheric options (relevant if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +callthermos (call thermosphere ?) +F +thermoswater (WATER: included without cycle only if water=F) +F +callconduct (call thermal conduction ? matter only if callthermos=T) +F +calleuv (call EUV heating ? matter only if callthermos=T) +F +callmolvis (call molecular viscosity ? matter only if callthermos=T) +F +callmoldiff (call molecular diffusion ? matter only if callthermos=T) +F +thermochem (call thermospheric photochemistry ? matter only if callthermos=T) +F +solarcondate (date for solar flux calculation: 1985 < date < 2002)) +1993.4 (Solar min=1996.4 ave=1993.4 max=1990.6) + + Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/dustopacity.def =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/dustopacity.def (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/dustopacity.def (revision 1291) @@ -0,0 +1,3 @@ +0.3 + +0.6 Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/input_coord =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/input_coord (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/input_coord (revision 1291) @@ -0,0 +1,6 @@ +0.0 lon +-60.0 lat +0.0 ls +09. lct +100. ztop +400. nz Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/input_more =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/input_more (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/input_more (revision 1291) @@ -0,0 +1,1 @@ + 1432.25 190.65 Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/input_sounding =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/input_sounding (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/input_sounding (revision 1291) @@ -0,0 +1,408 @@ + 5.07 199.89 0.00 + 1432.25 199.89 0.00 0.00 0.00 + 1433.25 197.46 0.00 0.49 1.32 + 1437.25 187.75 0.00 0.67 1.78 + 1442.25 187.61 0.00 0.79 1.99 + 1452.25 187.50 0.00 0.99 2.37 + 1482.25 187.82 0.00 1.06 2.91 + 1532.25 188.41 0.00 0.93 3.42 + 1632.25 189.97 0.00 0.94 3.55 + 1882.38 194.04 0.00 1.59 2.77 + 2132.50 196.97 0.00 2.37 2.51 + 2382.63 199.67 0.00 3.13 2.32 + 2632.75 201.34 0.00 3.75 2.27 + 2882.88 203.02 0.00 4.37 2.23 + 3133.00 204.71 0.00 4.99 2.18 + 3383.13 205.86 0.00 5.26 2.18 + 3633.25 206.90 0.00 5.44 2.20 + 3883.38 207.95 0.00 5.63 2.21 + 4133.50 209.00 0.00 5.82 2.23 + 4383.63 210.05 0.00 6.01 2.24 + 4633.75 211.11 0.00 6.19 2.26 + 4883.88 212.12 0.00 6.43 2.28 + 5134.00 213.07 0.00 6.72 2.30 + 5384.13 214.02 0.00 7.01 2.33 + 5634.26 214.98 0.00 7.30 2.35 + 5884.38 215.93 0.00 7.58 2.38 + 6134.51 216.90 0.00 7.87 2.40 + 6384.63 217.86 0.00 8.16 2.43 + 6634.76 218.83 0.00 8.45 2.46 + 6884.88 219.80 0.00 8.73 2.48 + 7135.01 220.89 0.00 9.14 2.43 + 7385.13 222.01 0.00 9.59 2.35 + 7635.26 223.14 0.00 10.03 2.27 + 7885.38 224.27 0.00 10.47 2.20 + 8135.51 225.41 0.00 10.92 2.12 + 8385.63 226.56 0.00 11.36 2.05 + 8635.76 227.71 0.00 11.80 1.97 + 8885.88 228.86 0.00 12.24 1.89 + 9136.01 230.02 0.00 12.69 1.82 + 9386.13 231.19 0.00 13.13 1.74 + 9636.26 232.36 0.00 13.57 1.66 + 9886.39 233.53 0.00 14.02 1.59 + 10136.51 234.79 0.00 14.48 1.52 + 10386.64 236.07 0.00 14.94 1.44 + 10636.76 237.35 0.00 15.40 1.37 + 10886.89 238.64 0.00 15.86 1.30 + 11137.01 239.93 0.00 16.32 1.23 + 11387.14 241.23 0.00 16.79 1.16 + 11637.25 242.53 0.00 17.25 1.09 + 11887.35 243.85 0.00 17.71 1.01 + 12137.55 245.16 0.00 18.17 0.94 + 12387.65 246.49 0.00 18.63 0.87 + 12637.75 247.82 0.00 19.10 0.80 + 12887.85 249.16 0.00 19.56 0.73 + 13138.05 250.50 0.00 20.02 0.66 + 13388.15 251.85 0.00 20.48 0.58 + 13638.25 253.21 0.00 20.94 0.51 + 13888.35 254.55 0.00 21.32 0.44 + 14138.55 255.88 0.00 21.67 0.38 + 14388.65 257.22 0.00 22.03 0.31 + 14638.75 258.57 0.00 22.38 0.24 + 14888.85 259.93 0.00 22.74 0.17 + 15139.05 261.29 0.00 23.09 0.11 + 15389.15 262.65 0.00 23.44 0.04 + 15639.25 264.02 0.00 23.80 -0.03 + 15889.35 265.40 0.00 24.15 -0.10 + 16139.55 266.79 0.00 24.51 -0.17 + 16389.65 268.18 0.00 24.86 -0.23 + 16639.75 269.57 0.00 25.21 -0.30 + 16889.85 270.97 0.00 25.57 -0.37 + 17140.05 272.38 0.00 25.92 -0.44 + 17390.15 273.79 0.00 26.28 -0.50 + 17640.25 275.21 0.00 26.63 -0.57 + 17890.35 276.64 0.00 26.99 -0.64 + 18140.55 278.07 0.00 27.34 -0.71 + 18390.65 279.53 0.00 27.73 -0.79 + 18640.75 281.00 0.00 28.11 -0.88 + 18890.85 282.48 0.00 28.49 -0.96 + 19141.05 283.96 0.00 28.88 -1.04 + 19391.15 285.44 0.00 29.26 -1.13 + 19641.25 286.93 0.00 29.64 -1.21 + 19891.35 288.43 0.00 30.03 -1.30 + 20141.55 289.94 0.00 30.41 -1.38 + 20391.65 291.45 0.00 30.80 -1.46 + 20641.75 292.96 0.00 31.18 -1.55 + 20891.85 294.48 0.00 31.56 -1.63 + 21142.05 296.01 0.00 31.95 -1.71 + 21392.15 297.55 0.00 32.33 -1.80 + 21642.25 299.09 0.00 32.71 -1.88 + 21892.45 300.63 0.00 33.10 -1.97 + 22142.55 302.19 0.00 33.48 -2.05 + 22392.65 303.75 0.00 33.87 -2.13 + 22642.75 305.31 0.00 34.25 -2.22 + 22892.95 306.89 0.00 34.63 -2.30 + 23143.05 308.82 0.00 35.07 -2.43 + 23393.15 310.76 0.00 35.50 -2.56 + 23643.25 312.71 0.00 35.94 -2.69 + 23893.45 314.67 0.00 36.37 -2.82 + 24143.55 316.64 0.00 36.81 -2.95 + 24393.65 318.63 0.00 37.24 -3.08 + 24643.75 320.62 0.00 37.68 -3.22 + 24893.95 322.63 0.00 38.11 -3.35 + 25144.05 324.65 0.00 38.55 -3.48 + 25394.15 326.68 0.00 38.98 -3.61 + 25644.25 328.73 0.00 39.41 -3.74 + 25894.45 330.78 0.00 39.85 -3.87 + 26144.55 332.85 0.00 40.28 -4.00 + 26394.65 334.93 0.00 40.72 -4.13 + 26644.75 337.02 0.00 41.15 -4.26 + 26894.95 339.12 0.00 41.59 -4.39 + 27145.05 341.23 0.00 42.02 -4.52 + 27395.15 343.36 0.00 42.46 -4.65 + 27645.25 345.50 0.00 42.89 -4.78 + 27895.45 347.94 0.00 43.33 -4.92 + 28145.55 350.67 0.00 43.76 -5.07 + 28395.65 353.43 0.00 44.19 -5.21 + 28645.75 356.20 0.00 44.63 -5.36 + 28895.95 359.00 0.00 45.06 -5.51 + 29146.05 361.81 0.00 45.50 -5.65 + 29396.15 364.65 0.00 45.93 -5.80 + 29646.25 367.51 0.00 46.37 -5.95 + 29896.45 370.40 0.00 46.80 -6.09 + 30146.55 373.31 0.00 47.24 -6.24 + 30396.65 376.24 0.00 47.67 -6.38 + 30646.75 379.19 0.00 48.11 -6.53 + 30896.95 382.17 0.00 48.54 -6.68 + 31147.05 385.16 0.00 48.98 -6.82 + 31397.15 388.19 0.00 49.41 -6.97 + 31647.25 391.23 0.00 49.85 -7.12 + 31897.45 394.30 0.00 50.28 -7.26 + 32147.55 397.40 0.00 50.72 -7.41 + 32397.65 400.52 0.00 51.15 -7.56 + 32647.75 403.66 0.00 51.59 -7.70 + 32897.95 407.32 0.00 52.07 -7.81 + 33148.05 411.07 0.00 52.57 -7.91 + 33398.15 414.86 0.00 53.06 -8.01 + 33648.25 418.67 0.00 53.56 -8.11 + 33898.45 422.53 0.00 54.05 -8.21 + 34148.55 426.42 0.00 54.55 -8.32 + 34398.65 430.34 0.00 55.04 -8.42 + 34648.75 434.30 0.00 55.53 -8.52 + 34898.95 438.30 0.00 56.03 -8.62 + 35149.05 442.33 0.00 56.52 -8.72 + 35399.15 446.39 0.00 57.02 -8.82 + 35649.25 450.50 0.00 57.51 -8.93 + 35899.45 454.64 0.00 58.01 -9.03 + 36149.55 458.82 0.00 58.50 -9.13 + 36399.65 463.04 0.00 59.00 -9.23 + 36649.75 467.29 0.00 59.49 -9.33 + 36899.95 471.59 0.00 59.98 -9.43 + 37150.05 475.92 0.00 60.48 -9.54 + 37400.15 480.29 0.00 60.97 -9.64 + 37650.25 484.70 0.00 61.47 -9.74 + 37900.45 489.27 0.00 62.03 -9.79 + 38150.55 493.99 0.00 62.65 -9.81 + 38400.65 498.76 0.00 63.27 -9.83 + 38650.75 503.56 0.00 63.89 -9.84 + 38900.95 508.41 0.00 64.51 -9.86 + 39151.05 513.31 0.00 65.14 -9.88 + 39401.15 518.25 0.00 65.76 -9.89 + 39651.25 523.24 0.00 66.38 -9.91 + 39901.45 528.27 0.00 67.00 -9.93 + 40151.55 533.35 0.00 67.62 -9.94 + 40401.65 538.48 0.00 68.24 -9.96 + 40651.75 543.65 0.00 68.86 -9.98 + 40901.95 548.87 0.00 69.48 -9.99 + 41152.05 554.14 0.00 70.10 -10.01 + 41402.15 559.46 0.00 70.72 -10.03 + 41652.35 564.83 0.00 71.34 -10.04 + 41902.45 570.25 0.00 71.97 -10.06 + 42152.55 575.72 0.00 72.59 -10.08 + 42402.65 581.23 0.00 73.21 -10.09 + 42652.85 586.80 0.00 73.83 -10.11 + 42902.95 592.42 0.00 74.45 -10.13 + 43153.05 598.10 0.00 75.07 -10.14 + 43403.15 603.73 0.00 75.65 -10.14 + 43653.35 608.99 0.00 76.07 -10.03 + 43903.45 614.29 0.00 76.48 -9.92 + 44153.55 619.64 0.00 76.89 -9.81 + 44403.65 625.03 0.00 77.31 -9.70 + 44653.85 630.47 0.00 77.72 -9.58 + 44903.95 635.95 0.00 78.14 -9.47 + 45154.05 641.48 0.00 78.55 -9.36 + 45404.15 647.06 0.00 78.96 -9.25 + 45654.35 652.69 0.00 79.38 -9.14 + 45904.45 658.36 0.00 79.79 -9.03 + 46154.55 664.09 0.00 80.20 -8.92 + 46404.65 669.86 0.00 80.62 -8.81 + 46654.85 675.68 0.00 81.03 -8.70 + 46904.95 681.55 0.00 81.45 -8.59 + 47155.05 687.47 0.00 81.86 -8.48 + 47405.15 693.44 0.00 82.27 -8.37 + 47655.35 699.46 0.00 82.69 -8.25 + 47905.45 705.53 0.00 83.10 -8.14 + 48155.55 711.65 0.00 83.51 -8.03 + 48405.65 717.83 0.00 83.93 -7.92 + 48655.85 724.06 0.00 84.34 -7.81 + 48905.95 730.34 0.00 84.75 -7.70 + 49156.05 736.67 0.00 85.17 -7.59 + 49406.15 743.06 0.00 85.58 -7.48 + 49656.35 748.44 0.00 85.61 -7.31 + 49906.45 753.78 0.00 85.62 -7.13 + 50156.55 759.16 0.00 85.63 -6.96 + 50406.65 764.58 0.00 85.64 -6.78 + 50656.85 770.05 0.00 85.65 -6.61 + 50906.95 775.55 0.00 85.66 -6.43 + 51157.05 781.09 0.00 85.67 -6.26 + 51407.15 786.66 0.00 85.68 -6.08 + 51657.35 792.28 0.00 85.69 -5.91 + 51907.45 797.94 0.00 85.70 -5.73 + 52157.55 803.64 0.00 85.71 -5.56 + 52407.65 809.38 0.00 85.72 -5.38 + 52657.85 815.16 0.00 85.73 -5.21 + 52907.95 820.98 0.00 85.74 -5.03 + 53158.05 826.84 0.00 85.75 -4.86 + 53408.15 832.75 0.00 85.76 -4.68 + 53658.35 838.70 0.00 85.77 -4.51 + 53908.45 844.69 0.00 85.78 -4.33 + 54158.55 850.72 0.00 85.79 -4.16 + 54408.65 856.80 0.00 85.80 -3.98 + 54658.85 862.92 0.00 85.81 -3.81 + 54908.95 869.08 0.00 85.82 -3.63 + 55159.05 875.28 0.00 85.83 -3.46 + 55409.15 881.53 0.00 85.84 -3.28 + 55659.35 887.51 0.00 85.79 -3.21 + 55909.45 892.76 0.00 85.58 -3.38 + 56159.55 898.05 0.00 85.38 -3.54 + 56409.65 903.36 0.00 85.18 -3.70 + 56659.85 908.71 0.00 84.97 -3.87 + 56909.95 914.08 0.00 84.77 -4.03 + 57160.05 919.49 0.00 84.56 -4.20 + 57410.15 924.92 0.00 84.36 -4.36 + 57660.35 930.38 0.00 84.16 -4.53 + 57910.45 935.88 0.00 83.95 -4.69 + 58160.55 941.40 0.00 83.75 -4.86 + 58410.65 946.96 0.00 83.55 -5.02 + 58660.85 952.54 0.00 83.34 -5.19 + 58910.95 958.16 0.00 83.14 -5.35 + 59161.05 963.80 0.00 82.93 -5.52 + 59411.15 969.48 0.00 82.73 -5.68 + 59661.35 975.19 0.00 82.53 -5.84 + 59911.45 980.93 0.00 82.32 -6.01 + 60161.55 986.70 0.00 82.12 -6.17 + 60411.65 992.50 0.00 81.92 -6.34 + 60661.85 998.34 0.00 81.71 -6.50 + 60911.95 1004.21 0.00 81.51 -6.67 + 61162.05 1010.11 0.00 81.30 -6.83 + 61412.15 1016.04 0.00 81.10 -7.00 + 61662.35 1021.97 0.00 80.89 -7.14 + 61912.45 1027.52 0.00 80.62 -7.01 + 62162.55 1033.09 0.00 80.35 -6.88 + 62412.75 1038.69 0.00 80.08 -6.74 + 62662.85 1044.31 0.00 79.81 -6.61 + 62912.95 1049.96 0.00 79.54 -6.48 + 63163.05 1055.63 0.00 79.27 -6.35 + 63413.25 1061.33 0.00 79.00 -6.21 + 63663.35 1067.05 0.00 78.73 -6.08 + 63913.45 1072.80 0.00 78.46 -5.95 + 64163.55 1078.57 0.00 78.19 -5.82 + 64413.75 1084.37 0.00 77.92 -5.69 + 64663.85 1090.19 0.00 77.65 -5.55 + 64913.95 1096.04 0.00 77.38 -5.42 + 65164.05 1101.91 0.00 77.11 -5.29 + 65414.25 1107.81 0.00 76.84 -5.16 + 65664.35 1113.73 0.00 76.57 -5.03 + 65914.45 1119.68 0.00 76.30 -4.89 + 66164.55 1125.65 0.00 76.03 -4.76 + 66414.75 1131.65 0.00 75.75 -4.63 + 66664.85 1137.67 0.00 75.48 -4.50 + 66914.95 1143.72 0.00 75.21 -4.36 + 67165.05 1149.79 0.00 74.94 -4.23 + 67415.25 1155.90 0.00 74.67 -4.10 + 67665.35 1163.47 0.00 74.04 -3.93 + 67915.45 1171.30 0.00 73.34 -3.75 + 68165.55 1179.19 0.00 72.65 -3.57 + 68415.75 1187.13 0.00 71.96 -3.39 + 68665.85 1195.12 0.00 71.26 -3.21 + 68915.95 1203.16 0.00 70.57 -3.03 + 69166.05 1211.25 0.00 69.88 -2.85 + 69416.25 1219.40 0.00 69.19 -2.67 + 69666.35 1227.60 0.00 68.49 -2.49 + 69916.45 1235.85 0.00 67.80 -2.31 + 70166.55 1244.15 0.00 67.11 -2.14 + 70416.75 1252.51 0.00 66.42 -1.96 + 70666.85 1260.92 0.00 65.72 -1.78 + 70916.95 1269.39 0.00 65.03 -1.60 + 71167.05 1277.91 0.00 64.34 -1.42 + 71417.25 1286.48 0.00 63.65 -1.24 + 71667.35 1295.11 0.00 62.95 -1.06 + 71917.45 1303.80 0.00 62.26 -0.88 + 72167.55 1312.54 0.00 61.57 -0.70 + 72417.75 1321.34 0.00 60.87 -0.52 + 72667.85 1330.19 0.00 60.18 -0.34 + 72917.95 1339.10 0.00 59.49 -0.16 + 73168.05 1348.54 0.00 58.61 -0.62 + 73418.25 1358.34 0.00 57.61 -1.48 + 73668.35 1368.20 0.00 56.62 -2.33 + 73918.45 1378.14 0.00 55.62 -3.19 + 74168.55 1388.15 0.00 54.63 -4.05 + 74418.75 1398.24 0.00 53.63 -4.90 + 74668.85 1408.39 0.00 52.63 -5.76 + 74918.95 1418.61 0.00 51.64 -6.61 + 75169.05 1428.91 0.00 50.64 -7.47 + 75419.25 1439.28 0.00 49.64 -8.33 + 75669.35 1449.72 0.00 48.65 -9.18 + 75919.45 1460.24 0.00 47.65 -10.04 + 76169.55 1470.83 0.00 46.66 -10.89 + 76419.75 1481.50 0.00 45.66 -11.75 + 76669.85 1492.25 0.00 44.66 -12.61 + 76919.95 1503.06 0.00 43.67 -13.46 + 77170.05 1513.96 0.00 42.67 -14.32 + 77420.25 1524.94 0.00 41.67 -15.17 + 77670.35 1535.99 0.00 40.68 -16.03 + 77920.45 1547.12 0.00 39.68 -16.89 + 78170.55 1558.32 0.00 38.69 -17.74 + 78420.75 1569.61 0.00 37.69 -18.60 + 78670.85 1580.35 0.00 37.07 -19.81 + 78920.95 1591.14 0.00 36.45 -21.04 + 79171.05 1602.00 0.00 35.84 -22.26 + 79421.25 1612.93 0.00 35.23 -23.49 + 79671.35 1623.92 0.00 34.62 -24.71 + 79921.45 1634.99 0.00 34.01 -25.94 + 80171.55 1646.13 0.00 33.40 -27.16 + 80421.75 1657.34 0.00 32.79 -28.39 + 80671.85 1668.62 0.00 32.18 -29.62 + 80921.95 1679.97 0.00 31.57 -30.84 + 81172.05 1691.39 0.00 30.96 -32.07 + 81422.25 1702.89 0.00 30.34 -33.29 + 81672.35 1714.46 0.00 29.73 -34.52 + 81922.45 1726.10 0.00 29.12 -35.74 + 82172.65 1737.82 0.00 28.51 -36.97 + 82422.75 1749.61 0.00 27.90 -38.19 + 82672.85 1761.47 0.00 27.29 -39.42 + 82922.95 1773.41 0.00 26.68 -40.64 + 83173.15 1785.43 0.00 26.07 -41.87 + 83423.25 1797.52 0.00 25.46 -43.10 + 83673.35 1809.69 0.00 24.84 -44.32 + 83923.45 1822.42 0.00 24.39 -45.55 + 84173.65 1835.31 0.00 23.97 -46.78 + 84423.75 1848.29 0.00 23.54 -48.01 + 84673.85 1861.36 0.00 23.12 -49.24 + 84923.95 1874.51 0.00 22.69 -50.46 + 85174.15 1887.75 0.00 22.27 -51.69 + 85424.25 1901.07 0.00 21.84 -52.92 + 85674.35 1914.49 0.00 21.42 -54.15 + 85924.45 1927.99 0.00 21.00 -55.38 + 86174.65 1941.59 0.00 20.57 -56.61 + 86424.75 1955.27 0.00 20.15 -57.84 + 86674.85 1969.04 0.00 19.72 -59.06 + 86924.95 1982.90 0.00 19.30 -60.29 + 87175.15 1996.86 0.00 18.87 -61.52 + 87425.25 2010.91 0.00 18.45 -62.75 + 87675.35 2025.05 0.00 18.02 -63.98 + 87925.45 2039.28 0.00 17.60 -65.21 + 88175.65 2053.61 0.00 17.17 -66.44 + 88425.75 2068.02 0.00 16.75 -67.66 + 88675.85 2082.54 0.00 16.32 -68.89 + 88925.95 2097.96 0.00 15.70 -69.93 + 89176.15 2114.62 0.00 14.80 -70.70 + 89426.25 2131.41 0.00 13.91 -71.47 + 89676.35 2148.32 0.00 13.01 -72.24 + 89926.45 2165.37 0.00 12.11 -73.01 + 90176.65 2182.56 0.00 11.22 -73.79 + 90426.75 2199.87 0.00 10.32 -74.56 + 90676.85 2217.32 0.00 9.42 -75.33 + 90926.95 2234.91 0.00 8.53 -76.10 + 91177.15 2252.64 0.00 7.63 -76.87 + 91427.25 2270.50 0.00 6.74 -77.64 + 91677.35 2288.50 0.00 5.84 -78.41 + 91927.45 2306.64 0.00 4.94 -79.18 + 92177.65 2324.93 0.00 4.05 -79.96 + 92427.75 2343.35 0.00 3.15 -80.73 + 92677.85 2361.92 0.00 2.25 -81.50 + 92927.95 2380.63 0.00 1.36 -82.27 + 93178.15 2399.49 0.00 0.46 -83.04 + 93428.25 2418.50 0.00 -0.43 -83.81 + 93678.35 2437.65 0.00 -1.33 -84.58 + 93928.45 2457.25 0.00 -2.45 -85.76 + 94178.65 2477.57 0.00 -3.96 -87.64 + 94428.75 2498.04 0.00 -5.46 -89.53 + 94678.85 2518.68 0.00 -6.97 -91.41 + 94928.95 2539.48 0.00 -8.48 -93.30 + 95179.15 2560.47 0.00 -9.99 -95.19 + 95429.25 2581.62 0.00 -11.50 -97.07 + 95679.35 2602.94 0.00 -13.01 -98.96 + 95929.45 2624.44 0.00 -14.52 -100.84 + 96179.65 2646.12 0.00 -16.03 -102.73 + 96429.75 2667.97 0.00 -17.54 -104.62 + 96679.85 2689.99 0.00 -19.05 -106.50 + 96929.95 2712.20 0.00 -20.55 -108.39 + 97180.15 2734.60 0.00 -22.06 -110.27 + 97430.25 2757.17 0.00 -23.57 -112.16 + 97680.35 2779.93 0.00 -25.08 -114.04 + 97930.45 2802.87 0.00 -26.59 -115.93 + 98180.65 2826.00 0.00 -28.10 -117.82 + 98430.75 2849.32 0.00 -29.61 -119.70 + 98680.85 2872.83 0.00 -31.12 -121.59 + 98930.95 2898.15 0.00 -32.26 -122.91 + 99181.15 2925.00 0.00 -33.12 -123.78 + 99431.25 2952.09 0.00 -33.97 -124.66 + 99681.35 2979.43 0.00 -34.83 -125.53 + 99931.45 3007.03 0.00 -35.68 -126.41 + 100181.60 3034.88 0.00 -36.54 -127.28 + 100431.80 3062.99 0.00 -37.39 -128.16 + 100681.90 3091.36 0.00 -38.25 -129.03 + 100932.00 3119.99 0.00 -39.10 -129.91 + 101182.10 3148.89 0.00 -39.96 -130.78 + 101432.30 3178.06 0.00 -40.81 -131.66 Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/input_therm =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/input_therm (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/input_therm (revision 1291) @@ -0,0 +1,407 @@ + 191.14 839.14 5.072664e+02 1.494174e-02 190.65 + 191.14 839.14 5.072107e+02 1.494010e-02 188.34 + 191.14 839.14 5.069883e+02 1.493355e-02 179.06 + 191.14 839.14 5.067104e+02 1.493098e-02 178.89 + 191.14 839.14 5.061549e+02 1.492278e-02 178.74 + 191.14 839.14 5.044931e+02 1.485232e-02 178.89 + 191.14 839.14 5.017410e+02 1.473463e-02 179.21 + 191.14 839.14 4.963062e+02 1.449251e-02 180.18 + 191.14 839.14 4.830946e+02 1.390212e-02 182.78 + 191.14 839.14 4.704177e+02 1.342560e-02 184.27 + 191.14 839.14 4.581017e+02 1.298002e-02 185.54 + 191.14 839.14 4.462144e+02 1.262198e-02 185.83 + 191.14 839.14 4.346354e+02 1.227323e-02 186.12 + 191.14 839.14 4.233572e+02 1.193353e-02 186.41 + 191.14 839.14 4.123535e+02 1.163460e-02 186.20 + 191.14 839.14 4.016326e+02 1.135003e-02 185.88 + 191.14 839.14 3.911900e+02 1.107286e-02 185.56 + 191.14 839.14 3.810193e+02 1.080289e-02 185.24 + 191.14 839.14 3.711127e+02 1.053994e-02 184.92 + 191.14 839.14 3.614640e+02 1.028383e-02 184.60 + 191.14 839.14 3.519988e+02 1.003410e-02 184.22 + 191.14 839.14 3.427144e+02 9.790657e-03 183.79 + 191.14 839.14 3.336746e+02 9.553624e-03 183.35 + 191.14 839.14 3.248731e+02 9.322844e-03 182.91 + 191.14 839.14 3.163042e+02 9.098161e-03 182.47 + 191.14 839.14 3.079609e+02 8.879395e-03 182.03 + 191.14 839.14 2.998380e+02 8.666408e-03 181.59 + 191.14 839.14 2.919291e+02 8.459031e-03 181.15 + 191.14 839.14 2.842291e+02 8.257132e-03 180.71 + 191.14 839.14 2.766171e+02 8.050816e-03 180.34 + 191.14 839.14 2.691763e+02 7.847213e-03 180.00 + 191.14 839.14 2.619354e+02 7.649080e-03 179.66 + 191.14 839.14 2.548895e+02 7.456284e-03 179.31 + 191.14 839.14 2.480329e+02 7.268666e-03 178.97 + 191.14 839.14 2.413610e+02 7.086103e-03 178.62 + 191.14 839.14 2.348683e+02 6.908444e-03 178.28 + 191.14 839.14 2.285505e+02 6.735570e-03 177.93 + 191.14 839.14 2.224024e+02 6.567340e-03 177.59 + 191.14 839.14 2.164199e+02 6.403642e-03 177.24 + 191.14 839.14 2.105982e+02 6.244341e-03 176.90 + 191.14 839.14 2.049330e+02 6.089326e-03 176.55 + 191.14 839.14 1.993070e+02 5.930856e-03 176.25 + 191.14 839.14 1.938262e+02 5.776122e-03 175.94 + 191.14 839.14 1.884964e+02 5.625649e-03 175.63 + 191.14 839.14 1.833129e+02 5.479309e-03 175.33 + 191.14 839.14 1.782722e+02 5.336998e-03 175.02 + 191.14 839.14 1.733699e+02 5.198595e-03 174.72 + 191.14 839.14 1.686027e+02 5.064008e-03 174.41 + 191.14 839.14 1.639669e+02 4.933127e-03 174.11 + 191.14 839.14 1.594567e+02 4.805795e-03 173.80 + 191.14 839.14 1.550723e+02 4.682014e-03 173.49 + 191.14 839.14 1.508085e+02 4.561637e-03 173.19 + 191.14 839.14 1.466619e+02 4.444569e-03 172.88 + 191.14 839.14 1.426277e+02 4.330676e-03 172.58 + 191.14 839.14 1.387060e+02 4.219959e-03 172.27 + 191.14 839.14 1.348922e+02 4.112286e-03 171.97 + 191.14 839.14 1.310953e+02 4.003646e-03 171.61 + 191.14 839.14 1.273811e+02 3.897009e-03 171.25 + 191.14 839.14 1.237737e+02 3.793434e-03 170.88 + 191.14 839.14 1.202683e+02 3.692793e-03 170.52 + 191.14 839.14 1.168623e+02 3.595002e-03 170.15 + 191.14 839.14 1.135514e+02 3.499942e-03 169.79 + 191.14 839.14 1.103356e+02 3.407612e-03 169.42 + 191.14 839.14 1.072109e+02 3.317897e-03 169.06 + 191.14 839.14 1.041747e+02 3.230723e-03 168.69 + 191.14 839.14 1.012232e+02 3.145984e-03 168.33 + 191.14 839.14 9.835655e+01 3.063678e-03 167.96 + 191.14 839.14 9.557108e+01 2.983704e-03 167.60 + 191.14 839.14 9.286446e+01 2.905994e-03 167.23 + 191.14 839.14 9.023347e+01 2.830455e-03 166.86 + 191.14 839.14 8.767804e+01 2.757085e-03 166.50 + 191.14 839.14 8.519498e+01 2.685793e-03 166.13 + 191.14 839.14 8.278223e+01 2.616521e-03 165.77 + 191.14 839.14 8.042753e+01 2.548678e-03 165.40 + 191.14 839.14 7.805846e+01 2.478316e-03 165.00 + 191.14 839.14 7.575916e+01 2.410028e-03 164.60 + 191.14 839.14 7.352759e+01 2.343750e-03 164.20 + 191.14 839.14 7.136089e+01 2.279400e-03 163.80 + 191.14 839.14 6.925889e+01 2.216970e-03 163.40 + 191.14 839.14 6.721879e+01 2.156380e-03 163.00 + 191.14 839.14 6.523878e+01 2.097574e-03 162.60 + 191.14 839.14 6.331634e+01 2.040478e-03 162.20 + 191.14 839.14 6.145129e+01 1.985086e-03 161.80 + 191.14 839.14 5.964117e+01 1.931326e-03 161.40 + 191.14 839.14 5.788438e+01 1.879149e-03 161.00 + 191.14 839.14 5.617865e+01 1.828489e-03 160.59 + 191.14 839.14 5.452385e+01 1.779342e-03 160.19 + 191.14 839.14 5.291780e+01 1.731642e-03 159.79 + 191.14 839.14 5.135844e+01 1.685329e-03 159.39 + 191.14 839.14 4.984560e+01 1.640399e-03 158.99 + 191.14 839.14 4.837735e+01 1.596792e-03 158.59 + 191.14 839.14 4.695235e+01 1.554469e-03 158.19 + 191.14 839.14 4.556747e+01 1.513262e-03 157.79 + 191.14 839.14 4.417502e+01 1.468897e-03 157.53 + 191.14 839.14 4.282513e+01 1.425888e-03 157.26 + 191.14 839.14 4.151649e+01 1.384193e-03 156.99 + 191.14 839.14 4.024734e+01 1.343756e-03 156.73 + 191.14 839.14 3.901746e+01 1.304571e-03 156.46 + 191.14 839.14 3.782518e+01 1.266583e-03 156.19 + 191.14 839.14 3.666933e+01 1.229757e-03 155.92 + 191.14 839.14 3.554835e+01 1.194041e-03 155.66 + 191.14 839.14 3.446207e+01 1.159431e-03 155.39 + 191.14 839.14 3.340898e+01 1.125878e-03 155.12 + 191.14 839.14 3.238808e+01 1.093351e-03 154.85 + 191.14 839.14 3.139799e+01 1.061806e-03 154.59 + 191.14 839.14 3.043852e+01 1.031236e-03 154.32 + 191.14 839.14 2.950839e+01 1.001601e-03 154.05 + 191.14 839.14 2.860668e+01 9.728714e-04 153.79 + 191.14 839.14 2.773218e+01 9.450089e-04 153.52 + 191.14 839.14 2.688474e+01 9.180085e-04 153.25 + 191.14 839.14 2.606321e+01 8.918334e-04 152.98 + 191.14 839.14 2.526677e+01 8.664582e-04 152.72 + 191.14 839.14 2.448795e+01 8.407207e-04 152.57 + 191.14 839.14 2.372769e+01 8.147671e-04 152.53 + 191.14 839.14 2.299104e+01 7.896196e-04 152.49 + 191.14 839.14 2.227726e+01 7.652527e-04 152.44 + 191.14 839.14 2.158537e+01 7.416331e-04 152.40 + 191.14 839.14 2.091522e+01 7.187558e-04 152.36 + 191.14 839.14 2.026589e+01 6.965889e-04 152.32 + 191.14 839.14 1.963671e+01 6.751104e-04 152.28 + 191.14 839.14 1.902683e+01 6.542904e-04 152.24 + 191.14 839.14 1.843612e+01 6.341249e-04 152.20 + 191.14 839.14 1.786375e+01 6.145854e-04 152.16 + 191.14 839.14 1.730915e+01 5.956527e-04 152.12 + 191.14 839.14 1.677156e+01 5.773006e-04 152.07 + 191.14 839.14 1.625087e+01 5.595252e-04 152.03 + 191.14 839.14 1.574635e+01 5.423019e-04 151.99 + 191.14 839.14 1.525748e+01 5.256133e-04 151.95 + 191.14 839.14 1.478362e+01 5.094364e-04 151.91 + 191.14 839.14 1.432464e+01 4.937680e-04 151.87 + 191.14 839.14 1.387991e+01 4.785862e-04 151.83 + 191.14 839.14 1.344900e+01 4.638757e-04 151.79 + 191.14 839.14 1.303591e+01 4.492990e-04 151.94 + 191.14 839.14 1.263620e+01 4.351378e-04 152.12 + 191.14 839.14 1.224875e+01 4.214111e-04 152.30 + 191.14 839.14 1.187318e+01 4.081050e-04 152.48 + 191.14 839.14 1.150898e+01 3.952019e-04 152.66 + 191.14 839.14 1.115609e+01 3.826995e-04 152.84 + 191.14 839.14 1.081403e+01 3.705806e-04 153.02 + 191.14 839.14 1.048244e+01 3.588331e-04 153.20 + 191.14 839.14 1.016091e+01 3.474414e-04 153.38 + 191.14 839.14 9.849350e+00 3.364034e-04 153.56 + 191.14 839.14 9.547352e+00 3.257040e-04 153.74 + 191.14 839.14 9.254609e+00 3.153325e-04 153.92 + 191.14 839.14 8.970732e+00 3.052751e-04 154.10 + 191.14 839.14 8.695670e+00 2.955300e-04 154.28 + 191.14 839.14 8.429046e+00 2.860839e-04 154.46 + 191.14 839.14 8.170593e+00 2.769272e-04 154.64 + 191.14 839.14 7.919968e+00 2.680479e-04 154.82 + 191.14 839.14 7.677124e+00 2.594443e-04 155.00 + 191.14 839.14 7.441730e+00 2.511046e-04 155.17 + 191.14 839.14 7.213550e+00 2.430205e-04 155.35 + 191.14 839.14 6.996545e+00 2.353695e-04 155.60 + 191.14 839.14 6.789512e+00 2.281004e-04 155.89 + 191.14 839.14 6.588607e+00 2.210466e-04 156.19 + 191.14 839.14 6.393645e+00 2.142013e-04 156.48 + 191.14 839.14 6.204379e+00 2.075561e-04 156.78 + 191.14 839.14 6.020786e+00 2.011101e-04 157.07 + 191.14 839.14 5.842629e+00 1.948549e-04 157.37 + 191.14 839.14 5.669741e+00 1.887846e-04 157.66 + 191.14 839.14 5.501904e+00 1.828918e-04 157.96 + 191.14 839.14 5.339098e+00 1.771756e-04 158.25 + 191.14 839.14 5.181112e+00 1.716286e-04 158.55 + 191.14 839.14 5.027799e+00 1.662457e-04 158.85 + 191.14 839.14 4.878964e+00 1.610201e-04 159.14 + 191.14 839.14 4.734592e+00 1.559511e-04 159.44 + 191.14 839.14 4.594494e+00 1.510321e-04 159.73 + 191.14 839.14 4.458485e+00 1.462568e-04 160.03 + 191.14 839.14 4.326556e+00 1.416247e-04 160.32 + 191.14 839.14 4.198530e+00 1.371296e-04 160.62 + 191.14 839.14 4.074294e+00 1.327676e-04 160.91 + 191.14 839.14 3.953684e+00 1.285330e-04 161.21 + 191.14 839.14 3.836693e+00 1.244253e-04 161.50 + 191.14 839.14 3.723162e+00 1.204392e-04 161.80 + 191.14 839.14 3.613731e+00 1.166319e-04 162.08 + 191.14 839.14 3.510807e+00 1.132119e-04 162.28 + 191.14 839.14 3.410856e+00 1.098907e-04 162.49 + 191.14 839.14 3.313749e+00 1.066640e-04 162.69 + 191.14 839.14 3.219409e+00 1.035292e-04 162.90 + 191.14 839.14 3.127716e+00 1.004824e-04 163.10 + 191.14 839.14 3.038671e+00 9.752356e-05 163.31 + 191.14 839.14 2.952161e+00 9.464895e-05 163.51 + 191.14 839.14 2.868114e+00 9.185622e-05 163.72 + 191.14 839.14 2.786427e+00 8.914187e-05 163.92 + 191.14 839.14 2.707098e+00 8.650593e-05 164.13 + 191.14 839.14 2.630028e+00 8.394499e-05 164.33 + 191.14 839.14 2.555152e+00 8.145699e-05 164.54 + 191.14 839.14 2.482378e+00 7.903881e-05 164.74 + 191.14 839.14 2.411706e+00 7.669050e-05 164.95 + 191.14 839.14 2.343045e+00 7.440901e-05 165.15 + 191.14 839.14 2.276340e+00 7.219250e-05 165.36 + 191.14 839.14 2.211507e+00 7.003820e-05 165.56 + 191.14 839.14 2.148546e+00 6.794612e-05 165.77 + 191.14 839.14 2.087377e+00 6.591358e-05 165.97 + 191.14 839.14 2.027951e+00 6.393893e-05 166.18 + 191.14 839.14 1.970192e+00 6.201971e-05 166.38 + 191.14 839.14 1.914102e+00 6.015591e-05 166.59 + 191.14 839.14 1.859608e+00 5.834515e-05 166.79 + 191.14 839.14 1.806666e+00 5.658597e-05 167.00 + 191.14 839.14 1.755997e+00 5.499950e-05 166.98 + 191.14 839.14 1.706817e+00 5.346544e-05 166.95 + 191.14 839.14 1.659013e+00 5.197432e-05 166.92 + 191.14 839.14 1.612549e+00 5.052498e-05 166.90 + 191.14 839.14 1.567367e+00 4.911565e-05 166.87 + 191.14 839.14 1.523469e+00 4.774638e-05 166.84 + 191.14 839.14 1.480801e+00 4.641544e-05 166.81 + 191.14 839.14 1.439328e+00 4.512179e-05 166.78 + 191.14 839.14 1.398999e+00 4.386386e-05 166.76 + 191.14 839.14 1.359817e+00 4.264167e-05 166.73 + 191.14 839.14 1.321732e+00 4.145370e-05 166.70 + 191.14 839.14 1.284714e+00 4.029902e-05 166.67 + 191.14 839.14 1.248718e+00 3.917621e-05 166.64 + 191.14 839.14 1.213745e+00 3.808532e-05 166.62 + 191.14 839.14 1.179751e+00 3.702495e-05 166.59 + 191.14 839.14 1.146709e+00 3.599432e-05 166.56 + 191.14 839.14 1.114580e+00 3.499212e-05 166.53 + 191.14 839.14 1.083364e+00 3.401841e-05 166.50 + 191.14 839.14 1.053021e+00 3.307195e-05 166.48 + 191.14 839.14 1.023529e+00 3.215202e-05 166.45 + 191.14 839.14 9.948509e-01 3.125748e-05 166.42 + 191.14 839.14 9.669880e-01 3.038836e-05 166.39 + 191.14 839.14 9.399049e-01 2.954358e-05 166.36 + 191.14 839.14 9.135808e-01 2.872247e-05 166.34 + 191.14 839.14 8.878531e-01 2.792712e-05 166.24 + 191.14 839.14 8.625717e-01 2.716162e-05 165.99 + 191.14 839.14 8.380096e-01 2.641790e-05 165.74 + 191.14 839.14 8.141472e-01 2.569537e-05 165.49 + 191.14 839.14 7.909548e-01 2.499312e-05 165.24 + 191.14 839.14 7.684324e-01 2.431116e-05 164.99 + 191.14 839.14 7.465510e-01 2.364861e-05 164.74 + 191.14 839.14 7.252931e-01 2.300494e-05 164.49 + 191.14 839.14 7.046318e-01 2.237934e-05 164.24 + 191.14 839.14 6.845675e-01 2.177180e-05 164.00 + 191.14 839.14 6.650742e-01 2.118156e-05 163.75 + 191.14 839.14 6.461362e-01 2.060814e-05 163.50 + 191.14 839.14 6.277299e-01 2.005081e-05 163.25 + 191.14 839.14 6.098553e-01 1.950959e-05 163.00 + 191.14 839.14 5.924895e-01 1.898376e-05 162.75 + 191.14 839.14 5.756184e-01 1.847292e-05 162.50 + 191.14 839.14 5.592208e-01 1.797642e-05 162.25 + 191.14 839.14 5.432971e-01 1.749426e-05 162.00 + 191.14 839.14 5.278264e-01 1.702582e-05 161.75 + 191.14 839.14 5.127966e-01 1.657073e-05 161.50 + 191.14 839.14 4.981887e-01 1.612842e-05 161.25 + 191.14 839.14 4.840029e-01 1.569888e-05 161.00 + 191.14 839.14 4.702206e-01 1.528157e-05 160.75 + 191.14 839.14 4.568312e-01 1.487614e-05 160.50 + 191.14 839.14 4.437747e-01 1.448030e-05 160.24 + 191.14 839.14 4.305899e-01 1.407449e-05 159.87 + 191.14 839.14 4.177965e-01 1.368072e-05 159.50 + 191.14 839.14 4.053785e-01 1.329851e-05 159.13 + 191.14 839.14 3.933342e-01 1.292781e-05 158.76 + 191.14 839.14 3.816480e-01 1.256812e-05 158.39 + 191.14 839.14 3.703088e-01 1.221911e-05 158.02 + 191.14 839.14 3.593023e-01 1.188035e-05 157.64 + 191.14 839.14 3.486269e-01 1.155178e-05 157.27 + 191.14 839.14 3.382690e-01 1.123297e-05 156.90 + 191.14 839.14 3.282186e-01 1.092364e-05 156.53 + 191.14 839.14 3.184631e-01 1.062337e-05 156.16 + 191.14 839.14 3.090012e-01 1.033215e-05 155.79 + 191.14 839.14 2.998206e-01 1.004958e-05 155.42 + 191.14 839.14 2.909126e-01 9.775404e-06 155.05 + 191.14 839.14 2.822659e-01 9.509270e-06 154.68 + 191.14 839.14 2.738794e-01 9.251145e-06 154.30 + 191.14 839.14 2.657421e-01 9.000691e-06 153.93 + 191.14 839.14 2.578469e-01 8.757685e-06 153.56 + 191.14 839.14 2.501829e-01 8.521797e-06 153.19 + 191.14 839.14 2.427496e-01 8.293012e-06 152.82 + 191.14 839.14 2.355373e-01 8.071024e-06 152.45 + 191.14 839.14 2.285394e-01 7.855640e-06 152.08 + 191.14 839.14 2.217465e-01 7.646563e-06 151.71 + 191.14 839.14 2.149170e-01 7.420948e-06 151.48 + 191.14 839.14 2.082628e-01 7.198933e-06 151.27 + 191.14 839.14 2.018148e-01 6.983797e-06 151.07 + 191.14 839.14 1.955637e-01 6.775233e-06 150.87 + 191.14 839.14 1.895088e-01 6.573209e-06 150.66 + 191.14 839.14 1.836413e-01 6.377441e-06 150.46 + 191.14 839.14 1.779556e-01 6.187739e-06 150.25 + 191.14 839.14 1.724435e-01 6.003832e-06 150.05 + 191.14 839.14 1.671044e-01 5.825693e-06 149.84 + 191.14 839.14 1.619305e-01 5.653068e-06 149.64 + 191.14 839.14 1.569171e-01 5.485795e-06 149.43 + 191.14 839.14 1.520567e-01 5.323629e-06 149.23 + 191.14 839.14 1.473487e-01 5.166550e-06 149.02 + 191.14 839.14 1.427866e-01 5.014334e-06 148.82 + 191.14 839.14 1.383658e-01 4.866836e-06 148.62 + 191.14 839.14 1.340800e-01 4.723843e-06 148.41 + 191.14 839.14 1.299287e-01 4.585333e-06 148.21 + 191.14 839.14 1.259059e-01 4.451113e-06 148.00 + 191.14 839.14 1.220077e-01 4.321053e-06 147.80 + 191.14 839.14 1.182286e-01 4.194964e-06 147.59 + 191.14 839.14 1.145681e-01 4.072831e-06 147.39 + 191.14 839.14 1.110208e-01 3.954478e-06 147.18 + 191.14 839.14 1.075253e-01 3.835007e-06 147.01 + 191.14 839.14 1.041030e-01 3.716316e-06 146.85 + 191.14 839.14 1.007909e-01 3.601448e-06 146.70 + 191.14 839.14 9.758420e-02 3.490235e-06 146.55 + 191.14 839.14 9.447963e-02 3.382564e-06 146.39 + 191.14 839.14 9.147253e-02 3.278274e-06 146.24 + 191.14 839.14 8.856231e-02 3.177342e-06 146.08 + 191.14 839.14 8.574467e-02 3.079623e-06 145.93 + 191.14 839.14 8.301676e-02 2.985015e-06 145.77 + 191.14 839.14 8.037451e-02 2.893377e-06 145.62 + 191.15 839.14 7.781737e-02 2.804692e-06 145.46 + 191.15 839.14 7.534159e-02 2.718828e-06 145.31 + 191.15 839.14 7.294464e-02 2.635698e-06 145.16 + 191.15 839.14 7.062296e-02 2.555179e-06 145.00 + 191.15 839.14 6.837606e-02 2.477253e-06 144.85 + 191.15 839.14 6.620066e-02 2.401807e-06 144.69 + 191.15 839.14 6.409454e-02 2.328763e-06 144.54 + 191.15 839.14 6.205453e-02 2.258013e-06 144.38 + 191.15 839.14 6.008025e-02 2.189542e-06 144.23 + 191.15 839.14 5.816878e-02 2.123249e-06 144.08 + 191.15 839.14 5.631818e-02 2.059067e-06 143.92 + 191.15 839.14 5.452568e-02 1.996901e-06 143.77 + 191.15 839.14 5.274843e-02 1.934156e-06 143.53 + 191.15 839.14 5.102769e-02 1.873369e-06 143.28 + 191.15 839.14 4.936314e-02 1.814568e-06 143.04 + 191.15 839.14 4.775220e-02 1.757660e-06 142.79 + 191.15 839.14 4.619445e-02 1.702632e-06 142.55 + 191.15 839.14 4.468752e-02 1.649398e-06 142.31 + 191.15 839.14 4.322979e-02 1.597903e-06 142.06 + 191.15 839.14 4.181901e-02 1.548066e-06 141.82 + 191.15 839.14 4.045480e-02 1.499874e-06 141.57 + 191.15 839.14 3.913511e-02 1.453255e-06 141.33 + 191.15 839.14 3.785850e-02 1.408158e-06 141.09 + 191.15 839.14 3.662301e-02 1.364513e-06 140.84 + 191.15 839.14 3.542831e-02 1.322310e-06 140.60 + 191.15 839.14 3.427259e-02 1.281483e-06 140.35 + 191.15 839.14 3.315415e-02 1.241973e-06 140.11 + 191.15 839.14 3.207261e-02 1.203767e-06 139.87 + 191.15 839.14 3.102635e-02 1.166808e-06 139.62 + 191.15 839.14 3.001423e-02 1.131053e-06 139.38 + 191.15 839.14 2.903476e-02 1.096453e-06 139.13 + 191.15 839.14 2.808760e-02 1.062994e-06 138.89 + 191.15 839.14 2.717134e-02 1.030626e-06 138.65 + 191.15 839.14 2.626059e-02 9.974335e-07 138.41 + 191.15 839.14 2.537629e-02 9.650491e-07 138.17 + 191.15 839.14 2.452209e-02 9.337668e-07 137.93 + 191.15 839.14 2.369664e-02 9.035374e-07 137.69 + 191.15 839.15 2.289897e-02 8.743256e-07 137.45 + 191.15 839.15 2.212787e-02 8.460868e-07 137.21 + 191.15 839.15 2.138301e-02 8.188088e-07 136.97 + 191.15 839.15 2.066323e-02 7.924491e-07 136.73 + 191.15 839.15 1.996767e-02 7.669767e-07 136.49 + 191.15 839.15 1.929528e-02 7.423527e-07 136.25 + 191.15 839.15 1.864577e-02 7.185666e-07 136.01 + 191.15 839.15 1.801812e-02 6.955812e-07 135.77 + 191.15 839.15 1.741160e-02 6.733695e-07 135.53 + 191.15 839.15 1.682529e-02 6.518976e-07 135.29 + 191.16 839.15 1.625892e-02 6.311565e-07 135.05 + 191.16 839.15 1.571162e-02 6.111134e-07 134.82 + 191.16 839.15 1.518274e-02 5.917451e-07 134.58 + 191.16 839.15 1.467148e-02 5.730217e-07 134.34 + 191.16 839.16 1.417761e-02 5.549356e-07 134.10 + 191.16 839.16 1.370037e-02 5.374583e-07 133.86 + 191.16 839.16 1.323558e-02 5.200246e-07 133.66 + 191.16 839.16 1.278160e-02 5.024430e-07 133.52 + 191.16 839.16 1.234335e-02 4.854705e-07 133.38 + 191.16 839.16 1.192013e-02 4.690801e-07 133.24 + 191.16 839.17 1.151141e-02 4.532517e-07 133.11 + 191.17 839.17 1.111657e-02 4.379603e-07 132.97 + 191.17 839.17 1.073541e-02 4.231989e-07 132.83 + 191.17 839.17 1.036732e-02 4.089436e-07 132.69 + 191.17 839.17 1.001185e-02 3.951771e-07 132.55 + 191.17 839.18 9.668438e-03 3.818777e-07 132.41 + 191.17 839.18 9.336930e-03 3.690392e-07 132.27 + 191.18 839.18 9.016789e-03 3.566410e-07 132.13 + 191.18 839.18 8.707624e-03 3.446678e-07 132.00 + 191.18 839.19 8.408951e-03 3.331009e-07 131.86 + 191.18 839.19 8.120627e-03 3.219348e-07 131.72 + 191.18 839.19 7.842191e-03 3.111517e-07 131.58 + 191.18 839.19 7.573301e-03 3.007382e-07 131.44 + 191.18 839.19 7.313535e-03 2.906781e-07 131.30 + 191.19 839.20 7.062771e-03 2.809666e-07 131.16 + 191.19 839.20 6.820606e-03 2.715881e-07 131.02 + 191.21 839.21 6.585661e-03 2.624410e-07 130.90 + 191.27 839.25 6.356860e-03 2.534489e-07 130.79 + 191.32 839.28 6.136089e-03 2.447724e-07 130.68 + 191.38 839.31 5.922986e-03 2.363972e-07 130.57 + 191.44 839.34 5.717284e-03 2.283129e-07 130.46 + 191.49 839.38 5.518652e-03 2.205065e-07 130.35 + 191.55 839.41 5.326992e-03 2.129741e-07 130.24 + 191.61 839.44 5.141987e-03 2.057032e-07 130.13 + 191.67 839.48 4.963409e-03 1.986849e-07 130.02 + 191.72 839.51 4.790969e-03 1.919078e-07 129.91 + 191.78 839.54 4.624581e-03 1.853686e-07 129.81 + 191.84 839.58 4.463972e-03 1.790565e-07 129.70 + 191.89 839.61 4.308939e-03 1.729636e-07 129.59 + 191.95 839.64 4.159237e-03 1.670801e-07 129.48 + 192.01 839.68 4.014789e-03 1.614031e-07 129.37 + 192.07 839.71 3.875357e-03 1.559233e-07 129.26 + 192.13 839.74 3.740768e-03 1.506338e-07 129.15 + 192.18 839.78 3.610805e-03 1.455262e-07 129.04 + 192.24 839.81 3.485403e-03 1.405977e-07 128.93 + 192.30 839.84 3.364357e-03 1.358405e-07 128.82 + 192.36 839.89 3.247567e-03 1.311597e-07 128.79 + 192.43 839.94 3.134828e-03 1.265704e-07 128.81 + 192.50 839.99 3.026044e-03 1.221421e-07 128.83 + 192.57 840.04 2.921034e-03 1.178674e-07 128.85 + 192.64 840.09 2.819669e-03 1.137411e-07 128.87 + 192.71 840.14 2.721785e-03 1.097565e-07 128.89 + 192.78 840.19 2.627333e-03 1.059117e-07 128.91 + 192.85 840.24 2.536160e-03 1.022002e-07 128.93 + 192.92 840.29 2.448150e-03 9.861758e-08 128.95 + 192.99 840.34 2.363164e-03 9.515800e-08 128.98 + 193.06 840.39 2.281157e-03 9.181973e-08 129.00 Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/levels =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/levels (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/levels (revision 1291) @@ -0,0 +1,1 @@ +link vert_level/levels Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/namelist.input =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/namelist.input (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/namelist.input (revision 1291) @@ -0,0 +1,65 @@ + &time_control + run_seconds = 14800 !! Duration of simulation [LES: !=0] + history_interval_s = 600 !! Frequency of outputs in seconds [LES: !=0] + frames_per_outfile = 6 !! Size of time dimension in files + start_year = 9999 !! Idealized simulation [LES: 9999] + end_year = 9999 !! Idealized simulation [LES: 9999] + restart = .false. !! (*) Output restart files ? + 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 + input_from_file = F, !! (n)(p2) Initialize a given domain with an input file + / + + &domains + ztop = 25000. !! Altitude above ground of model top [LES: !=0] + time_step = 10 !! Dynamical timestep + time_step_fract_num = 0 !! Additional fraction to time_step: numerator + time_step_fract_den = 1 !! Additional fraction to time_step: denominator + dx = 10000 !! (p2) Horizontal resolution + dy = 10000 !! (p2) Horizontal resolution (should be equal to dx) + e_we = 61 !! (r)(p2) Number of longitude grid points + e_sn = 61 !! (r)(p2) Number of latitude grid points + e_vert = 61 !! (r)(p2) Number of vertical levels + / + + &physics + init_TI = 250. !! Define constant thermal inertia value [LES: !=0] + init_AL = 0.25 !! Define constant albedo value [LES: !=0] + !! + !! OPTIONAL + !! + radt = 1 !! 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.001 !! (p3) Multiply ini & bdy zonal wind by init_U + init_MV = 0.001 !! (p3) Multiply ini & bdy meridional wind by init_V + init_U = 15. + init_V = 0. + / + + &dynamics + pert_coriolis = T !! Coriolis only acts on wind perturbation (idealized) + / + + &bdy_control + specified = F, !! (n)(p3) Boundary conditions specified by GCM + open_xs = T, !! (p3) Open boundary conditions @ western boundary + open_xe = T, !! (p3) Open boundary conditions @ eastern boundary + open_ys = T, !! (p3) Open boundary conditions @ southern boundary + open_ye = T, !! (p3) Open boundary conditions @ northern boundary + / + + &grib2 + / + + &fdda + / + + &namelist_quilt !! (*) + nio_tasks_per_group = 0, !! (*) + nio_groups = 1, !! (*) + / !! (*) Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/vert_level/levels =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/vert_level/levels (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/vert_level/levels (revision 1291) @@ -0,0 +1,21 @@ + 1.00000 + 0.995184 + 0.980846 + 0.957303 + 0.925073 + 0.884856 + 0.837508 + 0.784009 + 0.725431 + 0.662902 + 0.597572 + 0.530576 + 0.463003 + 0.395868 + 0.330092 + 0.266478 + 0.205706 + 0.148322 + 0.0947408 + 0.0452458 + 0.00000 Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/vert_level/levspe.pro =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/vert_level/levspe.pro (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/vert_level/levspe.pro (revision 1291) @@ -0,0 +1,95 @@ +pro levspe + + +; +; TWEAK PARAM +; +@param.idl + + +;; +;; HARD PARAM +;; +;!p.multi=[0,2,2] +set_plot, 'ps' +device, file='plot.ps';, /landscape +tinv=5 ;; intervalle plot +hache=10. ;; scale height +psurf=610. ;; pression surface -- ne change pas les eta levels + +; +; PRELIM +; +param = 'e_vert='+string(nlev,'(I0)')+' ztop='+string(altmax,'(I0)')+' e='+string(epsilon,'(I0)')+' c='+string(elong_cos,'(F4.2)') +; +epsilon = epsilon / (nlev-1) / 100. +; +exposant = !pi/nlev/elong_cos +; +; alpha (plus grand ecart en km) est determine pour que max(altitudes)=altmax +; +alpha = altmax / ( (nlev-1)/2. - sin(2.*exposant*(nlev-1))/4./exposant + epsilon*(nlev-1)^2/2. ) +print, 'alpha (km)', alpha +; +x=[findgen(nlev)] + + +; +; CALC +; +ecart = alpha - alpha*cos(exposant*x)^2 + epsilon*alpha*x +altitudes = alpha*x/2. - alpha*sin(2.*exposant*x)/4./exposant + epsilon*alpha*x^2/2. +logpressions = alog(psurf) - altitudes/hache +pressions=exp(logpressions) +ptop=psurf*exp(-altmax/hache) & print, ptop +etas=(pressions-ptop)/(psurf-ptop) +etas(nlev-1)=0. + + +; +; PRINT +; +print, 'eta levels' +print, etas +; +print, 'pressure' +pi=etas*(psurf-ptop)+ptop +print, pi +print, 'ptop', ptop +; +print, 'pseudo-altitude (meters)' +pseudo=10.*alog(610./pi) +print, 1000.*10.*alog(610./pi) + + +; +; PLOT +; +plot, x, etas, title='ETA '+param, xtickinterval=tinv +!psym=7 & oplot, x, etas +!psym=0 +plot, x, pseudo, title='ALTITUDE (km)', xtickinterval=tinv +!psym=7 & oplot, x, pseudo +!psym=0 +plot, x, pi, title='PRESSURE (Pa)', ylog=1, xtickinterval=tinv +!psym=7 & oplot, x, pi +!psym=0 +; +diff=pseudo - shift(pseudo,-1) & diff=-diff(0:nlev-2) +; +plot, x, ecart, title='ECART (km, L=theor x=calc)', xtickinterval=tinv +!psym=7 & oplot, x, diff +!psym=0 + + +; +; PRINT +; +;openw, 1, 'levels' +for i=0, nlev-1 do begin +; printf, 1, etas(i);,',' +print, etas(i);,',' +endfor +close, 1 + +end Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/vert_level/param.idl =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/vert_level/param.idl (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/vert_level/param.idl (revision 1291) @@ -0,0 +1,99 @@ +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +nlev=61 +altmax=120. ;; altitude maximale en km +;; +epsilon=30. ;; augmentation totale en pourcentage de l'ecart max (modifie inflexion) + ;; -- epsilon = 0 : point d'inflexion parfaitement plat + ;; -- epsilon = 100 : pas de point d'inflexion +elong_cos=2./3. ;; plus petit rapproche l'inflexion du sol ;; 2/3 parfait pour un nlev divisible par 3 +elong_cos=0.64 +;; plus petit ajoute des inflexions... +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;nlev=61 +;altmax=150. +;epsilon=30. +;elong_cos=2./3. +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +nlev=61 +altmax=150. +epsilon=100. +elong_cos=2. +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +nlev=201 +altmax=150. +epsilon=0. +elong_cos=2. +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +nlev=101 +altmax=80. ;;100. +epsilon=100. +elong_cos=60. +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +nlev=201 +altmax=150. +epsilon=0. +elong_cos=2. ;; pas inflexion +elong_cos=2.5 +elong_cos=3. +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;nlev=201 +;altmax=220. +;epsilon=0. +;elong_cos=2. ;; pas inflexion +;elong_cos=2.5 +;elong_cos=3. +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +nlev=61 +altmax=60. +epsilon=0. +elong_cos=2. ;; pas inflexion +elong_cos=2.5 +elong_cos=3. +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; + +epsilon=100. +elong_cos=60. + +;altmax=100. +altmax=130. +altmax=150. + +;;; ici OK avec simu + + +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +nlev=201 +altmax=150. +epsilon=0. +elong_cos=3. +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +nlev=151 +altmax=150. +epsilon=100. +elong_cos=60. +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;nlev=61 +nlev = 181 +altmax = 180. + + + +nlev=61 +altmax=25. +epsilon=100 +elong_cos=100. + +nlev=21 +altmax=15. +epsilon=100 +elong_cos=100. + +;nlev=21 +;altmax=40. +;epsilon=100 +;elong_cos=100. + Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/vert_level/simple.pro =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/vert_level/simple.pro (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/vert_level/simple.pro (revision 1291) @@ -0,0 +1,94 @@ +pro simple, z_scale + + +nlev=101 +; +; z-scale doit etre suffisamment bas pour que l avant-dernier niveau soit assez haut +; mais pas trop sinon instable (il faut avoir des niveaux espaces en haut semble-t-il) +; +z_scale = .25 +;z_scale = .15 marche pas (niveaux pas assez espaces en haut ?) +ptop = 0.2 +npbl = 5 + ;z_scale = .15 + ;ptop = 0.05 ;marche pas + +z_scale = .20 & ptop = 0.1 & npbl = 6 + + +psurf=610. + +pbl = [1.0000,0.9995,0.9980,0.9950,0.9850,0.9700,0.9400,0.9000] +x=[findgen(nlev)] +nlev = nlev - npbl + 1 + +etas=[findgen(nlev)] +FOR k=0, nlev-1 DO BEGIN + etas(k) = (exp(-float(k)/float(nlev-1)/z_scale) - exp(-1./z_scale)) / (1.-exp(-1./z_scale)) +ENDFOR +etas(nlev-1)=0. + +etas = [pbl(0:npbl-1),etas(1:nlev-1)] +nlev = nlev + npbl - 1 +print, nlev + +;pressions = ptop + etas * (psurf-ptop) +;zzz = 10. * alog ( psurf / pressions ) +;print, zzz + +; +; HARD PARAM +; +!p.multi=[0,2,2] +;set_plot, 'ps' +;device, file='plot.ps', /landscape + + + +; +; PRINT +; +print, 'eta levels' +print, etas +; +print, 'pressure' +pi=etas*(psurf-ptop)+ptop +print, pi +print, 'ptop', ptop +; +print, 'pseudo-altitude' +pseudo=10.*alog(610./pi) +print, 1000.*10.*alog(610./pi) + + +; +; PLOT +; +plot, x, etas, title='ETA', xtickinterval=tinv +!psym=7 & oplot, x, etas +!psym=0 +plot, x, pseudo, title='ALTITUDE (km)', xtickinterval=tinv +!psym=7 & oplot, x, pseudo +!psym=0 +plot, x, pi, title='PRESSURE (Pa)', ylog=1, xtickinterval=tinv +!psym=7 & oplot, x, pi +!psym=0 +; +diff=pseudo - shift(pseudo,-1) & diff=-diff(0:nlev-2) +; +plot, x, diff, title='ECART (km, L=theor x=calc)', xtickinterval=tinv +!psym=7 & oplot, x, diff +!psym=0 + +;stop + +; +; PRINT +; +openw, 1, 'levels' +for i=0, nlev-1 do begin + printf, 1, etas(i);,',' +endfor +close, 1 + +end Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/ze_hill =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/ze_hill (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_idealcrater/ze_hill (revision 1291) @@ -0,0 +1,47 @@ +0.0 +10.0 + +-4000. +2.0 + +0.0 +10.0 + +-4000. +2.0 + +1000. +4.0 + +4000. +2.0 + +1000. +4.0 + +6000. +4.0 + +4000. +8.0 + +50. +8.0 + +4000. +8.0 + +500. +5.0 + +4000. +8.0 + +100. +4.0 + +1000. +2.0 + +2000. +6.0 Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/callphys.def.MCD5 =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/callphys.def.MCD5 (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/callphys.def.MCD5 (revision 1291) @@ -0,0 +1,1 @@ +link ../../../../../LMDZ.MARS/deftank/callphys.def.MCD5 Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/callphys.def.mars0 =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/callphys.def.mars0 (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/callphys.def.mars0 (revision 1291) @@ -0,0 +1,191 @@ +##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 = /home/spiga/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.1 +# 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. +# Thermals +calltherm = .true. +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 = 6 +# 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/tests_newphys_THARSISREF/callphys.def.mars11_radact =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/callphys.def.mars11_radact (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/callphys.def.mars11_radact (revision 1291) @@ -0,0 +1,185 @@ +##General options +##~~~~~~~~~~~~~~~ +#Run with or without tracer transport ? +tracer=.true. + +#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 = /home/spiga/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.1 +# 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. +# Thermals +calltherm = .true. +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 = 6 +# 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 +# DUST: Radiatively active dust ? (matters if dustbin>0) +active = .true. +# DUST: use mass and number mixing ratios to predict dust size ? +doubleq = .true. +# 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. +# WATERICE: Radiatively active transported atmospheric water ice ? +activice = .true. +#activice = .false. +# WATER: Compute water cycle +water = .true. +# 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. +# 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 = 2 + + +## 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/tests_newphys_THARSISREF/callphys.def.mars12_radact_scav =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/callphys.def.mars12_radact_scav (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/callphys.def.mars12_radact_scav (revision 1291) @@ -0,0 +1,183 @@ +##General options +##~~~~~~~~~~~~~~~ +#Run with or without tracer transport ? +tracer=.true. + +#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 = /home/spiga/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.1 +# 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. +# Thermals +calltherm = .true. +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 = 6 +# 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 +# DUST: Radiatively active dust ? (matters if dustbin>0) +active = .true. +# DUST: use mass and number mixing ratios to predict dust size ? +doubleq = .true. +# 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. +# DUST/WATERICE: Gravitationnal sedimentation ? +sedimentation = .true. +# WATERICE: Radiatively active transported atmospheric water ice ? +activice = .true. +#activice = .false. +# WATER: Compute water cycle +water = .true. +# WATER: Microphysical scheme for water-ice clouds? +microphys = .true. +# 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. +# 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 = 2 + + +## 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/tests_newphys_THARSISREF/callphys.def.mars1_radact =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/callphys.def.mars1_radact (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/callphys.def.mars1_radact (revision 1291) @@ -0,0 +1,188 @@ +##General options +##~~~~~~~~~~~~~~~ +#Run with or without tracer transport ? +tracer=.true. + +#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 = /home/spiga/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.1 +# 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. +# Thermals +calltherm = .true. +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 = 6 +# 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. +# WATERICE: Radiatively active transported atmospheric water ice ? +activice = .true. +#activice = .false. +# WATER: Compute water cycle +water = .true. +# 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. +# 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 = 2 + + +## 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/tests_newphys_THARSISREF/callphys.def.mars2 =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/callphys.def.mars2 (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/callphys.def.mars2 (revision 1291) @@ -0,0 +1,190 @@ +##General options +##~~~~~~~~~~~~~~~ +#Run with or without tracer transport ? +tracer=.true. + +#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 = /home/spiga/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.1 +# 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. +# Thermals +calltherm = .true. +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 = 6 +# 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 = 1 +# 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. +# 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/tests_newphys_THARSISREF/callphys.def.mars3 =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/callphys.def.mars3 (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/callphys.def.mars3 (revision 1291) @@ -0,0 +1,188 @@ +freedust = .true. + +##General options +##~~~~~~~~~~~~~~~ +#Run with or without tracer transport ? +tracer=.true. + +#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 = /home/spiga/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.1 +# 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. +# Thermals +calltherm = .true. +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 = 6 +# 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 +# DUST: Radiatively active dust ? (matters if dustbin>0) +active = .true. +# DUST: use mass and number mixing ratios to predict dust size ? +doubleq = .true. +# DUST: lifted by GCM surface winds ? +lifting = .true. +# DUST: lifted by dust devils ? +callddevil = .false. +# DUST: Scavenging by H2O/CO2 snowfall ? +scavenging = .true. +scavenging = .false. +# DUST/WATERICE: Gravitationnal sedimentation ? +sedimentation = .true. +# 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/tests_newphys_THARSISREF/launch.ciclad =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/launch.ciclad (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/launch.ciclad (revision 1291) @@ -0,0 +1,14 @@ +#PBS -S /bin/bash +#PBS -j oe +#PBS -q infini +#PBS -N MARS_ATTACKS +#PBS -l "nodes=1:ppn=24" + +####### prep run +cd $PBS_O_WORKDIR +\rm *rsl.* 2> /dev/null +\rm wrfo* 2> /dev/null + +####### run +ulimit -s unlimited +mpirun wrf.exe Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/launch.gnome =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/launch.gnome (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/launch.gnome (revision 1291) @@ -0,0 +1,63 @@ +####################################################################### +####################################################################### nom +# @ job_name = THARSIS_REF.LAST.nodiff6 +####################################################################### MPICH / serial +# @ job_type = MPICH +####################################################################### 8 / 16 / 20 / 32 / 64 +# @ total_tasks = 12 +####################################################################### A / B / C / AP / BP / CP +# @ class = AP +####################################################################### 24gb (cf. llps pendant fonctionnement) +# @ resources = ConsumableCpus(1)ConsumableMemory(1gb) +####################################################################### real ou ideal ??? +# @ environment = real=1 ; tmpdir=$(job_name).$(jobid) ; para=$(job_type) +####################################################################### + +####################################################################### +####################################################################### +# @ output = $(job_name).$(jobid)/LOG +# @ error = $(output) +# @ blocking = unlimited +# @ notification = never +# @ queue +grep "^# @" $0 +set -vx +####################################################################### +mkdir $LOADL_STEP_INITDIR/${tmpdir} +cd $LOADL_STEP_INITDIR +\rm zeLAST +ln -sf ${tmpdir} zeLAST +cd zeLAST +####################################################################### +cp $LOADL_STEP_INITDIR/launch.gnome . +cp $LOADL_STEP_INITDIR/*.def . +cp $LOADL_STEP_INITDIR/*.exe . +cp $LOADL_STEP_INITDIR/namelist* . +#cp $LOADL_STEP_INITDIR/wrfi* . ### executer sur la frontale +# ### au prealable !!! +ln -sf $LOADL_STEP_INITDIR/wrfi* . +####################################################################### +if [[ "${real}" == "0" ]] +then + cp $LOADL_STEP_INITDIR/input_* . + cp $LOADL_STEP_INITDIR/ze_hill . ## topo + cp $LOADL_STEP_INITDIR/levels . ## custom vert levels + ##mpirun ideal.exe ##marche pas + ##mkdir ideal + ##mv rsl.* ideal/ +else +# cp $LOADL_STEP_INITDIR/wrfb* . +ln -sf $LOADL_STEP_INITDIR/wrfb* . +fi +####################################################################### +#### avant wrf.exe, necessaire ??? +#declare -x MP_STACK_SIZE=640000000 +#declare -x MP_STACK_SIZE=64000000 +####################################################################### +if [[ "${para}" == "MPICH" ]] +then + mpirun wrf.exe +else + ./wrf.exe +fi +####################################################################### Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/makemeso_input =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/makemeso_input (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/makemeso_input (revision 1291) @@ -0,0 +1,8 @@ +5 +12 +121 +121 +61 +1 +4 +2 Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/namelist.input =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/namelist.input (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/namelist.input (revision 1291) @@ -0,0 +1,63 @@ + &time_control + start_year = 2024, !! (p1) Start Martian Year (20XX for MY XX) + start_month = 05, !! (p1) Start Martian Month + start_day = 01, !! (p1) Start Martian Day + start_hour = 06, !! (p1) Start Martian Hour (at longitude 0) + end_year = 2024, !! (p1) End Martian Year (20XX for MY XX) + end_month = 05, !! (p1) End Martian Month + end_day = 05, !! (p1) End Martian Day + end_hour = 06, !! (p1) End Martian Hour (at longitude 0) + history_interval = 37, !! Frequency of outputs (37 --> 3700s = 1 Martian hour) + frames_per_outfile = 24, !! Size of time dimension in files + restart = .false. !! (*) Output restart files ? + 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 + time_step = 20 !! Dynamical timestep + dx = 30000, !! (p2) Horizontal resolution + dy = 30000, !! (p2) Horizontal resolution (should be equal to dx) + e_we = 121, !! (r)(p2) Number of longitude grid points + e_sn = 121, !! (r)(p2) Number of latitude grid points + e_vert = 61, !! (r)(p2) Number of vertical levels + p_top_requested = 5 !! (p3) Chosen value of pressure at the top of the model + !! + !! OPTIONAL + !! + num_metgrid_levels = 30 !! (p1) number of vertical levels in GCM inputs (+1) + / + + &physics + mars = 12 + num_soil_layers = 18 !!! NOUVELLE PHYSIQUE + / + + &dynamics + !! + !! OPTIONAL + !! + !tiso = 100. !! -- temperature for isotherm level (mesosphere) in init + !diff_6th_opt = 0 !! -- swith off correction filter for WRF integrations. normally works when switched on (default). + !khdif = 10., !! ** a more straightforward diffusion for tests. horizontal diffusion constant (m^2/s) + !kvdif = 10., !! ** a more straightforward diffusion for tests. vertical diffusion constant (m^2/s) + !km_opt = 1, !! ** a more straightforward diffusion for tests. eddy coefficient option + / + + &bdy_control + / + + &grib2 + / + + &fdda + / + + &namelist_quilt !! (*) + nio_tasks_per_group = 0, !! (*) + nio_groups = 1, !! (*) + / !! (*) Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/namelist.wps =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/namelist.wps (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/namelist.wps (revision 1291) @@ -0,0 +1,37 @@ +&share + wrf_core = 'ARW', !! [do not modify: choice of dynamical core] + max_dom = 1 !! (r)(n)(p2) Total number of domains + start_date = '2024-05-01_06:00:00' !! YYYY-MM-DD_HH:mm:ss start date + end_date = '2024-05-05_06:00:00' !! YYYY-MM-DD_HH:mm:ss end date + interval_seconds = 3700 !! (p2) Frequency of large-scale fields update (s) + io_form_geogrid = 2, !! [do not modify: choice of NETCDF outputs] + debug_level = 0, !! verbose level of the programs + opt_output_from_geogrid_path='./' !! location of the geogrid outputs +/ + +&geogrid + parent_id = 0, !! (n)(p2) Associated parent domain + parent_grid_ratio = 1, !! (n)(p2) Ratio of horizontal resolution parent/nest + i_parent_start = 0, !! (n)(p2) x-position of the bottom-left nest corner + j_parent_start = 0, !! (n)(p2) y-position of the bottom-left nest corner + e_we = 121, !! (r)(p2) Number of longitude grid points + e_sn = 121, !! (r)(p2) Number of latitude grid points + dx = 30000 + dy = 30000 + geog_data_res = '64ppd_noHRti', !! choice of static data sources + !! NB: possible: '64ppd', '32ppd', ... + !! NB: please glance at geogrid/GEOGRID.TBL + map_proj = 'mercator', !! map projection: 'mercator', 'lambert' or 'polar' + ref_lat = 10., !! north latitude of the center of the domain + ref_lon = -125., !! east longitude of the center of the domain + truelat1 = 10., !! (lambert or polar) lat position of projection cone + truelat2 = 0.0, !! [do not modify] + stand_lon = -125., !! (lambert or polar) lon position of projection cone + geog_data_path = './WPS_GEOG', !! [do not modify: symbolic link in the WPS folder] +/ + +&metgrid + fg_name = './WPSFEED/LMD' !! [do not modify: symbolic link in the WPS folder] + io_form_metgrid = 2, !! [do not modify: choice of NETCDF outputs] + opt_output_from_metgrid_path='./WRFFEED/current' !! [do not modify: symbolic link] +/ Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/run.def =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/run.def (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphys_THARSISREF/run.def (revision 1291) @@ -0,0 +1,2 @@ +# some definitions for the physics, in file 'callphys.def' +INCLUDEDEF=callphys.def Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphysfast/callphys.def =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphysfast/callphys.def (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphysfast/callphys.def (revision 1291) @@ -0,0 +1,118 @@ +##General options +##~~~~~~~~~~~~~~~ +#Directory where external input files are +datadir=/u/forget/WWW/datagcm/datafile/ + +#Run with or without tracer transport ? +tracer=.false. + +#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 =.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 (same as =24 for now) +# =24 Mars Year 24 from TES assimilation (ie: MCD reference case) +# =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm) +# =26 Mars Year 26 from TES assimilation +iaervar = 1 +# Dust opacity at 700 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 = .false. +# call CO2 NIR absorption ? matters only if callrad=T +callnirco2 = .true. +# call turbulent vertical diffusion ? +calldifv = .true. +# call convective adjustment ? +calladj = .true. +# call CO2 condensation ? +callcond =.true. +# call thermal conduction in the soil ? +callsoil = .true. +# call Lott's gravity wave/subgrid topography scheme ? +calllott = .false. + +## Radiative transfer options : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# the rad.transfer is computed every "iradia" physical timestep +iradia = 10 +# 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 = 0 +# DUST: Radiatively active dust ? (matters if dustbin>0) +active = .false. +# DUST: use mass and number mixing ratios to predict dust size ? +# (must also have dustbin=1) +doubleq = .false. +# DUST: lifted by GCM surface winds ? +lifting = .false. +# DUST: lifted by dust devils ? +callddevil = .false. +# DUST: Scavenging by CO2 snowfall ? +scavenging = .false. +# DUST/WATERICE: Gravitationnal sedimentation ? +sedimentation = .false. +# WATERICE: Radiatively active transported atmospheric water ice ? +activice = .false. +# WATER: Compute water cycle +water = .false. +# 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 = .false. +# PHOTOCHEMISTRY: include chemical species +photochem = .false. + +## Thermospheric options (relevant if tracer=T) : +##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call thermosphere ? +callthermos = .false. +# WATER: included without cycle (only if water=.false.) +thermoswater = .false. +# call thermal conduction ? (only if callthermos=.true.) +callconduct = .false. +# call EUV heating ? (only if callthermos=.true.) +calleuv=.false. +# call molecular viscosity ? (only if callthermos=.true.) +callmolvis = .false. +# call molecular diffusion ? (only if callthermos=.true.) +callmoldiff = .false. +# call thermospheric photochemistry ? (only if callthermos=.true.) +thermochem = .false. +# date for solar flux calculation: (1985 < date < 2002) +## (Solar min=1996.4 ave=1993.4 max=1990.6) +solarcondate = 1993.4 + + + Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphysfast/namelist.input =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphysfast/namelist.input (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphysfast/namelist.input (revision 1291) @@ -0,0 +1,71 @@ + &time_control + start_year = 2024, !! (p1) Start Martian Year (20XX for MY XX) + start_month = 01, !! (p1) Start Martian Month + start_day = 17, !! (p1) Start Martian Day + start_hour = 02, !! (p1) Start Martian Hour (at longitude 0) + end_year = 2024, !! (p1) End Martian Year (20XX for MY XX) + end_month = 01, !! (p1) End Martian Month + end_day = 17, !! (p1) End Martian Day + end_hour = 12, !! (p1) End Martian Hour (at longitude 0) + history_interval = 37, !! Frequency of outputs (37 --> 3700s = 1 Martian hour) + frames_per_outfile = 24, !! Size of time dimension in files + restart = .false. !! (*) Output restart files ? + 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 + time_step = 50 !! Dynamical timestep + dx = 20000, !! (p2) Horizontal resolution + dy = 20000, !! (p2) Horizontal resolution (should be equal to dx) + e_we = 61, !! (r)(p2) Number of longitude grid points + e_sn = 61, !! (r)(p2) Number of latitude grid points + e_vert = 61, !! (r)(p2) Number of vertical levels + p_top_requested = 5 !! (p3) Chosen value of pressure at the top of the model + !! + !! OPTIONAL + !! + num_metgrid_levels = 33 !! (p1) number of vertical levels in GCM inputs (+1) + / + + &physics + !! + !! OPTIONAL + !! + num_soil_layers = 18 !!! NOUVELLE PHYSIQUE + init_TI = 100., !! Define constant thermal inertia value + init_AL = 0.27, !! Define constant albedo value + / + + &dynamics + !! + !! OPTIONAL + !! + non_hydrostatic = F, !! Integrate in non-hydrostatic/hydrostatic mode + tiso = 100. + / + + &bdy_control + !! + !! OPTIONAL + !! + spec_bdy_width = 3 !! (p3) Width of transition zone with specified=T + !! NB: spec_bdy_width must equal relax_zone+1 + relax_zone = 2 !! (p3) Width of relaxation zone with specified=T + / + + &grib2 + / + + &fdda + / + + &namelist_quilt !! (*) + nio_tasks_per_group = 0, !! (*) + nio_groups = 1, !! (*) + / !! (*) Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphysfast/run.def =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphysfast/run.def (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphysfast/run.def (revision 1291) @@ -0,0 +1,2 @@ +# some definitions for the physics, in file 'callphys.def' +INCLUDEDEF=callphys.def Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/callphys.def =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/callphys.def (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/callphys.def (revision 1291) @@ -0,0 +1,116 @@ +##General options +##~~~~~~~~~~~~~~~ +#Run with or without tracer transport ? +tracer=.false. + +#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" ? +stats =.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 (same as =24 for now) +# =24 Mars Year 24 from TES assimilation (ie: MCD reference case) +# =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm) +# =26 Mars Year 26 from TES assimilation +iaervar = 1 +# Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case) +tauvis=0.5 +# 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 = .false. +# call CO2 NIR absorption ? matters only if callrad=T +callnirco2 = .false. +# call turbulent vertical diffusion ? +calldifv = .true. +# call convective adjustment ? +calladj = .false. +# call CO2 condensation ? +callcond =.false. +# call thermal conduction in the soil ? +callsoil = .true. +# call Lott's gravity wave/subgrid topography scheme ? +calllott = .false. + +## 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 = 0 +# DUST: Radiatively active dust ? (matters if dustbin>0) +active = .false. +# DUST: use mass and number mixing ratios to predict dust size ? +# (must also have dustbin=1) +doubleq = .false. +# DUST: lifted by GCM surface winds ? +lifting = .false. +# DUST: lifted by dust devils ? +callddevil = .false. +# DUST: Scavenging by CO2 snowfall ? +scavenging = .false. +# DUST/WATERICE: Gravitationnal sedimentation ? +sedimentation = .false. +# WATERICE: Radiatively active transported atmospheric water ice ? +activice = .false. +# WATER: Compute water cycle +water = .false. +# 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 = .false. +#caps = .true. +# PHOTOCHEMISTRY: include chemical species +photochem = .false. + +## Thermospheric options (relevant if tracer=T) : +##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call thermosphere ? +callthermos = .false. +# WATER: included without cycle (only if water=.false.) +thermoswater = .false. +# call thermal conduction ? (only if callthermos=.true.) +callconduct = .false. +# call EUV heating ? (only if callthermos=.true.) +calleuv=.false. +# call molecular viscosity ? (only if callthermos=.true.) +callmolvis = .false. +# call molecular diffusion ? (only if callthermos=.true.) +callmoldiff = .false. +# call thermospheric photochemistry ? (only if callthermos=.true.) +thermochem = .false. +# date for solar flux calculation: (1985 < date < 2002) +## (Solar min=1996.4 ave=1993.4 max=1990.6) +solarcondate = 1993.4 + + + Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/dustopacity.def =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/dustopacity.def (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/dustopacity.def (revision 1291) @@ -0,0 +1,1 @@ +0.5 Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/input_coord =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/input_coord (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/input_coord (revision 1291) @@ -0,0 +1,7 @@ +353.87 lon +-1.88 lat +244. ls +07. lct +20. ztop +100. nz + Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/input_more =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/input_more (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/input_more (revision 1291) @@ -0,0 +1,1 @@ + -1473.35 211.31 Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/input_sounding =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/input_sounding (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/input_sounding (revision 1291) @@ -0,0 +1,205 @@ + 7.17 203.68 0.00 + -1468.35 198.70 0.00 -1.91 1.26 + -1463.35 198.29 0.00 -2.23 1.42 + -1453.35 199.28 0.00 -2.87 1.73 + -1423.35 200.52 0.00 -3.52 1.54 + -1373.35 202.08 0.00 -3.72 0.05 + -1273.35 205.08 0.00 -3.96 -2.09 + -1173.35 207.71 0.00 -4.44 -3.74 + -1073.35 209.33 0.00 -5.70 -4.86 + -973.35 210.96 0.00 -6.96 -5.99 + -873.35 212.30 0.00 -7.88 -6.84 + -773.35 213.23 0.00 -8.30 -7.28 + -673.35 214.16 0.00 -8.72 -7.73 + -573.35 215.09 0.00 -9.14 -8.17 + -473.35 216.03 0.00 -9.56 -8.62 + -373.35 216.97 0.00 -9.99 -9.06 + -273.35 217.53 0.00 -9.56 -8.89 + -173.35 218.09 0.00 -9.12 -8.72 + -73.35 218.65 0.00 -8.69 -8.54 + 26.65 219.21 0.00 -8.26 -8.37 + 126.65 219.77 0.00 -7.83 -8.20 + 226.65 220.33 0.00 -7.40 -8.03 + 326.65 220.90 0.00 -6.97 -7.86 + 426.65 221.46 0.00 -6.54 -7.68 + 526.65 222.03 0.00 -6.10 -7.51 + 626.65 222.58 0.00 -5.74 -7.33 + 726.65 223.00 0.00 -5.78 -7.12 + 826.65 223.43 0.00 -5.81 -6.92 + 926.65 223.86 0.00 -5.85 -6.71 + 1026.65 224.29 0.00 -5.88 -6.50 + 1126.65 224.72 0.00 -5.92 -6.29 + 1226.65 225.15 0.00 -5.95 -6.08 + 1326.65 225.58 0.00 -5.99 -5.87 + 1426.65 226.01 0.00 -6.02 -5.66 + 1526.65 226.44 0.00 -6.06 -5.46 + 1626.65 226.88 0.00 -6.09 -5.25 + 1726.65 227.31 0.00 -6.13 -5.04 + 1826.65 227.75 0.00 -6.16 -4.83 + 1926.65 228.18 0.00 -6.20 -4.62 + 2026.65 228.62 0.00 -6.23 -4.41 + 2126.65 229.06 0.00 -6.27 -4.20 + 2226.65 229.49 0.00 -6.30 -4.00 + 2326.65 229.91 0.00 -6.28 -3.64 + 2426.65 230.31 0.00 -6.21 -3.11 + 2526.65 230.71 0.00 -6.13 -2.59 + 2626.65 231.11 0.00 -6.06 -2.07 + 2726.65 231.50 0.00 -5.99 -1.54 + 2826.65 231.90 0.00 -5.91 -1.02 + 2926.65 232.30 0.00 -5.84 -0.49 + 3026.65 232.70 0.00 -5.76 0.03 + 3126.65 233.11 0.00 -5.69 0.55 + 3226.65 233.51 0.00 -5.62 1.08 + 3326.65 233.91 0.00 -5.54 1.60 + 3426.65 234.31 0.00 -5.47 2.12 + 3526.65 234.72 0.00 -5.39 2.65 + 3626.65 235.12 0.00 -5.32 3.17 + 3726.65 235.53 0.00 -5.25 3.69 + 3826.65 235.93 0.00 -5.17 4.22 + 3926.65 236.34 0.00 -5.10 4.74 + 4026.65 236.75 0.00 -5.02 5.26 + 4126.65 237.16 0.00 -4.95 5.79 + 4226.65 237.57 0.00 -4.88 6.31 + 4326.65 237.98 0.00 -4.80 6.83 + 4426.65 238.39 0.00 -4.73 7.36 + 4526.65 238.80 0.00 -4.65 7.88 + 4626.65 239.21 0.00 -4.58 8.41 + 4726.65 239.62 0.00 -4.51 8.93 + 4826.65 240.03 0.00 -4.43 9.45 + 4926.65 240.47 0.00 -4.45 9.60 + 5026.65 240.92 0.00 -4.55 9.45 + 5126.65 241.37 0.00 -4.64 9.31 + 5226.65 241.82 0.00 -4.73 9.16 + 5326.65 242.27 0.00 -4.82 9.02 + 5426.65 242.72 0.00 -4.92 8.87 + 5526.65 243.17 0.00 -5.01 8.73 + 5626.65 243.63 0.00 -5.10 8.58 + 5726.65 244.08 0.00 -5.20 8.44 + 5826.65 244.54 0.00 -5.29 8.29 + 5926.65 244.99 0.00 -5.38 8.15 + 6026.65 245.45 0.00 -5.48 8.00 + 6126.65 245.91 0.00 -5.57 7.85 + 6226.65 246.37 0.00 -5.66 7.71 + 6326.65 246.83 0.00 -5.75 7.56 + 6426.65 247.29 0.00 -5.85 7.42 + 6526.65 247.75 0.00 -5.94 7.27 + 6626.65 248.21 0.00 -6.03 7.13 + 6726.65 248.68 0.00 -6.13 6.98 + 6826.65 249.14 0.00 -6.22 6.84 + 6926.65 249.60 0.00 -6.31 6.69 + 7026.65 250.07 0.00 -6.41 6.55 + 7126.65 250.54 0.00 -6.50 6.40 + 7226.65 251.00 0.00 -6.59 6.26 + 7326.65 251.47 0.00 -6.69 6.11 + 7426.65 251.94 0.00 -6.78 5.97 + 7526.65 252.41 0.00 -6.87 5.82 + 7626.65 252.88 0.00 -6.96 5.68 + 7726.65 253.35 0.00 -7.06 5.53 + 7826.65 253.83 0.00 -7.15 5.39 + 7926.65 254.30 0.00 -7.24 5.24 + 8026.65 254.78 0.00 -7.34 5.10 + 8126.65 255.25 0.00 -7.43 4.95 + 8226.65 255.73 0.00 -7.52 4.80 + 8326.65 256.20 0.00 -7.62 4.66 + 8426.65 256.68 0.00 -7.71 4.51 + 8526.65 257.16 0.00 -7.80 4.37 + 8626.65 257.56 0.00 -7.53 4.24 + 8726.65 257.92 0.00 -7.10 4.12 + 8826.65 258.29 0.00 -6.66 4.00 + 8926.65 258.66 0.00 -6.23 3.87 + 9026.65 259.02 0.00 -5.80 3.75 + 9126.65 259.39 0.00 -5.37 3.63 + 9226.65 259.75 0.00 -4.94 3.51 + 9326.65 260.12 0.00 -4.51 3.39 + 9426.65 260.49 0.00 -4.07 3.27 + 9526.65 260.86 0.00 -3.64 3.14 + 9626.65 261.23 0.00 -3.21 3.02 + 9726.65 261.59 0.00 -2.78 2.90 + 9826.65 261.96 0.00 -2.35 2.78 + 9926.65 262.33 0.00 -1.92 2.66 + 10026.65 262.70 0.00 -1.49 2.53 + 10126.65 263.07 0.00 -1.05 2.41 + 10226.65 263.44 0.00 -0.62 2.29 + 10326.65 263.82 0.00 -0.19 2.17 + 10426.65 264.19 0.00 0.24 2.05 + 10526.65 264.56 0.00 0.67 1.92 + 10626.65 264.93 0.00 1.10 1.80 + 10726.65 265.31 0.00 1.54 1.68 + 10826.65 265.68 0.00 1.97 1.56 + 10926.65 266.05 0.00 2.40 1.44 + 11026.65 266.43 0.00 2.83 1.32 + 11126.65 266.80 0.00 3.26 1.19 + 11226.65 267.18 0.00 3.69 1.07 + 11326.65 267.55 0.00 4.13 0.95 + 11426.65 267.93 0.00 4.56 0.83 + 11526.65 268.30 0.00 4.99 0.71 + 11626.65 268.68 0.00 5.42 0.58 + 11726.65 269.06 0.00 5.85 0.46 + 11826.65 269.43 0.00 6.28 0.34 + 11926.65 269.81 0.00 6.72 0.22 + 12026.65 270.19 0.00 7.15 0.10 + 12126.65 270.57 0.00 7.58 -0.02 + 12226.65 270.95 0.00 8.01 -0.15 + 12326.65 271.33 0.00 8.44 -0.27 + 12426.65 271.71 0.00 8.87 -0.39 + 12526.65 272.09 0.00 9.31 -0.51 + 12626.65 272.47 0.00 9.74 -0.63 + 12726.65 272.85 0.00 10.17 -0.76 + 12826.65 273.23 0.00 10.60 -0.88 + 12926.65 273.61 0.00 11.03 -1.00 + 13026.65 274.00 0.00 11.46 -1.12 + 13126.65 274.38 0.00 11.90 -1.24 + 13226.65 274.76 0.00 12.33 -1.37 + 13326.65 275.15 0.00 12.73 -1.48 + 13426.65 275.50 0.00 12.58 -1.41 + 13526.65 275.85 0.00 12.42 -1.33 + 13626.65 276.21 0.00 12.27 -1.26 + 13726.65 276.56 0.00 12.12 -1.19 + 13826.65 276.92 0.00 11.96 -1.11 + 13926.65 277.27 0.00 11.81 -1.04 + 14026.65 277.63 0.00 11.65 -0.97 + 14126.65 277.99 0.00 11.50 -0.90 + 14226.65 278.34 0.00 11.34 -0.82 + 14326.65 278.70 0.00 11.19 -0.75 + 14426.65 279.06 0.00 11.04 -0.68 + 14526.65 279.41 0.00 10.88 -0.61 + 14626.65 279.77 0.00 10.73 -0.53 + 14726.65 280.13 0.00 10.57 -0.46 + 14826.65 280.48 0.00 10.42 -0.39 + 14926.65 280.84 0.00 10.26 -0.31 + 15026.65 281.20 0.00 10.11 -0.24 + 15126.65 281.56 0.00 9.96 -0.17 + 15226.65 281.92 0.00 9.80 -0.10 + 15326.65 282.28 0.00 9.65 -0.02 + 15426.65 282.64 0.00 9.49 0.05 + 15526.65 283.00 0.00 9.34 0.12 + 15626.65 283.36 0.00 9.19 0.19 + 15726.65 283.72 0.00 9.03 0.27 + 15826.65 284.08 0.00 8.88 0.34 + 15926.65 284.44 0.00 8.72 0.41 + 16026.65 284.80 0.00 8.57 0.48 + 16126.65 285.16 0.00 8.41 0.56 + 16226.65 285.52 0.00 8.26 0.63 + 16326.65 285.88 0.00 8.11 0.70 + 16426.65 286.24 0.00 7.95 0.78 + 16526.65 286.60 0.00 7.80 0.85 + 16626.65 286.97 0.00 7.64 0.92 + 16726.65 287.33 0.00 7.49 0.99 + 16826.65 287.69 0.00 7.33 1.07 + 16926.65 288.05 0.00 7.18 1.14 + 17026.65 288.42 0.00 7.03 1.21 + 17126.65 288.78 0.00 6.87 1.28 + 17226.65 289.14 0.00 6.72 1.36 + 17326.65 289.51 0.00 6.56 1.43 + 17426.65 289.87 0.00 6.41 1.50 + 17526.65 290.24 0.00 6.26 1.58 + 17626.65 290.60 0.00 6.10 1.65 + 17726.65 290.97 0.00 5.95 1.72 + 17826.65 291.33 0.00 5.79 1.79 + 17926.65 291.70 0.00 5.64 1.87 + 18026.65 292.06 0.00 5.48 1.94 + 18126.65 292.43 0.00 5.33 2.01 + 18226.65 292.79 0.00 5.18 2.08 + 18326.65 293.16 0.00 5.02 2.16 + 18426.65 293.53 0.00 4.87 2.23 + 18526.65 293.89 0.00 4.71 2.30 Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/input_therm =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/input_therm (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/input_therm (revision 1291) @@ -0,0 +1,204 @@ + 191.14 839.14 7.166367E+02 1.827481E-02 2.061309E+02 + 191.14 839.14 7.162972E+02 1.822730E-02 2.056778E+02 + 191.14 839.14 7.156193E+02 1.811681E-02 2.066672E+02 + 191.14 839.14 7.136020E+02 1.796227E-02 2.078135E+02 + 191.14 839.14 7.102769E+02 1.775943E-02 2.092069E+02 + 191.14 839.14 7.037448E+02 1.737619E-02 2.118732E+02 + 191.14 839.14 6.973246E+02 1.703415E-02 2.141400E+02 + 191.14 839.14 6.910441E+02 1.678595E-02 2.153654E+02 + 191.14 839.14 6.848203E+02 1.653999E-02 2.165908E+02 + 191.14 839.14 6.786887E+02 1.632090E-02 2.175228E+02 + 191.14 839.14 6.726641E+02 1.613921E-02 2.180310E+02 + 191.14 839.14 6.666930E+02 1.595913E-02 2.185392E+02 + 191.14 839.14 6.607749E+02 1.578065E-02 2.190475E+02 + 191.14 839.14 6.549093E+02 1.560375E-02 2.195557E+02 + 191.14 839.14 6.490958E+02 1.542842E-02 2.200639E+02 + 191.14 839.14 6.433869E+02 1.528464E-02 2.201856E+02 + 191.14 839.14 6.377283E+02 1.514217E-02 2.203066E+02 + 191.14 839.14 6.321195E+02 1.500095E-02 2.204276E+02 + 191.14 839.14 6.265601E+02 1.486098E-02 2.205486E+02 + 191.14 839.14 6.210494E+02 1.472223E-02 2.206697E+02 + 191.14 839.14 6.155873E+02 1.458471E-02 2.207907E+02 + 191.14 839.14 6.101732E+02 1.444840E-02 2.209117E+02 + 191.14 839.14 6.048068E+02 1.431328E-02 2.210327E+02 + 191.14 839.14 5.994875E+02 1.417936E-02 2.211537E+02 + 191.14 839.14 5.942170E+02 1.404808E-02 2.212530E+02 + 191.14 839.14 5.890038E+02 1.392606E-02 2.212321E+02 + 191.14 839.14 5.838364E+02 1.380511E-02 2.212112E+02 + 191.14 839.14 5.787144E+02 1.368522E-02 2.211903E+02 + 191.14 839.14 5.736372E+02 1.356638E-02 2.211694E+02 + 191.14 839.14 5.686046E+02 1.344858E-02 2.211485E+02 + 191.14 839.14 5.636161E+02 1.333182E-02 2.211275E+02 + 191.14 839.14 5.586714E+02 1.321608E-02 2.211066E+02 + 191.14 839.14 5.537701E+02 1.310136E-02 2.210857E+02 + 191.14 839.14 5.489118E+02 1.298765E-02 2.210648E+02 + 191.14 839.14 5.440961E+02 1.287493E-02 2.210439E+02 + 191.14 839.14 5.393227E+02 1.276320E-02 2.210230E+02 + 191.14 839.14 5.345911E+02 1.265245E-02 2.210021E+02 + 191.14 839.14 5.299010E+02 1.254267E-02 2.209812E+02 + 191.14 839.14 5.252522E+02 1.243386E-02 2.209602E+02 + 191.14 839.14 5.206441E+02 1.232600E-02 2.209393E+02 + 191.14 839.14 5.160764E+02 1.221909E-02 2.209184E+02 + 191.14 839.14 5.115410E+02 1.211385E-02 2.208768E+02 + 191.14 839.14 5.070372E+02 1.201034E-02 2.208127E+02 + 191.14 839.14 5.025731E+02 1.190774E-02 2.207486E+02 + 191.14 839.14 4.981482E+02 1.180605E-02 2.206845E+02 + 191.14 839.14 4.937623E+02 1.170525E-02 2.206204E+02 + 191.14 839.14 4.894151E+02 1.160534E-02 2.205563E+02 + 191.14 839.14 4.851060E+02 1.150631E-02 2.204922E+02 + 191.14 839.14 4.808350E+02 1.140815E-02 2.204281E+02 + 191.14 839.14 4.766015E+02 1.131086E-02 2.203640E+02 + 191.14 839.14 4.724054E+02 1.121442E-02 2.202999E+02 + 191.14 839.14 4.682461E+02 1.111883E-02 2.202358E+02 + 191.14 839.14 4.641235E+02 1.102408E-02 2.201717E+02 + 191.14 839.14 4.600372E+02 1.093017E-02 2.201076E+02 + 191.14 839.14 4.559868E+02 1.083708E-02 2.200435E+02 + 191.14 839.14 4.519722E+02 1.074481E-02 2.199794E+02 + 191.14 839.14 4.479928E+02 1.065336E-02 2.199153E+02 + 191.14 839.14 4.440485E+02 1.056271E-02 2.198512E+02 + 191.14 839.14 4.401389E+02 1.047286E-02 2.197870E+02 + 191.14 839.14 4.362638E+02 1.038380E-02 2.197229E+02 + 191.14 839.14 4.324228E+02 1.029552E-02 2.196588E+02 + 191.14 839.14 4.286155E+02 1.020802E-02 2.195947E+02 + 191.14 839.14 4.248418E+02 1.012129E-02 2.195306E+02 + 191.14 839.14 4.211013E+02 1.003533E-02 2.194665E+02 + 191.14 839.14 4.173938E+02 9.950121E-03 2.194024E+02 + 191.14 839.14 4.137189E+02 9.865662E-03 2.193383E+02 + 191.14 839.14 4.100764E+02 9.781948E-03 2.192742E+02 + 191.14 839.14 4.064553E+02 9.697684E-03 2.192265E+02 + 191.14 839.14 4.028580E+02 9.613179E-03 2.191914E+02 + 191.14 839.14 3.992926E+02 9.529422E-03 2.191563E+02 + 191.14 839.14 3.957588E+02 9.446407E-03 2.191213E+02 + 191.14 839.14 3.922562E+02 9.364125E-03 2.190862E+02 + 191.14 839.14 3.887847E+02 9.282573E-03 2.190511E+02 + 191.14 839.14 3.853438E+02 9.201742E-03 2.190160E+02 + 191.14 839.14 3.819334E+02 9.121627E-03 2.189809E+02 + 191.14 839.14 3.785532E+02 9.042220E-03 2.189458E+02 + 191.14 839.14 3.752029E+02 8.963516E-03 2.189108E+02 + 191.14 839.14 3.718823E+02 8.885508E-03 2.188757E+02 + 191.14 839.14 3.685910E+02 8.808192E-03 2.188406E+02 + 191.14 839.14 3.653289E+02 8.731560E-03 2.188055E+02 + 191.14 839.14 3.620956E+02 8.655605E-03 2.187704E+02 + 191.14 839.14 3.588910E+02 8.580323E-03 2.187354E+02 + 191.14 839.14 3.557147E+02 8.505708E-03 2.187003E+02 + 191.14 839.14 3.525665E+02 8.431751E-03 2.186652E+02 + 191.14 839.14 3.494463E+02 8.358451E-03 2.186301E+02 + 191.14 839.14 3.463535E+02 8.285798E-03 2.185950E+02 + 191.14 839.14 3.432882E+02 8.213789E-03 2.185600E+02 + 191.14 839.14 3.402500E+02 8.142417E-03 2.185249E+02 + 191.14 839.14 3.372387E+02 8.071678E-03 2.184898E+02 + 191.14 839.14 3.342541E+02 8.001563E-03 2.184547E+02 + 191.14 839.14 3.312959E+02 7.932070E-03 2.184196E+02 + 191.14 839.14 3.283638E+02 7.863191E-03 2.183846E+02 + 191.14 839.14 3.254577E+02 7.794921E-03 2.183495E+02 + 191.14 839.14 3.225773E+02 7.727257E-03 2.183144E+02 + 191.14 839.14 3.197224E+02 7.660191E-03 2.182793E+02 + 191.14 839.14 3.168928E+02 7.593718E-03 2.182442E+02 + 191.14 839.14 3.140882E+02 7.527835E-03 2.182091E+02 + 191.14 839.14 3.113084E+02 7.462533E-03 2.181741E+02 + 191.14 839.14 3.085533E+02 7.397810E-03 2.181390E+02 + 191.14 839.14 3.058225E+02 7.333660E-03 2.181039E+02 + 191.14 839.14 3.031159E+02 7.270077E-03 2.180688E+02 + 191.14 839.14 3.004332E+02 7.207057E-03 2.180337E+02 + 191.14 839.14 2.977744E+02 7.144595E-03 2.179987E+02 + 191.14 839.14 2.951389E+02 7.082686E-03 2.179636E+02 + 191.14 839.14 2.924969E+02 7.022295E-03 2.178549E+02 + 191.14 839.14 2.898659E+02 6.962860E-03 2.177149E+02 + 191.14 839.14 2.872585E+02 6.903958E-03 2.175748E+02 + 191.14 839.14 2.846747E+02 6.845587E-03 2.174347E+02 + 191.14 839.14 2.821140E+02 6.787741E-03 2.172946E+02 + 191.14 839.14 2.795763E+02 6.730415E-03 2.171546E+02 + 191.14 839.14 2.770615E+02 6.673604E-03 2.170145E+02 + 191.14 839.14 2.745693E+02 6.617305E-03 2.168744E+02 + 191.14 839.14 2.720995E+02 6.561512E-03 2.167343E+02 + 191.14 839.14 2.696520E+02 6.506221E-03 2.165942E+02 + 191.14 839.14 2.672265E+02 6.451428E-03 2.164542E+02 + 191.14 839.14 2.648227E+02 6.397127E-03 2.163141E+02 + 191.14 839.14 2.624406E+02 6.343314E-03 2.161740E+02 + 191.14 839.14 2.600800E+02 6.289986E-03 2.160339E+02 + 191.14 839.14 2.577405E+02 6.237138E-03 2.158938E+02 + 191.14 839.14 2.554221E+02 6.184764E-03 2.157538E+02 + 191.14 839.14 2.531246E+02 6.132862E-03 2.156137E+02 + 191.14 839.14 2.508477E+02 6.081427E-03 2.154736E+02 + 191.14 839.14 2.485913E+02 6.030454E-03 2.153335E+02 + 191.14 839.14 2.463552E+02 5.979940E-03 2.151935E+02 + 191.14 839.14 2.441392E+02 5.929881E-03 2.150534E+02 + 191.14 839.14 2.419432E+02 5.880271E-03 2.149133E+02 + 191.14 839.14 2.397669E+02 5.831108E-03 2.147732E+02 + 191.14 839.14 2.376101E+02 5.782387E-03 2.146331E+02 + 191.14 839.14 2.354728E+02 5.734105E-03 2.144931E+02 + 191.14 839.14 2.333547E+02 5.686256E-03 2.143530E+02 + 191.14 839.14 2.312557E+02 5.638838E-03 2.142129E+02 + 191.14 839.14 2.291755E+02 5.591847E-03 2.140728E+02 + 191.14 839.14 2.271141E+02 5.545278E-03 2.139328E+02 + 191.14 839.14 2.250712E+02 5.499128E-03 2.137927E+02 + 191.14 839.14 2.230466E+02 5.453393E-03 2.136526E+02 + 191.14 839.14 2.210403E+02 5.408070E-03 2.135125E+02 + 191.14 839.14 2.190520E+02 5.363154E-03 2.133724E+02 + 191.14 839.14 2.170816E+02 5.318642E-03 2.132324E+02 + 191.14 839.14 2.151290E+02 5.274531E-03 2.130923E+02 + 191.14 839.14 2.131938E+02 5.230817E-03 2.129522E+02 + 191.14 839.14 2.112761E+02 5.187495E-03 2.128121E+02 + 191.14 839.14 2.093757E+02 5.144564E-03 2.126720E+02 + 191.14 839.14 2.074924E+02 5.102018E-03 2.125320E+02 + 191.14 839.14 2.056260E+02 5.059855E-03 2.123919E+02 + 191.14 839.14 2.037763E+02 5.018072E-03 2.122518E+02 + 191.14 839.14 2.019433E+02 4.976665E-03 2.121117E+02 + 191.14 839.14 2.001268E+02 4.935629E-03 2.119716E+02 + 191.14 839.14 1.983267E+02 4.894963E-03 2.118316E+02 + 191.14 839.14 1.965427E+02 4.854663E-03 2.116915E+02 + 191.14 839.14 1.947748E+02 4.814725E-03 2.115514E+02 + 191.14 839.14 1.930228E+02 4.775146E-03 2.114113E+02 + 191.14 839.14 1.912836E+02 4.735827E-03 2.112695E+02 + 191.14 839.14 1.894983E+02 4.694828E-03 2.110901E+02 + 191.14 839.14 1.877297E+02 4.654211E-03 2.109107E+02 + 191.14 839.14 1.859776E+02 4.613973E-03 2.107313E+02 + 191.14 839.14 1.842419E+02 4.574111E-03 2.105519E+02 + 191.14 839.14 1.825223E+02 4.534621E-03 2.103726E+02 + 191.14 839.14 1.808188E+02 4.495499E-03 2.101932E+02 + 191.14 839.14 1.791312E+02 4.456743E-03 2.100138E+02 + 191.14 839.14 1.774594E+02 4.418348E-03 2.098344E+02 + 191.14 839.14 1.758031E+02 4.380311E-03 2.096550E+02 + 191.14 839.14 1.741624E+02 4.342630E-03 2.094756E+02 + 191.14 839.14 1.725369E+02 4.305300E-03 2.092962E+02 + 191.14 839.14 1.709266E+02 4.268319E-03 2.091169E+02 + 191.14 839.14 1.693313E+02 4.231683E-03 2.089375E+02 + 191.14 839.14 1.677509E+02 4.195388E-03 2.087581E+02 + 191.14 839.14 1.661853E+02 4.159433E-03 2.085787E+02 + 191.14 839.14 1.646343E+02 4.123813E-03 2.083993E+02 + 191.14 839.14 1.630977E+02 4.088526E-03 2.082199E+02 + 191.14 839.14 1.615755E+02 4.053568E-03 2.080406E+02 + 191.14 839.14 1.600675E+02 4.018935E-03 2.078612E+02 + 191.14 839.14 1.585736E+02 3.984627E-03 2.076818E+02 + 191.14 839.14 1.570936E+02 3.950639E-03 2.075024E+02 + 191.14 839.14 1.556275E+02 3.916967E-03 2.073230E+02 + 191.14 839.14 1.541750E+02 3.883610E-03 2.071436E+02 + 191.14 839.14 1.527361E+02 3.850565E-03 2.069642E+02 + 191.14 839.14 1.513106E+02 3.817827E-03 2.067849E+02 + 191.14 839.14 1.498984E+02 3.785396E-03 2.066055E+02 + 191.14 839.14 1.484993E+02 3.753267E-03 2.064261E+02 + 191.14 839.14 1.471134E+02 3.721437E-03 2.062467E+02 + 191.14 839.14 1.457404E+02 3.689905E-03 2.060673E+02 + 191.14 839.14 1.443802E+02 3.658668E-03 2.058879E+02 + 191.14 839.14 1.430326E+02 3.627721E-03 2.057086E+02 + 191.14 839.14 1.416977E+02 3.597064E-03 2.055292E+02 + 191.14 839.14 1.403752E+02 3.566693E-03 2.053498E+02 + 191.14 839.14 1.390651E+02 3.536605E-03 2.051704E+02 + 191.14 839.14 1.377672E+02 3.506798E-03 2.049910E+02 + 191.14 839.14 1.364814E+02 3.477269E-03 2.048116E+02 + 191.14 839.14 1.352076E+02 3.448016E-03 2.046322E+02 + 191.14 839.14 1.339457E+02 3.419036E-03 2.044529E+02 + 191.14 839.14 1.326956E+02 3.390326E-03 2.042735E+02 + 191.14 839.14 1.314571E+02 3.361884E-03 2.040941E+02 + 191.14 839.14 1.302302E+02 3.333708E-03 2.039147E+02 + 191.14 839.14 1.290148E+02 3.305794E-03 2.037353E+02 + 191.14 839.14 1.278107E+02 3.278141E-03 2.035559E+02 + 191.14 839.14 1.266178E+02 3.250747E-03 2.033766E+02 + 191.14 839.14 1.254361E+02 3.223608E-03 2.031972E+02 + 191.14 839.14 1.242654E+02 3.196722E-03 2.030178E+02 + 191.14 839.14 1.231056E+02 3.170087E-03 2.028384E+02 + 191.14 839.14 1.219566E+02 3.143701E-03 2.026590E+02 + 191.14 839.14 1.208184E+02 3.117560E-03 2.024796E+02 + 191.14 839.14 1.196908E+02 3.091665E-03 2.023002E+02 + 191.14 839.14 1.185737E+02 3.066010E-03 2.021209E+02 + 191.14 839.14 1.174671E+02 3.040595E-03 2.019415E+02 Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/launch =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/launch (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/launch (revision 1291) @@ -0,0 +1,42 @@ +#PBS -S /bin/bash +#PBS -j oe +#PBS -q infini +#PBS -N EXOMARS +#PBS -l "nodes=1:ppn=8" +##PBS -l "nodes=2:ppn=4" +##PBS -l "nodes=1:ppn=12" + +##PBS -l walltime=10:30,mem=320kb + +######## openMP threads +#cat $PBS_NODEFILE | sort | uniq > /tmp/coucou +#cat /tmp/coucou +#declare -x OMP_NUM_THREADS=4 + +#declare -x MP_STACK_SIZE=64000000 + +####### prep run +#cd $PBS_O_WORKDIR +#\rm *rsl.* 2> /dev/null +#\rm wrfo* 2> /dev/null +#\rm wrfi* 2> /dev/null +#./ideal.exe +#mv rsl.out.0000 ideal_rsl.out.0000 + +####### run +#mpirun yeye +mpirun -np 4 wrf.exe & + + +################ + +###PBS -l "nodes=2:ppn=8" + +#declare -x OMP_NUM_THREADS=4 +#declare -x MP_STACK_SIZE=64000000 + +#mpirun wrf.exe + +####PBS -m ae +####PBS -M spiga@lmd.jussieu.fr + Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/namelist.input =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/namelist.input (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/namelist.input (revision 1291) @@ -0,0 +1,77 @@ + &time_control + run_seconds = 40700 !! Duration of simulation [LES: !=0] + history_interval_s = 100 !! Frequency of outputs in seconds [LES: !=0] + frames_per_outfile = 37 !! Size of time dimension in files + start_year = 9999 !! Idealized simulation [LES: 9999] + end_year = 9999 !! Idealized simulation [LES: 9999] + restart = .false. !! (*) Output restart files ? + 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 = 10000. !! Altitude above ground of model top [LES: !=0] + time_step = 1 !! Dynamical timestep + time_step_fract_num = 1 !! Additional fraction to time_step: numerator + time_step_fract_den = 2 !! Additional fraction to time_step: denominator + dx = 100 !! (p2) Horizontal resolution + dy = 100 !! (p2) Horizontal resolution (should be equal to dx) + e_we = 45 !! (r)(p2) Number of longitude grid points + e_sn = 45 !! (r)(p2) Number of latitude grid points + e_vert = 71 !! (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 = 199. !! Define constant thermal inertia value [LES: !=0] + init_AL = 0.153 !! Define constant albedo value [LES: !=0] + num_soil_layers = 18 !!! NOUVELLE PHYSIQUE + !! + !! OPTIONAL + !! + radt = 45 !! 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 + / + + &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.10 !! 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/tests_newphyslesfast/run.def =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/run.def (revision 1291) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/tests_newphyslesfast/run.def (revision 1291) @@ -0,0 +1,2 @@ +# some definitions for the physics, in file 'callphys.def' +INCLUDEDEF=callphys.def