Changeset 3276 for trunk

Timestamp:

Mar 20, 2024, 4:29:20 PM (8 months ago)

Author:

jleconte

Message:

moist_adj_generic cleaned, completed and commented

File:

• trunk/LMDZ.GENERIC/libf/phystd/moistadj_generic.F90 (modified) (14 diffs)

trunk/LMDZ.GENERIC/libf/phystd/moistadj_generic.F90

@@ -1 +1 @@
 subroutine moistadj_generic(ngrid, nlayer, nq, pt, pq, pdq, pplev, pplay, pdtmana, pdqmana, ptimestep, rneb)
 
-   !use watercommon_h, only: T_h2O_ice_liq, RLVTT, RCPD, RCPV, Psat_water, Lcpdqsat_water
-   !USE tracer_h, only: igcm_h2o_vap, igcm_h2o_ice
    use generic_cloud_common_h
    use generic_tracer_index_mod, only: generic_tracer_index
@@ -21 +19 @@
 !     -------
 !     Adapted from the moistadj.F90 routine
-!     for generic tracers (condensable species)
+!     for generic condensable species (GCS) tracers
 !     by Noe CLEMENT (2023)
 !     
@@ -75 +73 @@
       DOUBLE PRECISION :: zdp, zdpm
 
-      real q_cri(ngrid,nlayer)
+      real q_cri(ngrid,nlayer) ! moist convection inhibition criterion
 
       REAL zsat ! super-saturation
@@ -99 +97 @@
 
          RCPD = cpp
-         ! RCPV   = 1.88e3 ! specific heat capacity of water vapor at 350K ! MUST BE CHANGED !!!
-         ! 1.709e3 for methane - should be added in table_tracers_condensable (datagcm)
 
          write(*,*) "value for metallicity? "
@@ -116 +112 @@
                RLVTT_generic = constants_RLVTT_generic(iq)
                RCPV_generic = constants_RCPV_generic(iq)
+               Tref = constants_Tref(iq)
+
+               write(*,*) noms(igcm_generic_vap),", q_cri at ", Tref, "K (in kg/kg): ", ( 1 / (1 - 1/epsi_generic)) * (r * mugaz/1000.) / delta_vapH * Tref
+
             endif
          enddo
@@ -158 +158 @@
             call Lcpdqsat_generic(v_t,v_p,zpsat(i,k),zqs(i,k),zdqs(i,k),zdlnpsat(i,k))
 
-            !call Psat_water(v_t,v_p,zpsat(i,k),zqs(i,k))
-            !call Lcpdqsat_water(v_t,v_p,zpsat(i,k),zqs(i,k),zdqs(i,k),zdlnpsat(i,k))
          ENDDO
       ENDDO
@@ -179 +177 @@
             gamcpdz(i,k) = ( (R/RCPD*v_cptt2*(1.- v_zqs) + RLVTT_generic*v_zqs) * (pplay(i,k-1)-pplay(i,k))/pplev(i,k) )  &
                / (((1.- v_zqs) + v_zqs * RCPV_generic/RCPD)*v_pratio + v_zqs  * v_dlnpsat)                
+            ! Note that gamcpdz is defined as positive, so -gamcpdz is the real moist-adiabatic gradient [dT/dz]_ad
          ENDDO
       ENDDO
       
-      !     calculates moist convection inhibition criterion
-      DO k = 1, nlayer
-         DO i = 1, ngrid
-            q_cri(i,k) = ( 1 / (1 - 1/epsi_generic)) * r * mugaz/1000. / delta_vapH * zt(i,k)
-         ENDDO
-      ENDDO
-
-      write(*,*) "q_cri at 80K : ", ( 1 / (1 - 1/epsi_generic)) * (r * mugaz/1000.) / delta_vapH * 80
-
-!------------------------------------ modification of unstable profile
+      ! calculate moist convection inhibition criterion
+      IF (epsi_generic .gt. 1) THEN ! GCS molecular weight is heavier than dry gas: 
+         ! inhibition of moist convection if vapor amount exceeds q_cri (Eq. 17 of Leconte et al. 2017)
+         DO k = 1, nlayer
+            DO i = 1, ngrid
+               q_cri(i,k) = ( 1 / (1 - 1/epsi_generic)) * r * mugaz/1000. / delta_vapH * zt(i,k)
+            ENDDO
+         ENDDO
+      ELSE ! GCS molecular weight is lighter than dry gas
+         DO k = 1, nlayer
+            DO i = 1, ngrid
+               q_cri(i,k) = 2. ! vapor amount will never exceed q_cri, q_cri call becomes transparent in the next lines
+            ENDDO
+         ENDDO
+      ENDIF
+
+
+!------------------------------------ defining the bottom (k1) and the top (k2) of the moist-convective column
       DO 9999 i = 1, ngrid
 
@@ -209 +216 @@
          +(local_q(i,k2)-zqs(i,k2))*(pplev(i,k2)-pplev(i,k2+1))
 
+      ! if: (gradient is not steeper than moist-adiabat) or (level is not saturated) or (moist convection is inhibited because criterion is satisfied)
+      ! then: no moist convection (GO TO 810)
       IF ( zflo.LE.0.0 .OR. zsat.LE.0.0 .OR. local_q(i,k2-1).GT.q_cri(i,k2-1)) GOTO 810
       k1 = k2 - 1
@@ -224 +233 @@
    821 CONTINUE
 
-!------------------------------------------------------ local adjustment
+!------------------------------------------------------ local adjustment of the moist-convective column between k1 and k2
    830 CONTINUE ! actual adjustment
    Do nn=1,niter
@@ -254 +263 @@
          call Psat_generic(v_t,v_p,metallicity,zpsat(i,k),zqs(i,k))
          call Lcpdqsat_generic(v_t,v_p,zpsat(i,k),zqs(i,k),zdqs(i,k),zdlnpsat(i,k))
-         
-         !call Psat_water(v_t,v_p,zpsat(i,k),zqs(i,k))
-         !call Lcpdqsat_water(v_t,v_p,zpsat(i,k),zqs(i,k),zdqs(i,k),zdlnpsat(i,k))
 
       ENDDO
@@ -293 +299 @@
 !     Test to see if we've reached the bottom
 
-      IF (k1 .EQ. 1) GOTO 841 ! yes we have!
+      IF (k1 .EQ. 1) GOTO 841 ! yes we have! the bottom of moist-convective column is the bottom of the model
+
       zflo = v_cptt(i,k1-1) - v_cptt(i,k1) - gamcpdz(i,k1)
       zsat=(local_q(i,k1-1)-zqs(i,k1-1))*(pplev(i,k1-1)-pplev(i,k1))   &
          + (local_q(i,k1)-zqs(i,k1))*(pplev(i,k1)-pplev(i,k1+1))
-      IF (zflo.LE.0.0 .OR. zsat.LE.0.0) GOTO 841 ! yes we have!
+
+      IF (zflo.LE.0.0 .OR. zsat.LE.0.0 .OR. local_q(i,k1-1).GT.q_cri(i,k1-1)) GOTO 841 ! yes we have!
+      ! the bottom of the moist-convective column is no more convective: GOTO 841
 
    840 CONTINUE
       k1 = k1 - 1
       IF (k1 .EQ. 1) GOTO 830 ! GOTO 820 (a tester, Z.X.Li, mars 1995)
+      ! bottom of moist-convective column is bottom of the model but is still convective: GOTO 830
       zsat = zsat + (local_q(i,k1-1)-zqs(i,k1-1))               &
                   *(pplev(i,k1-1)-pplev(i,k1))
       zflo = v_cptt(i,k1-1) - v_cptt(i,k1) - gamcpdz(i,k1)
-      IF (zflo.GT.0.0 .AND. zsat.GT.0.0) THEN
+      IF (zflo.GT.0.0 .AND. zsat.GT.0.0 .AND. local_q(i,k1-1).LT.q_cri(i,k1-1)) THEN
+         ! the bottom of the moist-convective column is still convective
+         ! we continue to search for the non-convective bottom
          GOTO 840
       ELSE
+         ! the bottom of the moist-convective column is no more convective now,
+         ! but, since the bottom of moist-convective column has changed, we must do again the adjustment
          GOTO 830 ! GOTO 820 (a tester, Z.X.Li, mars 1995)
       ENDIF
@@ -313 +327 @@
 
       GOTO 810 ! look for other layers higher up
+      ! because there could be other moist-convective columns higher up (but separated from the one we have just calculated before)
 
 9999 CONTINUE ! loop over all the points
 
 !-----------------------------------------------------------------------
-! Determine the cloud fraction (hypothese: la nebulosite a lieu
-! a l'endroit ou la vapeur d'eau est diminuee par l'ajustement):
+! Determine the cloud fraction (hypothesis: nebulosity occurs
+! where GCS vapor is reduced by adjustment):
 
       DO k = 1, nlayer
@@ -329 +344 @@
       ENDDO
 
-! Distribuer l'eau condensee en eau liquide nuageuse (hypothese:
-! l'eau liquide est distribuee aux endroits ou la vapeur d'eau
-! diminue et d'une maniere proportionnelle a cet diminution):
+! Distribute GCS condensates into cloudy liquid/solid condensates (hypothesis:
+! liquid/solid condensates are distributed to areas where GCS vapor
+! decreases and are distributed in proportion to this decrease):
       DO i = 1, ngrid
          IF (itest(i)) THEN
@@ -368 +383 @@
       ENDDO
 
-!     now subroutine -----> GCM variables
+      ! now subroutine -----> GCM variables
       DO k = 1, nlayer
          DO i = 1, ngrid