Changeset 1257 for palm/trunk/SOURCE/boundary_conds.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/boundary_conds.f90 (modified) (4 diffs)

palm/trunk/SOURCE/boundary_conds.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! loop independent clauses added
 !  
 ! Former revisions:
@@ -171 +171 @@
     IF ( ibc_pt_b == 0 )  THEN
        !$acc kernels present( nzb_s_inner, pt, pt_p )
+       !$acc loop independent
        DO  i = nxlg, nxrg
+          !$acc loop independent
           DO  j = nysg, nyng
              pt_p(nzb_s_inner(j,i),j,i) = pt(nzb_s_inner(j,i),j,i)
@@ -179 +181 @@
     ELSEIF ( ibc_pt_b == 1 )  THEN
        !$acc kernels present( nzb_s_inner, pt_p )
+       !$acc loop independent
        DO  i = nxlg, nxrg
+          !$acc loop independent
           DO  j = nysg, nyng
              pt_p(nzb_s_inner(j,i),j,i) = pt_p(nzb_s_inner(j,i)+1,j,i)
@@ -208 +212 @@
     IF ( .NOT. constant_diffusion )  THEN
        !$acc kernels present( e_p, nzb_s_inner )
+       !$acc loop independent
        DO  i = nxlg, nxrg
+          !$acc loop independent
           DO  j = nysg, nyng
              e_p(nzb_s_inner(j,i),j,i) = e_p(nzb_s_inner(j,i)+1,j,i)