Changeset 4366 for palm/trunk/SOURCE/surface_mod.f90

Timestamp:

Jan 9, 2020 8:12:43 AM (5 years ago)

Author:

raasch

Message:

code vectorization for NEC Aurora: vectorized version of Temperton FFT, vectorization of Newtor iteration for calculating the Obukhov length

File:

• palm/trunk/SOURCE/surface_mod.f90 (modified) (2 diffs)

palm/trunk/SOURCE/surface_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! workaround implemented to avoid vectorization bug on NEC Aurora
+! 
+! 4360 2020-01-07 11:25:50Z suehring
 ! Fix also remaining message calls.
 ! 
@@ -2375 +2378 @@
                                                  num_def_h_kji(0) - 1
              start_index_def_h(0)           = surf_def_h(0)%end_index(j,i) + 1
+!
+!--          ATTENTION:
+!--          workaround to prevent vectorization bug on NEC Aurora
+             IF ( start_index_def_h(0) < -99999 )  THEN
+                PRINT*, 'i=', i, ' j=',j, ' s=',surf_def_h(0)%start_index(j,i),                    &
+                        ' e=', surf_def_h(0)%end_index(j,i)
+             ENDIF
 !
 !--          Downward-facing surfaces, except model top