Changeset 5200 for LMDZ6/trunk/libf/dyn3d_common

Timestamp:

Sep 19, 2024, 4:11:02 PM (5 months ago)

Author:

dcugnet

Message:

Replace the hard-coded "tnat1" (TRUE by default) by a variable that can be changed using "tnateq1" in *.def files.
"iso_verif_dyn" now uses also this variable (fix)

File:

• LMDZ6/trunk/libf/dyn3d_common/iso_verif_dyn.F (modified) (2 diffs)

LMDZ6/trunk/libf/dyn3d_common/iso_verif_dyn.F

@@ -64 +64 @@
         function iso_verif_aberrant_nostop
      :           (x,iso,q,err_msg)
+#ifdef CPP_IOIPSL
+        USE IOIPSL, ONLY: getin
+#else
+        USE ioipsl_getincom, ONLY: getin
+#endif
         USE infotrac, ONLY: isoName, getKey
         implicit none
@@ -77 +82 @@
         parameter (qmin=1e-11)
         parameter (deltaDmax=200.0,deltaDmin=-999.9)
+        LOGICAL, SAVE :: ltnat1
+        LOGICAL, SAVE :: lFirst=.TRUE.
 
         ! output
         integer iso_verif_aberrant_nostop
 
+        IF(lFirst) THEN
+           ltnat1 = .TRUE.; CALL getin('tnateq1', ltnat1)
+           lFirst = .FALSE.
+        END IF
         iso_verif_aberrant_nostop=0
 
         ! verifier que HDO est raisonable
          if (q.gt.qmin) then
-             IF(getKey('tnat', tnat, isoName(iso))) THEN
+             IF(ltnat1) THEN
+                tnat = 1.0
+             ELSE IF(getKey('tnat', tnat, isoName(iso))) THEN
                   err_msg = 'Missing isotopic parameter "tnat"'
                   iso_verif_aberrant_nostop=1