Changeset 4878

Timestamp:

Mar 27, 2024, 5:50:13 PM (6 weeks ago)

Author:

evignon

Message:

correction de l'humidite a saturation pour la sublimation
dans poprecip
Lea Raillard

File:

• LMDZ6/trunk/libf/phylmd/lmdz_lscp_poprecip.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/lmdz_lscp_poprecip.F90

@@ -163 +163 @@
     !--which does not need a barrier on snowclr, because included in the formula
     dsnowsub = precipfracclr(i) * ( MAX(0., &
-             - coef_sub * ( 1. - expo_sub ) * (1. - qvap(i) / qsatl(i)) * dz(i) &
+             - coef_sub * ( 1. - expo_sub ) * (1. - qvap(i) / qsati(i)) * dz(i) &
              + ( snowclr(i) / MAX(thresh_precip_frac, precipfracclr(i)) )**( 1. - expo_sub ) &
              ) )**( 1. / ( 1. - expo_sub ) ) - snowclr(i)