Changeset 276 for trunk/LMDZ.MARS/libf

Timestamp:

Sep 1, 2011, 3:11:19 PM (13 years ago)

Author:

acolaitis

Message:

Minor modification in vdif_cd.F. Important modification to thermals scheme : added specific set of parameter values for use in mesoscale model. Thermals will be usable in Mesoscale after an other verification is made about subrgrid gustiness.

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• calltherm_mars.F90 (modified) (1 diff)
• vdif_cd.F (modified) (9 diffs)

trunk/LMDZ.MARS/libf/phymars/calltherm_mars.F90

@@ -87 +87 @@
 !  ===================
  
-         r_aspect_thermals=2.   ! ultimately conrols the amount of mass going through the thermals
+!         r_aspect_thermals     ! ultimately conrols the amount of mass going through the thermals
                                 ! decreasing it increases the thermals effect. Tests at gcm resolution
                                 ! shows that too low values destabilize the model
                                 ! when changing this value, one should check that the surface layer model
                                 ! outputs the correct Cd*u and Ch*u through changing the gustiness coefficient beta
+
+
+#ifdef MESOSCALE
+         nsplit_thermals=2
+         r_aspect_thermals=0.7
+#else
          nsplit_thermals=40
+         r_aspect_thermals=2.
+#endif
+
          call getin("nsplit_thermals",nsplit_thermals)
 

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

@@ -61 +61 @@
 c    Local(2):
 c    ---------
+      REAL z1,zcd0
 
       REAL rib(ngrid)  ! Bulk Richardson number
@@ -81 +82 @@
       REAL ite
       REAL residual
-      REAL zu2
+      REAL zu2(ngrid)
 c-----------------------------------------------------------------------
 c   couche de surface:
@@ -90 +91 @@
 c Original formulation :
 
-!      DO ig=1,ngrid
-!         z1=1.E+0 + pz(ig,1)/pz0(ig)
-!         zcd0=karman/log(z1)
-!         zcd0=zcd0*zcd0
-!         pcdv(ig)=zcd0
-!         pcdh(ig)=zcd0
-!      ENDDO
-
+c      DO ig=1,ngrid
+c         z1=1.E+0 + pz(ig,1)/pz0(ig)
+c         zcd0=karman/log(z1)
+c         zcd0=zcd0*zcd0
+c         pcdv(ig)=zcd0
+c         pcdh(ig)=zcd0
+c      ENDDO
+     
 !      print*,'old : cd,ch; ',pcdv,pcdh
 
@@ -109 +110 @@
       pz0tcomp(:) = 0.1*pz0(:)
       rib(:)=0.
+
       pcdv(:)=0.
       pcdh(:)=0.
@@ -143 +145 @@
 !         zu2=pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1)
 !         zu2=MAX(pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1),wmax(ig)**2)
-         zu2= pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) + (1.*wmax(ig))**2
+       zu2(ig)= pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) + (1.*wmax(ig))**2
 
                ! we add the wmax to simulate
@@ -157 +159 @@
      &      *(((log(pz(ig,1)/pz0(ig)))**2)/(log(pz(ig,1)/pz0t)))
      &      *(ph(ig,1)-pts(ig))
-     &  /zu2
+     &  /zu2(ig)
 
          else
@@ -195 +197 @@
       endif
 
-       reynolds(ig)=karman*sqrt(fm(ig))*sqrt(pu(ig,1)**2 + pv(ig,1)**2)
+       reynolds(ig)=karman*sqrt(fm(ig))
+     &      *sqrt(zu2(ig))
+c     &      *sqrt(pu(ig,1)**2 + pv(ig,1)**2)
      &       *pz0(ig)/(log(z1z0)*nu)
        pz0tcomp(ig)=pz0(ig)*exp(-karman*7.3*
@@ -229 +233 @@
 !     &              'Prandtl nb','no units',
 !     &                         0,prandtl)
-!        call WRITEDIAGFI(ngrid,'Re',
+!       call WRITEDIAGFI(ngrid,'Re',
 !     &              'Reynolds nb','no units',
 !     &                         0,reynolds)
@@ -236 +240 @@
 !     &                         2,pz0tcomp)
 
+
       RETURN
       END