From 0094f5636f57e107576649c49d94b3571e479dc8 Mon Sep 17 00:00:00 2001
From: ESCOBAR Juan <escj@nuwa>
Date: Thu, 10 Mar 2022 15:55:18 +0100
Subject: [PATCH] Juan 10/03/2022:ZSOLVER/zsolver_inv.f90, add mppdb_check
---
src/ZSOLVER/zsolver_inv.f90 | 32 ++++++++++++++++++++++++++++++--
1 file changed, 30 insertions(+), 2 deletions(-)
diff --git a/src/ZSOLVER/zsolver_inv.f90 b/src/ZSOLVER/zsolver_inv.f90
index 588cd940c..48dc2b7e1 100644
--- a/src/ZSOLVER/zsolver_inv.f90
+++ b/src/ZSOLVER/zsolver_inv.f90
@@ -148,6 +148,8 @@ SUBROUTINE ZSOLVER_INV(HLBCX,HLBCY,PDXHATM,PDYHATM,PRHOM,PAF,PBF,PCF, &
#ifdef MNH_OPENACC
USE MODE_MNH_ZWORK, ONLY: MNH_MEM_GET, MNH_MEM_POSITION_PIN, MNH_MEM_RELEASE
#endif
+ !
+ USE MODE_MPPDB
!
IMPLICIT NONE
!
@@ -206,6 +208,24 @@ SUBROUTINE ZSOLVER_INV(HLBCX,HLBCY,PDXHATM,PDYHATM,PRHOM,PAF,PBF,PCF, &
!* 1. COMPUTE LOOP BOUNDS
! -------------------
!
+ if ( mppdb_initialized ) then
+ !Check all in arrays
+ call Mppdb_check( PRHOM, "Zsolver_inv beg:PRHOM" )
+ call Mppdb_check( PY, "Zsolver_inv beg:PY" )
+ call Mppdb_check( PF_1_Y, "Zsolver_inv beg:PF_1_Y" )
+ call Mppdb_check( PAF_ZS, "Zsolver_inv beg:PAF_ZS" )
+ call Mppdb_check( PBF_ZS, "Zsolver_inv beg:PBF_ZS" )
+ call Mppdb_check( PCF_ZS, "Zsolver_inv beg:PCF_ZS" )
+ call Mppdb_check( PDXATH_ZS,"Zsolver_inv beg:PDXATH_ZS" )
+ call Mppdb_check( PDYATH_ZS,"Zsolver_inv beg:PDYATH_ZS" )
+ call Mppdb_check( PRHO_ZS, "Zsolver_inv beg:PRHO_ZS" )
+ call Mppdb_check( PBFB, "Zsolver_inv beg:PBFB" )
+ call Mppdb_check( A_K, "Zsolver_inv beg:A_K" )
+ call Mppdb_check( B_K, "Zsolver_inv beg:B_K" )
+ call Mppdb_check( C_K, "Zsolver_inv beg:C_K" )
+ call Mppdb_check( D_K, "Zsolver_inv beg: D_K" )
+ end if
+ !
!CALL GET_PHYSICAL_ll(IIB,IJB,IIE,IJE)
CALL GET_INDICE_ll(IIB,IJB,IIE,IJE)
CALL GET_DIM_EXT_ll('B',IIU,IJU)
@@ -297,7 +317,11 @@ SUBROUTINE ZSOLVER_INV(HLBCX,HLBCY,PDXHATM,PDYHATM,PRHOM,PAF,PBF,PCF, &
IF (L2D) THEN
STOP ' NO 2D MULTIGRID YET '
END IF
-
+ if ( mppdb_initialized ) then
+ !Check all in arrays
+ call Mppdb_check( ZY, "Zsolver_inv before MG :ZY" )
+ call Mppdb_check( PF_1_Y, "Zsolver_inv before MG :PF_1_Y" )
+ end if
call mg_main_initialise_rhs_mnh(IIB,IIE,IIU,IJB,IJE,IJU,IKU,ZY)
call mg_main_initialise_u_mnh(IIB,IIE,IIU,IJB,IJE,IJU,IKU,PF_1_Y)
@@ -314,7 +338,11 @@ SUBROUTINE ZSOLVER_INV(HLBCX,HLBCY,PDXHATM,PDYHATM,PRHOM,PAF,PBF,PCF, &
CALL GET_HALO_DDC(PF_1_Y)
#endif
!
- CALL PF_1_Y_BOUND(PF_1_Y)
+ CALL PF_1_Y_BOUND(PF_1_Y)
+ if ( mppdb_initialized ) then
+ !Check all in arrays
+ call Mppdb_check( PF_1_Y, "Zsolver_inv end:PF_1_Y" )
+ end if
!-------------------------------------------------------------------------------
#ifndef MNH_OPENACC
DEALLOCATE(ZY)
--
GitLab