From b104adba42769425e026f7571aa0252e8168f15c Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Date: Fri, 9 Jul 2021 09:14:24 +0200
Subject: [PATCH] Philippe 09/07/2021: bugfix: altitude levels are on the
correct grid position (mass point)
---
src/MNH/ini_lesn.f90 | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/src/MNH/ini_lesn.f90 b/src/MNH/ini_lesn.f90
index a0776d1fe..45d7b9f9a 100644
--- a/src/MNH/ini_lesn.f90
+++ b/src/MNH/ini_lesn.f90
@@ -39,6 +39,7 @@
! P. Wautelet 12/08/2020: bugfix: use NUNDEF instead of XUNDEF for integer variables
! P. Wautelet 04/01/2021: bugfix: nles_k was used instead of nspectra_k for a loop index
! P. Wautelet 30/03/2021: budgets: LES cartesian subdomain limits are defined in the physical domain
+! P. Wautelet 09/07/2021: bugfix: altitude levels are on the correct grid position (mass point)
! --------------------------------------------------------------------------
!
!* 0. DECLARATIONS
@@ -72,6 +73,7 @@ USE MODD_CONDSAMP, ONLY : LCONDSAMP
!
USE MODI_INI_LES_CART_MASKn
USE MODI_COEF_VER_INTERP_LIN
+USE MODI_SHUMAN
!
IMPLICIT NONE
!
@@ -365,7 +367,7 @@ IF (ANY(XLES_ALTITUDES(:)/=XUNDEF)) THEN
END DO
END DO
END DO
- CALL COEF_VER_INTERP_LIN(XZZ,ZZ_LES,NKLIN_LES,XCOEFLIN_LES)
+ CALL COEF_VER_INTERP_LIN(MZF(XZZ),ZZ_LES,NKLIN_LES,XCOEFLIN_LES)
!
DEALLOCATE(ZZ_LES)
END IF
--
GitLab