From 3e67238099faf9d583a2af24c74a3aa0e893d883 Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Date: Tue, 16 May 2023 16:12:28 +0200
Subject: [PATCH] Philippe 16/05/2023: bugfix: flat_inv: wrong array was passed
to FFT991
---
src/MNH/flat_inv.f90 | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/src/MNH/flat_inv.f90 b/src/MNH/flat_inv.f90
index e1408a5b4..ab0bf1965 100644
--- a/src/MNH/flat_inv.f90
+++ b/src/MNH/flat_inv.f90
@@ -496,7 +496,7 @@ IF (.NOT. L2D) THEN
IF (HLBCY(1) == 'CYCL') THEN
ZBAND1D( 1 : SIZE(ZBAND_YRT,1)*SIZE(ZBAND_YRT,2)*SIZE(ZBAND_YRT(:,:,IKB-1:),3) ) => ZBAND_YRT(:,:,IKB-1:)
ZWORK1D( 1 : SIZE(ZWORKY) ) => ZWORKY(:,:,:)
- CALL FFT991( ZBAND1D, ZBAND_YRT, PTRIGSY, KIFAXY, INC2Y, IJMAX,ILOTY, +1, &
+ CALL FFT991( ZBAND1D, ZWORK1D, PTRIGSY, KIFAXY, INC2Y, IJMAX,ILOTY, +1, &
SIZE(ZBAND1D), SIZE(ZWORK1D), SIZE(PTRIGSY) )
ELSE
CALL FFT55( ZBAND_YRT(1:,1:,IKB-1:), ZWORKY, PTRIGSY, KIFAXY, INC2Y, IJMAX, ILOTY, +1 )
--
GitLab