Changeset 256 for trunk/LMDZ.MARS/libf/phymars/meso_inc

Timestamp:

Aug 3, 2011, 4:49:20 PM (13 years ago)

Author:

acolaitis

Message:

--- AC 03/08/2011 ---

M 255 libf/phymars/physiq.F
<> Modified to interface new SL parametrization

M 255 libf/phymars/vdif_cd.F
<> New SL parametrization based on a bulk Richardson Monin-Obukhov theory formulation

Stability functions are taken from D.E. England et al. (95)
Similarity functions coefficients based on Dyer and Hicks (70).
Includes thermal roughness length computation, heat and momentum drag coefficient computation
Can be used to output hydrodynamic-related SL quantities (bulk Richardson, turbulent Prandtl number estimation, Reynolds number)

M 255 libf/phymars/thermcell_dqupdown.F90
<> Minor modification to suit picky compilers

M 255 libf/phymars/vdifc.F
<> Now takes into account sub-grid gustiness, evaluated from thermals activity (it's proxy being the maximum vertical velocity)

M 255 libf/phymars/meso_inc/meso_inc_les.F
<> Minor modification : u* is now taken from the new vdifc and not recomputed from a simple law

---

File:

• trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_les.F (modified) (1 diff)

trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_les.F

@@ -3 +3 @@
           hfx(ig) = zflubid(ig)-capcal(ig)*zdtsdif(ig)
           !! u star in similarity theory in m/s
-          ust(ig) = 0.4
-     .               * sqrt( pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) )
-     .               / log( 1.E+0 + zzlay(ig,1)/z0_default )
+!          ust(ig) = 0.4
+!     .               * sqrt( pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) )
+!     .               / log( 1.E+0 + zzlay(ig,1)/z0_default )
+
+! New SL parametrization, ust is more accurately computed in vdif_cd :
+          ust(ig) = zcdv(ig)
+
          ENDDO   
-
 !         write (*,*) 'PHYS HFX cp zdts', hfx(100), zflubid(100), 
 !     .       capcal(100),