Changeset 5756 for LMDZ6/trunk/libf/phylmd

Timestamp:

Jul 3, 2025, 3:25:59 PM (3 days ago)

Author:

dcugnet

Message:

Add "isoFamilies", the list of defined isotopes families (==H2O? for now).

File:

• LMDZ6/trunk/libf/phylmd/infotrac_phy.F90 (modified) (11 diffs)

LMDZ6/trunk/libf/phylmd/infotrac_phy.F90

@@ -3 +3 @@
 MODULE infotrac_phy
 
-   USE       strings_mod, ONLY: msg, maxlen, cat, dispTable, num2str, strStack, strParse, strCount, strIdx
+   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,         &
-          niso,   ntiso,   nzone,   nphas,   maxTableWidth,  iqIsoPha, iqWIsoPha, ixIso, new2oldH2O
+          niso,   ntiso,   nzone,   nphas, isoF=>isoFamilies,iqIsoPha, iqWIsoPha, ixIso, new2oldH2O
    IMPLICIT NONE
 
@@ -24 +24 @@
    !=== 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 family selection tool + selected family index
    !=== FOR ISOTOPES: Specific to water
    PUBLIC :: iH2O                                          !--- Value of "ixIso" for "H2O" isotopes class
@@ -84 +84 @@
 !  | 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                |
 !  +-------------+------------------------------------------------------+-------------+------------------------+
 !
@@ -106 +106 @@
 !$OMP THREADPRIVATE(ivap, iliq, isol, ibs, icf, irvc)
 
-   !=== 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)
@@ -118 +118 @@
    INTEGER, SAVE, ALLOCATABLE :: conv_flg(:), pbl_flg(:)        !--- Convection / boundary layer activation (nbtr)
 !$OMP THREADPRIVATE(conv_flg, pbl_flg)
+
+   !=== LIST OF DEFINED ISOTOPES FAMILIES
+   CHARACTER(LEN=maxlen), SAVE, ALLOCATABLE :: isoFamilies(:)   !--- Generation 0 tracer name for each isotopes family (nbIso)
+!$OMP THREADPRIVATE(isoFamilies)
 
   !=== SPECIFIC TO STRATOSPHERIC AEROSOLS (CK/OB)
@@ -167 +171 @@
    INTEGER :: nqtrue                                                 !--- Tracers nb from tracer.def (no higher order moments)
    INTEGER :: iad                                                    !--- Advection scheme number
-   INTEGER :: iq, jq, nt, im, nm                                     !--- Indexes and temporary variables
+   INTEGER :: iq, jq, it, nt, im, nm                                 !--- Indexes and temporary variables
    LOGICAL :: lerr, lInit
    TYPE(trac_type), ALLOCATABLE, TARGET :: ttr(:)
@@ -294 +298 @@
       IF(iad == -1) CALL abort_physic(modname, msg1, 1)
 
-      !--- SET FIELDS longName, isInPhysics
+      !--- SET FIELDS longName AND isInPhysics
       t1%longName   = t1%name; IF(iad > 0) t1%longName=TRIM(t1%name)//descrq(iad)
       t1%isInPhysics= iad >= 0 .AND. (delPhase(t1%gen0Name) /= 'H2O' .OR. t1%component /= 'lmdz')
@@ -315 +319 @@
    CALL MOVE_ALLOC(FROM=ttr, TO=tracers)
 
-   !--- SET FIELDS iqParent, iqDescen, nqDescen, nqChildren, iGeneration
+   !--- SET FIELDS iqParent, iqDescen, nqDescen, nqChildren
    IF(indexUpdate(tracers)) CALL abort_physic(modname, 'problem with tracers indices update', 1)
 
-   !=== READ PHYSICAL PARAMETERS FOR ISOTOPES
-   niso = 0; nzone = 0; nphas = nqo; ntiso = 0; isoCheck = .FALSE.
+   !=== DETERMINE ISOTOPES RELATED PARAMETERS ; DEFINE THE EXPLICIT KEYS iso_i*
    IF(processIsotopes()) CALL abort_physic(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_p('ok_iso_verif', isotope%check)
+      END IF
+   END IF
 
 !##############################################################################################################################
@@ -326 +337 @@
 !##############################################################################################################################
    DO iq = 1, nqtrue
-      t1 => tracers(iq)
       IF(hadv(iq)     ==    vadv(iq)    ) iad = hadv(iq)
       IF(hadv(iq)==10 .AND. vadv(iq)==16) iad = 11
-      tracers(iq)%isInPhysics= iad >= 0 .AND. (delPhase(t1%gen0Name) /= 'H2O' .OR. t1%component /= 'lmdz')
+      tracers(iq)%isInPhysics= iad >= 0 .AND. (delPhase(tracers(iq)%gen0Name) /= 'H2O' .OR. tracers(iq)%component /= 'lmdz')
    END DO
-!##############################################################################################################################
-   END IF
-!##############################################################################################################################
+   tracers(nqtrue+1:nqtot)%isInPhysics = .FALSE.
+!##############################################################################################################################
+   END IF
+!##############################################################################################################################
+   isoFamilies = isoF(:)
 
    !--- Convection / boundary layer activation for all tracers
@@ -377 +389 @@
    IF(dispTable('issssssssiiiiiiii', ['iq  ', 'name', 'lNam', 'g0Nm', 'prnt', 'type', 'phas', 'comp',      &
                               'isPh', 'iGen', 'iqPr', 'nqDe', 'nqCh', 'iGrp', 'iNam', 'iZon', 'iPha'],     &
-      cat(t%name, t%longName, t%gen0Name, t%parent, t%type, t%phase, t%component, num2str(t%isInPhysics)), &
+      cat(t%name, t%longName,  t%gen0Name, t%parent, t%type, t%phase, t%component, num2str(t%isInPhysics)),&
       cat([(iq, iq=1, nqtot)], t%iGeneration, t%iqParent, t%nqDescen, t%nqChildren, t%iso_iGroup,          &
                   t%iso_iName, t%iso_iZone, t%iso_iPhase), nColMax=maxTableWidth, nHead=2, sub=modname))   &
@@ -383 +395 @@
    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.)
 
    IF(CPPKEY_STRATAER .AND. type_trac == 'coag') THEN