Changeset 578

Timestamp:

10/13/17 13:59:36 (7 years ago)

Author:

dubos

Message:

devel : OpenMP fixes for NH kernels

Location:

codes/icosagcm/devel/src/kernels

Files:

• caldyn_slow_NH.k90 (modified) (1 diff)
• caldyn_solver.k90 (modified) (2 diffs)
• caldyn_vert_NH.k90 (modified) (1 diff)

codes/icosagcm/devel/src/kernels/caldyn_slow_NH.k90

@@ -145 +145 @@
       END DO
    END DO
+   ! We need a barrier here because we compute gradPhi2, F_el and w_il above and do a vertical average below
+   !$OMP BARRIER
    DO l = ll_begin, ll_end
       !DIR$ SIMD

codes/icosagcm/devel/src/kernels/caldyn_solver.k90

@@ -55 +55 @@
       END DO
    END DO
+   ! We need a barrier here because we compute pres above and do a vertical difference below
    !$OMP BARRIER
    IF (ll_begin==1) THEN
@@ -82 +83 @@
       END DO
    END IF
+   ! We need a barrier here because we update W above and do a vertical average below
+   !$OMP BARRIER
    DO l = ll_begin, ll_end
       !DIR$ SIMD

codes/icosagcm/devel/src/kernels/caldyn_vert_NH.k90

@@ -48 +48 @@
       END DO
    END DO
-   ! We need a barrier here because previously computed W_etadot is needed at 'CELL' and 'DOWN(CELL)'
+   ! We need a barrier here because we compute W_etadot above and do a vertical difference below
    !$OMP BARRIER
    IF (ll_begin==1) THEN