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

Timestamp:

Oct 28, 2023, 4:47:32 PM (8 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/physiq_mod.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

@@ -197 +197 @@
        d_q_dyn2d,d_ql_dyn2d,d_qs_dyn2d,d_qbs_dyn2d, &
        ! Physic tendencies
-       d_t_con,d_q_con,d_u_con,d_v_con, &
+       d_t_con,d_q_con,d_q_con_zmasse,d_u_con,d_v_con, &
        d_tr, &                              !! to be removed?? (jyg)
        d_t_wake,d_q_wake, &
@@ -1332 +1332 @@
     real xtql1(ntraciso,klon),xtql2(ntraciso,klon),corrxtql(ntraciso)
 #endif
-    
+
+    !--OB flag to activate better conservation of water tendency when convection is not called every timestep
+    LOGICAL, PARAMETER :: ok_conserv_d_q_con=.FALSE.
 
     REAL, dimension(klon,klev) :: t_env,q_env
@@ -4042 +4044 @@
     ENDIF
 
+    !--saving d_q_con * zmass for next timestep if convection is not called every timestep
+    IF (ok_conserv_d_q_con) THEN
+      d_q_con_zmasse(:,:) = d_q_con(:,:) * zmasse(:,:)
+    ENDIF
+
     !     CALL homogene(paprs, q_seri, d_q_con, u_seri,v_seri,
     !    .              d_u_con, d_v_con)
@@ -4071 +4078 @@
 !!
 !!
+
+    !--recompute d_q_con with zmasse from new timestep
+    IF (ok_conserv_d_q_con) THEN
+      d_q_con(:,:)=d_q_con_zmasse(:,:)/zmasse(:,:)
+    ENDIF
+
     CALL add_phys_tend(d_u_con, d_v_con, d_t_con, d_q_con, dql0, dqi0, dqbs0, paprs, &
          'convection',abortphy,flag_inhib_tend,itap,0 &