Index: LMDZ4/branches/LMDZ4_AR5/libf/phylmd/clesphys.h
===================================================================
--- LMDZ4/branches/LMDZ4_AR5/libf/phylmd/clesphys.h	(revision 1483)
+++ LMDZ4/branches/LMDZ4_AR5/libf/phylmd/clesphys.h	(revision 1500)
@@ -69,4 +69,5 @@
        LOGICAL :: ok_strato
        LOGICAL :: ok_hines
+       LOGICAL :: ok_precip_fonte
 
        COMMON/clesphys/cycle_diurne, soil_model, new_oliq,              &
@@ -88,5 +89,5 @@
      &     , ok_lic_melt, cvl_corr, aer_type                            &
      &     , qsol0, iflag_rrtm, ok_strato,ok_hines,ecrit_LES            &
-     &     , co2_ppm0
+     &     , co2_ppm0, ok_precip_fonte
      
 !$OMP THREADPRIVATE(/clesphys/)
Index: LMDZ4/branches/LMDZ4_AR5/libf/phylmd/conf_phys.F90
===================================================================
--- LMDZ4/branches/LMDZ4_AR5/libf/phylmd/conf_phys.F90	(revision 1483)
+++ LMDZ4/branches/LMDZ4_AR5/libf/phylmd/conf_phys.F90	(revision 1500)
@@ -160,4 +160,5 @@
   LOGICAL,SAVE      :: carbon_cycle_tr_omp
   LOGICAL,SAVE      :: carbon_cycle_cpl_omp
+  LOGICAL,SAVE :: ok_precip_fonte_omp
 
   integer, intent(out):: read_climoz ! read ozone climatology, OpenMP shared
@@ -1180,4 +1181,14 @@
   ok_hfCOSP_omp = .false.
   call getin('ok_hfCOSP',ok_hfCOSP_omp)
+
+!
+!Config Key  = ok_hfCOSP
+!Config Desc =
+!Config Def  = .false.
+!Config Help =
+!
+  ok_precip_fonte_omp = .false.
+  call getin('ok_precip_fonte',ok_precip_fonte_omp)
+
 
 !
@@ -1538,4 +1549,5 @@
     ok_journeCOSP = ok_journeCOSP_omp
     ok_hfCOSP = ok_hfCOSP_omp
+    ok_precip_fonte = ok_precip_fonte_omp
     seuil_inversion=seuil_inversion_omp
     lonmin_ins = lonmin_ins_omp
@@ -1691,4 +1703,5 @@
   write(numout,*)' ok_journeCOSP = ',ok_journeCOSP
   write(numout,*)' ok_hfCOSP =',ok_hfCOSP
+  write(numout,*)' ok_precip_fonte =',ok_precip_fonte
   write(numout,*)' solarlong0 = ', solarlong0
   write(numout,*)' qsol0 = ', qsol0
Index: LMDZ4/branches/LMDZ4_AR5/libf/phylmd/fonte_neige_mod.F90
===================================================================
--- LMDZ4/branches/LMDZ4_AR5/libf/phylmd/fonte_neige_mod.F90	(revision 1483)
+++ LMDZ4/branches/LMDZ4_AR5/libf/phylmd/fonte_neige_mod.F90	(revision 1500)
@@ -267,5 +267,9 @@
                (1. - coeff_rel) * run_off_lic_0(j)
           run_off_lic_0(j) = run_off_lic(i)
-          run_off_lic(i)   = run_off_lic(i) + fqfonte(i)/dtime
+          if (ok_precip_fonte) then
+            run_off_lic(i)   = run_off_lic(i) + fqfonte(i) + precip_rain(i)
+          else
+            run_off_lic(i)   = run_off_lic(i) + fqfonte(i)/dtime
+          endif
        END DO
     ENDIF
