Changeset 4742 for LMDZ6/trunk/libf/phylmdiso/phys_local_var_mod.F90

Timestamp:

Oct 28, 2023, 4:47:32 PM (9 months ago)

Author:

oboucher

Message:

Introduction of a flag ok_conserv_d_q_con to conserve water when convection is called every second timestep.
The idea is to save d_q_con_zmasse from a timestep to the next rather than d_q_con. This is only as good as concvl for conserving water.

File:

• LMDZ6/trunk/libf/phylmdiso/phys_local_var_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmdiso/phys_local_var_mod.F90

@@ -36 +36 @@
       REAL, SAVE, ALLOCATABLE :: d_tr_dyn(:,:,:)
       !$OMP THREADPRIVATE(d_tr_dyn)
-      REAL, SAVE, ALLOCATABLE :: d_t_con(:,:),d_q_con(:,:)
-      !$OMP THREADPRIVATE(d_t_con,d_q_con)
+      REAL, SAVE, ALLOCATABLE :: d_t_con(:,:),d_q_con(:,:),d_q_con_zmasse(:,:)
+      !$OMP THREADPRIVATE(d_t_con,d_q_con,d_q_con_zmasse)
       REAL, SAVE, ALLOCATABLE :: d_u_con(:,:),d_v_con(:,:)
       !$OMP THREADPRIVATE(d_u_con,d_v_con)
@@ -737 +737 @@
       ALLOCATE(d_u_dyn(klon,klev),d_v_dyn(klon,klev))
       ALLOCATE(d_tr_dyn(klon,klev,nbtr))                   !RomP
-      ALLOCATE(d_t_con(klon,klev),d_q_con(klon,klev))
+      ALLOCATE(d_t_con(klon,klev),d_q_con(klon,klev),d_q_con_zmasse(klon,klev))
       ALLOCATE(d_u_con(klon,klev),d_v_con(klon,klev))
       ALLOCATE(d_t_wake(klon,klev),d_q_wake(klon,klev))
@@ -1139 +1139 @@
       DEALLOCATE(d_u_dyn,d_v_dyn)
       DEALLOCATE(d_tr_dyn)                      !RomP
-      DEALLOCATE(d_t_con,d_q_con)
+      DEALLOCATE(d_t_con,d_q_con,d_q_con_zmasse)
       DEALLOCATE(d_u_con,d_v_con)
       DEALLOCATE(d_t_wake,d_q_wake)