Juan 14/09/2021:ZSOLVER/turb_hor_dyn_corr.f90, use GET_HALO_D for GPU optimisation

src/ZSOLVER/turb_hor_dyn_corr.f90

@@ -177,6 +177,7 @@ USE MODI_BITREP
 #ifdef MNH_OPENACC
 USE MODE_MNH_ZWORK , ONLY : MNH_ALLOCATE_ZT3D , MNH_REL_ZT3D, MNH_ALLOCATE_ZT3DP , MNH_ALLOCATE_ZT2D
 #endif
+USE MODI_GET_HALO
 !
 IMPLICIT NONE
 !
@@ -659,9 +660,10 @@ ZFLX(:,:,IKB-1) = 2. * ZFLX(:,:,IKB-1) -  ZFLX(:,:,IKB)
 !!! to be absolutely sure, we do a wait
 !$acc wait
 !
-!$acc update self(ZFLX)
-CALL UPDATE_HALO_ll(TZFIELDS_ll, IINFO_ll)
-!$acc update device(ZFLX) async(10)
+!!$!$acc update self(ZFLX)
+!!$CALL UPDATE_HALO_ll(TZFIELDS_ll, IINFO_ll)
+!!$!$acc update device(ZFLX) async(10)
+CALL GET_HALO_D(ZFLX)
 !
 IF ( OCLOSE_OUT .AND. OTURB_FLX ) THEN
   ! stores <U U>