Changeset 2318 for LMDZ5/trunk/libf/phylmd/physiq.F90

Timestamp:

Jun 26, 2015, 9:41:13 AM (9 years ago)

Author:

jyg

Message:

Bug fixing concerning the adiabatic adjustment of
the wake region (svn2309).

File:

• LMDZ5/trunk/libf/phylmd/physiq.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -2164 +2164 @@
      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)
+          IF (wake_s(i) .GT. 1.e-3) THEN
+            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)
+          ENDIF
         ENDDO
      ENDDO