Changeset 1574

Timestamp:

Sep 20, 2011, 2:09:47 PM (13 years ago)

Author:

jghattas

Message:

Corrected calculation of the variables wbeff and plfc (diagnostics).

Location:

LMDZ5/trunk/libf/phylmd

Files:

• concvl.F (modified) (1 diff)
• cv3p1_closure.F (modified) (1 diff)
• phys_output_mod.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/concvl.F

@@ -248 +248 @@
           cbmf(i) = 0.
           plcl(i) = 0.
-          plfc(i) = 0.
-          wbeff(i) = 0.
           sigd(i) = 0.
          ENDDO
       ENDIF   !(ifrst .EQ. 0)
+
+c Initialisation a chaque pas de temps
+      plfc(:)  = 0.
+      wbeff(:) = 100.
 
       DO k = 1, klev+1

LMDZ5/trunk/libf/phylmd/cv3p1_closure.F

@@ -516 +516 @@
         wb2(il) = sqrt(2.*max(Ale(il)+cin(il),0.))
       ENDDO
-c
-      IF (flag_wb==0) THEN
-        wbeff(:) = wbmax
-      ELSE IF (flag_wb==1) THEN
-        wbeff(1:ncum) = wbmax/(1.+500./(ph(1:ncum,1)-plfc(1:ncum)))
-      ELSE IF (flag_wb==2) THEN
-        wbeff(1:ncum) = wbmax*(0.01*(ph(1:ncum,1)-plfc(1:ncum)))**2
-      ENDIF
-c
+
+      DO il = 1, ncum
+         IF (plfc(il) .lt. 100.) THEN
+c        This is an irealistic value for plfc => no calculation of wbeff
+            wbeff(il) = 100.1
+         ELSE
+c        Calculate wbeff
+            IF (flag_wb==0) THEN
+               wbeff(il) = wbmax
+            ELSE IF (flag_wb==1) THEN
+               wbeff(il) = wbmax/(1.+500./(ph(il,1)-plfc(il)))
+            ELSE IF (flag_wb==2) THEN
+               wbeff(il) = wbmax*(0.01*(ph(il,1)-plfc(il)))**2
+            ENDIF
+         END IF
+      END DO
+
+
       DO il = 1,ncum
 cjyg    Modification du coef de wb*wb pour conformite avec papier Wake

LMDZ5/trunk/libf/phylmd/phys_output_mod.F90

@@ -1213 +1213 @@
                      o_plfc%flag,o_plfc%name, "Level of Free Convection", "hPa")
                 CALL histdef2d(iff,clef_stations(iff), &
-                     o_wbeff%flag,o_wbeff%name, "Conv. updraft velocity at LFC", "m/s")
+                     o_wbeff%flag,o_wbeff%name, "Conv. updraft velocity at LFC (<100)", "m/s")
              end if
              CALL histdef2d(iff,clef_stations(iff), &