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

Timestamp:

Jan 24, 2012, 6:31:51 PM (13 years ago)

Author:

acolaitis

Message:

Mesoscale/Les?. Added improvments to mode 21. Corrected variable names for new subgrid gustiness definitions.

File:

• trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_les.F (modified) (2 diffs)

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

@@ -20 +20 @@
             sensheat(ig) = (pplay(ig,1)/(r*pt(ig,1)))*cpp
      &        *sqrt(pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1)
-     &        + (1.2*wmax_th(ig))**2)
+     &        + (log(1.+0.7*wstar(ig) + 2.3*wstar(ig)**2)**2))
      &        *zcdh(ig)*(tsurf(ig)-zh(ig,1))
 
 ! New SL parametrization, ustar is more accurately computed in vdif_cd :
+
             ustar(ig) = sqrt(zcdv(ig)*
-     &   (pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) + (0.7*wmax_th(ig))**2)
+     &   (pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) 
+     &     + (log(1.+0.7*wstar(ig) + 2.3*wstar(ig)**2)**2))
      &                     )
 
@@ -49 +51 @@
 
          DO ig=1,ngrid
-          wmax_th(ig)=0.  !! no additional gustiness needed in surface layer (see vdifc.F)
+          wstar(ig)=0.  !! no additional gustiness needed in surface layer (see vdifc.F)
           DO l=1,nlayer
             zdvdif(ig,l) = 0.