From 61ea3847e6be7b1fd02308a0d782107efaad95d7 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Beno=C3=AEt=20Vi=C3=A9?= <benoit.vie@meteo.fr>
Date: Fri, 18 Nov 2022 09:50:48 +0100
Subject: [PATCH] bugfix radiations for cases with any of NMOM_C, _R, _I = 1
---
src/MNH/radiations.f90 | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/src/MNH/radiations.f90 b/src/MNH/radiations.f90
index ecf0f36fe..67dba74ce 100644
--- a/src/MNH/radiations.f90
+++ b/src/MNH/radiations.f90
@@ -766,9 +766,9 @@ IF( CCLOUD == 'LIMA' ) THEN
DO JJ=IJB,IJE
DO JI=IIB,IIE
IIJ = 1 + (JI-IIB) + (IIE-IIB+1)*(JJ-IJB)
- IF (LWARM) ZCCT_LIMA(IIJ,JKRAD) = MAX(0.,PSVT(JI,JJ,JK,NSV_LIMA_NC))
- IF (LWARM .AND. LRAIN) ZCRT_LIMA(IIJ,JKRAD) = MAX(0.,PSVT(JI,JJ,JK,NSV_LIMA_NR))
- IF (LCOLD) ZCIT_LIMA(IIJ,JKRAD) = MAX(0.,PSVT(JI,JJ,JK,NSV_LIMA_NI))
+ IF (LWARM .AND. NSV_LIMA_NC.GT.0) ZCCT_LIMA(IIJ,JKRAD) = MAX(0.,PSVT(JI,JJ,JK,NSV_LIMA_NC))
+ IF (LWARM .AND. LRAIN .AND. NSV_LIMA_NR.GT.0) ZCRT_LIMA(IIJ,JKRAD) = MAX(0.,PSVT(JI,JJ,JK,NSV_LIMA_NR))
+ IF (LCOLD .AND. NSV_LIMA_NI.GT.0) ZCIT_LIMA(IIJ,JKRAD) = MAX(0.,PSVT(JI,JJ,JK,NSV_LIMA_NI))
END DO
END DO
END DO
--
GitLab