Changeset 6053 for LMDZ6/trunk/libf/phylmd/conf_phys_m.f90

Timestamp:

Jan 23, 2026, 1:11:53 PM (2 days ago)

Author:

evignon

Message:

corrections dans la routine simplehydrol et ajout de diagnostiques
pour quantifier les flux d'eau liés à la fonte de la neige sur les surfaces
englacees

File:

• LMDZ6/trunk/libf/phylmd/conf_phys_m.f90 (modified) (7 diffs)

LMDZ6/trunk/libf/phylmd/conf_phys_m.f90

@@ -176 +176 @@
     REAL,SAVE :: inertie_sol_omp,inertie_sno_omp,inertie_sic_omp
     REAL,SAVE :: inertie_lic_omp
+    REAL,SAVE :: chasno_tun_omp
+    LOGICAL,SAVE :: forc_ts_melt_omp
     REAL,SAVE :: qsol0_omp
     REAL,SAVE :: evap0_omp
@@ -1623 +1625 @@
 
 
+    !Config Key = forc_ts_melt
+    !Config Desc = force surface temperature when snow melts
+    !Config Def  = .TRUE.
+    !Config Help = set to .FALSE. to avoid unrealistically forcing ts=273.15 during snow melt
+    forc_ts_melt = .TRUE.
+    CALL getin('forc_ts_melt', forc_ts_melt)
+    
+
     !Config Key = ok_lic_cond
     !Config Desc = Prise en compte depot de vapeur d'eau sur la calotte dans le bilan d'eau
@@ -2060 +2070 @@
     z0min_omp = 0.000015
     CALL getin('z0min',z0min_omp)
-
+ 
+    ! PARAMETERS FOR SNOW AND ICE MELTING
+    chasno_tun_omp=0.15
+    CALL getin('chasno_tun',chasno_tun_omp)
 
     ! PARAMETERS FOR CONVECTIVE INHIBITION BY TROPOS. DRYNESS
@@ -2640 +2653 @@
     cvl_sig2feed = cvl_sig2feed_omp
     cvl_corr = cvl_corr_omp
+    forc_ts_melt = forc_ts_melt_omp
     ok_lic_melt = ok_lic_melt_omp
     ok_lic_cond = ok_lic_cond_omp
@@ -2658 +2672 @@
     ratio_z0hz0m_landice=ratio_z0hz0m_landice_omp
 
+    chasno_tun=chasno_tun_omp
     f_rugoro=f_rugoro_omp
 
@@ -2882 +2897 @@
     WRITE(lunout,*) ' ok_lic_melt=', ok_lic_melt
     WRITE(lunout,*) ' ok_lic_cond=', ok_lic_cond
+    WRITE(lunout,*) ' forc_ts_melt=', forc_ts_melt
     WRITE(lunout,*) ' iflag_cycle_diurne=',iflag_cycle_diurne
     WRITE(lunout,*) ' soil_model=',soil_model
@@ -2978 +2994 @@
     WRITE(lunout,*) ' iflag_sic = ', iflag_sic
     WRITE(lunout,*) ' iflag_inertie = ', iflag_inertie
+    WRITE(lunout,*) ' chasno_tun = ', chasno_tun
     WRITE(lunout,*) ' inertie_sol = ', inertie_sol
     WRITE(lunout,*) ' inertie_sic = ', inertie_sic