From 99ade97f5ddc60226271bbfded1e8c375eacc2ba Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Date: Wed, 20 Jan 2021 11:10:42 +0100
Subject: [PATCH] Philippe 20/01/2021: budgets: bugfix: ztrajz was not reset
 when necessary

---
 src/MNH/mode_les_diachro.f90 | 10 +++++-----
 1 file changed, 5 insertions(+), 5 deletions(-)

diff --git a/src/MNH/mode_les_diachro.f90 b/src/MNH/mode_les_diachro.f90
index 9031c723c..d694de82c 100644
--- a/src/MNH/mode_les_diachro.f90
+++ b/src/MNH/mode_les_diachro.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2020 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2021 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
 !MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
@@ -325,7 +325,6 @@ REAL    :: ZMAX_NORM_M
 !* normalization height (usually maximum BL height)
 !
 IF (OAVG) THEN
-
   ITEMP_MEAN_START = COUNT( xles_times(:)<=XLES_TEMP_MEAN_START ) + 1
   ITEMP_MEAN_END   = COUNT( xles_times(:)<=XLES_TEMP_MEAN_END   )
 
@@ -912,9 +911,6 @@ ikh = nles_levels(iles_k)
 
 ztrajx(:, :, :) = ( iil + iih ) / 2
 ztrajy(:, :, :) = ( ijl + ijh ) / 2
-do jk = 1, iles_k
-  ztrajz(jk, :, :) = xles_current_z(jk)
-end do
 
 if ( Present( hsuffixes ) ) then
   if ( Size( hsuffixes ) /= Size( pfield, 3) ) &
@@ -959,6 +955,10 @@ Allocate( tzdates( nles_current_times ) )
 
 tzdates(:) = tles_dates(:)
 
+do jk = 1, iles_k
+  ztrajz(jk, :, :) = xles_current_z(jk)
+end do
+
 !Copy all fields from tpfield
 tzfields(:) = tpfield
 
-- 
GitLab