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

Timestamp:

Jul 2, 2025, 4:26:06 PM (41 hours ago)

Author:

dcugnet

Message:

Make the ioipsl_getin_p routine more flexible: optional default value
and optional flag to print or not the effective value of the key
=> get rid of the get_in routine

File:

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

LMDZ6/trunk/libf/phylmdiso/isotopes_mod.F90

@@ -6 +6 @@
    USE infotrac_phy, ONLY: isoName, niso, ntiso
    USE iso_params_mod
+   USE ioipsl_getin_p_mod, ONLY : getin_p
    IMPLICIT NONE
-   INTERFACE get_in; MODULE PROCEDURE getinp_s, getinp_i, getinp_r, getinp_l;  END INTERFACE get_in
 
   !--- Contains all isotopic variables + their initialization
@@ -187 +187 @@
 
          !--- Initialisation: reading the isotopic parameters.
-         CALL get_in('lambda',     lambda_sursat, 0.004); IF(ok_nocinsat) lambda_sursat = 0.
-         CALL get_in('thumxt1',    thumxt1,       0.75*1.2)
-         CALL get_in('ntot',       ntot,          20,  .FALSE.)
-         CALL get_in('h_land_ice', h_land_ice,    20., .FALSE.)
-         CALL get_in('P_veg',      P_veg,         1.0, .FALSE.)
-         CALL get_in('bidouille_anti_divergence', bidouille_anti_divergence, .FALSE.)
-         CALL get_in('essai_convergence',         essai_convergence,         .FALSE.)
-         CALL get_in('initialisation_iso',        initialisation_iso,        0)
+         CALL getin_p('lambda',     lambda_sursat, 0.004); IF(ok_nocinsat) lambda_sursat = 0.
+         CALL getin_p('thumxt1',    thumxt1,       0.75*1.2)
+         CALL getin_p('ntot',       ntot,          20,  lDisp=.FALSE.)
+         CALL getin_p('h_land_ice', h_land_ice,    20., lDisp=.FALSE.)
+         CALL getin_p('P_veg',      P_veg,         1.0, lDisp=.FALSE.)
+         CALL getin_p('bidouille_anti_divergence', bidouille_anti_divergence, .FALSE.)
+         CALL getin_p('essai_convergence',         essai_convergence,         .FALSE.)
+         CALL getin_p('initialisation_iso',        initialisation_iso,        0)
 
 !        IF(nzone>0 .AND. initialisation_iso==0) &
-!           CALL get_in('initialisation_isotrac',initialisation_isotrac)
-         CALL get_in('modif_sst',      modif_sst,         0)
-         CALL get_in('deltaTtest',     deltaTtest,      0.0)     !--- For modif_sst>=1
-         CALL get_in('deltaTtestpoles',deltaTtestpoles, 0.0)     !--- For modif_sst>=2
-         CALL get_in( 'sstlatcrit',    sstlatcrit,     30.0)     !--- For modif_sst>=3
-         CALL get_in('dsstlatcrit',   dsstlatcrit,      0.0)     !--- For modif_sst>=3
+!           CALL getin_p('initialisation_isotrac',initialisation_isotrac)
+         CALL getin_p('modif_sst',      modif_sst,         0)
+         CALL getin_p('deltaTtest',     deltaTtest,      0.0)     !--- For modif_sst>=1
+         CALL getin_p('deltaTtestpoles',deltaTtestpoles, 0.0)     !--- For modif_sst>=2
+         CALL getin_p( 'sstlatcrit',    sstlatcrit,     30.0)     !--- For modif_sst>=3
+         CALL getin_p('dsstlatcrit',   dsstlatcrit,      0.0)     !--- For modif_sst>=3
 #ifdef ISOVERIF
          CALL msg('iso_init 270:  sstlatcrit='//num2str( sstlatcrit), modname, sstlatcrit < 0.0) !--- For modif_sst>=2
@@ -209 +209 @@
 #endif              
 
-         CALL get_in('modif_sic', modif_sic,  0)
+         CALL getin_p('modif_sic', modif_sic,  0)
          IF(modif_sic >= 1) &
-         CALL get_in('deltasic',  deltasic, 0.1)
-
-         CALL get_in('albedo_prescrit', albedo_prescrit, 0)
+         CALL getin_p('deltasic',  deltasic, 0.1)
+
+         CALL getin_p('albedo_prescrit', albedo_prescrit, 0)
          IF(albedo_prescrit == 1) THEN
-            CALL get_in('lon_min_albedo', lon_min_albedo, -200.)
-            CALL get_in('lon_max_albedo', lon_max_albedo,  200.)
-            CALL get_in('lat_min_albedo', lat_min_albedo, -100.)
-            CALL get_in('lat_max_albedo', lat_max_albedo,  100.)
-         END IF
-         CALL get_in('deltaO18_oce',        deltaO18_oce,   0.0)
-         CALL get_in('deltaP_BL',           deltaP_BL,     10.0)
-         CALL get_in('ruissellement_pluie', ruissellement_pluie, 0)
-         CALL get_in('alphak_stewart',      alphak_stewart,      1)
-         CALL get_in('tdifexp_sol',         tdifexp_sol,      0.67)
-         CALL get_in('calendrier_guide',    calendrier_guide,    0)
-         CALL get_in('cste_surf_cond',      cste_surf_cond,      0)
-         CALL get_in('mixlen',              mixlen,           35.0)
-         CALL get_in('evap_cont_cste',      evap_cont_cste,      0)
-         CALL get_in('deltaO18_evap_cont',  deltaO18_evap_cont,0.0)
-         CALL get_in('d_evap_cont',         d_evap_cont,       0.0)
-         CALL get_in('nudge_qsol',          nudge_qsol,          0)
-         CALL get_in('region_nudge_qsol',   region_nudge_qsol,   1)
+            CALL getin_p('lon_min_albedo', lon_min_albedo, -200.)
+            CALL getin_p('lon_max_albedo', lon_max_albedo,  200.)
+            CALL getin_p('lat_min_albedo', lat_min_albedo, -100.)
+            CALL getin_p('lat_max_albedo', lat_max_albedo,  100.)
+         END IF
+         CALL getin_p('deltaO18_oce',        deltaO18_oce,   0.0)
+         CALL getin_p('deltaP_BL',           deltaP_BL,     10.0)
+         CALL getin_p('ruissellement_pluie', ruissellement_pluie, 0)
+         CALL getin_p('alphak_stewart',      alphak_stewart,      1)
+         CALL getin_p('tdifexp_sol',         tdifexp_sol,      0.67)
+         CALL getin_p('calendrier_guide',    calendrier_guide,    0)
+         CALL getin_p('cste_surf_cond',      cste_surf_cond,      0)
+         CALL getin_p('mixlen',              mixlen,           35.0)
+         CALL getin_p('evap_cont_cste',      evap_cont_cste,      0)
+         CALL getin_p('deltaO18_evap_cont',  deltaO18_evap_cont,0.0)
+         CALL getin_p('d_evap_cont',         d_evap_cont,       0.0)
+         CALL getin_p('nudge_qsol',          nudge_qsol,          0)
+         CALL getin_p('region_nudge_qsol',   region_nudge_qsol,   1)
          nlevmaxO17 = 50
          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)
+         CALL getin_p('no_pce',   no_pce,     0)
+         CALL getin_p('A_satlim', A_satlim, 1.0)
+         CALL getin_p('ok_restrict_A_satlim', ok_restrict_A_satlim, 0)
 #ifdef ISOVERIF
-         CALL msg(' 315: A_satlim='//num2str(A_satlim), modname, A_satlim > 1.0)
+         CALL msg(' 315: A_satlim='//TRIM(num2str(A_satlim)), modname, A_satlim > 1.0)
          IF(A_satlim > 1.0) STOP
 #endif
-!        CALL get_in('slope_limiterxy',   slope_limiterxy,  2.0) 
-!        CALL get_in('slope_limiterz',    slope_limiterz,   2.0) 
-         CALL get_in('modif_ratqs',       modif_ratqs,        0)
-         CALL get_in('Pcrit_ratqs',       Pcrit_ratqs,    500.0)
-         CALL get_in('ratqsbasnew',       ratqsbasnew,     0.05)
-         CALL get_in('fac_modif_evaoce',  fac_modif_evaoce, 1.0)
-         CALL get_in('ok_bidouille_wake', ok_bidouille_wake,  0)
+!        CALL getin_p('slope_limiterxy',   slope_limiterxy,  2.0) 
+!        CALL getin_p('slope_limiterz',    slope_limiterz,   2.0) 
+         CALL getin_p('modif_ratqs',       modif_ratqs,        0)
+         CALL getin_p('Pcrit_ratqs',       Pcrit_ratqs,    500.0)
+         CALL getin_p('ratqsbasnew',       ratqsbasnew,     0.05)
+         CALL getin_p('fac_modif_evaoce',  fac_modif_evaoce, 1.0)
+         CALL getin_p('ok_bidouille_wake', ok_bidouille_wake,  0)
          ! si oui, la temperature de cond est celle de l'environnement, pour eviter
          ! bugs quand temperature dans ascendances convs est mal calculee
-         CALL get_in('cond_temp_env',        cond_temp_env,        .FALSE.)
+         CALL getin_p('cond_temp_env',        cond_temp_env,        .FALSE.)
          IF(ANY(isoName == 'HTO')) &
-         CALL get_in('ok_prod_nucl_tritium', ok_prod_nucl_tritium, .FALSE., .FALSE.)
-         CALL get_in('tnateq1', ltnat1, .TRUE.)
+         CALL getin_p('ok_prod_nucl_tritium', ok_prod_nucl_tritium, .FALSE., lDisp=.FALSE.)
+         CALL getin_p('tnateq1', ltnat1, .TRUE.)
 
          CALL msg('iso_O18, iso_HDO, iso_eau = '//TRIM(strStack(num2str([iso_O18, iso_HDO, iso_eau]))), modname)