Changeset 5229 for LMDZ6/branches/Amaury_dev/libf/dyn3dmem/dynetat0_loc.f90

Timestamp:

Sep 25, 2024, 12:03:08 PM (5 weeks ago)

Author:

abarral

Message:

Merge r5214

File:

• LMDZ6/branches/Amaury_dev/libf/dyn3dmem/dynetat0_loc.f90 (modified) (3 diffs)

LMDZ6/branches/Amaury_dev/libf/dyn3dmem/dynetat0_loc.f90

@@ -8 +8 @@
   USE parallel_lmdz
   USE lmdz_infotrac, ONLY: nqtot, tracers, niso, iqIsoPha, iH2O, isoName, &
-          new2oldH2O, newHNO3, oldHNO3, getKey
+          new2oldH2O, newHNO3, oldHNO3
   USE lmdz_strings, ONLY: maxlen, msg, strStack, real2str, int2str, strIdx
   USE netcdf, ONLY: nf90_open, nf90_nowrite, nf90_inquire_dimension, nf90_inq_varid, &
@@ -23 +23 @@
   USE lmdz_iniprint, ONLY: lunout, prt_level
   USE lmdz_comgeom
+  USE iso_params_mod   ! tnat_* and alpha_ideal_*
 
   USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
@@ -188 +189 @@
           CALL msg(' !!!  Beware: alpha_ideal put to 1  !!!', modname)
         ELSE
-          IF(getKey('tnat', tnat, isoName(iName)) .OR. getKey('alpha', alpha_ideal, isoName(iName))) &
-                  CALL abort_gcm(TRIM(modname), 'missing isotopic parameters', 1)
+          SELECT CASE(isoName(iName))
+            CASE('H216O'); tnat = tnat_H216O; alpha_ideal = alpha_ideal_H216O
+            CASE('H217O'); tnat = tnat_H217O; alpha_ideal = alpha_ideal_H217O
+            CASE('H218O'); tnat = tnat_H218O; alpha_ideal = alpha_ideal_H218O
+            CASE('HDO');   tnat = tnat_HDO;   alpha_ideal = alpha_ideal_HDO
+            CASE('HTO');   tnat = tnat_HTO;   alpha_ideal = alpha_ideal_HTO
+            CASE DEFAULT; CALL abort_gcm(TRIM(modname), 'unknown isotope "' // TRIM(isoName(iName)) // '" ; check tracer.def file', 1)
+          END SELECT
         END IF
         CALL msg('Missing tracer <' // TRIM(var) // '> => initialized with a simplified Rayleigh distillation law.', modname)