Changeset 6158

Timestamp:

Apr 3, 2026, 5:56:32 PM (7 hours ago)

Author:

Sebastien Nguyen

Message:

modifications of isotopes_mod.F90 for snow surface isotopic exchanges

File:

• LMDZ6/branches/LMDZISOdevSN/libf/phylmdiso/isotopes_mod.F90 (modified) (2 diffs)

LMDZ6/branches/LMDZISOdevSN/libf/phylmdiso/isotopes_mod.F90

@@ -14 +14 @@
    !--- Negligible lower thresholds: no need to check for absurd values under these lower limits
    REAL, PARAMETER :: &
-      ridicule      = 1e-12,              & ! For mixing ratios
-      ridicule_rain = 1e-8,               & ! For rain fluxes (rain, zrfl...) in kg/s <-> 1e-3 mm/day
-      ridicule_evap = ridicule_rain*1e-2, & ! For evaporations                in kg/s <-> 1e-3 mm/day
-      ridicule_qsol = ridicule_rain,      & ! For qsol                        in kg <-> 1e-8 kg
-      ridicule_snow = ridicule_qsol         ! For snow                        in kg <-> 1e-8 kg
+      ridicule      = 1e-12,              & ! For mixing ratios               in kg/kg
+      ridicule_rain = 1e-8,               & ! For rain fluxes (rain, zrfl...) in kg/m2/s <-> 1e-3 mm/day
+      ridicule_evap = ridicule_rain*1e-2, & ! For evaporations                in kg/m2/s <-> 1e-3 mm/day
+      ridicule_qsol = 1e-8,               & ! For qsol                        in kg/m2 <-> 1e-8 kg/m2 = 1e-8 mm
+      ridicule_snow = ridicule_qsol         ! For snow                        in kg/m2 <-> 1e-8 kg/m2 = 1e-8 mm
    REAL, PARAMETER :: expb_max = 30.0
 
@@ -222 +222 @@
          !--- Initialisation: reading the isotopic parameters.
          CALL getin_p('lambda',     lambda_sursat, 0.004); IF(ok_nocinsat) lambda_sursat = 0.
+         CALL get_in('surf_snow_corr_d18O', surf_snow_corr_d18O, 0.0)
+         CALL get_in('surf_snow_corr_dD', surf_snow_corr_dD, 0.0)
+         CALL get_in('snow_max', snow_max, 3000.)
          CALL getin_p('thumxt1',    thumxt1,       0.75*1.2)
          CALL getin_p('ntot',       ntot,          20,  lDisp=.FALSE.)
+         ! mass_surfsnow : default = 3 kg m-2 ~ 1 cm of snow at density 350 kg m-3 (Dome C)
+         CALL get_in('mass_surfsnow', mass_surfsnow, 3.0)
+         ! min_surfsnow : default = 0.03 kg m-2 ~ 0.1 mm of snow at density 350 kg m-3 (Dome C)
+         CALL get_in('min_surfsnow', min_surfsnow, 0.03)
          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 get_in('ok_turbfrac_sublim', ok_turbfrac_sublim, .FALSE.)
+         CALL get_in('ok_turbfrac_condens', ok_turbfrac_condens, .FALSE.)
+         CALL get_in('ok_evap_surfsnow', ok_evap_surfsnow, .FALSE.)
+         CALL get_in('ddelta_evap_max', ddelta_evap_max, 0.1)
+         CALL get_in('h1m1_min_abs', h1m1_min_abs, 0.05)
+
          CALL getin_p('essai_convergence',         essai_convergence,         .FALSE.)
          CALL getin_p('initialisation_iso',        initialisation_iso,        0)