From 260280dbb5ae05725d2ab2d47975f05f79e65156 Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Date: Tue, 5 Jan 2021 14:19:41 +0100
Subject: [PATCH] Philippe 05/01/2021: bugfix: CSPECTRA_LEVEL_TYPE='Z'
computation was wrong
(cherry picked from commit 025cea89e874ea1246be65d942451fcaaabe7ee0)
---
src/MNH/spec_ver_int.f90 | 11 ++++++-----
1 file changed, 6 insertions(+), 5 deletions(-)
diff --git a/src/MNH/spec_ver_int.f90 b/src/MNH/spec_ver_int.f90
index abd3ed1e0..8cfc2223a 100644
--- a/src/MNH/spec_ver_int.f90
+++ b/src/MNH/spec_ver_int.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2000-2018 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2000-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.
@@ -50,8 +50,9 @@ END MODULE MODI_SPEC_VER_INT
!! -------------
!! Original 07/02/00
!! Philippe Wautelet: 05/2016-04/2018: new data structures and calls for I/O
-!!
-!! --------------------------------------------------------------------------
+! P. Wautelet 10/04/2019: replace ABORT and STOP calls by Print_msg
+! P. Wautelet 05/01/2021: bugfix: CSPECTRA_LEVEL_TYPE='Z' computation was wrong
+! --------------------------------------------------------------------------
!
!* 0. DECLARATIONS
! ------------
@@ -96,10 +97,10 @@ ELSE IF (CSPECTRA_LEVEL_TYPE=='K') THEN
ZA(:,:,JK) = PA_MNH(:,:,NSPECTRA_LEVELS(JK))
END DO
ELSE IF (CSPECTRA_LEVEL_TYPE=='Z') THEN
- PA_SPEC = VER_INTERP_LIN(PA_MNH,NKLIN_CURRENT_SPEC,XCOEFLIN_CURRENT_SPEC)
+ ZA(:,:,:) = VER_INTERP_LIN(PA_MNH,NKLIN_CURRENT_SPEC,XCOEFLIN_CURRENT_SPEC)
!
WHERE(NKLIN_CURRENT_SPEC<2)
- PA_SPEC = XUNDEF
+ ZA(:,:,:) = XUNDEF
END WHERE
ELSE
PRINT*, '-------> STOP in SPEC_VER_INT <----------'
--
GitLab