Changeset 3634 for palm/trunk/SOURCE/pres.f90

Timestamp:

Dec 18, 2018 12:31:28 PM (5 years ago)

Author:

knoop

Message:

OpenACC port for SPEC

File:

• palm/trunk/SOURCE/pres.f90 (modified) (15 diffs)

palm/trunk/SOURCE/pres.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! OpenACC port for SPEC
+! 
+! 3347 2018-10-15 14:21:08Z suehring
 ! Bugfixes in offline nesting.
 ! Add comment. 
@@ -232 +235 @@
 
     CALL cpu_log( log_point(8), 'pres', 'start' )
+
+!$ACC DATA COPYOUT(d, tend) &
+!$ACC COPY(u, v, w, p) &
+!$ACC COPYIN(rho_air, rho_air_zw, ddzu, ddzw, wall_flags_0)
+
+!$ACC DATA &
+!$ACC COPYIN(bc_h(0:1)) &
+!$ACC COPYIN(bc_h(0)%i(1:bc_h(0)%ns)) &
+!$ACC COPYIN(bc_h(0)%j(1:bc_h(0)%ns)) &
+!$ACC COPYIN(bc_h(0)%k(1:bc_h(0)%ns)) &
+!$ACC COPYIN(bc_h(1)%i(1:bc_h(1)%ns)) &
+!$ACC COPYIN(bc_h(1)%j(1:bc_h(1)%ns)) &
+!$ACC COPYIN(bc_h(1)%k(1:bc_h(1)%ns))
 
 
@@ -442 +458 @@
     ELSE
        !$OMP PARALLEL DO SCHEDULE( STATIC ) PRIVATE (i,j,k)
+       !$ACC PARALLEL LOOP COLLAPSE(3) PRIVATE(i, j, k) &
+       !$ACC PRESENT(d)
        DO  i = nxl, nxr
           DO  j = nys, nyn
@@ -490 +508 @@
     !$OMP PARALLEL PRIVATE (i,j,k)
     !$OMP DO SCHEDULE( STATIC )
+    !$ACC PARALLEL LOOP COLLAPSE(3) PRIVATE(i, j, k) &
+    !$ACC PRESENT(u, v, w, rho_air, rho_air_zw, ddzw, wall_flags_0) &
+    !$ACC PRESENT(d)
     DO  i = nxl, nxr
        DO  j = nys, nyn
@@ -513 +534 @@
        !$OMP PARALLEL PRIVATE (i,j,k) FIRSTPRIVATE(threadsum) REDUCTION(+:localsum)
        !$OMP DO SCHEDULE( STATIC )
+       !$ACC PARALLEL LOOP COLLAPSE(3) PRIVATE(i,j,k) &
+       !$ACC REDUCTION(+:threadsum) COPY(threadsum) &
+       !$ACC PRESENT(d)
        DO  i = nxl, nxr
           DO  j = nys, nyn
@@ -547 +571 @@
 !--    z-direction
        !$OMP PARALLEL DO PRIVATE (i,j,k)
+       !$ACC PARALLEL LOOP COLLAPSE(3) PRIVATE(i, j, k) &
+       !$ACC PRESENT(d, tend)
        DO  i = nxl, nxr
           DO  j = nys, nyn
@@ -565 +591 @@
 !--       Upward facing
           !$OMP PARALLEL DO PRIVATE( i, j, k )
+          !$ACC PARALLEL LOOP PRIVATE(i, j, k) &
+          !$ACC PRESENT(bc_h, tend)
           DO  m = 1, bc_h(0)%ns
              i = bc_h(0)%i(m)
@@ -574 +602 @@
 !--       Downward facing
           !$OMP PARALLEL DO PRIVATE( i, j, k )
+          !$ACC PARALLEL LOOP PRIVATE(i, j, k) &
+          !$ACC PRESENT(bc_h, tend)
           DO  m = 1, bc_h(1)%ns
              i = bc_h(1)%i(m)
@@ -621 +651 @@
 !--       Dirichlet
           !$OMP PARALLEL DO PRIVATE (i,j,k)
+          !$ACC PARALLEL LOOP COLLAPSE(2) PRIVATE(i, j) &
+          !$ACC PRESENT(tend)
           DO  i = nxlg, nxrg
              DO  j = nysg, nyng
@@ -688 +720 @@
        !$OMP PARALLEL PRIVATE (i,j,k)
        !$OMP DO
+       !$ACC PARALLEL LOOP COLLAPSE(3) PRIVATE(i, j, k) &
+       !$ACC PRESENT(p, tend)
        DO  i = nxl-1, nxr+1
           DO  j = nys-1, nyn+1
@@ -701 +735 @@
        !$OMP PARALLEL PRIVATE (i,j,k)
        !$OMP DO
+       !$ACC PARALLEL LOOP COLLAPSE(3) PRIVATE(i, j, k) &
+       !$ACC PRESENT(p, tend)
        DO  i = nxl-1, nxr+1
           DO  j = nys-1, nyn+1
@@ -733 +769 @@
     !$OMP PARALLEL PRIVATE (i,j,k)
     !$OMP DO
+    !$ACC PARALLEL LOOP COLLAPSE(2) PRIVATE(i, j, k) &
+    !$ACC PRESENT(u, v, w, tend, ddzu, wall_flags_0)
     DO  i = nxl, nxr   
        DO  j = nys, nyn
@@ -900 +938 @@
 #else
        !$OMP DO SCHEDULE( STATIC )
+       !$ACC PARALLEL LOOP COLLAPSE(3) PRIVATE(i, j, k) &
+       !$ACC PRESENT(u, v, w, rho_air, rho_air_zw, ddzw, wall_flags_0) &
+       !$ACC PRESENT(d)
        DO  i = nxl, nxr
           DO  j = nys, nyn
@@ -916 +957 @@
 !--    Compute possible PE-sum of divergences for flow_statistics
        !$OMP DO SCHEDULE( STATIC )
+       !$ACC PARALLEL LOOP COLLAPSE(3) PRIVATE(i, j, k) &
+       !$ACC REDUCTION(+:threadsum) COPY(threadsum) &
+       !$ACC PRESENT(d)
        DO  i = nxl, nxr
           DO  j = nys, nyn
@@ -939 +983 @@
     CALL cpu_log( log_point(8), 'pres', 'stop' )
 
+!$ACC END DATA
+!$ACC END DATA
 
  END SUBROUTINE pres