From 722c7f987b2907234f21671b84e5fa1108c7d3fb Mon Sep 17 00:00:00 2001
From: Juan Escobar <escj@aero.obs-mip.fr>
Date: Tue, 27 Sep 2022 15:58:53 +0200
Subject: [PATCH] Juan 27/09/2022:prandtl.f90, Bug, missing (:,:,:) in array
explansion for no reprod side
---
src/MNH/prandtl.f90 | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/MNH/prandtl.f90 b/src/MNH/prandtl.f90
index def839a36..892901e71 100644
--- a/src/MNH/prandtl.f90
+++ b/src/MNH/prandtl.f90
@@ -985,7 +985,7 @@ ELSE IF (L2D) THEN ! 3D case in a 2D model
!$acc kernels
!$mnh_expand_array(JI=1:JIU,JJ=1:JJU,JK=1:JKU)
#if !defined(MNH_BITREP) && !defined(MNH_BITREP_OMP)
- ZTMP1_DEVICE(:,:,:) = (XG / PTHVREF * PLM * PLEPS / PTKEM)**2
+ ZTMP1_DEVICE(:,:,:) = (XG / PTHVREF(:,:,:) * PLM(:,:,:) * PLEPS(:,:,:) / PTKEM(:,:,:))**2
#else
ZTMP1_DEVICE(:,:,:) = BR_P2(XG / PTHVREF(:,:,:) * PLM(:,:,:) * PLEPS(:,:,:) / PTKEM(:,:,:))
#endif
@@ -1066,7 +1066,7 @@ ELSE ! 3D case in a 3D model
!$acc kernels
!$mnh_expand_array(JI=1:JIU,JJ=1:JJU,JK=1:JKU)
#if !defined(MNH_BITREP) && !defined(MNH_BITREP_OMP)
- ZTMP1_DEVICE(:,:,:) = (XG / PTHVREF * PLM * PLEPS / PTKEM)**2
+ ZTMP1_DEVICE(:,:,:) = (XG / PTHVREF(:,:,:) * PLM(:,:,:) * PLEPS(:,:,:) / PTKEM(:,:,:))**2
#else
ZTMP1_DEVICE(:,:,:) = BR_P2(XG / PTHVREF(:,:,:) * PLM(:,:,:) * PLEPS(:,:,:) / PTKEM(:,:,:))
#endif
--
GitLab