Index: LMDZ6/trunk/libf/phylmd/lmdz_lscp_poprecip.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/lmdz_lscp_poprecip.F90	(revision 4876)
+++ LMDZ6/trunk/libf/phylmd/lmdz_lscp_poprecip.F90	(revision 4878)
@@ -163,5 +163,5 @@
     !--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)