Changeset 5609 for LMDZ6/branches/contrails/libf/misc

Timestamp:

Apr 13, 2025, 7:10:19 PM (4 months ago)

Author:

aborella

Message:

• changed treatment of prognostic variables for prognostic clouds
• adapted sedimentation and autoconversion for prognostic cirrus clouds
• cloud mixing, ice sedimentation and ISSR diagnosis are now consistent with the water vapor PDF
• simplified assumptions for ice crystals deposition / sublimation
• first version of the coupling between prognostic cirrus clouds and deep convection
• added persistent contrail cirrus clouds in radiative diagnostics

File:

• LMDZ6/branches/contrails/libf/misc/readTracFiles_mod.f90 (modified) (2 diffs)

LMDZ6/branches/contrails/libf/misc/readTracFiles_mod.f90

@@ -63 +63 @@
     LOGICAL               :: isAdvected  = .FALSE.              !--- "true" tracers: iadv > 0.   COUNT(isAdvected )=nqtrue
     LOGICAL               :: isInPhysics = .TRUE.               !--- "true" tracers: in tr_seri. COUNT(isInPhysics)=nqtottr
-    INTEGER               :: iso_iGroup  = 0                    !--- Isotopes group index in isotopes(:)
+    INTEGER               :: iso_iGroup  = -1                   !--- Isotopes group index in isotopes(:)
     INTEGER               :: iso_iName   = 0                    !--- Isotope  name  index in isotopes(iso_iGroup)%trac(:)
     INTEGER               :: iso_iZone   = 0                    !--- Isotope  zone  index in isotopes(iso_iGroup)%zone(:)
@@ -122 +122 @@
   !--- SOME PARAMETERS THAT ARE NOT LIKELY TO CHANGE OFTEN
   CHARACTER(LEN=maxlen), SAVE      :: tran0        = 'air'      !--- Default transporting fluid
-  CHARACTER(LEN=maxlen), PARAMETER :: old_phases   = 'vlibfcaq'    !--- Old phases for water (no separator)
-  CHARACTER(LEN=maxlen), PARAMETER :: known_phases = 'glsbfcaq'    !--- Known phases initials
+  CHARACTER(LEN=maxlen), PARAMETER :: old_phases   = 'vlib'     !--- Old phases for water (no separator)
+  CHARACTER(LEN=maxlen), PARAMETER :: known_phases = 'glsb'     !--- Known phases initials
   INTEGER, PARAMETER :: nphases = LEN_TRIM(known_phases)        !--- Number of phases
   CHARACTER(LEN=maxlen), SAVE      :: phases_names(nphases) &   !--- Known phases names
-                                = ['gaseous  ', 'liquid   ', 'solid    ','blownSnow', 'fracCloud', 'cldVapor ', 'aviContra', 'qtContra ']
+                                = ['gaseous  ', 'liquid   ', 'solid    ','blownSnow']
   CHARACTER(LEN=1),      SAVE :: phases_sep  =  '_'             !--- Phase separator
   CHARACTER(LEN=maxlen), SAVE :: isoFile = 'isotopes_params.def'!--- Name of the isotopes parameters file