Changeset 5791 for LMDZ6/branches/contrails/libf/phylmdiso/isotrac_mod.F90

Timestamp:

Jul 28, 2025, 7:23:15 PM (7 days ago)

Author:

aborella

Message:

Merge with trunk r5789

Location:

LMDZ6/branches/contrails

Files:

• . (modified) (1 prop)
• libf/phylmdiso/isotrac_mod.F90 (modified) (15 diffs)

LMDZ6/branches/contrails/libf/phylmdiso/isotrac_mod.F90

@@ -4 +4 @@
 MODULE isotrac_mod
   USE infotrac_phy, ONLY: niso, ntiso, nzone, delPhase
-  USE isotopes_mod, ONLY: ridicule, get_in
-
+  USE isotopes_mod, ONLY: ridicule
+  USE ioipsl_getin_p_mod, ONLY : getin_p
   IMPLICIT NONE
   SAVE
@@ -148 +148 @@
    USE isotopes_mod, ONLY: iso_eau, ntracisoOR, initialisation_iso
    USE dimphy,       ONLY: klon, klev
-   USE  strings_mod, ONLY: int2str, strStack, strTail, strHead, strIdx, fmsg
+   USE  strings_mod, ONLY: num2str, strStack, strTail, strHead, strIdx, maxlen
 
    IMPLICIT NONE
@@ -156 +156 @@
    INTEGER :: izone_pres, izone_lat
    INTEGER :: nzone_opt
-
-   IF(fmsg("traceurs_init 18: isotrac ne marche que si on met l'eau comme isotope", 'iso_traceurs_init', iso_eau==0)) STOP
+   LOGICAL :: lerr
+   CHARACTER(LEN=maxlen) :: modname
+
+   modname = 'iso_traceurs_init'
+   lerr = iso_eau == 0
+   IF(lerr) CALL abort_physics(TRIM(modname)//' 18', 'isotrac does not work without H216O isotope', 1)
 
    !--- Initialize
@@ -171 +175 @@
    ALLOCATE( boite_map(klon,klev))
 
-   IF(initialisation_iso == 0) CALL get_in('initialisation_isotrac', initialisation_isotrac)
+   IF(initialisation_iso == 0) CALL getin_p('initialisation_isotrac', initialisation_isotrac)
 
    !--- Read tracing option
-   CALL get_in('option_traceurs', option_traceurs)
+   CALL getin_p('option_traceurs', option_traceurs)
 
    !--- Genral case: no traceurs in ORCHIDEE
@@ -214 +218 @@
       !========================================================================================================================
          ! lire les use_bassin
-         CALL get_in('use_bassin_Atlantic',   use_bassin_Atlantic)
-         CALL get_in('use_bassin_Medit',      use_bassin_Medit)
-         CALL get_in('use_bassin_Indian',     use_bassin_Indian)
-         CALL get_in('use_bassin_Austral',    use_bassin_Austral)
-         CALL get_in('use_bassin_Pacific',    use_bassin_Pacific)
-         CALL get_in('use_bassin_MerArabie',  use_bassin_MerArabie)
-         CALL get_in('use_bassin_BengalGolf', use_bassin_BengalGolf)
-         CALL get_in('use_bassin_SouthIndian',use_bassin_SouthIndian)
-         CALL get_in('use_bassin_Tropics',    use_bassin_Tropics)
-         CALL get_in('use_bassin_Midlats',    use_bassin_Midlats)
-         CALL get_in('use_bassin_HighLats',   use_bassin_HighLats)
+         CALL getin_p('use_bassin_Atlantic',   use_bassin_Atlantic)
+         CALL getin_p('use_bassin_Medit',      use_bassin_Medit)
+         CALL getin_p('use_bassin_Indian',     use_bassin_Indian)
+         CALL getin_p('use_bassin_Austral',    use_bassin_Austral)
+         CALL getin_p('use_bassin_Pacific',    use_bassin_Pacific)
+         CALL getin_p('use_bassin_MerArabie',  use_bassin_MerArabie)
+         CALL getin_p('use_bassin_BengalGolf', use_bassin_BengalGolf)
+         CALL getin_p('use_bassin_SouthIndian',use_bassin_SouthIndian)
+         CALL getin_p('use_bassin_Tropics',    use_bassin_Tropics)
+         CALL getin_p('use_bassin_Midlats',    use_bassin_Midlats)
+         CALL getin_p('use_bassin_HighLats',   use_bassin_HighLats)
          nzone_opt  =  2  +  COUNT([use_bassin_Atlantic, use_bassin_Medit,     use_bassin_Indian,     &
             use_bassin_Austral,     use_bassin_Pacific,  use_bassin_MerArabie, use_bassin_BengalGolf, &
@@ -291 +295 @@
          option_cond=0
          DO izone=1,nzone_temp
-            strtrac(izone) = 't'//TRIM(int2str(izone))
+            strtrac(izone) = 't'//TRIM(num2str(izone))
             WRITE(*,*) 'izone, strtrac=', izone, strtrac(izone)
          END DO
@@ -352 +356 @@
          ! 1 par defaut pour colorier a la fois condensat LS et condensat convectif.
          ! Mais on peut mettre 2 si on ne veut que colorier que le condensat convectif.
-         CALL get_in('option_cond',option_cond)
+         CALL getin_p('option_cond',option_cond)
          strtrac(izone_poubelle)='res'
          strtrac(izone_cond)='con'
@@ -435 +439 @@
          option_cond=0
          DO izone=1,nzone_temp
-            strtrac(izone) = 't'//TRIM(int2str(izone))
+            strtrac(izone) = 't'//TRIM(num2str(izone))
             WRITE(*,*) 'izone, strtrac = ', izone, strtrac(izone)
          END DO
@@ -470 +474 @@
             DO izone_lat=1,nzone_lat
                izone=izone_lat+(izone_pres-1)*nzone_lat
-               strtrac(izone) = 't'//TRIM(int2str(izone_pres))//TRIM(int2str(izone_lat))
+               strtrac(izone) = 't'//TRIM(num2str(izone_pres))//TRIM(num2str(izone_lat))
                write(*,*) 'izone_pres, izone_lat, izone, strtrac = ',izone_pres, izone_lat, izone, strtrac(izone)
             END DO
@@ -530 +534 @@
          izone_cond=nzone_temp+2
          DO izone=1,nzone_temp
-            strtrac(izone) = 't'//TRIM(int2str(izone))
+            strtrac(izone) = 't'//TRIM(num2str(izone))
             WRITE(*,*) 'izone, strtrac = ', izone, strtrac(izone)
          END DO !do izone=1,nzone_temp
@@ -576 +580 @@
          option_cond=1
          DO izone=1,nzone_temp
-            strtrac(izone) = 't'//TRIM(int2str(izone))
+            strtrac(izone) = 't'//TRIM(num2str(izone))
             WRITE(*,*) 'izone, strtrac = ', izone, strtrac(izone)
          END DO
@@ -586 +590 @@
       CASE(20)     !=== TRACING TROPICAL/EXTRATROPICAL/EXTRATROPICAL RECYCLING TO STUDY HUMIDITY AND SUBTROPICAL ISOTOPES CONTROL
       !========================================================================================================================
-         CALL get_in('lim_tag20', lim_tag20, 35.0)
+         CALL getin_p('lim_tag20', lim_tag20, 35.0)
          nzone_opt=3
          izone_cont=1
@@ -621 +625 @@
       CASE(22)     !=== TRACING WATER VAPOUR PROCESSED IN THE 3-LEVELS SCONVECTION ZONES BT, MT AND UT
       !========================================================================================================================
-         CALL get_in('lim_precip_tag22', lim_precip_tag22, 20.0)
+         CALL getin_p('lim_precip_tag22', lim_precip_tag22, 20.0)
          nzone_opt=3
          izone_cont=1
@@ -649 +653 @@
    IF(option_tmin == 1) THEN
       seuil_tag_tmin = 0.01
-      CALL get_in('seuil_tag_tmin',        seuil_tag_tmin,        0.01)
-      CALL get_in('seuil_tag_tmin_ls',     seuil_tag_tmin_ls,     seuil_tag_tmin)
-      CALL get_in('option_seuil_tag_tmin', option_seuil_tag_tmin, 1)
+      CALL getin_p('seuil_tag_tmin',        seuil_tag_tmin,        0.01)
+      CALL getin_p('seuil_tag_tmin_ls',     seuil_tag_tmin_ls,     seuil_tag_tmin)
+      CALL getin_p('option_seuil_tag_tmin', option_seuil_tag_tmin, 1)
    END IF
 
@@ -659 +663 @@
 #ifdef ISOVERIF
    WRITE(*,*) 'traceurs_init 65: bilan de l''init:'
-   WRITE(*,*) 'index_zone = '//TRIM(strStack(int2str(index_zone(1:ntiso))))
-   WRITE(*,*) 'index_iso  = '//TRIM(strStack(int2str(index_iso (1:ntiso))))
+   WRITE(*,*) 'index_zone = '//TRIM(strStack(num2str(index_zone(1:ntiso))))
+   WRITE(*,*) 'index_iso  = '//TRIM(strStack(num2str(index_iso (1:ntiso))))
    DO izone=1,nzone
-      WRITE(*,*)'itZonIso('//TRIM(int2str(izone))//',:) = '//strStack(int2str(itZonIso(izone,:)))
+      WRITE(*,*)'itZonIso('//TRIM(num2str(izone))//',:) = '//strStack(num2str(itZonIso(izone,:)))
    END DO
    DO izone=1,nzone
-      WRITE(*,*)'strtrac('//TRIM(int2str(izone))//',:) = '//TRIM(strtrac(izone))
+      WRITE(*,*)'strtrac('//TRIM(num2str(izone))//',:) = '//TRIM(strtrac(izone))
    END DO
    WRITE(*,*) 'ntracisoOR=',ntracisoOR