From af5ce936219a79b7915076549a14b507455119e5 Mon Sep 17 00:00:00 2001
From: Juan ESCOBAR <juan.escobar@aero.obs-mip.fr>
Date: Wed, 21 Sep 2022 19:15:58 +0200
Subject: [PATCH] Juan 21/09/2022:rain_ice_fast_rg.f90, GPU nvhpc optimization
 , mv mnh_expand over tab indirection I1(1:IGDRY) <-> not parallelizing
 without "acc independent"

---
 src/MNH/rain_ice_fast_rg.f90 | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/MNH/rain_ice_fast_rg.f90 b/src/MNH/rain_ice_fast_rg.f90
index 9d3fd8ae8..f512efac6 100644
--- a/src/MNH/rain_ice_fast_rg.f90
+++ b/src/MNH/rain_ice_fast_rg.f90
@@ -342,6 +342,7 @@ END IF
 !*       6.2.3  select the (PLBDAG,PLBDAS) couplet
 !
 !$acc kernels
+   !$mnh_expand_where(JL=1:IGDRY)
     ZVECLBDAG(1:IGDRY) = PLBDAG(I1(1:IGDRY))
     ZVECLBDAS(1:IGDRY) = PLBDAS(I1(1:IGDRY))
 !
@@ -349,7 +350,6 @@ END IF
 !               in the geometrical set of (Lbda_g,Lbda_s) couplet use to
 !               tabulate the SDRYG-kernel
 !
-   !$mnh_expand_where(JL=1:IGDRY)
     ZVEC1(1:IGDRY) = MAX( 1.00001, MIN( REAL(NDRYLBDAG)-0.00001,           &
 #if !defined(MNH_BITREP) && !defined(MNH_BITREP_OMP)
                           XDRYINTP1G * LOG( ZVECLBDAG(1:IGDRY) ) + XDRYINTP2G ) )
@@ -464,6 +464,7 @@ END IF
 !*       6.2.8  select the (PLBDAG,PLBDAR) couplet
 !
 !$acc kernels
+ !$mnh_expand_where(JL=1:IGDRY)    
     ZVECLBDAG(1:IGDRY) = PLBDAG(I1(1:IGDRY))
     ZVECLBDAR(1:IGDRY) = PLBDAR(I1(1:IGDRY))
 !
@@ -471,7 +472,6 @@ END IF
 !               in the geometrical set of (Lbda_g,Lbda_r) couplet use to
 !               tabulate the RDRYG-kernel
 !
- !$mnh_expand_where(JL=1:IGDRY)    
     ZVEC1(1:IGDRY) = MAX( 1.00001, MIN( REAL(NDRYLBDAG)-0.00001,           &
 #if !defined(MNH_BITREP) && !defined(MNH_BITREP_OMP)
                           XDRYINTP1G * LOG( ZVECLBDAG(1:IGDRY) ) + XDRYINTP2G ) )
-- 
GitLab