From 43775537e49295948405602eb01732320f653990 Mon Sep 17 00:00:00 2001
From: Quentin Rodier <quentin.rodier@meteo.fr>
Date: Fri, 1 Mar 2019 15:20:43 +0100
Subject: [PATCH] Quentin 01/03/2019: backup from commit 7d17f6da9: the TINY
function added on PURS lead to non-zero values to other variables such as
convective precip. and snow variables such as T%TSNOW_ROOF%WSNOW in
urban_snow_evol.F90 which lead to activate the snow cover 1layer model even
if snow is not present
---
src/MNH/convect_updraft.f90 | 3 ---
1 file changed, 3 deletions(-)
diff --git a/src/MNH/convect_updraft.f90 b/src/MNH/convect_updraft.f90
index 13b65f82a..d860b6240 100644
--- a/src/MNH/convect_updraft.f90
+++ b/src/MNH/convect_updraft.f90
@@ -144,7 +144,6 @@ END MODULE MODI_CONVECT_UPDRAFT
!! Original 07/11/95
!! Last modified 10/12/97
!! V.Masson, C.Lac, Sept. 2010 : Correction of a loop for reproducibility
-! P. Wautelet 12/02/2019: bugfix: PURR/PURS were sometimes =-0. -> problems later
!-------------------------------------------------------------------------------
!
!* 0. DECLARATIONS
@@ -380,9 +379,7 @@ DO JK = MAX( IKB + 1, JKMIN ), IKE - 1
PUTPR(:) = PUTPR(:) + PUPR(:,JKP) ! total precipitation rate
ZWORK2(:) = PURR(:,JKP) / MAX( 1.E-8, PURC(:,JKP) + PURI(:,JKP) )
PURR(:,JKP) = ZWORK2(:) * PURC(:,JKP) ! liquid precipitation
- PURR(:,JKP)=MAX(TINY(0.),PURR(:,JKP)) !To prevent problems when =-0. (occur sometimes)
PURS(:,JKP) = ZWORK2(:) * PURI(:,JKP) ! solid precipitation
- PURS(:,JKP)=MAX(TINY(0.),PURS(:,JKP)) !To prevent problems when =-0. (occur sometimes)
!
!
!* 7. Update r_c, r_i, enthalpy, r_w for precipitation
--
GitLab