Changeset 3279 for trunk

Timestamp:

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

Author:

jleconte

Message:

flag moist_convection_inhibition (Noe C)

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• inifis_mod.F90 (modified) (1 diff)
• moistadj_generic.F90 (modified) (2 diffs)

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

@@ -1029 +1029 @@
      call getin_p("moistadjustment_generic",moistadjustment_generic)
      if (is_master) write(*,*)trim(rname)//": moistadjustment_generic = ", moistadjustment_generic
+
+     if (is_master) write(*,*)trim(rname)//": Moist convection inhibition for GCS ?"
+     moist_convection_inhibition=.false. ! default value
+     call getin_p("moist_convection_inhibition",moist_convection_inhibition)
+     if (is_master) write(*,*)trim(rname)//": moist_convection_inhibition = ", moist_convection_inhibition
      
      if (is_master) write(*,*)trim(rname)//": Virtual correction ?"

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

@@ -6 +6 @@
    use ioipsl_getin_p_mod, only: getin_p !-> to get the metallicity
    use comcstfi_mod, only: r, cpp, mugaz
+   use callkeys_mod, only: moist_convection_inhibition
 
    implicit none
@@ -182 +183 @@
       
       ! calculate moist convection inhibition criterion
-      IF (epsi_generic .gt. 1) THEN ! GCS molecular weight is heavier than dry gas: 
+      IF (moist_convection_inhibition .and. (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)
+         write(*,*) 'inhibition of moist convection if vapor amount exceeds a critical mixing ratio (see Leconte et al. 2017)'
          DO k = 1, nlayer
             DO i = 1, ngrid