Index: LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_precip.f90
===================================================================
--- LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_precip.f90	(revision 5580)
+++ LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_precip.f90	(revision 5581)
@@ -1090,7 +1090,7 @@
       IF ( (snowclr_tmp(i) .GT. 0.) .AND. (precipfracclr_tmp(i) .GT. eps) &
               .AND. (rainclr_tmp(i) .LT. eps) ) THEN
+        dsnowres = snowclr_tmp(i) * EXP( - snowclr_tmp(i) / precipfracclr_tmp(i) &
+                                           / snow_thresh_resuspension )
         IF ( dsnowres .GT. 0. ) THEN
-          dsnowres = snowclr_tmp(i) * EXP( - snowclr_tmp(i) / precipfracclr_tmp(i) &
-                                             / snow_thresh_resuspension )
           dqsres(i) = dqsres(i) - dsnowres / dhum_to_dflux(i)
           !--The following line determines the in-cloud ice water content of the newly formed cloud
@@ -1111,7 +1111,7 @@
       IF ( (snowcld_tmp(i) .GT. 0.) .AND. (precipfraccld_tmp(i) .GT. 0.) &
               .AND. (raincld_tmp(i) .LT. eps) ) THEN
+        dsnowres = snowcld_tmp(i) * EXP( - snowcld_tmp(i) / precipfraccld_tmp(i) &
+                                           / snow_thresh_resuspension )
         IF ( dsnowres .GT. 0. ) THEN
-          dsnowres = snowcld_tmp(i) * EXP( - snowcld_tmp(i) / precipfraccld_tmp(i) &
-                                             / snow_thresh_resuspension )
           dqsres(i) = dqsres(i) - dsnowres / dhum_to_dflux(i)
           !--Add tendencies