diff --git a/src/PHYEX/turb/modd_turbn.f90 b/src/PHYEX/turb/modd_turbn.f90
index a3f2c6e9df4054b8b84a257cb31d1d5b27f8be5c..ef3ad4cbe1da3b08d6e7aca87a58f73b0079519e 100644
--- a/src/PHYEX/turb/modd_turbn.f90
+++ b/src/PHYEX/turb/modd_turbn.f90
@@ -122,6 +122,7 @@ REAL               :: XCEI_MAX  !< maximum threshold for the instability index C
                                      !(beginning of the saturation of the amplification)
 REAL, DIMENSION(:,:,:), POINTER  :: XCEI !< Cloud Entrainment instability index to emphasize localy 
                                          ! turbulent fluxes
+  LOGICAL            :: LTURB_PRECIP ! switch to apply turbulence to precipitating hydrometeor mixing ratios
 
 !  
 END TYPE TURB_t
@@ -180,6 +181,7 @@ REAL, POINTER :: XCOEF_AMPL_SAT=>NULL()
 REAL, POINTER :: XCEI_MIN=>NULL()
 REAL, POINTER :: XCEI_MAX =>NULL()
 REAL, DIMENSION(:,:,:), POINTER  :: XCEI=>NULL()
+LOGICAL, POINTER :: LTURB_PRECIP=>NULL()
 !
 NAMELIST/NAM_TURBn/XIMPL,CTURBLEN,CTURBDIM,LTURB_FLX,LTURB_DIAG,  &
                    LSIG_CONV,LRMC01,CTOM,&
@@ -188,7 +190,7 @@ NAMELIST/NAM_TURBn/XIMPL,CTURBLEN,CTURBDIM,LTURB_FLX,LTURB_DIAG,  &
                    XALTHGRAD, XCLDTHOLD, XLINI, LHARAT, &
                    LPROJQITURB, LSMOOTH_PRANDTL, XMINSIGS, NTURBSPLIT, &
                    LCLOUDMODIFLM, CTURBLEN_CLOUD, &
-                   XCOEF_AMPL_SAT, XCEI_MIN, XCEI_MAX
+                   XCOEF_AMPL_SAT, XCEI_MIN, XCEI_MAX, LTURB_PRECIP
 !
 !-------------------------------------------------------------------------------
 !
@@ -276,6 +278,7 @@ XCOEF_AMPL_SAT=>TURB_MODEL(KTO)%XCOEF_AMPL_SAT
 XCEI_MIN=>TURB_MODEL(KTO)%XCEI_MIN
 XCEI_MAX =>TURB_MODEL(KTO)%XCEI_MAX
 XCEI=>TURB_MODEL(KTO)%XCEI
+LTURB_PRECIP=>TURB_MODEL(KTO)%LTURB_PRECIP
 !
 ENDIF
 !
@@ -386,6 +389,7 @@ IF(LLDEFAULTVAL) THEN
   XCOEF_AMPL_SAT = 5.
   XCEI_MIN = 0.001E-06
   XCEI_MAX = 0.01E-06
+  LTURB_PRECIP =.FALSE.
   !
   IF(HPROGRAM=='AROME') THEN
     XTKEMIN=1.E-6
diff --git a/src/PHYEX/turb/turb.f90 b/src/PHYEX/turb/turb.f90
index 79eb73fe6a3d5c0d31702c58b15384470af34e5d..e333635ea5d10e4a6c6af8485e19e52de5f1681e 100644
--- a/src/PHYEX/turb/turb.f90
+++ b/src/PHYEX/turb/turb.f90
@@ -1346,11 +1346,13 @@ IF ( KRRL >= 1 ) THEN
 END IF!
 !
 PRSVS(:,:,:)=ZWORKS(:,:,1:KSV)
-IF (KRR.GE.3) PRRS(:,:,3)=ZWORKS(:,:,KSV+3)
-IF (KRR.GE.5) PRRS(:,:,5)=ZWORKS(:,:,KSV+5)
-IF (KRR.GE.6) PRRS(:,:,6)=ZWORKS(:,:,KSV+6)
-IF (KRR.GE.7) PRRS(:,:,7)=ZWORKS(:,:,KSV+7)
-IF (OFLYER)   PWSV(:,:,:)=ZWORKWSV(:,:,1:KSV)
+IF (TURBN%LTURB_PRECIP) THEN
+   IF (KRR.GE.3) PRRS(:,:,3)=ZWORKS(:,:,KSV+3)
+   IF (KRR.GE.5) PRRS(:,:,5)=ZWORKS(:,:,KSV+5)
+   IF (KRR.GE.6) PRRS(:,:,6)=ZWORKS(:,:,KSV+6)
+   IF (KRR.GE.7) PRRS(:,:,7)=ZWORKS(:,:,KSV+7)
+END IF
+   IF (OFLYER)   PWSV(:,:,:)=ZWORKWSV(:,:,1:KSV)
 !
 ! Remove non-physical negative values (unnecessary in a perfect world) + corresponding budgets
 CALL SOURCES_NEG_CORRECT_PHY(D,KSV,HCLOUD,HELEC,'NETUR',KRR,PTSTEP,PPABST,PTHLT,PRT,PRTHLS,PRRS,PRSVS)