Changeset 5791 for LMDZ6/branches/contrails/libf/dyn3d_common/infotrac.f90

Timestamp:

Jul 28, 2025, 7:23:15 PM (6 days ago)

Author:

aborella

Message:

Merge with trunk r5789

Location:

LMDZ6/branches/contrails

Files:

• . (modified) (1 prop)
• libf/dyn3d_common/infotrac.f90 (modified) (13 diffs)

LMDZ6/branches/contrails/libf/dyn3d_common/infotrac.f90

@@ -3 +3 @@
 MODULE infotrac
 
-   USE       strings_mod, ONLY: msg, maxlen, cat, dispTable, int2str, bool2str, strStack, strParse, strCount, strIdx
-   USE readTracFiles_mod, ONLY: trac_type, readTracersFiles, tracers,  addPhase,  addKey, iH2O,  &
-       isoSelect,  indexUpdate, isot_type, testTracersFiles, isotope,  delPhase,  getKey, tran0, &
-       isoKeys, isoName, isoZone, isoPhas, processIsotopes,  isoCheck, itZonIso,  nbIso,         &
-          niso,   ntiso,   nzone,   nphas,   maxTableWidth,  iqIsoPha, iqWIsoPha, ixIso, new2oldH2O, newHNO3, oldHNO3
+   USE       strings_mod, ONLY: msg, maxlen, cat, dispTable, num2str, strStack, strParse, strCount, strIdx, maxTableWidth
+   USE readTracFiles_mod, ONLY: trac_type, readTracersFiles, tracers,  addPhase,  addKey, iH2O,    &
+       isoSelect,  indexUpdate, isot_type, testTracersFiles, isotope,  delPhase,  getKey, tran0,   &
+       isoKeys, isoName, isoZone, isoPhas, processIsotopes,  isoCheck, itZonIso,  nbIso, newHNO3,  &
+          niso,   ntiso,   nzone,   nphas, isoF=>isoFamilies,iqIsoPha, iqWIsoPha, ixIso, oldHNO3, new2oldH2O
    IMPLICIT NONE
 
@@ -22 +22 @@
    !=== FOR ISOTOPES: General
    PUBLIC :: isot_type, nbIso                              !--- Derived type, full isotopes families database + nb of families
-   PUBLIC :: isoSelect, ixIso                              !--- Isotopes family selection tool + selected family index
+   PUBLIC :: isoSelect, ixIso, isoFamilies                 !--- Isotopes families selection tool + selected family index
    !=== FOR ISOTOPES: Specific to water
    PUBLIC :: iH2O                                          !--- Value of "ixIso" for "H2O" isotopes class
@@ -81 +81 @@
 !  | iso_iName   | Isotope  name  index in isotopes(iso_iGroup)%trac(:) | iso_indnum  | 1:niso                 |
 !  | iso_iZone   | Isotope  zone  index in isotopes(iso_iGroup)%zone(:) | zone_num    | 1:nzone                |
-!  | iso_iPhas   | Isotope  phase index in isotopes(iso_iGroup)%phas(:) | phase_num   | 1:nphas                |
+!  | iso_iPhase  | Isotope  phase index in isotopes(iso_iGroup)%phas(:) | phase_num   | 1:nphas                |
 !  +-------------+------------------------------------------------------+-------------+------------------------+
 !
@@ -98 +98 @@
 !  +-----------------+--------------------------------------------------+--------------------+-----------------+
 
+   !=== THRESHOLDS FOR WATER
    REAL, PARAMETER :: min_qParent = 1.e-30, min_qMass = 1.e-18, min_ratio = 1.e-16 ! MVals et CRisi
 
-   !=== DIMENSIONS OF THE TRACERS TABLES AND OTHER SCALAR VARIABLES
+   !=== DIMENSIONS OF THE TRACERS TABLES, TRACERS TYPE(S)
    INTEGER, SAVE :: nqtot                                       !--- Tracers nb in dynamics (incl. higher moments + H2O)
    INTEGER, SAVE :: nbtr                                        !--- Tracers nb in physics  (excl. higher moments + H2O)
@@ -111 +112 @@
    INTEGER, SAVE, ALLOCATABLE :: conv_flg(:), pbl_flg(:)        !--- Convection / boundary layer activation (nbtr)
 
+   !=== LIST OF DEFINED ISOTOPES FAMILIES
+   CHARACTER(LEN=maxlen), SAVE, ALLOCATABLE :: isoFamilies(:)   !--- Generation 0 tracer name for each isotopes family (nbIso)
+
 CONTAINS
 
@@ -116 +120 @@
    USE iniprint_mod_h
    USE control_mod, ONLY: planet_type
+   USE ioipsl, ONLY: getin
    USE lmdz_reprobus_wrappers, ONLY: Init_chem_rep_trac
    USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA, CPPKEY_REPROBUS, CPPKEY_STRATAER
    USE dimensions_mod, ONLY: iim, jjm, llm, ndm
-IMPLICIT NONE
+   IMPLICIT NONE
 !==============================================================================================================================
 !
@@ -152 +157 @@
    INTEGER :: nqtrue                                                 !--- Tracers nb from tracer.def (no higher order moments)
    INTEGER :: iad                                                    !--- Advection scheme number
-   INTEGER :: iq, jq, nt, im, nm, ig                                 !--- Indexes and temporary variables
+   INTEGER :: iq, jq, it, nt, im, nm, ig                             !--- Indexes and temporary variables
    LOGICAL :: lerr
    TYPE(trac_type), ALLOCATABLE, TARGET :: ttr(:)
@@ -256 +261 @@
    IF( nqtot /= nqtrue ) THEN
       CALL msg('The choice of advection scheme for one or more tracers makes it necessary to add tracers')
-      CALL msg('The number of true tracers is '//TRIM(int2str(nqtrue)))
-      CALL msg('The total number of tracers needed is '//TRIM(int2str(nqtot)))
+      CALL msg('The number of true tracers is '//TRIM(num2str(nqtrue)))
+      CALL msg('The total number of tracers needed is '//TRIM(num2str(nqtot)))
    END IF
 
@@ -290 +295 @@
       IF(iad == -1) CALL abort_gcm(modname, msg1, 1)
 
-      !--- SET FIELDS longName and iadv
+      !--- SET FIELDS longName AND iadv
       t1%longName   = t1%name; IF(iad > 0) t1%longName=TRIM(t1%name)//descrq(iad)
       t1%iadv       = iad
@@ -319 +324 @@
       ig  = t1%iGeneration
       nam = t1%name
-      val = 'iadv='//TRIM(int2str(iad))
+      val = 'iadv='//TRIM(num2str(iad))
 
       !--- ONLY TESTED VALUES FOR TRACERS FOR NOW:               iadv = 14, 10 (and 0 for non-transported tracers)
@@ -335 +340 @@
    END DO
 
-   !=== READ PHYSICAL PARAMETERS FOR ISOTOPES ; DONE HERE BECAUSE dynetat0 AND iniacademic NEED "tnat" AND "alpha_ideal"
-   niso = 0; nzone = 0; nphas = nqo; ntiso = 0; isoCheck = .FALSE.
+   !=== DETERMINE ISOTOPES RELATED PARAMETERS ; DEFINE THE EXPLICIT KEYS iso_i*
    IF(processIsotopes()) CALL abort_gcm(modname, 'problem when processing isotopes parameters', 1)
+   iH2O = -1
+   IF(nbIso /= 0) THEN
+      IF(isoSelect('H2O', .TRUE.)) THEN
+         IF(isoSelect(1, .TRUE.)) CALL abort_physic(modname, "Can't select the first isotopes family", 1)
+      ELSE
+         iH2O = ixIso; CALL getin('ok_iso_verif', isotope%check)
+      END IF
+   END IF
+   isoFamilies = isoF(:)
 
    !--- Convection / boundary layer activation for all tracers
@@ -356 +369 @@
    !=== DISPLAY THE RESULTS
    IF(.NOT..TRUE.) RETURN
-   CALL msg('nqo    = '//TRIM(int2str(nqo)),    modname)
-   CALL msg('nbtr   = '//TRIM(int2str(nbtr)),   modname)
-   CALL msg('nqtrue = '//TRIM(int2str(nqtrue)), modname)
-   CALL msg('nqtot  = '//TRIM(int2str(nqtot)),  modname)
-   CALL msg('niso   = '//TRIM(int2str(niso)),   modname)
-   CALL msg('ntiso  = '//TRIM(int2str(ntiso)),  modname)
-   CALL msg('nqCO2  = '//TRIM(int2str(nqCO2)),  modname, CPPKEY_INCA)
-   CALL msg('nqINCA = '//TRIM(int2str(nqINCA)), modname, CPPKEY_INCA)
+   CALL msg('nqo    = '//TRIM(num2str(nqo)),    modname)
+   CALL msg('nbtr   = '//TRIM(num2str(nbtr)),   modname)
+   CALL msg('nqtrue = '//TRIM(num2str(nqtrue)), modname)
+   CALL msg('nqtot  = '//TRIM(num2str(nqtot)),  modname)
+   CALL msg('niso   = '//TRIM(num2str(niso)),   modname)
+   CALL msg('ntiso  = '//TRIM(num2str(ntiso)),  modname)
+   CALL msg('nqCO2  = '//TRIM(num2str(nqCO2)),  modname, CPPKEY_INCA)
+   CALL msg('nqINCA = '//TRIM(num2str(nqINCA)), modname, CPPKEY_INCA)
    t => tracers
    CALL msg('Information stored in '//TRIM(modname)//': ', modname)
@@ -374 +387 @@
    CALL msg('No isotopes identified.', modname, nbIso == 0)
    IF(nbIso == 0) RETURN
-   CALL msg('For isotopes family "H2O":', modname)
-   CALL msg('  isoKeys%name = '//strStack(isoKeys%name), modname)
-   CALL msg('  isoName = '//strStack(isoName),      modname)
-   CALL msg('  isoZone = '//strStack(isoZone),      modname)
-   CALL msg('  isoPhas = '//TRIM(isoPhas),          modname)
+   DO it = 1, nbIso
+      IF(isoSelect(it, .TRUE.)) CALL abort_physic(modname, 'Problem when selecting isotopes class', 1)
+      CALL msg('For isotopes family "'//TRIM(isoFamilies(it))//'":', modname)
+      CALL msg('  isoKeys%name = '//strStack(isoKeys%name), modname)
+      CALL msg('  isoName = '//strStack(isoName), modname)
+      CALL msg('  isoZone = '//strStack(isoZone), modname)
+      CALL msg('  isoPhas = '//TRIM(isoPhas),     modname)
+   END DO
+   IF(isoSelect(iH2O, .TRUE.)) lerr = isoSelect(1, .TRUE.)
 
 END SUBROUTINE init_infotrac