From 6932f716865e406f5982b0dc305e64b71ae14986 Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Date: Mon, 30 Jan 2023 13:54:27 +0100
Subject: [PATCH] Philippe 30/01/2023: bugfix: set PCIT vertical values only
for ICEx
---
src/MNH/resolved_cloud.f90 | 9 +++++----
1 file changed, 5 insertions(+), 4 deletions(-)
diff --git a/src/MNH/resolved_cloud.f90 b/src/MNH/resolved_cloud.f90
index c2ec86f73..275735ad4 100644
--- a/src/MNH/resolved_cloud.f90
+++ b/src/MNH/resolved_cloud.f90
@@ -1084,12 +1084,13 @@ SELECT CASE ( HCLOUD )
ENDIF
!
END SELECT
-! TODO: code a generic routine to update vertical lower and upper levels to 0, a
-! specific value or to IKB or IKE and apply it to every output prognostic variable of physics
-PCIT(:,:,1)=0.
-PCIT(:,:,IKE+1)=0.
!
IF(HCLOUD=='ICE3' .OR. HCLOUD=='ICE4' ) THEN
+! TODO: code a generic routine to update vertical lower and upper levels to 0, a
+! specific value or to IKB or IKE and apply it to every output prognostic variable of physics
+ PCIT(:,:,1) = 0.
+ PCIT(:,:,IKE+1) = 0.
+
PINPRC3D=ZFPR(:,:,:,2) / CST%XRHOLW
PINPRR3D=ZFPR(:,:,:,3) / CST%XRHOLW
PINPRS3D=ZFPR(:,:,:,5) / CST%XRHOLW
--
GitLab