Changeset 5214 for LMDZ6/trunk/libf/dyn3d/check_isotopes.F90

Timestamp:

Sep 22, 2024, 10:07:56 PM (4 weeks ago)

Author:

dcugnet

Message:

The fortran parameters file "iso_params_mod.F90" is introduced so that "tnat" and "alpha_ideal" are defined in a single place but used in several.
The "getKey" routine is only used in "infotrac" and "infotrac_phy" routines, but could be used outside this scope to get tracers parameters (read from "tracer.def") or isotopic parameters (read from "isotopes_params.def" - disabled for now).

File:

• LMDZ6/trunk/libf/dyn3d/check_isotopes.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/dyn3d/check_isotopes.F90

@@ -2 +2 @@
    USE strings_mod, ONLY: maxlen, msg, strIdx, strStack, int2str, real2str
    USE infotrac,    ONLY: nqtot, niso, nphas, isotope, isoCheck, iqIsoPha, isoSelect, &
-                          ntiso, iH2O, nzone, tracers, isoName,  itZonIso, getKey
+                          ntiso, iH2O, nzone, tracers, isoName,  itZonIso
+   USE iso_params_mod,  ONLY: tnat_H216O, tnat_H217O, tnat_H218O, tnat_HDO, tnat_HTO
 #ifdef CPP_IOIPSL
    USE ioipsl,          ONLY: getin
@@ -15 +16 @@
    CHARACTER(LEN=maxlen) :: modname, msg1, nm(2)
    INTEGER :: ixt, ipha, k, i, iq, iiso, izon, ieau, iqeau, iqpar
-   INTEGER, ALLOCATABLE ::   ix(:)
+   INTEGER, ALLOCATABLE       ::   ix(:)
    REAL,    ALLOCATABLE, SAVE :: tnat(:)
-   REAL    :: xtractot, xiiso, deltaD, q1, q2
+   REAL :: xtractot, xiiso, deltaD, q1, q2
    REAL, PARAMETER :: borne     = 1e19,  &
                       errmax    = 1e-8,  &       !--- Max. absolute error
@@ -26 +27 @@
                       ridicule  = 1e-12
    INTEGER, SAVE :: iso_eau, iso_O17, iso_O18, iso_HDO, iso_HTO
-   LOGICAL, SAVE :: ltnat1, first=.TRUE.
+   LOGICAL       :: ltnat1
+   LOGICAL, SAVE :: first=.TRUE.
 
    modname='check_isotopes'
@@ -35 +37 @@
       ltnat1 = .TRUE.; CALL getin('tnateq1', ltnat1)
       ALLOCATE(tnat(niso))
-      iso_eau = strIdx(isoName,'H216O')
-      iso_O17 = strIdx(isoName,'H217O')
-      iso_O18 = strIdx(isoName,'H218O')
-      iso_HDO = strIdx(isoName,'HDO')
-      iso_HTO = strIdx(isoName,'HTO')
-      IF(ltnat1) THEN
-         tnat(:)=1.0
-      ELSE
-         IF(getKey('tnat', tnat)) CALL abort_gcm(modname, 'missing isotopic parameter', 1)
-      END IF
+      iso_eau = strIdx(isoName,'H216O'); IF(iso_eau /= 0) tnat(iso_eau) = tnat_H216O
+      iso_O17 = strIdx(isoName,'H217O'); IF(iso_O17 /= 0) tnat(iso_O17) = tnat_H217O
+      iso_O18 = strIdx(isoName,'H218O'); IF(iso_O18 /= 0) tnat(iso_O18) = tnat_H218O
+      iso_HDO = strIdx(isoName,'HDO');   IF(iso_HDO /= 0) tnat(iso_HDO) = tnat_HDO
+      iso_HTO = strIdx(isoName,'HTO');   IF(iso_HTO /= 0) tnat(iso_HTO) = tnat_HTO
+      IF(ltnat1) tnat(:) = 1.
       first = .FALSE.
    END IF
@@ -56 +54 @@
          DO k = 1, llm
             DO i = 1, ip1jmp1
-               IF(ABS(q(i,k,iq)) <= borne) CYCLE
+               IF(ABS(q(i,k,iq)) < borne) CYCLE
                WRITE(msg1,'(s,"(",i0,",",i0,",",i0,") = ",ES12.4)')TRIM(isoName(ixt)),i,k,iq,q(i,k,iq)
                CALL msg(msg1, modname)