Changeset 765 for trunk/LMDZ.MARS/libf

Timestamp:

Aug 31, 2012, 9:52:22 AM (12 years ago)

Author:

acolaitis

Message:

Added a low threshold subgrid gustiness in richardson computations for the LES, after observations of unreallistically high HFX values for very low winds conditions.

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• vdif_cd.F (modified) (3 diffs)
• vdifc.F (modified) (1 diff)

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

@@ -1 +1 @@
       SUBROUTINE vdif_cd(ngrid,nlay,pz0,
-     &           pg,pz,pu,pv,wstar,pts,ph,pcdv,pcdh)
+     &           pg,pz,pu,pv,wstar,pts,ph,pcdv,pcdh
+#ifdef MESOSCALE
+     &                ,flag_LES
+#endif
+     &                  )
       IMPLICIT NONE
 c=======================================================================
@@ -48 +52 @@
       REAL, INTENT(IN) :: wstar(ngrid)
       REAL, INTENT(OUT) :: pcdv(ngrid),pcdh(ngrid) ! momentum and heat drag coefficient
+
+#ifdef MESOSCALE
+      LOGICAL, INTENT(IN) :: flag_LES     !! pour LES avec isfflx!=0
+#endif
 
 c   Local:
@@ -152 +160 @@
           zu2(ig)=pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) 
      &     + (log(1.+0.7*wstar(ig) + 2.3*wstar(ig)**2))**2 
-
+#ifdef MESOSCALE
+          if(flag_LES) then
+             zu2(ig)=MAX(zu2(ig),1.)
+          endif
+#endif
 !       zu2(ig)=pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) + (0.5*wstar(ig))**2
 

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

@@ -334 +334 @@
 
       CALL vdif_cd(ngrid,nlay,pz0,g,pzlay,pu,pv,wstar,ptsrf,ph
-     &             ,zcdv_true,zcdh_true)
+     &             ,zcdv_true,zcdh_true
+#ifdef MESOSCALE
+     &                ,flag_LES
+#endif
+     &           )
 
         zu2(:)=pu(:,1)*pu(:,1)+pv(:,1)*pv(:,1)