From a4584cbdb1df0bd984570c2678c2459bcf01efd4 Mon Sep 17 00:00:00 2001
From: Juan Escobar <escj@aero.obs-mip.fr>
Date: Wed, 17 Apr 2019 18:24:42 +0200
Subject: [PATCH] Juan 17/04/2019:communication.f90 , correction for
scalarprod_mnh in Transpose called by main CG
---
tensorproductmultigrid_Source/communication.f90 | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/tensorproductmultigrid_Source/communication.f90 b/tensorproductmultigrid_Source/communication.f90
index e549f3a09..2cbb3ee82 100644
--- a/tensorproductmultigrid_Source/communication.f90
+++ b/tensorproductmultigrid_Source/communication.f90
@@ -478,8 +478,8 @@ contains
if (LUseT) then
local_sumt = 0.0_rl
do iz=0,nz+1
- do iy=a%iy_min,a%iy_max
- do ix=a%ix_min,a%ix_max
+ do iy=a%icompy_min,a%icompy_max
+ do ix=a%icompx_min,a%icompx_max
local_sumt = local_sumt &
+ a%st(ix,iy,iz)*b%st(ix,iy,iz)
end do
--
GitLab