| 62 |
|
|
| 63 |
\subsubsection{One-Dimensional Arrays} |
\subsubsection{One-Dimensional Arrays} |
| 64 |
|
|
| 65 |
The following arrays are indexed to the tile number, and the indices are |
The following arrays are of size {\em NTILES}, are indexed to the tile number, |
| 66 |
omitted in their descriptions. |
and the indices are omitted in their descriptions. |
| 67 |
|
|
| 68 |
The arrays {\em exch2\_tnx} and {\em exch2\_tny} |
The arrays {\em exch2\_tnx} and {\em exch2\_tny} |
| 69 |
express the x and y dimensions of each tile. At present for each tile |
express the x and y dimensions of each tile. At present for each tile |
| 72 |
MITgcm are to allow varying tile sizes. |
MITgcm are to allow varying tile sizes. |
| 73 |
|
|
| 74 |
The location of the tiles' Cartesian origin within a subdomain are determined |
The location of the tiles' Cartesian origin within a subdomain are determined |
| 75 |
by the arrays {\em exch2\_tbasex} and {\em exch2\_tbasey}. These |
by the arrays {\em exch2\_tbasex} and {\em exch2\_tbasey}. These variables |
| 76 |
|
are used to relate the location of the edges of the tiles to each other. As |
| 77 |
|
an example, in the default six-tile topology (the degenerate case) |
| 78 |
|
each index in these arrays are |
| 79 |
|
set to 0. The twenty-four, 32x32 cube face case discussed above will have |
| 80 |
|
values of 0 or 16, depending on the quadrant the tile falls within the |
| 81 |
|
subdomain. {\em exch2\_myFace} contains the number of the |
| 82 |
|
cubeface/subdomain of each tile, numbered 1-6 in the case of the standard |
| 83 |
|
cube topology. |
| 84 |
|
|
| 85 |
|
The arrays {\em exch2\_txglobalo} and {\em exch2\_txglobalo} are similar to |
| 86 |
|
{\em exch2\_tbasex} and {\em exch2\_tbasey}, but locate the tiles within |
| 87 |
|
the global address space, similar to that used by global files. |
| 88 |
|
|
| 89 |
|
The arrays {\em exch2\_isWedge}, {\em exch2\_isEedge}, {\em exch2\_isSedge}, |
| 90 |
|
and {\em exch2\_isNedge} are set to 1 if the indexed tile lies on the edge |
| 91 |
|
of a subdomain, 0 if not. The values are used within the topology generator |
| 92 |
|
to determine the orientation of neighboring tiles and to indicate whether |
| 93 |
|
a tile lies on the corner of a subdomain. The latter case indicates |
| 94 |
|
special exchange and numerical handling for the singularities at the eight |
| 95 |
|
corners of the cube. {\em exch2\_isNedge} contains a count of how many |
| 96 |
|
neighboring tiles each tile has, and is used for setting bounds for looping |
| 97 |
|
over neighboring tiles. {\em exch2\_tProc} holds the process rank of each tile, |
| 98 |
|
and is used in interprocess communication. |
| 99 |
|
|
| 100 |
\subsubsection{Two-Dimensional Arrays} |
\subsubsection{Two-Dimensional Arrays} |
| 101 |
|
|
| 102 |
|
The following arrays are all of size {\em MAX\_NEIGHBOURS}x{\em NTILES} and |
| 103 |
|
describe the orientations between the the tiles. |
| 104 |
|
|
| 105 |
|
The array {\em exch2\_neighbourId(a,T)} holds the tile number $T_{n}$ for each tile |
| 106 |
|
{\em T}'s neighbor tile {\em a}, and {\em exch2\_opposingSend\_record(a,T)} holds |
| 107 |
|
the index c in {\em exch2\_neighbourId(b,$T_{n}$)} that holds the tile number T. |
| 108 |
|
In other words, |
| 109 |
|
|
| 110 |
|
\begin{verbatim} |
| 111 |
|
exch2_neighbourId( exch2_opposingSend_record(a,T), exch2_neighbourId(a,T) ) = T |
| 112 |
|
\end{verbatim} |
| 113 |
|
|
| 114 |
|
{\em exch2\_neighbourId(exch2\_opposingSend\_record(a,T),exch2\_neighbourId(a,T))=T}. |
| 115 |
|
This is to provide a backreference from the neighbor tiles. |
| 116 |
|
|
| 117 |
|
|
| 118 |
// |
// |
| 119 |
|
|
| 120 |
\begin{verbatim} |
\begin{verbatim} |
| 121 |
C NTILES :: Number of tiles in this topology |
|
| 122 |
C MAX_NEIGHBOURS :: Maximum number of neighbours any tile has. |
|
|
C exch2_domain_nxt :: Total domain length in tiles. |
|
|
C exch2_domain_nyt :: Maximum domain height in tiles. |
|
|
C exch2_tnx :: Size in X for each tile. |
|
|
C exch2_tny :: Size in Y for each tile. |
|
|
C exch2_tbasex :: Tile offset in X within its sub-domain (cube face) |
|
|
C exch2_tbasey :: Tile offset in Y within its sub-domain (cube face) |
|
|
C exch2_tglobalxlo :: Tile base X index within global index space. |
|
|
C exch2_tglobalylo :: Tile base Y index within global index space. |
|
|
C exch2_isWedge :: 0 if West not at domain edge, 1 if it is. |
|
|
C exch2_isNedge :: 0 if North not at domain edge, 1 if it is. |
|
|
C exch2_isEedge :: 0 if East not at domain edge, 1 if it is. |
|
|
C exch2_isSedge :: 0 if South not at domain edge, 1 if it is. |
|
|
C exch2_myFace :: Cube face number used for I/O. |
|
|
C exch2_nNeighbours :: Tile neighbour entries count. |
|
|
C exch2_tProc :: Rank of process owning tile |
|
|
C :: (filled at run time). |
|
| 123 |
C exch2_neighbourId :: Tile number for each neighbour entry. |
C exch2_neighbourId :: Tile number for each neighbour entry. |
| 124 |
C exch2_opposingSend_record :: Record for entry in target tile send |
C exch2_opposingSend_record :: Record for entry in target tile send |
| 125 |
C :: list that has this tile and face |
C :: list that has this tile and face |