eesupp/src/ini_jam.F

C $Header: /u/gcmpack/models/MITgcmUV/eesupp/src/ini_jam.F,v 1.3 2000/04/05 14:52:27 adcroft Exp $
C $Name: $

#include "CPP_EEOPTIONS.h"

#define USE_JAM_INIT

      SUBROUTINE INI_JAM

C     Initialise JAM based communication
      IMPLICIT NONE

#ifdef LETS_MAKE_JAM

C     == Global data ==
#include "SIZE.h"
#include "EEPARAMS.h"
#include "EESUPPORT.h"

#include "JAM_INFO.h"
#include "MPI_INFO.h"

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

C     == Local variables ==
      Real*8  dummyVal
      INTEGER rc
      CHARACTER*(MAX_LEN_FNAM) fnam
#ifdef JAM_WITH_TWO_PROCS_PER_NODE
      INTEGER myTwoProcRank
#endif

C     Names used in JAM exch testing for MPI based exchange
      mpi_northId = mpiPidN
      mpi_southId = mpiPidS

#ifdef USE_JAM_INIT
C     JAM initialisation. This should work with or without
C     MPI. If we dont use MPI we have to start procs. by hand
C     though!
#ifndef JAM_WITH_TWO_PROCS_PER_NODE
      CALL JAM_collective_init
#else  /* ! JAM_WITH_TWO_PROCS_PER_NODE */
      myTwoProcRank = MOD(myProcId,2)
      IF ( myTwoProcRank .EQ. 0 ) THEN
       myTwoProcRank = 1
      ELSE
       myTwoProcRank = 0
      ENDIF
      CALL JAM_collective_init( myTwoProcRank )
#endif /* JAM_WITH_TWO_PROCS_PER_NODE */

      CALL GET_JAM_self_address(   jam_pid )
      CALL GET_JAM_partition_size( jam_np  )
      jam_exchKey = 100
      dummyVal    = 0.
      CALL JAM_barrier_start( dummyVal )
      CALL JAM_barrier_done(  dummyVal )

C     Set up connectivity
      jam_northId = jam_pid+1
      IF ( jam_northId .EQ. jam_np ) jam_northId = 0
      jam_southId = jam_pid-1
      IF ( jam_southId .LT. 0      ) jam_southId = jam_np-1

      numberOfProcs = jam_np
      myProcId      = jam_pid

#ifdef JAM_WITH_TWO_PROCS_PER_NODE
C     The following fixes the instance where MPI and
C     James JAM arrive at a different id for this process
C     *note* this assumes nPx=1 (no parallization in X)
      myPid         = jam_pid
      pidN          = mod(jam_pid+1,jam_np)
      pidS          = mod(jam_pid-1+jam_np,jam_np)
      pidW          = jam_pid
      pidE          = jam_pid
      myYGlobalLo   = sNy*nSy*jam_pid+1
      myXGlobalLo   = 1
#endif

#endif
#endif /* LETS_MAKE_JAM */

      RETURN
      END
1	cnh	1.4	C $Header: /u/gcmpack/models/MITgcmUV/eesupp/src/ini_jam.F,v 1.3 2000/04/05 14:52:27 adcroft Exp $
2			C $Name: $
3	adcroft	1.1
4			#include "CPP_EEOPTIONS.h"
5
6			#define USE_JAM_INIT
7
8			SUBROUTINE INI_JAM
9
10			C Initialise JAM based communication
11			IMPLICIT NONE
12
13			#ifdef LETS_MAKE_JAM
14
15			C == Global data ==
16			#include "SIZE.h"
17			#include "EEPARAMS.h"
18			#include "EESUPPORT.h"
19
20			#include "JAM_INFO.h"
21			#include "MPI_INFO.h"
22
23			#ifdef ALLOW_MPI
24			#include "mpif.h"
25			#endif
26
27			C == Local variables ==
28	adcroft	1.2	Real*8 dummyVal
29	adcroft	1.1	INTEGER rc
30			CHARACTER*(MAX_LEN_FNAM) fnam
31			#ifdef JAM_WITH_TWO_PROCS_PER_NODE
32			INTEGER myTwoProcRank
33			#endif
34
35			C Names used in JAM exch testing for MPI based exchange
36			mpi_northId = mpiPidN
37			mpi_southId = mpiPidS
38
39			#ifdef USE_JAM_INIT
40			C JAM initialisation. This should work with or without
41			C MPI. If we dont use MPI we have to start procs. by hand
42			C though!
43			#ifndef JAM_WITH_TWO_PROCS_PER_NODE
44			CALL JAM_collective_init
45			#else /* ! JAM_WITH_TWO_PROCS_PER_NODE */
46			myTwoProcRank = MOD(myProcId,2)
47			IF ( myTwoProcRank .EQ. 0 ) THEN
48			myTwoProcRank = 1
49			ELSE
50			myTwoProcRank = 0
51			ENDIF
52			CALL JAM_collective_init( myTwoProcRank )
53			#endif /* JAM_WITH_TWO_PROCS_PER_NODE */
54
55			CALL GET_JAM_self_address( jam_pid )
56			CALL GET_JAM_partition_size( jam_np )
57			jam_exchKey = 100
58			dummyVal = 0.
59			CALL JAM_barrier_start( dummyVal )
60			CALL JAM_barrier_done( dummyVal )
61
62			C Set up connectivity
63			jam_northId = jam_pid+1
64			IF ( jam_northId .EQ. jam_np ) jam_northId = 0
65			jam_southId = jam_pid-1
66			IF ( jam_southId .LT. 0 ) jam_southId = jam_np-1
67
68			numberOfProcs = jam_np
69			myProcId = jam_pid
70	adcroft	1.3
71			#ifdef JAM_WITH_TWO_PROCS_PER_NODE
72			C The following fixes the instance where MPI and
73			C James JAM arrive at a different id for this process
74			C note this assumes nPx=1 (no parallization in X)
75			myPid = jam_pid
76			pidN = mod(jam_pid+1,jam_np)
77			pidS = mod(jam_pid-1+jam_np,jam_np)
78			pidW = jam_pid
79			pidE = jam_pid
80			myYGlobalLo = sNynSyjam_pid+1
81			myXGlobalLo = 1
82			#endif
83	adcroft	1.1
84			#endif
85			#endif /* LETS_MAKE_JAM */
86
87			RETURN
88			END