Changeset 5384

Timestamp:

Dec 5, 2024, 12:25:15 PM (8 days ago)

Author:

evignon

Message:

ajout de diagnostiques des tendances de precip, Audran Borella

Location:

LMDZ6/trunk/libf/phylmd

Files:

• lmdz_lscp.f90 (modified) (4 diffs)
• phys_output_write_mod.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/lmdz_lscp.f90

@@ -670 +670 @@
                 ENDIF
 
+                !--Diagnostics
+                dqreva(i,k) = - zqev / dtime
+                dqssub(i,k) = - zqevti / dtime
 
                 ! New solid and liquid precipitation fluxes after evap and sublimation
@@ -676 +679 @@
                 zifln(i) = Max(0.,zifl(i) - zqevi*(paprs(i,k)-paprs(i,k+1)) &
                          /RG/dtime)
-
+                
 
                 ! vapor, temperature, precip fluxes update
@@ -1383 +1386 @@
                 ENDIF
 
+
                 IF (iflag_autoconversion .GE. 1) THEN
                    ! debugged version with phase conservation through the autoconversion process
@@ -1394 +1398 @@
                    zoliq(i)  = MAX(zoliq(i)-zprecip , 0.0)
                 ENDIF
+
+                !--Diagnostics
+                dqrauto(i,k) = dqrauto(i,k) + zrain / dtime
+                dqsauto(i,k) = dqsauto(i,k) + zsnow / dtime
+
 
                 ! c_iso: call isotope_conversion (for readibility) or duplicate 

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -2086 +2086 @@
            CALL histwrite_phy(o_dqsfreez, dqsfreez)
            CALL histwrite_phy(o_dqsrim, dqsrim)
+           ELSE
+            CALL histwrite_phy(o_dqreva, dqreva)
+            CALL histwrite_phy(o_dqssub, dqssub)
+            CALL histwrite_phy(o_dqrauto, dqrauto)
+            CALL histwrite_phy(o_dqsauto, dqsauto)
            ENDIF
        ENDIF