diff --git a/src/ZSOLVER/gdiv.f90 b/src/ZSOLVER/gdiv.f90
index 5bfd218bf519a7bf8bbbcae521aeffcf7e33ff0e..7072d22c6427b62af01e8b8caab48baab2e7b0bb 100644
--- a/src/ZSOLVER/gdiv.f90
+++ b/src/ZSOLVER/gdiv.f90
@@ -19,8 +19,8 @@ CHARACTER (LEN=4), DIMENSION(2), INTENT(IN) :: HLBCY    ! y-direction LBC type
                                                  ! Metric coefficients:       
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PDXX      ! d*xx 
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PDYY      ! d*yy 
-REAL, DIMENSION(:,:,:), INTENT(IN)  :: PDZX      ! d*zx 
-REAL, DIMENSION(:,:,:), INTENT(IN)  :: PDZY      ! d*zy 
+REAL, DIMENSION(:,:,:), INTENT(INOUT)  :: PDZX      ! d*zx 
+REAL, DIMENSION(:,:,:), INTENT(INOUT)  :: PDZY      ! d*zy 
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PDZZ      ! d*zz
 !
                                                  ! Field components
@@ -120,8 +120,8 @@ CHARACTER (LEN=4), DIMENSION(2), INTENT(IN) :: HLBCY    ! y-direction LBC type
                                                  ! Metric coefficients:       
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PDXX      ! d*xx 
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PDYY      ! d*yy 
-REAL, DIMENSION(:,:,:), INTENT(IN)  :: PDZX      ! d*zx 
-REAL, DIMENSION(:,:,:), INTENT(IN)  :: PDZY      ! d*zy 
+REAL, DIMENSION(:,:,:), INTENT(INOUT)  :: PDZX      ! d*zx 
+REAL, DIMENSION(:,:,:), INTENT(INOUT)  :: PDZY      ! d*zy 
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PDZZ      ! d*zz
 !
                                                  ! Field components
@@ -150,8 +150,8 @@ INTEGER :: IKE          ! indice K for the last inner mass point along z
 !
 INTEGER :: JI,JJ,JK                         ! loop indexes
 !
-#ifdef MNH_OPENACC
 INTEGER :: IIU,IJU,IKU
+#ifdef MNH_OPENACC
 !
 REAL, DIMENSION(:,:,:) , POINTER , CONTIGUOUS :: ZTMP1,ZTMP2
 INTEGER :: IZTMP1,IZTMP2
@@ -176,7 +176,7 @@ IKU=SIZE(PU,3)
 !
 #ifndef MNH_OPENACC
 ALLOCATE(ZUC(IIU,IJU,IKU),ZVC(IIU,IJU,IKU),ZWC(IIU,IJU,IKU))
-ALLOCATE(Z1(IIU,IJU,IKU),Z2(IIU,IJU,IKU),Z2(IIU,IJU,IKU))
+ALLOCATE(Z1(IIU,IJU,IKU),Z2(IIU,IJU,IKU),Z3(IIU,IJU,IKU))
 #else
 IZUC = MNH_ALLOCATE_ZT3D(ZUC,IIU,IJU,IKU  )
 IZVC = MNH_ALLOCATE_ZT3D(ZVC,IIU,IJU,IKU  )