Changeset 4819 for LMDZ6/trunk/libf/phylmd/lmdz_lscp.F90

Timestamp:

Feb 14, 2024, 8:55:10 PM (4 months ago)

Author:

evignon

Message:

modifications du commit precedent a la suite de l'atelier nuages

File:

• LMDZ6/trunk/libf/phylmd/lmdz_lscp.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_lscp.F90

@@ -19 +19 @@
      Tcontr, qcontr, qcontr2, fcontrN, fcontrP,         &
      cloudth_sth,cloudth_senv,cloudth_sigmath,cloudth_sigmaenv, &
-     dqreva,dqssub,dqrauto,dqrcol,dqrmelt,dqrfreez,dqsauto, &
-     dqsagg,dqsrim,dqsmelt,dqsfreez)
+     qrain, qsnow, dqreva, dqssub, dqrauto, dqrcol,     &
+     dqrmelt, dqrfreez, dqsauto, dqsagg, dqsrim,        &
+     dqsmelt, dqsfreez)
 
 !++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
@@ -203 +204 @@
   ! for POPRECIP
 
-  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: dqreva         !-- rain tendendy due to evaporation [kg/kg/s]
-  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: dqssub         !-- snow tendency due to sublimation [kg/kg/s]
-  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: dqrcol         !-- rain tendendy due to collection by rain of liquid cloud droplets [kg/kg/s]
-  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: dqsagg         !-- snow tendency due to collection of lcoud ice by aggregation [kg/kg/s]
-  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: dqrauto        !-- rain tendency due to autoconversion of cloud liquid [kg/kg/s]
-  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: dqsauto        !-- snow tendency due to autoconversion of cloud ice [kg/kg/s]
-  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: dqsrim         !-- snow tendency due to riming [kg/kg/s]
-  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: dqsmelt        !-- snow tendency due to melting [kg/kg/s]
-  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: dqrmelt        !-- rain tendency due to melting [kg/kg/s]
-  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: dqsfreez       !-- snow tendency due to freezing [kg/kg/s]
-  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: dqrfreez       !-- rain tendency due to freezing [kg/kg/s]
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: qrain          !--specific rain content [kg/kg]
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: qsnow          !--specific snow content [kg/kg]
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: dqreva         !--rain tendendy due to evaporation [kg/kg/s]
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: dqssub         !--snow tendency due to sublimation [kg/kg/s]
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: dqrcol         !--rain tendendy due to collection by rain of liquid cloud droplets [kg/kg/s]
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: dqsagg         !--snow tendency due to collection of lcoud ice by aggregation [kg/kg/s]
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: dqrauto        !--rain tendency due to autoconversion of cloud liquid [kg/kg/s]
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: dqsauto        !--snow tendency due to autoconversion of cloud ice [kg/kg/s]
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: dqsrim         !--snow tendency due to riming [kg/kg/s]
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: dqsmelt        !--snow tendency due to melting [kg/kg/s]
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: dqrmelt        !--rain tendency due to melting [kg/kg/s]
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: dqsfreez       !--snow tendency due to freezing [kg/kg/s]
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: dqrfreez       !--rain tendency due to freezing [kg/kg/s]
 
 
@@ -371 +374 @@
 temp_cltop(:,:)=0.
 !-- poprecip
-dqreva(:,:)=0.0
-dqrauto(:,:)=0.0
-dqrmelt(:,:)=0.0
-dqrfreez(:,:)=0.0
-dqrcol(:,:)=0.0
-dqssub(:,:)=0.0
-dqsauto(:,:)=0.0
-dqsrim(:,:)=0.0
-dqsagg(:,:)=0.0
-dqsfreez(:,:)=0.0
-dqsmelt(:,:)=0.0
+qrain(:,:)    = 0.
+qsnow(:,:)    = 0.
+dqreva(:,:)   = 0.
+dqrauto(:,:)  = 0.
+dqrmelt(:,:)  = 0.
+dqrfreez(:,:) = 0.
+dqrcol(:,:)   = 0.
+dqssub(:,:)   = 0.
+dqsauto(:,:)  = 0.
+dqsrim(:,:)   = 0.
+dqsagg(:,:)   = 0.
+dqsfreez(:,:) = 0.
+dqsmelt(:,:)  = 0.
 
 
@@ -992 +997 @@
                             zrfl, zrflclr, zrflcld, &
                             zifl, ziflclr, ziflcld, &
-                            dqrauto(:,k),dqrcol(:,k),dqrmelt(:,k),dqrfreez(:,k), &
-                            dqsauto(:,k),dqsagg(:,k),dqsrim(:,k),dqsmelt(:,k),dqsfreez(:,k) &
+                            qrain(:,k), qsnow(:,k), dqrauto(:,k), &
+                            dqrcol(:,k), dqrmelt(:,k), dqrfreez(:,k), &
+                            dqsauto(:,k), dqsagg(:,k), dqsrim(:,k), &
+                            dqsmelt(:,k), dqsfreez(:,k) &
                             )
 
@@ -1003 +1010 @@
                 zfice(i) = 0.0
         ENDIF
-        ! when poprecip activated, radiation does not see any precipitation content
-        radocond(i,k) = zoliq(i)
-        radocondl(i,k)= radocond(i,k)*(1.-zfice(i))
-        radocondi(i,k)= radocond(i,k)*zfice(i)
+
+        ! calculation of specific content of condensates seen by radiative scheme 
+        IF (ok_radocond_snow) THEN
+           radocond(i,k) = zoliq(i)
+           radocondl(i,k)= radocond(i,k)*(1.-zfice(i))
+           radocondi(i,k)= radocond(i,k)*zfice(i)+qsnow(i,k)
+        ELSE
+           radocond(i,k) = zoliq(i)
+           radocondl(i,k)= radocond(i,k)*(1.-zfice(i))
+           radocondi(i,k)= radocond(i,k)*zfice(i)
+        ENDIF
       ENDDO