Changeset 4152

Timestamp:

May 16, 2022, 11:43:50 AM (2 years ago)

Author:

Ehouarn Millour

Message:

Update Openacc directives for driver.
EM

File:

• LMDZ6/branches/Portage_acc/libf/dyn3d/driver.F90 (modified) (7 diffs)

LMDZ6/branches/Portage_acc/libf/dyn3d/driver.F90

@@ -55 +55 @@
   REAL,ALLOCATABLE :: pr(:,:) ! vorticity
 
-  !$acc data create(pplev, pplay, pphi, pphis, pt, pu, pv, pq, pr)
 
 ! 0. Preliminary stuff: read in parameters from run.def
@@ -169 +168 @@
   ALLOCATE(pv(ngrid,llm)) ! meridional wind at full levels
   ALLOCATE(pq(ngrid,llm,nqtot)) ! tracers
+
+  !$acc data create(pplev, pplay, pphi, pphis, pt, pu, pv, pq, pr)
 
   CALL init_temperature_pressure_geopot(ps,masse,teta,phis, &
@@ -224 +225 @@
    REAL :: pdq(ngrid,llm,nqtot)
    REAL :: pdps(ngrid)
-   !$acc data create(pdu,pdv,pdt,pdq,pdps)
    
    INTEGER :: iday,istep
@@ -232 +232 @@
    REAL :: flxw(ngrid,llm) ! vertical mass flux ! set to zero here.
    REAL :: phi_top(ngrid)
-   !$acc data create(phi_top,flxw)
+   !$acc data present(pplev,pplay,pphi,pphis,pu,pv,pt,pq) &
+   !$acc & create(phi_top,flxw,pdu,pdv,pdt,pdq,pdps)
    
    flxw(1:ngrid,1:llm)=0.
@@ -254 +255 @@
        pv(:,:)=pv(:,:)+dtphys*pdv(:,:)
        pq(:,:,:)=pq(:,:,:)+dtphys*pdq(:,:,:)
-       !$acc end kernels
        
        ! Recomputation of relative geopotential:
@@ -273 +273 @@
    ENDDO ! of DO iday=0,nday-1
    
+   !$acc end data
+   
    END SUBROUTINE timeloop
 
@@ -409 +411 @@
     USE exner_hyb_m, ONLY: exner_hyb
     USE comconst_mod, ONLY: cpp
-    USE geopot_mod, ONLY: geopot
 
     IMPLICIT NONE