Changeset 3208 for LMDZ6/trunk/libf/phylmd/add_wake_tend.F90

Timestamp:

Feb 16, 2018, 12:42:18 PM (6 years ago)

Author:

jyg

Message:

Implementation of a first crude model of the
dynamic of wake population.

File:

• LMDZ6/trunk/libf/phylmd/add_wake_tend.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/add_wake_tend.F90

@@ -1 @@
-SUBROUTINE add_wake_tend(zddeltat, zddeltaq, zds, zddens, zoccur, text, abortphy)
+SUBROUTINE add_wake_tend(zddeltat, zddeltaq, zds, zddensaw, zddensw, zoccur, text, abortphy)
 !===================================================================
 ! Ajoute les tendances liées aux diverses parametrisations physiques aux
@@ -9 +9 @@
 
 USE dimphy, ONLY: klon, klev
-USE phys_state_var_mod, ONLY: wake_deltat, wake_deltaq, wake_s, wake_dens
+USE phys_state_var_mod, ONLY: wake_deltat, wake_deltaq, wake_s,  &
+                              awake_dens, wake_dens
 
 USE print_control_mod, ONLY: prt_level
@@ -17 +18 @@
 !------------
   REAL, DIMENSION(klon, klev),   INTENT (IN)         :: zddeltat, zddeltaq
-  REAL, DIMENSION(klon),         INTENT (IN)         :: zds, zddens
+  REAL, DIMENSION(klon),         INTENT (IN)         :: zds, zddensaw, zddensw
   INTEGER, DIMENSION(klon),      INTENT (IN)         :: zoccur
   CHARACTER*(*),                 INTENT (IN)         :: text
@@ -53 +54 @@
          DO i = 1, klon
            IF (zoccur(i) .GE. 1) THEN
-             wake_s(i)    = wake_s(i)    + zds(i)
-             wake_dens(i) = wake_dens(i) + zddens(i)
+             wake_s(i)     = wake_s(i)    + zds(i)
+             awake_dens(i) = awake_dens(i) + zddensaw(i)
+             wake_dens(i)  = wake_dens(i) + zddensw(i)
            ELSE
-             wake_s(i)    = 0.
-             wake_dens(i) = 0.
+             wake_s(i)     = 0.
+             awake_dens(i) = 0.
+             wake_dens(i)  = 0.
            ENDIF   ! (zoccur(i) .GE. 1)
          END DO