Index: /LMDZ6/branches/LMDZISOdevSN/libf/phylmdiso/isotopes_mod.F90
===================================================================
--- /LMDZ6/branches/LMDZISOdevSN/libf/phylmdiso/isotopes_mod.F90	(revision 6157)
+++ /LMDZ6/branches/LMDZISOdevSN/libf/phylmdiso/isotopes_mod.F90	(revision 6158)
@@ -14,9 +14,9 @@
    !--- 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,9 +222,22 @@
          !--- 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)