Philippe 15/01/2024: OpenACC: anel_balance_n: exchange data correctly for pressure solver

src/MNH/anel_balancen.f90

@@ -317,6 +317,9 @@ IRRL = 0             !      |                                              |
 IRRI = 0             !       ==============================================
 !
 IMI = GET_CURRENT_MODEL_INDEX()
+!$acc update device( XRHODJ, ZRHOM, ZAF, ZBFY, ZCF, &
+!$acc&                ZTRIGSX, ZTRIGSY, IIFAXX, IIFAXY, &
+!$acc&                ZTH, ZRR, ZRU, ZRV, ZRW, ZPABST, ZBFB, ZBF_SXP2_YP1_Z )
 CALL PRESSUREZ(CLBCX,CLBCY,CPRESOPT,NITR,LITRADJ,ITCOUNT,XRELAX,IMI,  &
                XRHODJ,XDXX,XDYY,XDZZ,XDZX,XDZY,ZDXHATM,ZDYHATM,ZRHOM, &
                ZAF,ZBFY,ZCF,ZTRIGSX,ZTRIGSY,IIFAXX,IIFAXY,            &
@@ -324,6 +327,7 @@ CALL PRESSUREZ(CLBCX,CLBCY,CPRESOPT,NITR,LITRADJ,ITCOUNT,XRELAX,IMI,  &
                ZTH,ZRR,XRHODREF,XTHVREF,XRVREF,XEXNREF, XLINMASS,     &
                ZRU,ZRV,ZRW,ZPABST,                                    &
                ZBFB,ZBF_SXP2_YP1_Z,PRESIDUAL                          )
+!$acc update self( ZRU, ZRV, ZRW, ZPABST )
 !
 CALL MPPDB_CHECK3D(XRHODJ,"anel_balancen3.2-after pressurez halo::XRHODJ",PRECISION)
 CALL MPPDB_CHECK3D(ZRU,"anel_balancen3.2-after pressurez::ZRU",PRECISION)