Changeset 1664 for LMDZ5/branches/testing/libf/phylmd/cv3p1_closure.F

Timestamp:

Oct 9, 2012, 3:29:15 PM (12 years ago)

Author:

Laurent Fairhead

Message:

Version testing basée sur la r1575

---

Testing release based on r1575

Location:

LMDZ5/branches/testing

Files:

• . (modified) (1 prop)
• libf/phylmd/cv3p1_closure.F (modified) (1 diff)

LMDZ5/branches/testing/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