MITgcm_contrib/cg2d_bench/gsum.F

      SUBROUTINE GSUM_R8( ans, phi )

#ifdef ALLOW_MPI
#include "mpif.h"
#endif

      Real*8 ans
      Real*8 phi

      INTEGER RC

      ans = phi

#ifdef USE_MPI_GSUM
      CALL MPI_Allreduce( phi,
     &                    ans,
     &                      1,
     &                    MPI_REAL8,
     &                    MPI_SUM,
     &                    MPI_COMM_WORLD,
     &                    rc
     &                  )
#endif

#ifdef USE_JAM_GSUM
C     JAM global sum
      CALL JAM_barrier_start( phi )
      CALL JAM_barrier_done(  ans )
#endif

C     WRITE(6,*) ' phi = ', phi, ' Sum = ', ans
C     CALL MPI_Finalize( rc )
C     STOP

      RETURN
      END
1	SUBROUTINE GSUM_R8( ans, phi )
2
3	#ifdef ALLOW_MPI
4	#include "mpif.h"
5	#endif
6
7	Real*8 ans
8	Real*8 phi
9
10	INTEGER RC
11
12	ans = phi
13
14	#ifdef USE_MPI_GSUM
15	CALL MPI_Allreduce( phi,
16	& ans,
17	& 1,
18	& MPI_REAL8,
19	& MPI_SUM,
20	& MPI_COMM_WORLD,
21	& rc
22	& )
23	#endif
24
25	#ifdef USE_JAM_GSUM
26	C JAM global sum
27	CALL JAM_barrier_start( phi )
28	CALL JAM_barrier_done( ans )
29	#endif
30
31	C WRITE(6,*) ' phi = ', phi, ' Sum = ', ans
32	C CALL MPI_Finalize( rc )
33	C STOP
34
35	RETURN
36	END