From 9714ce05333445322284434e2de9051ff3d7997b Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?S=C3=A9bastien=20Riette?= <sebastien.riette@meteo.fr>
Date: Mon, 13 Jun 2022 15:59:25 +0200
Subject: [PATCH] S. Riette 13/6/2022 Bugfix for unphysical points
---
src/common/micro/mode_ice4_nucleation.F90 | 6 +++++-
1 file changed, 5 insertions(+), 1 deletion(-)
diff --git a/src/common/micro/mode_ice4_nucleation.F90 b/src/common/micro/mode_ice4_nucleation.F90
index 1ba3335b9..8cf213442 100644
--- a/src/common/micro/mode_ice4_nucleation.F90
+++ b/src/common/micro/mode_ice4_nucleation.F90
@@ -69,7 +69,11 @@ INTEGER :: JI
IF (LHOOK) CALL DR_HOOK('ICE4_NUCLEATION', 0, ZHOOK_HANDLE)!
!
!$mnh_expand_where(JI=1:KSIZE)
-GNEGT(:)=PT(:)<CST%XTT .AND. PRVT(:)>ICED%XRTMIN(1) .AND. ODCOMPUTE(:)
+WHERE(ODCOMPUTE(:))
+ GNEGT(:)=PT(:)<CST%XTT .AND. PRVT(:)>ICED%XRTMIN(1)
+ELSEWHERE
+ GNEGT(:)=.FALSE.
+ENDWHERE
!$mnh_end_expand_where(JI=1:KSIZE)
ZUSW(:)=0.
--
GitLab