Changeset 5927 for LMDZ6/trunk/libf/phylmd/infotrac_phy.F90

Timestamp:

Dec 12, 2025, 5:26:10 PM (7 weeks ago)

Author:

Sebastien Nguyen

Message:

Changes to compile LMDZ-OR-ISO and wrtie output variables Rsol isotopes soil ratio) xtevap xtcoastal xtrivflu. Changes from CA and ND to write output variables xtprw (precipitatble water) uxt and vxt (meridional and zonal advected humidity) Rlandice and xtsnow.

File:

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

LMDZ6/trunk/libf/phylmd/infotrac_phy.F90

@@ -252 +252 @@
 !==============================================================================================================================
    nqtrue = SIZE(tracers)                                                                               !--- "true" tracers
-   nqo    =      COUNT(tracers(:)%component == 'lmdz' .AND. delPhase(tracers(:)%name)     == 'H2O')     !-- Water phases
+   nqo    =      COUNT(tracers(:)%component == 'lmdz' .AND. delPhase(tracers(:)%name)     == 'H2O')     !--- Water phases
    nqtke  =      COUNT(tracers(:)%component == 'lmdz' .AND. tracers(:)%name               == 'TKE')     !--- TKE tracers
    nbtr = nqtrue-COUNT(tracers(:)%component == 'lmdz' .AND. delPhase(tracers(:)%gen0Name) == 'H2O')     !--- Passed to phytrac
@@ -322 +322 @@
    CALL MOVE_ALLOC(FROM=ttr, TO=tracers)
 
-   !--- SET FIELDS iqParent, iqDescen, nqDescen, nqChildren
+   !--- SET FIELDS iqParent, iqDescen, nqDescen, nqChildren, iGeneration
    IF(indexUpdate(tracers)) CALL abort_physic(modname, 'problem with tracers indices update', 1)
 
    !=== DETERMINE ISOTOPES RELATED PARAMETERS ; DEFINE THE EXPLICIT KEYS iso_i*
    IF(processIsotopes()) CALL abort_physic(modname, 'Problem when processing isotopes parameters', 1)
+   !!SN do we really need to call isoSelect('H2O',.TRUE.) in processIsotopes and here ?
    iH2O = -1
    IF(nbIso /= 0) THEN
@@ -415 +416 @@
       CALL msg('  isoPhas = '//TRIM(isoPhas),     modname)
    END DO
+   !!SN do we need to call isoSelect(iH2O, .TRUE.) three times ?
    IF(isoSelect(iH2O, .TRUE.)) lerr = isoSelect(1, .TRUE.)