Changeset 4742 for LMDZ6

Timestamp:

Oct 28, 2023, 4:47:32 PM (14 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.

Location:

LMDZ6/trunk/libf

Files:

• phylmd/phys_local_var_mod.F90 (modified) (4 diffs)
• phylmd/physiq_mod.F90 (modified) (4 diffs)
• phylmdiso/phys_local_var_mod.F90 (modified) (3 diffs)
• phylmdiso/physiq_mod.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90

@@ -38 +38 @@
       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)
@@ -640 +640 @@
       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))
@@ -767 +767 @@
       ALLOCATE(toplwad0_aero_s2(klon), sollwad0_aero_s2(klon))
 
-
 ! FH Ajout de celles necessaires au phys_output_write_mod
 
@@ -987 +986 @@
       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)

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -162 +162 @@
        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, &
@@ -1256 +1256 @@
     LOGICAL, PARAMETER :: mass_fixer=.FALSE.
     REAL qql1(klon),qql2(klon),corrqql
+
+    !--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
@@ -3377 +3380 @@
     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)
@@ -3406 +3414 @@
 !!
 !!
+
+    !--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)

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)

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 &