From 1d76957ef933f012484435f5f7b9a42ce6d84d80 Mon Sep 17 00:00:00 2001
From: Christine Lac <christine.lac@meteo.fr>
Date: Fri, 12 May 2017 14:14:59 +0200
Subject: [PATCH] C.Lac: Correction on droplet deposition
---
src/MNH/rain_c2r2_khko.f90 | 14 ++++++++++----
1 file changed, 10 insertions(+), 4 deletions(-)
diff --git a/src/MNH/rain_c2r2_khko.f90 b/src/MNH/rain_c2r2_khko.f90
index 834f410cb..77ed5a590 100644
--- a/src/MNH/rain_c2r2_khko.f90
+++ b/src/MNH/rain_c2r2_khko.f90
@@ -218,6 +218,7 @@ END MODULE MODI_RAIN_C2R2_KHKO
!! activation by cooling (OACTIT : mis en commentaires)
!! M.Mazoyer : 04/2016 : Add supersaturation diagnostics
!! C.Lac : 07/2016 : Add droplet deposition
+!! C.Lac : 01/2017 : Correction on droplet deposition
!-------------------------------------------------------------------------------
!
!* 0. DECLARATIONS
@@ -1741,6 +1742,9 @@ INTEGER :: JL ! and PACK intrinsics
! optimization by looking for locations where
! the precipitating fields are larger than a minimal value only !!!
!
+IF (OSEDC) PINPRC (:,:) = 0.
+IF (LDEPOC) PINDEP (:,:) = 0.
+!
DO JN = 1 , KSPLITR
GSEDIM(:,:,:) = .FALSE.
IF( OSEDC ) THEN
@@ -1933,10 +1937,12 @@ IF (LDEPOC) THEN
GDEP(IIB:IIE,IJB:IJE) = PRCS(IIB:IIE,IJB:IJE,2) >0 .AND. &
PCCS(IIB:IIE,IJB:IJE,2) >0
WHERE (GDEP)
- PRCS(:,:,2) = PRCS(:,:,2) - XVDEPOC * PRCT(:,:,2) * PRHODJ(:,:,2)
- PCCS(:,:,2) = PCCS(:,:,2) - XVDEPOC * PCCT(:,:,2) * PRHODJ(:,:,2)
- PINPRC(:,:) = PINPRC(:,:) + XVDEPOC * PRCT(:,:,2) * PRHODJ(:,:,2) /XRHOLW
- PINDEP(:,:) = XVDEPOC * PRCT(:,:,2) * PRHODJ(:,:,2) /XRHOLW
+ PRCS(:,:,2) = PRCS(:,:,2) - XVDEPOC * PRCS(:,:,2)*PTSTEP / &
+ ( PZZ(:,:,3) - PZZ(:,:,2))
+ PCCS(:,:,2) = PCCS(:,:,2) - XVDEPOC * PCCS(:,:,2)*PTSTEP / &
+ ( PZZ(:,:,3) - PZZ(:,:,2))
+ PINPRC(:,:) = PINPRC(:,:) + XVDEPOC * PRCS(:,:,2) * PTSTEP * PRHODREF(:,:,2) /XRHOLW
+ PINDEP(:,:) = XVDEPOC * PRCS(:,:,2) * PTSTEP * PRHODREF(:,:,2) /XRHOLW
END WHERE
END IF
!
--
GitLab