Juan 25/08/2014: version reprod on GPU until next bug on advecuvw_4th

MNH/ppm_rhodj.f90

@@ -73,7 +73,11 @@ END MODULE MODI_PPM_RHODJ
 !*       0.    DECLARATIONS
 !              ------------
 !
-USE MODI_PPM             
+USE MODI_PPM
+     
+USE&
+    OPENACC    
+   
 !
 !
 !
@@ -87,11 +91,14 @@ CHARACTER (LEN=4), DIMENSION(2), INTENT(IN) :: HLBCY ! Y direction LBC type
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PCRU  ! contravariant
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PCRV  !  components
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PCRW  ! of momentum
+!$acc declare pcopyin(PCRU,PCRV,PCRW) 
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PRHODJ ! density
+!$acc declare present(PRHODJ)
 ! Temporary advected rhodj
 REAL, DIMENSION(:,:,:),   INTENT(OUT)   :: PRHOX1,PRHOX2
 REAL, DIMENSION(:,:,:),   INTENT(OUT)   :: PRHOY1,PRHOY2
 REAL, DIMENSION(:,:,:),   INTENT(OUT)   :: PRHOZ1,PRHOZ2
+!$acc declare pcopyout(PRHOX1,PRHOX2, PRHOY1,PRHOY2, PRHOZ1,PRHOZ2 )
 !
 REAL,                     INTENT(IN)    :: PTSTEP ! Time step 
 !
@@ -100,19 +107,32 @@ REAL,                     INTENT(IN)    :: PTSTEP ! Time step
 INTEGER :: IGRID ! localisation on the model grid
 !
 REAL, DIMENSION(SIZE(PCRU,1),SIZE(PCRU,2),SIZE(PCRU,3)) :: ZUNIT
+!BUG PGI ? acc declare create(ZUNIT)
 !
 !-------------------------------------------------------------------------------
 !
+
+!!$if ( acc_is_present(PRHODJ) )  print*,"PRHODJ JPRESENT"
+
 !
 IGRID = 1
 !
+!$acc data create (ZUNIT) 
+
+!$acc kernels
 ZUNIT = 1.0
+!$acc end kernels
+
+!$acc update device(PRHODJ)
+
 CALL PPM_S0_X(HLBCX, IGRID, ZUNIT, PCRU, PRHODJ, PTSTEP,PRHOX1)
 CALL PPM_S0_Y(HLBCY, IGRID, ZUNIT, PCRV, PRHOX1, PTSTEP,PRHOY1)
 CALL PPM_S0_Z(IGRID, ZUNIT, PCRW, PRHOY1, PTSTEP,PRHOZ1)
 CALL PPM_S0_Z(IGRID, ZUNIT, PCRW, PRHODJ, PTSTEP,PRHOZ2)
 CALL PPM_S0_Y(HLBCY, IGRID, ZUNIT, PCRV, PRHOZ2, PTSTEP,PRHOY2)
 CALL PPM_S0_X(HLBCX, IGRID, ZUNIT, PCRU, PRHOY2, PTSTEP,PRHOX2)
+
+!$acc end data
 !
 !
 END SUBROUTINE PPM_RHODJ