Changeset 5748 for LMDZ6/trunk/libf/dyn3d/check_isotopes.f90

Timestamp:

Jul 2, 2025, 12:00:08 PM (3 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/dyn3d/check_isotopes.f90 (modified) (5 diffs)

LMDZ6/trunk/libf/dyn3d/check_isotopes.f90

@@ -1 +1 @@
 SUBROUTINE check_isotopes_seq(q, ip1jmp1, err_msg)
-   USE strings_mod, ONLY: maxlen, msg, strIdx, strStack, int2str, real2str
+   USE strings_mod, ONLY: maxlen, msg, strIdx, strStack, num2str
    USE infotrac,    ONLY: nqtot, niso, nphas, isotope, isoCheck, iqIsoPha, isoSelect, &
                           ntiso, iH2O, nzone, tracers, isoName,  itZonIso
@@ -80 +80 @@
                   CYCLE
                END IF
-               CALL msg('ixt, iq = '//TRIM(strStack(int2str([ixt,iq]))), modname)
-               msg1 = '('//TRIM(strStack(int2str([i,k])))//')'
-               CALL msg(TRIM(tracers(iqpar)%name)//TRIM(msg1)//' = '//TRIM(real2str(q1)), modname)
-               CALL msg(TRIM(tracers(iq   )%name)//TRIM(msg1)//' = '//TRIM(real2str(q2)), modname)
+               CALL msg('ixt, iq = '//TRIM(strStack(num2str([ixt,iq]))), modname)
+               msg1 = '('//TRIM(strStack(num2str([i,k])))//')'
+               CALL msg(TRIM(tracers(iqpar)%name)//TRIM(msg1)//' = '//TRIM(num2str(q1)), modname)
+               CALL msg(TRIM(tracers(iq   )%name)//TRIM(msg1)//' = '//TRIM(num2str(q2)), modname)
                CALL abort_gcm(modname, 'Error with isotopes: '//TRIM(err_msg), 1)
             END DO
@@ -111 +111 @@
                deltaD = (q2/q1/tnat(ixt)-1.)*1000.
                IF(deltaD <= deltaDmax .AND. deltaD >= deltaDmin) CYCLE
-               CALL msg('ixt, iq = '//TRIM(strStack(int2str([ixt,iq]))), modname)
-               msg1 = '('//TRIM(strStack(int2str([i,k])))//')'
-               CALL msg(TRIM(tracers(iqpar)%name)//TRIM(msg1)//' = '//TRIM(real2str(q1)), modname)
-               CALL msg(TRIM(tracers(iq   )%name)//TRIM(msg1)//' = '//TRIM(real2str(q2)), modname)
-               CALL msg(TRIM(nm(iiso))//TRIM(real2str(deltaD)), modname)
+               CALL msg('ixt, iq = '//TRIM(strStack(num2str([ixt,iq]))), modname)
+               msg1 = '('//TRIM(strStack(num2str([i,k])))//')'
+               CALL msg(TRIM(tracers(iqpar)%name)//TRIM(msg1)//' = '//TRIM(num2str(q1)), modname)
+               CALL msg(TRIM(tracers(iq   )%name)//TRIM(msg1)//' = '//TRIM(num2str(q2)), modname)
+               CALL msg(TRIM(nm(iiso))//TRIM(num2str(deltaD)), modname)
                CALL abort_gcm(modname, 'Error with isotopes: '//TRIM(err_msg), 1)
             END DO
@@ -140 +140 @@
                   deltaD = (q2/q1/tnat(iso_HDO)-1.)*1000.
                   IF(deltaD <= deltaDmax .AND. deltaD >= deltaDmin) CYCLE
-                  CALL msg('izon, ipha = '//TRIM(strStack(int2str([izon, ipha]))), modname)
-                  CALL msg( 'ixt, ieau = '//TRIM(strStack(int2str([ ixt, ieau]))), modname)
-                  msg1 = '('//TRIM(strStack(int2str([i,k])))//')'
-                  CALL msg(TRIM(tracers(iqeau)%name)//TRIM(msg1)//' = '//TRIM(real2str(q1)), modname)
-                  CALL msg(TRIM(tracers(iq   )%name)//TRIM(msg1)//' = '//TRIM(real2str(q2)), modname)
-                  CALL msg('deltaD = '//TRIM(real2str(deltaD)), modname)
+                  CALL msg('izon, ipha = '//TRIM(strStack(num2str([izon, ipha]))), modname)
+                  CALL msg( 'ixt, ieau = '//TRIM(strStack(num2str([ ixt, ieau]))), modname)
+                  msg1 = '('//TRIM(strStack(num2str([i,k])))//')'
+                  CALL msg(TRIM(tracers(iqeau)%name)//TRIM(msg1)//' = '//TRIM(num2str(q1)), modname)
+                  CALL msg(TRIM(tracers(iq   )%name)//TRIM(msg1)//' = '//TRIM(num2str(q2)), modname)
+                  CALL msg('deltaD = '//TRIM(num2str(deltaD)), modname)
                   CALL abort_gcm(modname, 'Error with isotopes: '//TRIM(err_msg), 1)
                END DO
@@ -163 +163 @@
                IF(ABS(xtractot-xiiso) > errmax .AND. ABS(xtractot-xiiso)/MAX(MAX(ABS(xtractot),ABS(xiiso)),1e-18) > errmaxrel) THEN
                   CALL msg('Error in iso_verif_aberrant trac: '//TRIM(err_msg))
-                  CALL msg('iiso, ipha = '//TRIM(strStack(int2str([iiso, ipha]))), modname)
-                  CALL msg('q('//TRIM(strStack(int2str([i,k])))//',:) = '//TRIM(strStack(real2str(q(i,k,:)))), modname)
+                  CALL msg('iiso, ipha = '//TRIM(strStack(num2str([iiso, ipha]))), modname)
+                  CALL msg('q('//TRIM(strStack(num2str([i,k])))//',:) = '//TRIM(strStack(num2str(q(i,k,:)))), modname)
                   CALL abort_gcm(modname, 'Error with isotopes: '//TRIM(err_msg), 1)
                END IF