Changeset 499 for trunk/LMDZ.MARS/libf/phymars/vdifc.F

Timestamp:

Jan 16, 2012, 10:16:46 AM (13 years ago)

Author:

acolaitis

Message:

M 498 libf/phymars/calltherm_interface.F90
^{------------------> Added computation of wstar. Wmax is no longer an output.}

M 498 libf/phymars/physiq.F
^{------------------> Minor changes.}

M 498 libf/phymars/vdif_cd.F
^{------------------> Changed computation of bulk Richardson number to allow more unstable configurations.}

M 498 libf/phymars/vdifc.F
^{------------------> Changed computation of surface conductance coefficients following a new approach. Comparisons with LES are much better.}

File:

• trunk/LMDZ.MARS/libf/phymars/vdifc.F (modified) (5 diffs)

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

@@ -5 +5 @@
      $                pdufi,pdvfi,pdhfi,pdqfi,pfluxsrf,
      $                pdudif,pdvdif,pdhdif,pdtsrf,pq2,
-     $                pdqdif,pdqsdif,wmax,zcdv_true,zcdh_true,hfmax)
+     $                pdqdif,pdqsdif,wstar,zcdv_true,zcdh_true,hfmax)
       IMPLICIT NONE
 
@@ -60 +60 @@
 c    Argument added to account for subgrid gustiness :
 
-      REAL wmax(ngridmx), hfmax(ngridmx)
+      REAL wstar(ngridmx), hfmax(ngridmx), zi(ngridmx)
 
 c    Traceurs :
@@ -313 +313 @@
 c       ---------------------
 
-      CALL vdif_cd(ngrid,nlay,pz0,g,pzlay,pu,pv,wmax,ptsrf,ph
+      CALL vdif_cd(ngrid,nlay,pz0,g,pzlay,pu,pv,wstar,ptsrf,ph
      &             ,zcdv_true,zcdh_true)
 
@@ -319 +319 @@
 
         IF (callrichsl) THEN
-       zcdv(:)=zcdv_true(:)*sqrt(zu2(:)+(6.*hfmax(:))**2)     ! subgrid gustiness added by quadratic interpolation with a coeff beta
-       zcdh(:)=zcdh_true(:)*sqrt(zu2(:)+(10.*hfmax(:))**2)     ! LES comparisons. This parameter is linked to the thermals settings)
-
-!       zcdv(:)=zcdv_true(:)*sqrt(zu2(:)+(0.7*wmax(:))**2)     ! subgrid gustiness added by quadratic interpolation with a coeff beta
-!       zcdh(:)=zcdh_true(:)*sqrt(zu2(:)+(1.2*wmax(:))**2)     ! LES comparisons. This parameter is linked to the thermals settings)
-
-        ust(:)=sqrt(zcdv_true(:))*sqrt(zu2(:)+(6.*hfmax(:))**2)
+       zcdv(:)=zcdv_true(:)*sqrt(zu2(:)+
+     &  (14.75*sqrt(hfmax(:)) - 10.*hfmax(:) + 1.7*hfmax(:)**2)**2)     ! subgrid gustiness added by quadratic interpolation with a coeff beta
+       zcdh(:)=zcdh_true(:)*sqrt(zu2(:)+
+     &  (20.6*hfmax(:) - 22.8*hfmax(:)**2 + 9.*hfmax(:)**4)**2)     ! LES comparisons. This parameter is linked to the thermals settings)
+
+        ust(:)=sqrt(zcdv_true(:))*sqrt(zu2(:)+
+     &  (14.75*sqrt(hfmax(:)) - 10.*hfmax(:) + 1.7*hfmax(:)**2)**2)
         tst(:)=(ph(:,1)-ptsrf(:))*zcdh(:)/ust(:)
-!        ust(:)=sqrt(zcdv_true(:))*sqrt(zu2(:)+(0.7*wmax(:))**2)
-!        tst(:)=(ph(:,1)-ptsrf(:))*zcdh(:)/ust(:)
 !        ust(:)=sqrt(zcdv_true(:))*sqrt(zu2(:))
 !        tst(:)=(ph(:,1)-ptsrf(:))*zcdh_true(:)/sqrt(zcdv_true(:))
@@ -334 +332 @@
         zcdv(:)=zcdv_true(:)*sqrt(zu2(:))     ! 1 / bulk aerodynamic momentum conductance 
         zcdh(:)=zcdh_true(:)*sqrt(zu2(:))     ! 1 / bulk aerodynamic heat conductance
-!        ust(:)=sqrt(zcdv_true(:))*sqrt(zu2(:))
-!        tst(:)=(ph(:,1)-ptsrf(:))*zcdh(:)/ust(:)
         ENDIF