Index: /trunk/LMDZ.MARS/README
===================================================================
--- /trunk/LMDZ.MARS/README	(revision 566)
+++ /trunk/LMDZ.MARS/README	(revision 567)
@@ -1444,2 +1444,122 @@
    value was asked to user, but twice! Once is enough.
 
+== 07/03/12 == AC
+********************************************************
+>> UPDATE ON PARAMETERS TO SET TO USE THE THERMALS MODEL 
+********************************************************
+
+>> IN CALLPHYS.DEF <<
+---------------------
+
+- calltherm = .true.       > activates the thermals
+- callrichsl = .true.      > activate the Richardson Monin Obukhov surface layer
+with subgrid gustiness, and the yamada4 diffusion model from earth (corrected version of vdif_kc)
+- leave calladj = .true.   > so that the convective adjustment may work ABOVE the PBL.
+
+>> IN RUN.DEF <<
+----------------
+
+- 96 physicial timesteps per day (the internal sub-timestep of thermals will adapt
+so that it costs only 10% more CPU time to go from "48pdt+thermals" to "96pdt+thermals".)
+> this has to be respected because of a coupling between the radiative transfer and thermals
+> if you really wish to run with 48 pdt/day, you may want to deactive the thermals...
+
+>> USEFULL DIAGNOSTICS <<
+-------------------------
+
+- activating the thermals will automatically add 3 2D-outputs to diagfi.nc files:
+> zmax_th : interpolated boundary layer height (m)
+> hfmax_th : maximum vertical eddy heat flux reached in thermals (K.m/s)
+> wstar : free convection velocity reached in thermals (m/s)
+
+A rule of thumb for determining the mean velocities in subgrid-scale structures is:
+
+Mean updraft velocity : ~ wstar
+Maximum updraft velocity : ~ 2 or 3.*wstar
+Mean downdraft velocity : ~ 0.66*wstar
+Maximum absolute value of downdraft velocity : ~ 2.*wstar
+
+>> AT COMPILATION TIME <<
+-------------------------
+
+- +1 vertical level during the compilation, and then using the updated z2sig.def
+
+>> UPDATED Z2SIG.DEF <<
+-----------------------
+
+- one level has been added in the PBL, so that to reach 150km you have to compile with
+33 vertical levels instead of 32.
+
+10.00000     H: atmospheric scale height (km) (used as a reference only)
+0.0020736
+0.015735
+0.0600151
+0.163344
+0.36358
+0.708007
+1.25334 
+2.06571
+3.22069
+4.80327
+6.90787
+9.63834
+13.108
+18.9666
+25.0626
+31.5527
+38.4369
+45.4369
+52.4369
+59.4369
+66.4369
+73.4369
+80.4369
+87.4369
+94.4369
+101.4369
+108.437
+115.437
+122.437
+129.437
+136.437
+143.437
+150.437
+157.437
+164.437
+171.437
+178.437
+185.437
+192.437
+199.437
+206.437
+213.437
+220.437
+227.437
+234.437
+241.437
+248.437
+255.437
+262.437
+269.437
+276.437
+283.437
+290.437
+297.437
+304.437
+311.437
+318.437
+325.437
+332.437
+339.437
+346.437
+353.437
+360.437
+367.437
+374.437
+381.437
+388.437
+395.437
+
+********************************************************
+>> END OF UPDATE ON PARAMETERS FOR THE THERMALS MODEL
+********************************************************