Changeset 1257 for palm/trunk/SOURCE/diffusion_v.f90

Timestamp:

Nov 8, 2013 3:18:40 PM (10 years ago)

Author:

raasch

Message:

New:
---

openACC porting of timestep calculation
(modules, timestep, time_integration)

Changed:

---

openACC loop directives and vector clauses removed (because they do not give any performance improvement with PGI
compiler versions > 13.6)
(advec_ws, buoyancy, coriolis, diffusion_e, diffusion_s, diffusion_u, diffusion_v, diffusion_w, diffusivities, exchange_horiz, fft_xy, pres, production_e, transpose, tridia_solver, wall_fluxes)

openACC loop independent clauses added
(boundary_conds, prandtl_fluxes, pres)

openACC declare create statements moved after FORTRAN declaration statement
(diffusion_u, diffusion_v, diffusion_w, fft_xy, poisfft, production_e, tridia_solver)

openACC end parallel replaced by end parallel loop
(flow_statistics, pres)

openACC "kernels do" replaced by "kernels loop"
(prandtl_fluxes)

output format for theta* changed to avoid output of *
(run_control)

Errors:

---

bugfix for calculation of advective timestep (old version may cause wrong timesteps in case of
vertixcally stretched grids)
Attention: standard run-control output has changed!
(timestep)

File:

• palm/trunk/SOURCE/diffusion_v.f90 (modified) (7 diffs)

palm/trunk/SOURCE/diffusion_v.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-!
+! openacc loop and loop vector clauses removed, declare create moved after
+! the FORTRAN declaration statement
 !
 ! Former revisions:
@@ -262 +263 @@
        REAL    ::  kmxm, kmxp, kmzm, kmzp
 
+       REAL, DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) ::  vsus
        !$acc declare create ( vsus )
-       REAL, DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) ::  vsus
 
 !
@@ -276 +277 @@
        !$acc         present ( ddzu, ddzw, fxm, fxp, wall_v )           &
        !$acc         present ( nzb_v_inner, nzb_v_outer, nzb_diff_v )
-       !$acc loop
        DO  i = i_left, i_right
           DO  j = j_south, j_north
 !
 !--          Compute horizontal diffusion
-             !$acc loop vector(32)
              DO  k = 1, nzt
                 IF ( k > nzb_v_outer(j,i) )  THEN
@@ -306 +305 @@
 !
 !--          Wall functions at the left and right walls, respectively
-             !$acc loop vector(32)
              DO  k = 1, nzt
                 IF( k > nzb_v_inner(j,i)  .AND.  k <= nzb_v_outer(j,i)  .AND. &
@@ -337 +335 @@
 !--          Compute vertical diffusion. In case of simulating a Prandtl
 !--          layer, index k starts at nzb_v_inner+2.
-             !$acc loop vector(32)
              DO  k = 1, nzt_diff
                 IF ( k >= nzb_diff_v(j,i) )  THEN
@@ -373 +370 @@
        IF ( use_surface_fluxes )  THEN
 
-          !$acc loop
           DO  i = i_left, i_right
-             !$acc loop vector(32)
              DO  j = j_south, j_north
           
@@ -404 +399 @@
           k = nzt
 
-          !$acc loop
           DO  i = i_left, i_right
-             !$acc loop vector(32)
              DO  j = j_south, j_north