Changeset 3571 for LMDZ6/trunk/libf/phylmd

Timestamp:

Oct 1, 2019, 6:10:13 AM (5 years ago)

Author:

fhourdin

Message:

iflag_wb devient un real pour pouvoir lui donner une valeur continue
quand on l'utilise comme valeur de wbmax en surface.
Pas très joli mais important pour le tuning automatique pour lequel on
genere des valeurs continues de wb_max.

Location:

LMDZ6/trunk/libf/phylmd

Files:

• cv3p1_closure.F90 (modified) (1 diff)
• cv3p2_closure.F90 (modified) (1 diff)
• cv3param.h (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/cv3p1_closure.F90

@@ -537 +537 @@
     ELSE
       ! Calculate wbeff
-      IF (flag_wb==0) THEN
+      IF (NINT(flag_wb)==0) THEN
         wbeff(il) = wbmax
-      ELSE IF (flag_wb==1) THEN
+      ELSE IF (NINT(flag_wb)==1) THEN
         wbeff(il) = wbmax/(1.+500./(ph(il,1)-plfc(il)))
-      ELSE IF (flag_wb==2) THEN
+      ELSE IF (NINT(flag_wb)==2) THEN
         wbeff(il) = wbmax*(0.01*(ph(il,1)-plfc(il)))**2
       ELSE ! Option provisoire ou le iflag_wb/10 est considere comme une vitesse

LMDZ6/trunk/libf/phylmd/cv3p2_closure.F90

@@ -609 +609 @@
     ELSE
       ! Calculate wbeff
-      IF (flag_wb==0) THEN
+      IF (NINT(flag_wb)==0) THEN
         wbeff(il) = wbmax
-      ELSE IF (flag_wb==1) THEN
+      ELSE IF (NINT(flag_wb)==1) THEN
         wbeff(il) = wbmax/(1.+500./(ph(il,1)-plfc(il)))
-      ELSE IF (flag_wb==2) THEN
+      ELSE IF (NINT(flag_wb)==2) THEN
         wbeff(il) = wbmax*(0.01*(ph(il,1)-plfc(il)))**2
       END IF

LMDZ6/trunk/libf/phylmd/cv3param.h

@@ -7 +7 @@
 !------------------------------------------------------------
 
-      integer flag_epKEorig, flag_wb
+      integer flag_epKEorig
+      real flag_wb
       integer cv_flag_feed
       integer noff, minorig, nl, nlp, nlm