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

Timestamp:

Jul 3, 2025, 3:25:59 PM (3 weeks ago)

Author:

dcugnet

Message:

Add "isoFamilies", the list of defined isotopes families (==H2O? for now).

File:

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

LMDZ6/trunk/libf/phylmdiso/isotopes_mod.F90

@@ -4 +4 @@
 MODULE isotopes_mod
    USE strings_mod,  ONLY: msg, num2str, maxlen, strIdx, strStack
-   USE infotrac_phy, ONLY: isoName, niso, ntiso
+   USE infotrac_phy, ONLY: isoName, niso, ntiso, nbIso, isoFamilies, isoSelect, isoCheck
    USE iso_params_mod
    USE ioipsl_getin_p_mod, ONLY : getin_p
@@ -150 +150 @@
 
    !=== Local variables:
-   INTEGER :: ixt, is
+   INTEGER :: ixt, is, ii
    LOGICAL :: ltnat1
    CHARACTER(LEN=maxlen) :: modname, sxt
@@ -175 +175 @@
    CALL msg('64: niso = '//TRIM(num2str(niso)), modname)
 
+   DO ii = 1, nbIso
+      CALL msg('Can''t select isotopes class "'//TRIM(isoFamilies(ii))//'"', modname, isoSelect(ii, lVerbose=.TRUE.))
+
+!==============================================================================================================================
+      IF(isoFamilies(ii) == 'H2O') THEN
+!==============================================================================================================================
          !--- Init de ntracisoOR: on ecrasera en cas de traceurs de tagging isotopiques
          !                     (nzone>0) si complications avec ORCHIDEE
@@ -203 +209 @@
          CALL getin_p( 'sstlatcrit',    sstlatcrit,     30.0)     !--- For modif_sst>=3
          CALL getin_p('dsstlatcrit',   dsstlatcrit,      0.0)     !--- For modif_sst>=3
-#ifdef ISOVERIF
+         IF(isoCheck) THEN
          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              
-
+         END IF
          CALL getin_p('modif_sic', modif_sic,  0)
          IF(modif_sic >= 1) &
@@ -238 +243 @@
          CALL getin_p('A_satlim', A_satlim, 1.0)
          CALL getin_p('ok_restrict_A_satlim', ok_restrict_A_satlim, 0)
-#ifdef ISOVERIF
+         IF(isoCheck) THEN
          CALL msg(' 315: A_satlim='//TRIM(num2str(A_satlim)), modname, A_satlim > 1.0)
          IF(A_satlim > 1.0) STOP
-#endif
+         END IF
 !        CALL getin_p('slope_limiterxy',   slope_limiterxy,  2.0) 
 !        CALL getin_p('slope_limiterz',    slope_limiterz,   2.0) 
@@ -374 +379 @@
          CALL msg('69: h_land_ice = '//TRIM(num2str(h_land_ice)),    modname)
          CALL msg('69:      P_veg = '//TRIM(num2str(P_veg)),         modname)
+!==============================================================================================================================
+      ELSE
+!==============================================================================================================================
+         CALL abort_physic('"isotopes_mod" is not set up yet for isotopes family "'//TRIM(isoFamilies(ii))//'"', modname, 1)
+!==============================================================================================================================
+      END IF
+!==============================================================================================================================
+   END DO
 
 END SUBROUTINE iso_init
-
-
-SUBROUTINE getinp_s(nam, val, def, lDisp)
-   USE ioipsl, ONLY: getin
-   USE mod_phys_lmdz_mpi_data, ONLY :  is_mpi_root
-   USE mod_phys_lmdz_omp_data, ONLY :  is_omp_root
-   USE mod_phys_lmdz_transfert_para, ONLY : bcast
-   CHARACTER(LEN=*),           INTENT(IN)    :: nam
-   CHARACTER(LEN=*),           INTENT(INOUT) :: val
-   CHARACTER(LEN=*), OPTIONAL, INTENT(IN)    :: def
-   LOGICAL,          OPTIONAL, INTENT(IN)    :: lDisp
-   LOGICAL :: lD
-!$OMP BARRIER
-   IF(is_mpi_root.AND.is_omp_root) THEN
-      IF(PRESENT(def)) val=def; CALL getin(nam,val)
-      lD=.TRUE.; IF(PRESENT(lDisp)) lD=lDisp
-      IF(lD) CALL msg(TRIM(nam)//' = '//TRIM(val))
-  END IF
-  CALL bcast(val)
-END SUBROUTINE getinp_s
-
-SUBROUTINE getinp_i(nam, val, def, lDisp)
-   USE ioipsl, ONLY: getin
-   USE mod_phys_lmdz_mpi_data, ONLY :  is_mpi_root
-   USE mod_phys_lmdz_omp_data, ONLY :  is_omp_root
-   USE mod_phys_lmdz_transfert_para, ONLY : bcast
-   CHARACTER(LEN=*),  INTENT(IN)    :: nam
-   INTEGER,           INTENT(INOUT) :: val
-   INTEGER, OPTIONAL, INTENT(IN)    :: def
-   LOGICAL, OPTIONAL, INTENT(IN)    :: lDisp
-   LOGICAL :: lD
-!$OMP BARRIER
-   IF(is_mpi_root.AND.is_omp_root) THEN
-      IF(PRESENT(def)) val=def; CALL getin(nam,val)
-      lD=.TRUE.; IF(PRESENT(lDisp)) lD=lDisp
-      IF(lD) CALL msg(TRIM(nam)//' = '//TRIM(num2str(val)))
-  END IF
-  CALL bcast(val)
-END SUBROUTINE getinp_i
-
-SUBROUTINE getinp_r(nam, val, def, lDisp)
-   USE ioipsl, ONLY: getin
-   USE mod_phys_lmdz_mpi_data, ONLY :  is_mpi_root
-   USE mod_phys_lmdz_omp_data, ONLY :  is_omp_root
-   USE mod_phys_lmdz_transfert_para, ONLY : bcast
-   CHARACTER(LEN=*),  INTENT(IN)    :: nam
-   REAL,              INTENT(INOUT) :: val
-   REAL,    OPTIONAL, INTENT(IN)    :: def
-   LOGICAL, OPTIONAL, INTENT(IN)    :: lDisp
-   LOGICAL :: lD
-!$OMP BARRIER
-   IF(is_mpi_root.AND.is_omp_root) THEN
-      IF(PRESENT(def)) val=def; CALL getin(nam,val)
-      lD=.TRUE.; IF(PRESENT(lDisp)) lD=lDisp
-      IF(lD) CALL msg(TRIM(nam)//' = '//TRIM(num2str(val)))
-  END IF
-  CALL bcast(val)
-END SUBROUTINE getinp_r
-
-SUBROUTINE getinp_l(nam, val, def, lDisp)
-   USE ioipsl, ONLY: getin
-   USE mod_phys_lmdz_mpi_data, ONLY :  is_mpi_root
-   USE mod_phys_lmdz_omp_data, ONLY :  is_omp_root
-   USE mod_phys_lmdz_transfert_para, ONLY : bcast
-   CHARACTER(LEN=*),  INTENT(IN)    :: nam
-   LOGICAL,           INTENT(INOUT) :: val
-   LOGICAL, OPTIONAL, INTENT(IN)    :: def
-   LOGICAL, OPTIONAL, INTENT(IN)    :: lDisp
-   LOGICAL :: lD
-!$OMP BARRIER
-   IF(is_mpi_root.AND.is_omp_root) THEN
-      IF(PRESENT(def)) val=def; CALL getin(nam,val)
-      lD=.TRUE.; IF(PRESENT(lDisp)) lD=lDisp
-      IF(lD) CALL msg(TRIM(nam)//' = '//TRIM(num2str(val)))
-  END IF
-  CALL bcast(val)
-END SUBROUTINE getinp_l
 
 END MODULE isotopes_mod