From 4005167e3a967dfd98205f51563a11f42f9df825 Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Date: Thu, 4 Jun 2020 17:44:56 +0200
Subject: [PATCH] Philippe 04/06/2020: bugfix: initialize ZCONCM/T

---
 src/MNH/ini_one_wayn.f90 | 3 ++-
 src/MNH/one_wayn.f90     | 3 ++-
 2 files changed, 4 insertions(+), 2 deletions(-)

diff --git a/src/MNH/ini_one_wayn.f90 b/src/MNH/ini_one_wayn.f90
index 1c0ad2b59..61626d6fe 100644
--- a/src/MNH/ini_one_wayn.f90
+++ b/src/MNH/ini_one_wayn.f90
@@ -140,7 +140,7 @@ SUBROUTINE INI_ONE_WAY_n(KDAD,HLUOUT,PTSTEP,KMI,KTCOUNT,                 &
 !!      J.Escobar : 18/12/2015 : Correction of bug in bound in // for NHALO <>1 
 !!      B.VIE   2016 : LIMA
 !  P. Wautelet 05/2016-04/2018: new data structures and calls for I/O
-!  P. Wautelet 04/06/2020: correct call to Set_conc_lima
+!  P. Wautelet 04/06/2020: correct call to Set_conc_lima + initialize ZCONCM
 !------------------------------------------------------------------------------
 !
 !*      0.   DECLARATIONS
@@ -402,6 +402,7 @@ ENDIF
 IF (HCLOUD=="LIMA"  ) THEN
    IF (CCLOUD/="LIMA") THEN
       ALLOCATE(ZCONCM(SIZE(XRHODJ,1),SIZE(XRHODJ,2),SIZE(XRHODJ,3),NSV_LIMA_A(KMI)))
+      ZCONCM(:, :, :, :) = 0.
       IF (CCLOUD == "REVE") THEN
          ZINIT_TYPE = "INI1"
       ELSE
diff --git a/src/MNH/one_wayn.f90 b/src/MNH/one_wayn.f90
index baa51ac83..e31ffce7b 100644
--- a/src/MNH/one_wayn.f90
+++ b/src/MNH/one_wayn.f90
@@ -187,7 +187,7 @@ SUBROUTINE ONE_WAY_n(KDAD,PTSTEP,KMI,KTCOUNT,                            &
 !!      J.Escobar : 18/12/2015 : Correction of bug in bound in // for NHALO <>1 
 !!      Modification    01/2016  (JP Pinty) Add LIMA
 !  P. Wautelet: 05/2016-04/2018: new data structures and calls for I/O
-!  P. Wautelet 04/06/2020: correct call to Set_conc_lima
+!  P. Wautelet 04/06/2020: correct call to Set_conc_lima + initialize ZCONCT
 !------------------------------------------------------------------------------
 !
 !*      0.   DECLARATIONS
@@ -483,6 +483,7 @@ ENDIF
 IF (HCLOUD=="LIMA"  ) THEN
    IF (CCLOUD/="LIMA") THEN
       ALLOCATE(ZCONCT(SIZE(XRHODJ,1),SIZE(XRHODJ,2),SIZE(XRHODJ,3),NSV_LIMA_A(KMI)))
+      ZCONCT(:, :, :, :) = 0.
       IF (CCLOUD == "REVE") THEN
          ZINIT_TYPE = "INI1"
       ELSE
-- 
GitLab