Changeset 265

Timestamp:

Aug 5, 2011, 12:54:31 PM (14 years ago)

Author:

acolaitis

Message:

Corrected wrong definition for u* in mesoscale model

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

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

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

@@ -8 +8 @@
 
 ! New SL parametrization, ust is more accurately computed in vdif_cd :
-          ust(ig) = zcdv(ig)
+        ust(ig) = sqrt(zcdv(ig)*(pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1)))
 
          ENDDO   

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

@@ -74 +74 @@
       REAL pz0t        ! initial thermal roughness length. (local)
       REAL ric         ! critical richardson number
+      REAL reynolds(ngrid)    ! reynolds number for UBL
       REAL prandtl(ngrid)     ! prandtl number for UBL
-      REAL reynolds(ngrid)    ! reynolds number for UBL
       REAL pz0tcomp(ngrid)     ! computed z0t
       REAL ite
@@ -101 +101 @@
 c phim = 1+betam*zeta   or   (1-bm*zeta)**am
 c phih = alphah + betah*zeta    or   alphah(1.-bh*zeta)**ah
-c am=-0.5, ah=-0.25
+c am=-0.25, ah=-0.5
 
       pz0t = 0.     ! for the sake of simplicity
@@ -199 +199 @@
        
       ENDDO
-
+!
 !      print*,'new : cd,ch; ',pcdv,pcdh
 

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

@@ -5 +5 @@
      $                pdufi,pdvfi,pdhfi,pdqfi,pfluxsrf,
      $                pdudif,pdvdif,pdhdif,pdtsrf,pq2,
-     $                pdqdif,pdqsdif,wmax,zcdv)
+     $                pdqdif,pdqsdif,wmax,zcdv_true)
       IMPLICIT NONE
 
@@ -253 +253 @@
 !     &              'heat drag','no units',
 !     &                         2,zcdh_true)
-!        call WRITEDIAGFI(ngridmx,'u*',
+!        call WRITEDIAGFI(ngridmx,'u*2/u',
 !     &              'friction velocity','m/s',
 !     &                         2,zcdv)
-!        call WRITEDIAGFI(ngridmx,'T*',
+!        call WRITEDIAGFI(ngridmx,'u*T*/(T0-T)',
 !     &              'friction temperature','m/s',
 !     &                         2,zcdh)