Changeset 4370 for palm/trunk/SOURCE/prognostic_equations.f90

Timestamp:

Jan 10, 2020 2:00:44 PM (5 years ago)

Author:

raasch

Message:

bugfixes for previous commit: unused variables removed, Temperton-fft usage on GPU, openacc porting of vector version of Obukhov length calculation, collective read switched off on NEC to avoid hanging; some vector directives added in prognostic equations to force vectorization on Intel19 compiler, configuration files for NEC Aurora added

File:

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

palm/trunk/SOURCE/prognostic_equations.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! vector directives added to force vectorization on Intel19 compiler
+! 
+! 4360 2020-01-07 11:25:50Z suehring
 ! Introduction of wall_flags_total_0, which currently sets bits based on static
 ! topography information used in wall_flags_static_0
@@ -936 +939 @@
     DO  i = nxlu, nxr
        DO  j = nys, nyn
+          !following directive is required to vectorize on Intel19
+          !DIR$ IVDEP
           DO  k = nzb+1, nzt
              u_p(k,j,i) = u(k,j,i) + ( dt_3d * ( tsc(2) * tend(k,j,i) +          &
@@ -1039 +1044 @@
     DO  i = nxl, nxr
        DO  j = nysv, nyn
+          !following directive is required to vectorize on Intel19
+          !DIR$ IVDEP
           DO  k = nzb+1, nzt
              v_p(k,j,i) = v(k,j,i) + ( dt_3d * ( tsc(2) * tend(k,j,i) +        &
@@ -1138 +1145 @@
     DO  i = nxl, nxr
        DO  j = nys, nyn
+          !following directive is required to vectorize on Intel19
+          !DIR$ IVDEP
           DO  k = nzb+1, nzt-1
              w_p(k,j,i) = w(k,j,i) + ( dt_3d * ( tsc(2) * tend(k,j,i) +        &
@@ -1277 +1286 @@
        DO  i = nxl, nxr
           DO  j = nys, nyn
+             !following directive is required to vectorize on Intel19
+             !DIR$ IVDEP
              DO  k = nzb+1, nzt
                 pt_p(k,j,i) = pt(k,j,i) + ( dt_3d * ( sbt * tend(k,j,i) +      &
@@ -1400 +1411 @@
        DO  i = nxl, nxr
           DO  j = nys, nyn
+             !following directive is required to vectorize on Intel19
+             !DIR$ IVDEP
              DO  k = nzb+1, nzt
                 q_p(k,j,i) = q(k,j,i) + ( dt_3d * ( sbt * tend(k,j,i) +        &
@@ -1520 +1533 @@
        DO  i = nxl, nxr
           DO  j = nys, nyn
+             !following directive is required to vectorize on Intel19
+             !DIR$ IVDEP
              DO  k = nzb+1, nzt
                 s_p(k,j,i) = s(k,j,i) + ( dt_3d * ( sbt * tend(k,j,i) +        &