Changeset 4952

Timestamp:

May 24, 2024, 6:48:07 PM (8 months ago)

Author:

oboucher

Message:

flag for mass fixer is now ok_water_mass_fixer and can be changed in namelist

Location:

LMDZ6/trunk/libf/phylmd

Files:

• clesphys.h (modified) (1 diff)
• physiq_mod.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/clesphys.h

@@ -110 +110 @@
        LOGICAL :: ok_3Deffect
 
+!OB flag to activate water mass fixer in physiq
+       LOGICAL :: ok_water_mass_fixer
+
        COMMON/clesphys/                                                 &
 ! REAL FIRST

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -1265 +1265 @@
 
     !--OB variables for mass fixer (hard coded for now)
-    LOGICAL, PARAMETER :: mass_fixer=.FALSE.
     REAL qql1(klon),qql2(klon),corrqql
 
@@ -1509 +1508 @@
        iflag_phytrac = 1 ! by default we do want to call phytrac
        CALL getin_p('iflag_phytrac',iflag_phytrac)
+
+       ok_water_mass_fixer=.FALSE.  ! OB: by default we do not apply the mass fixer
+       CALL getin_p('ok_water_mass_fixer',ok_water_mass_fixer)
 #ifdef CPP_Dust
        IF (iflag_phytrac.EQ.0) THEN 
@@ -1533 +1535 @@
        WRITE(lunout,*) 'fl_cor_ebil=',        fl_cor_ebil
        WRITE(lunout,*) 'iflag_phytrac=',      iflag_phytrac
+       WRITE(lunout,*) 'ok_water_mass_fixer=',ok_water_mass_fixer
        WRITE(lunout,*) 'NVM=',                nvm_lmdz
 
@@ -2458 +2461 @@
     ENDDO
     !
-    !--OB mass fixer 
-    IF (mass_fixer) THEN
+    !--OB water mass fixer 
+    IF (ok_water_mass_fixer) THEN
     !--store initial water burden
     qql1(:)=0.0
@@ -5475 +5478 @@
     !--currently flag is turned off
     !==================================================================
-    IF (mass_fixer) THEN
+    IF (ok_water_mass_fixer) THEN
     qql2(:)=0.0
     DO k = 1, klev