Changeset 5763 for LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

Timestamp:

Jul 8, 2025, 11:09:14 AM (8 days ago)

Author:

jyg

Message:

Sequel of commit 5761: two versions of subroutine wake are implemented in the module
lmdz_wake.f90, one with the old computations for the vertical differential advection
(called "wake") and the other with the new computations introduced in commit 5761
(called wake2).

The switch between the two versions is controled by the flag iflag_wake.
iflag_wake is now made of two digits: the first one points to the version of wake

to be used (0 for wake or 2 for wake2); the second one defines the differential convective
heating.

+ iflag_wake= 01, 02, 03 ==> as before
+ iflag_wake= 21, 22, 23 ==> use wake2 instead of wake

File:

• LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -2527 +2527 @@
           ENDIF
           CALL histwrite_phy(o_tntc, zx_tmp_fi3d)
-       ELSE IF(iflag_thermals.GE.1.AND.iflag_wake.EQ.1) THEN
+       ELSE IF(iflag_thermals.GE.1.AND.iflag_wake.GE.1) THEN
           IF (vars_defined) THEN
              zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys + &
@@ -2547 +2547 @@
           IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
           CALL histwrite_phy(o_tnhusc, zx_tmp_fi3d)
-       ELSE IF (iflag_thermals.GE.1.AND.iflag_wake.EQ.1) THEN
+       ELSE IF (iflag_thermals.GE.1.AND.iflag_wake.GE.1) THEN
           IF (vars_defined) THEN
              zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys + &