Changeset 503

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.

Location:

trunk

Files:

• LMDZ.MARS/libf/phymars/meso_inc/meso_inc_les.F (modified) (2 diffs)
• MESOSCALE/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM (modified) (1 diff)
• MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F (modified) (1 diff)

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.

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM

@@ -107 +107 @@
 state  real  TAU_ICE   ij   misc  1  -  rd   "TAU_ICE"   "CLOUD OD at 825 cm-1 TES"        ""        #SAVEMARS2 tauTES
 state  real  PDTZ      ikj  misc  1  -  rd   "PDT"       "TEMP TENDENCY"                   "K s-1"   #SAVEMARS3 pdt
-state  real  LMAX_TH   ij   misc  1  -  rd   "LMAX_TH"   "MAXIMUM LEVEL REACHED IN TH"     ""        #SAVEMARS2 lmax_th_out
-state  real  HFMAX_TH  ij   misc  1  -  rd   "HFMAX_TH"  "MAXIMUM TH HEAT FLUX"            "m.K/s"   #SAVEMARS2 hfmax_th
-state  real  WMAX_TH   ij   misc  1  -  rd   "WMAX_TH"   "MAXIMUM TH VERTICAL VELOCITY"    "m/s"     #SAVEMARS2 wmax_th
+state  real  LMAX_TH   ij   misc  1  -  rhd   "LMAX_TH"   "MAXIMUM LEVEL REACHED IN TH"     ""        #SAVEMARS2 lmax_th_out
+state  real  HFMAX_TH  ij   misc  1  -  rhd   "HFMAX_TH"  "MAXIMUM TH HEAT FLUX"            "m.K/s"   #SAVEMARS2 hfmax_th
+state  real  WSTAR     ij   misc  1  -  rhd  "WSTAR"   "FREE CONVECTION VELOCITY FROM TH"    "m/s"     #SAVEMARS2 wstar
 state  real  Z0SET     ij   misc  1  -  rd   "Z0SET"     "SET SURFACE ROUGHNESS"           "m"       #SAVEMARS2 z0
 

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F

@@ -1303 +1303 @@
       ENDDO
       IF(lctindex .eq. 0) print*, 'probleme index'
-        ziindex=zipbl(lctindex)
+        IF ((zilct(lctindex) .gt. 12.) .and.( zilct(lctindex) .lt. 18)) THEN
+            ziindex=MAX(0.,zipbl(lctindex)-7.)
+        ELSE
+            ziindex=0.
+        ENDIF
 ENDIF