Changeset 5994

Timestamp:

Jan 5, 2026, 5:08:57 PM (7 days ago)

Author:

evignon

Message:

debut harmonisation des routines lscp: update des commentaires dans lscp_main

File:

• LMDZ6/trunk/libf/phylmd/lmdz_lscp_main.f90 (modified) (12 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_lscp_main.f90

@@ -35 +35 @@
      dqsmelt, dqsfreez)
 
-!++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
-! Authors: Z.X. Li (LMD), J-L Dufresne (LMD), C. Rio (LMD), J-Y Grandpeix (LMD)
-!          A. JAM (LMD), J-B Madeleine (LMD), E. Vignon (LMD), L. Touzze-Peiffert (LMD)
+!++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+! Contacts: JB Madeleine jbm@lmd.ipsl.fr ;  E Vignon etienne.vignon@lmd.ipsl.fr
 !--------------------------------------------------------------------------------
-! Date: 01/2021
+! First version: 01/2021
 !--------------------------------------------------------------------------------
 ! Aim: Large Scale Clouds and Precipitation (LSCP)
@@ -47 +46 @@
 ! and 'ilp' (il pleut, L. Li)
 !
-! Compared to the original fisrtilp code, lscp
-! -> assumes thermcep = .TRUE. all the time (fisrtilp inconsistent when .FALSE.)
-! -> consider always precipitation thermalisation (fl_cor_ebil>0)
-! -> option iflag_fisrtilp_qsat<0 no longer possible (qsat does not evolve with T) 
-! -> option "oldbug" by JYG has been removed
-! -> iflag_t_glace >0 always
-! -> the 'all or nothing' cloud approach is no longer available (cpartiel=T always)
-! -> rectangular distribution from L. Li no longer available
-! -> We always account for the Wegener-Findeisen-Bergeron process (iflag_bergeron = 2 in fisrt) 
 !--------------------------------------------------------------------------------
 ! References: 
 !
 ! - Bony, S., & Emanuel, K. A. 2001, JAS, doi: 10.1175/1520-0469(2001)058<3158:APOTCA>2.0.CO;2
-! - Hourdin et al. 2013, Clim Dyn, doi:10.1007/s00382-012-1343-y
 ! - Jam et al. 2013, Boundary-Layer Meteorol, doi:10.1007/s10546-012-9789-3
-! - Jouhaud, et al. 2018. JAMES, doi:10.1029/2018MS001379
 ! - Madeleine et al. 2020, JAMES, doi:10.1029/2020MS002046
-! - Touzze-Peifert Ludo, PhD thesis, p117-124
+! - Borella et al. 2025, JAMES, doi:10.1029/2024MS004918
+! - Raillard et al. 2026, JAMES, doi:10.22541/essoar.175096287.71557703/v1 
+! - Vignon et al. 2026, ACP, doi:0.5194/egusphere-2025-4641 
 ! -------------------------------------------------------------------------------
 ! Code structure:
 !
-! P0>     Thermalization of the precipitation coming from the overlying layer
-! P1>     Evaporation of the precipitation (falling from the k+1 level)
-! P2>     Cloud formation (at the k level) 
-! P2.A.1> With the PDFs, calculation of cloud properties using the inital
-!         values of T and Q
-! P2.A.2> Coupling between condensed water and temperature
-! P2.A.3> Calculation of final quantities associated with cloud formation
-! P2.B>   Release of Latent heat after cloud formation
-! P3>     Autoconversion to precipitation (k-level)
-! P4>     Wet scavenging
-!------------------------------------------------------------------------------
-! Some preliminary comments (JBM) :
+! the formation of cloud and precipitation is treated as a top-down process,
+! starting from the top of the atmosphere down to the ground
 !
-! The cloud water that the radiation scheme sees is not the same that the cloud
-! water used in the physics and the dynamics 
-! 
-! During the autoconversion to precipitation (P3 step), radocond (cloud water used
-! by the radiation scheme) is calculated as an average of the water that remains
-! in the cloud during the precipitation and not the water remaining at the end
-! of the time step. The latter is used in the rest of the physics and advected
-! by the dynamics. 
-!
-! In summary:
-!
-! Radiation:
-! xflwc(newmicro)+xfiwc(newmicro) =
-!   radocond=lwcon(nc)+iwcon(nc)
-!
-! Notetheless, be aware of:
-!
-! radocond .NE. ocond(nc)
-! i.e.:
-! lwcon(nc)+iwcon(nc) .NE. ocond(nc)
-! but oliq+(ocond-oliq) .EQ. ocond
-! (which is not trivial)
+! P1>     Precipitation processes before cloud formation
+! P2>     Cloud formation and cloud phase determination  
+! P3>     Precipitation formation and interaction with clouds 
+! P4>     Computation of cloud water seen by radiation
+! P5>     Wet scavenging
+! P6>     Diagnostics and outputs
 !++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 !
@@ -290 +254 @@
   REAL, DIMENSION(klon) :: zdqsdT_raw
   REAL, DIMENSION(klon) :: gammasat,dgammasatdt                   ! coefficient to make cold condensation at the correct RH and derivative wrt T
-  REAL, DIMENSION(klon) :: Tbef,Tbefth,Tbefthm1,qlibef,DT                  ! temperature, humidity and temp. variation during condensation iteration
+  REAL, DIMENSION(klon) :: Tbef,Tbefth,Tbefthm1,qlibef,DT         ! temperature, humidity and temp. variation during condensation iteration
   REAL :: num,denom
   REAL :: cste
@@ -513 +477 @@
    
     !================================================================
-    ! Flag for the new and more microphysical treatment of precipitation from Atelier Nuage (R)
+    ! Flag for the new and more microphysical treatment of precipitation (poprecip) from Atelier Nuage (R)
     IF ( ok_poprecip ) THEN
 
@@ -565 +529 @@
     ! -------------------------------------------------------------------
 
-        ! P2.1> With the PDFs (log-normal, bigaussian) 
-        ! cloud properties calculation with the initial values of t and q
-        ! ----------------------------------------------------------------
-
         ! initialise gammasat and stratiform_or_all_distrib
         stratiform_or_all_distrib(:)=.TRUE.
         gammasat(:)=1.
 
-        ! part of the code that is supposed to become obsolete soon
-        !++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+        !!!! part of the code that is obsolete and should be deleted soon !!!
+        !++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
         IF (.NOT. ok_lscp_mergecond) THEN
           IF (iflag_cld_th.GE.5) THEN
@@ -684 +644 @@
                         ! we calculate a new iteration
                         keepgoing(i) = .TRUE.
-
-                        ! P2.2.1> cloud fraction and condensed water mass calculation
-                        ! Calculated variables: 
-                        ! rneb : cloud fraction
-                        ! zqn : total water within the cloud
-                        ! zcond: mean condensed water within the mesh
-                        ! rhcl: clear-sky relative humidity
-                        !---------------------------------------------------------------
 
                         ! new temperature that only serves in the iteration process:
@@ -759 +711 @@
                     ENDDO
 
-                    !---------------------------------------------
-                    !--   CONDENSATION AND ICE SUPERSATURATION  --
-                    !---------------------------------------------
+                    ! P2.1> CONDENSATION AND ICE SUPERSATURATION 
+                    !-------------------------------------------
 
                     CALL condensation_ice_supersat( &
@@ -802 +753 @@
 
 
-                  ! Cloud phase determination
-
+                  ! P2.2> Cloud phase determination
+                  !--------------------------------
 
                   IF (iflag_icefrac .LE. 1) THEN
@@ -891 +842 @@
                       IF (keepgoing(i)) THEN
 
-                       ! P2.2.2> Approximative calculation of temperature variation DT 
-                       !           due to condensation.
-                       ! Calculated variables: 
-                       ! dT : temperature change due to condensation
-                       !---------------------------------------------------------------
+                       ! Approximative calculation of temperature variation DT due to condensation.
+                       !----------------------------------------------------------------------------
 
                 
@@ -932 +880 @@
         ENDDO         ! iter=1,iflag_fisrtilp_qsat+1
 
-        ! P2.2> Final quantities calculation 
+        ! P2.3> Final quantities calculation 
         ! Calculated variables: 
         !   rneb : cloud fraction
         !   zcond: mean condensed water in the mesh
         !   zqn  : mean water vapor in the mesh
-        !   zfice: ice fraction in clouds
         !   zt   : temperature
         !   rhcl : clear-sky relative humidity
@@ -1014 +961 @@
 
     !================================================================
-    ! Flag for the new and more microphysical treatment of precipitation from Atelier Nuage (R)
+    ! Flag for the new and more microphysical treatment of precipitation (poprecip) from Atelier Nuage (R)
     IF (ok_poprecip) THEN
 
@@ -1064 +1011 @@
         ENDIF
       ELSE
+        ! During the autoconversion to precipitation with the historical precipitation scheme (P3 step),
+        ! radocond (cloud water used by the radiation scheme) is calculated as an average 
+        ! of the water that remains in the cloud during the precipitation 
+        ! and not the water remaining at the end of the time step
         radocond(i,k) = zradocond(i)
       ENDIF