Changeset 5748 for LMDZ6/trunk/libf/phylmdiso/isotopes_mod.F90

Timestamp:

Jul 2, 2025, 12:00:08 PM (2 days ago)

Author:

dcugnet

Message:

• Use REAL(KIND=REAL32) and REAL(KIND=REAL64) Iinstead of REAL and DOUBLE PRECISION

to avoid ambiguity problems in generic procedure when reals are promoted to doubles.

• generic "num2str" replaces "str2int", "str2real", "str2dble" and "str2bool" functions.

File:

• LMDZ6/trunk/libf/phylmdiso/isotopes_mod.F90 (modified) (10 diffs)

LMDZ6/trunk/libf/phylmdiso/isotopes_mod.F90

@@ -3 +3 @@
 
 MODULE isotopes_mod
-   USE strings_mod,  ONLY: msg, real2str, int2str, bool2str, maxlen, strIdx, strStack
+   USE strings_mod,  ONLY: msg, num2str, maxlen, strIdx, strStack
    USE infotrac_phy, ONLY: isoName, niso, ntiso
    USE iso_params_mod
@@ -173 +173 @@
 
    !--- Check number of isotopes
-   CALL msg('64: niso = '//TRIM(int2str(niso)), modname)
+   CALL msg('64: niso = '//TRIM(num2str(niso)), modname)
 
          !--- Init de ntracisoOR: on ecrasera en cas de traceurs de tagging isotopiques
@@ -180 +180 @@
 
          !--- Type of water isotopes:
-         iso_eau = strIdx(isoName, 'H216O'); CALL msg('iso_eau='//int2str(iso_eau), modname)
-         iso_HDO = strIdx(isoName, 'HDO');   CALL msg('iso_HDO='//int2str(iso_HDO), modname)
-         iso_O18 = strIdx(isoName, 'H218O'); CALL msg('iso_O18='//int2str(iso_O18), modname)
-         iso_O17 = strIdx(isoName, 'H217O'); CALL msg('iso_O17='//int2str(iso_O17), modname)
-         iso_HTO = strIdx(isoName, 'HTO');   CALL msg('iso_HTO='//int2str(iso_HTO), modname)
+         iso_eau = strIdx(isoName, 'H216O'); CALL msg('iso_eau='//num2str(iso_eau), modname)
+         iso_HDO = strIdx(isoName, 'HDO');   CALL msg('iso_HDO='//num2str(iso_HDO), modname)
+         iso_O18 = strIdx(isoName, 'H218O'); CALL msg('iso_O18='//num2str(iso_O18), modname)
+         iso_O17 = strIdx(isoName, 'H217O'); CALL msg('iso_O17='//num2str(iso_O17), modname)
+         iso_HTO = strIdx(isoName, 'HTO');   CALL msg('iso_HTO='//num2str(iso_HTO), modname)
 
          !--- Initialisation: reading the isotopic parameters.
@@ -204 +204 @@
          CALL get_in('dsstlatcrit',   dsstlatcrit,      0.0)     !--- For modif_sst>=3
 #ifdef ISOVERIF
-         CALL msg('iso_init 270:  sstlatcrit='//real2str( sstlatcrit), modname, sstlatcrit < 0.0) !--- For modif_sst>=2
-         CALL msg('iso_init 279: dsstlatcrit='//real2str(dsstlatcrit), modname, sstlatcrit < 0.0) !--- For modif_sst>=3
+         CALL msg('iso_init 270:  sstlatcrit='//num2str( sstlatcrit), modname, sstlatcrit < 0.0) !--- For modif_sst>=2
+         CALL msg('iso_init 279: dsstlatcrit='//num2str(dsstlatcrit), modname, sstlatcrit < 0.0) !--- For modif_sst>=3
          IF(modif_sst >= 2 .AND. sstlatcrit < 0.0) STOP
 #endif              
@@ -234 +234 @@
          CALL get_in('region_nudge_qsol',   region_nudge_qsol,   1)
          nlevmaxO17 = 50
-         CALL msg('nlevmaxO17='//TRIM(int2str(nlevmaxO17)))
+         CALL msg('nlevmaxO17='//TRIM(num2str(nlevmaxO17)))
          CALL get_in('no_pce',   no_pce,     0)
          CALL get_in('A_satlim', A_satlim, 1.0)
          CALL get_in('ok_restrict_A_satlim', ok_restrict_A_satlim, 0)
 #ifdef ISOVERIF
-         CALL msg(' 315: A_satlim='//real2str(A_satlim), modname, A_satlim > 1.0)
+         CALL msg(' 315: A_satlim='//num2str(A_satlim), modname, A_satlim > 1.0)
          IF(A_satlim > 1.0) STOP
 #endif
@@ -256 +256 @@
          CALL get_in('tnateq1', ltnat1, .TRUE.)
 
-         CALL msg('iso_O18, iso_HDO, iso_eau = '//TRIM(strStack(int2str([iso_O18, iso_HDO, iso_eau]))), modname)
+         CALL msg('iso_O18, iso_HDO, iso_eau = '//TRIM(strStack(num2str([iso_O18, iso_HDO, iso_eau]))), modname)
 
          !--------------------------------------------------------------
@@ -363 +363 @@
          CALL msg('285: verif initialisation:', modname)
          DO ixt=1,niso
-            sxt=int2str(ixt)
+            sxt=num2str(ixt)
             CALL msg(' * isoName('//TRIM(sxt)//') = <'//TRIM(isoName(ixt))//'>',  modname)
-            CALL msg(  '    tnat('//TRIM(sxt)//') = '//TRIM(real2str(tnat(ixt))), modname)
-!           CALL msg('    alpha_liq_sol('//TRIM(sxt)//') = '//TRIM(real2str(alpha_liq_sol(ixt))), modname)
-!           CALL msg(        '   tkcin0('//TRIM(sxt)//') = '//TRIM(real2str(tkcin0(ixt))),        modname)
-!           CALL msg(       '   tdifrel('//TRIM(sxt)//') = '//TRIM(real2str(tdifrel(ixt))),       modname)
+            CALL msg(  '    tnat('//TRIM(sxt)//') = '//TRIM(num2str(tnat(ixt))), modname)
+!           CALL msg('    alpha_liq_sol('//TRIM(sxt)//') = '//TRIM(num2str(alpha_liq_sol(ixt))), modname)
+!           CALL msg(        '   tkcin0('//TRIM(sxt)//') = '//TRIM(num2str(tkcin0(ixt))),        modname)
+!           CALL msg(       '   tdifrel('//TRIM(sxt)//') = '//TRIM(num2str(tdifrel(ixt))),       modname)
          END DO
-         CALL msg('69:     lambda = '//TRIM(real2str(lambda_sursat)), modname)
-         CALL msg('69:    thumxt1 = '//TRIM(real2str(thumxt1)),       modname)
-         CALL msg('69: h_land_ice = '//TRIM(real2str(h_land_ice)),    modname)
-         CALL msg('69:      P_veg = '//TRIM(real2str(P_veg)),         modname)
+         CALL msg('69:     lambda = '//TRIM(num2str(lambda_sursat)), modname)
+         CALL msg('69:    thumxt1 = '//TRIM(num2str(thumxt1)),       modname)
+         CALL msg('69: h_land_ice = '//TRIM(num2str(h_land_ice)),    modname)
+         CALL msg('69:      P_veg = '//TRIM(num2str(P_veg)),         modname)
 
 END SUBROUTINE iso_init
@@ -411 +411 @@
       IF(PRESENT(def)) val=def; CALL getin(nam,val)
       lD=.TRUE.; IF(PRESENT(lDisp)) lD=lDisp
-      IF(lD) CALL msg(TRIM(nam)//' = '//TRIM(int2str(val)))
+      IF(lD) CALL msg(TRIM(nam)//' = '//TRIM(num2str(val)))
   END IF
   CALL bcast(val)
@@ -430 +430 @@
       IF(PRESENT(def)) val=def; CALL getin(nam,val)
       lD=.TRUE.; IF(PRESENT(lDisp)) lD=lDisp
-      IF(lD) CALL msg(TRIM(nam)//' = '//TRIM(real2str(val)))
+      IF(lD) CALL msg(TRIM(nam)//' = '//TRIM(num2str(val)))
   END IF
   CALL bcast(val)
@@ -449 +449 @@
       IF(PRESENT(def)) val=def; CALL getin(nam,val)
       lD=.TRUE.; IF(PRESENT(lDisp)) lD=lDisp
-      IF(lD) CALL msg(TRIM(nam)//' = '//TRIM(bool2str(val)))
+      IF(lD) CALL msg(TRIM(nam)//' = '//TRIM(num2str(val)))
   END IF
   CALL bcast(val)