| 1 |
C $Header: /u/gcmpack/MITgcm/eesupp/inc/BARRIER.h,v 1.6 2005/11/07 18:16:07 cnh Exp $ |
| 2 |
C $Name: $ |
| 3 |
CBOP |
| 4 |
C !ROUTINE: BARRIER.h |
| 5 |
C !INTERFACE: |
| 6 |
C include "BARRIER.h" |
| 7 |
C !DESCRIPTION: |
| 8 |
C *==========================================================* |
| 9 |
C | BARRIER.h |
| 10 |
C | o Globals used by Fortran barrier routine. |
| 11 |
C *==========================================================* |
| 12 |
C | These variables support a simple Fortran shared memory |
| 13 |
C | barrier routine. They do busy waiting, that is the |
| 14 |
C | thread that is waiting sits in a tight loop reading |
| 15 |
C | some memory location and waiting for all other threads. |
| 16 |
C | On some systems this is sometimes a good method to use. |
| 17 |
C | On the T3E memory is not shared so the routine should |
| 18 |
C | use the T3E "eureka" barriers. On CRAY and NEC there are |
| 19 |
C | hardware barriers that are accessed through a compiler |
| 20 |
C | directives. Finally proper multi-threading compilers |
| 21 |
C | support barrier compile directives - sometimes these |
| 22 |
C | are good, sometimes they are lousy. |
| 23 |
C | The barrier mechanism is used as follows |
| 24 |
C | 1. In the single-threaded part of the code |
| 25 |
C | CALL FBAR\_INIT |
| 26 |
C | on CRAY, NEC this routine does nothing |
| 27 |
C | on T3E there is no single-threaded code |
| 28 |
C | but there may be barrier initialisation - |
| 29 |
C | need to check. |
| 30 |
C | 2. When we need to synchronize everybody just |
| 31 |
C | CALL FBAR( myThid ) |
| 32 |
C | where myThid is myThreadId |
| 33 |
C | on CRAY, NEC FBAR will just do C\$DIR BARRIER |
| 34 |
C | or the like. |
| 35 |
C | on T3E FBAR does CALL BARRIER(...) or something |
| 36 |
C | need to check this. |
| 37 |
C *==========================================================* |
| 38 |
CEOP |
| 39 |
COMMON / BARRIER_COMMON / key1, key2, key3, |
| 40 |
& door1, door2, door3, |
| 41 |
& bCount, masterSet |
| 42 |
INTEGER key1(lShare4,MAX_NO_THREADS) |
| 43 |
INTEGER key2(lShare4,MAX_NO_THREADS) |
| 44 |
INTEGER key3(lShare4,MAX_NO_THREADS) |
| 45 |
INTEGER door1 |
| 46 |
INTEGER door2 |
| 47 |
INTEGER door3 |
| 48 |
INTEGER VALID |
| 49 |
PARAMETER ( VALID = 1 ) |
| 50 |
INTEGER INVALID |
| 51 |
PARAMETER ( INVALID = 0 ) |
| 52 |
INTEGER OPEN |
| 53 |
PARAMETER ( OPEN = 1 ) |
| 54 |
INTEGER SHUT |
| 55 |
PARAMETER ( SHUT = 0 ) |
| 56 |
INTEGER bCount(MAX_NO_THREADS) |
| 57 |
INTEGER masterSet(MAX_NO_THREADS) |