Philippe 04/07/2016: compute compute SUM by hand instead of intrinsic Fortran...

Philippe 04/07/2016: compute compute SUM by hand instead of intrinsic Fortran in SUM_1DFIELD_ll to allow bit reproductibility with PGI compiler (16.4)

src/LIB/SURCOUCHE/src/mode_sum_ll.f90

@@ -1209,6 +1209,8 @@ REAL, DIMENSION(:,:), ALLOCATABLE :: ZSUM_ll
 !!    MODIFICATIONS
 !!    -------------
 !     Original 16/09/98
+!     2016/07/04 Philippe Wautelet: compute SUM by hand instead of intrinsic Fortran
+!                                   to allow bit reproductibility with PGI compiler (16.4)
 !
 !-------------------------------------------------------------------------------
 !
@@ -1248,6 +1250,7 @@ REAL, DIMENSION(:,:), ALLOCATABLE :: ZSUM_ll
   INTEGER :: IB, IE
   INTEGER :: IGB, IGE
   INTEGER :: IWEST, IEAST, INORTH, ISOUTH
+  INTEGER :: JI
 !
 !-------------------------------------------------------------------------------
 !
@@ -1345,8 +1348,12 @@ REAL, DIMENSION(:,:), ALLOCATABLE :: ZSUM_ll
 !*        3.   CALCULATE THE SUM
 !              -----------------
 !
-!OCL SCALAR
-  ZSUM = SUM(ZGLOBFIELD(1:IGE-IGB+1))
+!!OCL SCALAR
+!  ZSUM = SUM(ZGLOBFIELD(1:IGE-IGB+1))
+  ZSUM = 0.
+  DO JI = 1, IGE-IGB+1
+    ZSUM = ZSUM + ZGLOBFIELD(JI)
+  END DO
 !
   DEALLOCATE(ZGLOBFIELD)
 !