Philippe 05/04/2022: OpenACC: ZSOLVER: modification to get correct result with ZRESI

src/ZSOLVER/gdiv.f90

@@ -181,7 +181,7 @@ CALL MNH_MEM_GET( ZWC, IIU, IJU, IKU )
 CALL MNH_MEM_GET( Z1,  IIU, IJU, IKU )
 CALL MNH_MEM_GET( Z2,  IIU, IJU, IKU )
 CALL MNH_MEM_GET( Z3,  IIU, IJU, IKU )
-!$acc data present(ZUC,ZVC,ZWC,Z1,Z2,Z3)
+!$acc data present( ZUC, ZVC, ZWC, Z1, Z2, Z3 ) copyout( PGDIV )
 #endif
 !
 GWEST  = ( HLBCX(1) /= 'CYCL' .AND. LWEST_ll() )

src/ZSOLVER/pressurez.f90

@@ -719,11 +719,13 @@ ELSE
     A_K,B_K,C_K,D_K)
 !
   CASE('ZRESI')     ! Conjugate Residual method
+!$acc update host( ZDV_SOURCE, ZPHIT, ZTHETAV )
     CALL CONRESOLZ(HLBCX,HLBCY,PDXX,PDYY,PDZX,PDZY,PDZZ,PRHODJ,ZTHETAV,       &
     PDXHATM,PDYHATM,PRHOT,PAF,PBF,PCF,PTRIGSX,PTRIGSY,                       &
     KIFAXX,KIFAXY,KITR,ZDV_SOURCE,ZPHIT,                                     &
     PBFB,&
     PBF_SXP2_YP1_Z) !JUAN Z_SPLITING
+!$acc update device( ZPHIT )
   END SELECT
 END IF
 !

src/ZSOLVER/qlap.f90

@@ -310,6 +310,7 @@ CALL ADD3DFIELD_ll( TZFIELDS_ll, ZW, 'QLAP::ZW' )
 CALL UPDATE_HALO_ll(TZFIELDS_ll,IINFO_ll)
 CALL CLEANLIST_ll(TZFIELDS_ll)
 !
+!$acc update device( ZU, ZV, ZW )
 CALL GDIV(HLBCX,HLBCY,PDXX,PDYY,PDZX,PDZY,PDZZ,ZU,ZV,ZW,PQLAP)
 !
 IF (LIBM)  THEN