Changeset 263 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Aug 5, 2011, 11:56:23 AM (13 years ago)

Author:

aslmd

Message:

MESOSCALE: Something was broken for MESOSCALE model in r259-r260 (flag_LES is not known in vdifc). Moved this to meso_inc/meso_inc_les.F. Also there was a problem with the new surface layer if calltherm is false because then wmax has no value... Fixed.

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• meso_inc/meso_inc_les.F (modified) (1 diff)
• physiq.F (modified) (2 diffs)
• vdifc.F (modified) (1 diff)

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

@@ -28 +28 @@
 
          DO ig=1,ngrid
+          wmax_th(ig)=0.  !! no additional gustiness needed in surface layer (see vdifc.F)
           DO l=1,nlayer
             zdvdif(ig,l) = 0.

trunk/LMDZ.MARS/libf/phymars/physiq.F

@@ -322 +322 @@
       REAL lmax_th_out(ngridmx),zmax_th(ngridmx)
       REAL, SAVE :: wmax_th(ngridmx)
-      REAL, SAVE :: hfmax_th(ngridmx)
+      REAL, SAVE :: hfmax_th(ngridmx)  !! AS: pourquoi en save ?
       REAL pdu_th(ngridmx,nlayermx),pdv_th(ngridmx,nlayermx)
       REAL pdt_th(ngridmx,nlayermx),pdq_th(ngridmx,nlayermx,nqmx)
@@ -784 +784 @@
         else   !of if calltherm
         lmax_th(:)=0
+        wmax_th(:)=0.
         end if
 

trunk/LMDZ.MARS/libf/phymars/vdifc.F

@@ -234 +234 @@
       CALL vdif_cd(ngrid,nlay,pz0,g,pzlay,pu,pv,ptsrf,ph
      &             ,zcdv_true,zcdh_true)
-
-#ifdef MESOSCALE
-      IF (flag_LES) THEN
-      wmax(:)=0.
-      ENDIF
-#endif
 
       DO ig=1,ngrid