Changeset 2309 for LMDZ5/trunk

Timestamp:

Jun 23, 2015, 6:08:54 PM (9 years ago)

Author:

jyg

Message:

Modification of physiq.F90: when wakes are
present, an adiabatic adjustment is applied to the
wake region before calling the Emanuel convection
scheme. The switch 'ok_adjwk', read from .def
files, controls this adjustment.

File:

• LMDZ5/trunk/libf/phylmd/physiq.F90 (modified) (4 diffs)

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -406 +406 @@
   REAL dt_a(klon,klev)
   REAL dq_a(klon,klev)
+  REAL d_t_adjwk(klon,klev)                !jyg
+  REAL d_q_adjwk(klon,klev)                !jyg
+  LOGICAL,SAVE :: ok_adjwk=.FALSE.
+  !$OMP THREADPRIVATE(ok_adjwk) 
   REAL, dimension(klon) :: www
   REAL, SAVE :: alp_offset
@@ -999 +1003 @@
      CALL suphel ! initialiser constantes et parametres phys.
      CALL getin_p('random_notrig_max',random_notrig_max)
+     CALL getin_p('ok_adjwk',ok_adjwk)
   ENDIF
 
@@ -2139 +2144 @@
         enddo
      enddo
+!
+!jyg<
+     ! Perform dry adiabatic adjustment on wake profile
+     ! The corresponding tendencies are added to the convective tendencies
+     ! after the call to the convective scheme.
+     IF (iflag_wake>=1) then
+      IF (ok_adjwk) THEN
+        limbas(:) = 1
+        CALL ajsec(paprs, pplay, t_wake, q_wake, limbas, &
+                  d_t_adjwk, d_q_adjwk)
+      ELSE
+      DO k=1,klev
+         DO i=1,klon
+           d_t_adjwk(i,k) = 0.
+           d_q_adjwk(i,k) = 0.
+         ENDDO
+      ENDDO
+     ENDIF
+!
+     DO k=1,klev
+        DO i=1,klon
+          t_wake(i,k) = t_wake(i,k) + d_t_adjwk(i,k)
+          q_wake(i,k) = q_wake(i,k) + d_q_adjwk(i,k)
+          wake_deltat(i,k) = wake_deltat(i,k) + d_t_adjwk(i,k)
+          wake_deltaq(i,k) = wake_deltaq(i,k) + d_q_adjwk(i,k)
+        ENDDO
+     ENDDO
+     ENDIF ! (iflag_wake>=1)
+!>jyg
+!
 
      ! Calcul de l'energie disponible ALE (J/kg) et de la puissance
@@ -2328 +2363 @@
            if (iflagctrl(i).le.1) itau_con(i)=itau_con(i)+1
         enddo
-
+!
+!jyg<
+!    Add the tendency due to the dry adjustment of the wake profile
+       DO k=1,klev
+          DO i=1,klon
+          ftd(i,k) = ftd(i,k) + wake_s(i)*d_t_adjwk(i,k)/dtime
+          fqd(i,k) = fqd(i,k) + wake_s(i)*d_q_adjwk(i,k)/dtime
+          d_t_con(i,k) = d_t_con(i,k) + wake_s(i)*d_t_adjwk(i,k)
+          d_q_con(i,k) = d_q_con(i,k) + wake_s(i)*d_q_adjwk(i,k)
+          enddo
+       enddo
+!>jyg
+!
      ELSE ! ok_cvl