Changeset 3000 for LMDZ5/trunk/libf/phylmd/calwake.F90

Timestamp:

Oct 2, 2017, 10:16:06 AM (7 years ago)

Author:

jyg

Message:

(i) Implementation of the conditionning of the Alp
provided by the wakes: when iflag_alp_wk_cond > 0,
Alp_wk is divided by the probability that there is
a gust front within the grid cell.

New subroutine alpale_wk.F90 called by physiq.

(ii) Some changes concerning the initialization of
the variable wake_dens.

File:

• LMDZ5/trunk/libf/phylmd/calwake.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/calwake.F90

@@ -72 +72 @@
   ! Variable internes
   ! -----------------
+  LOGICAL, SAVE                                      :: first = .TRUE.
+  !$OMP THREADPRIVATE(first)
   INTEGER                                            :: i, l
   INTEGER, DIMENSION(klon)                           :: znatsurf    ! 0 if pctsrf(is_oce)>0.1; 1 else.
@@ -177 +179 @@
 
   DO i = 1, klon
-    wdens(i) = wake_dens(i)
+    wdens(i) = max(0., wake_dens(i))
   END DO
 
@@ -298 +300 @@
       wake_dens(i) = wdens(i)
     END DO
-  ENDIF
+  ENDIF  ! (iflag_wake_tend .EQ. 0)
+!
+  IF (first) THEN
+    DO i = 1,klon
+      IF (wake_dens(i) < -1.) THEN
+        wake_dens(i) = wdens(i)
+      ENDIF
+    ENDDO
+    first=.false.
+  ENDIF  ! (first)
 !>jyg