Philippe 24/01/2022: bugfix: wrong gradient subroutines (bug introduced in commit 2a75fe0d)

src/ZSOLVER/turb_hor_thermo_flux.f90

@@ -1394,7 +1394,7 @@ IF (KSPLT==1 .AND. LLES_CALL) THEN
   CALL MZF_DEVICE( ZTMP4_DEVICE, ZTMP1_DEVICE )
   CALL LES_MEAN_SUBGRID( ZTMP1_DEVICE, X_LES_RES_ddxa_W_SBG_UaThl , .TRUE. )
   !
-  CALL GX_M_M_DEVICE(PTHLM,PDYY,PDZZ,PDZY,ZTMP1_DEVICE)
+  CALL GY_M_M_DEVICE(PTHLM,PDYY,PDZZ,PDZY,ZTMP1_DEVICE)
   CALL MYF_DEVICE(ZFLX,ZTMP2_DEVICE)
   !$acc kernels
   ZTMP3_DEVICE(:,:,:) = ZTMP1_DEVICE(:,:,:) * ZTMP2_DEVICE(:,:,:)
@@ -1405,7 +1405,7 @@ IF (KSPLT==1 .AND. LLES_CALL) THEN
   !
   IF (KRR>=1) THEN
 !$acc data copy(X_LES_RES_ddxa_Rt_SBG_UaThl)
-    CALL GX_M_M_DEVICE(PRM(:,:,:,1),PDYY,PDZZ,PDZY,ZTMP1_DEVICE)
+    CALL GY_M_M_DEVICE(PRM(:,:,:,1),PDYY,PDZZ,PDZY,ZTMP1_DEVICE)
     CALL MYF_DEVICE(ZFLX,ZTMP2_DEVICE)
     !$acc kernels
     ZTMP3_DEVICE(:,:,:) = ZTMP1_DEVICE(:,:,:) * ZTMP2_DEVICE(:,:,:)
@@ -1764,14 +1764,14 @@ IF (KRR/=0) THEN
     CALL MZF_DEVICE( ZTMP4_DEVICE, ZTMP1_DEVICE )
     CALL LES_MEAN_SUBGRID( ZTMP1_DEVICE,X_LES_RES_ddxa_W_SBG_UaRt , .TRUE. )
     !
-    CALL GX_M_M_DEVICE(PTHLM,PDYY,PDZZ,PDZY,ZTMP1_DEVICE)
+    CALL GY_M_M_DEVICE(PTHLM,PDYY,PDZZ,PDZY,ZTMP1_DEVICE)
     CALL MYF_DEVICE(ZFLX,ZTMP2_DEVICE)
     !$acc kernels
     ZTMP3_DEVICE(:,:,:) = ZTMP1_DEVICE(:,:,:) * ZTMP2_DEVICE(:,:,:)
     !$acc end kernels
     CALL LES_MEAN_SUBGRID( ZTMP3_DEVICE, X_LES_RES_ddxa_Thl_SBG_UaRt , .TRUE. )
     !
-    CALL GX_M_M_DEVICE(PRM(:,:,:,1),PDYY,PDZZ,PDZY,ZTMP1_DEVICE)
+    CALL GY_M_M_DEVICE(PRM(:,:,:,1),PDYY,PDZZ,PDZY,ZTMP1_DEVICE)
     CALL MYF_DEVICE(ZFLX,ZTMP2_DEVICE)
     !$acc kernels
     ZTMP3_DEVICE(:,:,:) = ZTMP1_DEVICE(:,:,:) * ZTMP2_DEVICE(:,:,:)