10 |
%% o automatically inserted at \section{Reference} |
%% o automatically inserted at \section{Reference} |
11 |
|
|
12 |
|
|
13 |
\section{Extended Cubed Sphere Exchange} |
\section{exch2: Extended Cubed Sphere \mbox{Topology}} |
14 |
\label{sec:exch2} |
\label{sec:exch2} |
15 |
|
|
16 |
|
|
24 |
individual tiles may be run on separate processors in different |
individual tiles may be run on separate processors in different |
25 |
combinations, and whether exchanges between particular tiles occur |
combinations, and whether exchanges between particular tiles occur |
26 |
between different processors is determined at runtime. This |
between different processors is determined at runtime. This |
27 |
flexibility provides for manual load balancing across a relatively |
flexibility provides for manual compile-time load balancing across a |
28 |
arbitrary number of processors. |
relatively arbitrary number of processors. \\ |
29 |
|
|
30 |
The exchange parameters are declared in |
The exchange parameters are declared in |
31 |
\filelink{pkg/exch2/W2\_EXCH2\_TOPOLOGY.h}{pkg-exch2-W2_EXCH2_TOPOLOGY.h} |
\filelink{pkg/exch2/W2\_EXCH2\_TOPOLOGY.h}{pkg-exch2-W2_EXCH2_TOPOLOGY.h} |
32 |
and assigned in |
and assigned in |
33 |
\filelink{pkg/exch2/w2\_e2setup.F}{pkg-exch2-w2_e2setup.F}. The |
\filelink{pkg/exch2/w2\_e2setup.F}{pkg-exch2-w2_e2setup.F}. The |
34 |
validity of the cube topology depends on the \texttt{SIZE.h} file as |
validity of the cube topology depends on the \texttt{SIZE.h} file as |
35 |
detailed below. Both files are generated by Matlab scripts in ?? |
detailed below. Both files are generated by Matlab scripts in |
36 |
check these in already! and should not be edited. The default files |
\texttt{utils/exch2/matlab-topology-generator}; see Section |
37 |
provided in the release configure a cubed sphere arrangement of six |
\ref{sec:topogen} for details on creating alternate topologies. The |
38 |
tiles, one per subdomain, each with 32$\times$32 grid points, all |
default files provided in the release configure a cubed sphere |
39 |
running on a single processor. Pregenerated examples of these files |
topology of six tiles, one per subdomain, each with 32$\times$32 grid |
40 |
with alternate topologies are provided in ??. |
points, all running on a single processor. Pregenerated examples of |
41 |
|
these files with alternate topologies are provided under |
42 |
|
\texttt{utils/exch2/code-mods} along with the appropriate |
43 |
|
\texttt{SIZE.h} file for single-processor execution. |
44 |
|
|
45 |
\subsection{Invoking exch2} |
\subsection{Invoking exch2} |
46 |
|
|
47 |
To use exch2 with the cubed sphere, the following conditions must be met: |
To use exch2 with the cubed sphere, the following conditions must be |
48 |
|
met: \\ |
49 |
|
|
50 |
- the exch2 package is included when \texttt{genmake2} is run. The |
$\bullet$ The exch2 package is included when \texttt{genmake2} is run. The |
51 |
easiest way to do this is to add the line \texttt{exch2} to the |
easiest way to do this is to add the line \texttt{exch2} to the |
52 |
\texttt{profile.conf} file -- see Section \ref{sect:buildingCode} |
\texttt{profile.conf} file -- see Section \ref{sect:buildingCode} |
53 |
for general details. \\ |
for general details. \\ |
54 |
|
|
55 |
- an example of \texttt{W2\_EXCH2\_TOPOLOGY.h} and |
$\bullet$ An example of \texttt{W2\_EXCH2\_TOPOLOGY.h} and |
56 |
\texttt{w2\_e2setup.F} must reside in a directory containing code |
\texttt{w2\_e2setup.F} must reside in a directory containing code |
57 |
linked when \texttt{genmake2} runs. The safest place to put these |
linked when \texttt{genmake2} runs. The safest place to put these |
58 |
is the directory indicated in the \texttt{-mods=DIR} command line |
is the directory indicated in the \texttt{-mods=DIR} command line |
59 |
modifier (typically \texttt{../code}), or the build directory. The |
modifier (typically \texttt{../code}), or the build directory. The |
60 |
default versions of these files reside in \texttt{pkg/exch2}, but |
default versions of these files reside in \texttt{pkg/exch2} and are |
61 |
they should be left untouched to avoid breaking configurations other |
linked automatically if no other versions exist elsewhere in the |
62 |
than the one you intend to modify.\\ |
link path, but they should be left untouched to avoid breaking |
63 |
|
configurations other than the one you intend to modify.\\ |
64 |
- files containing grid parameters, named |
|
65 |
\texttt{tile}xxx\texttt{.mitgrid} where xxx is \texttt{001} through |
$\bullet$ Files containing grid parameters, named |
66 |
\texttt{006}, must be in the working directory when the MITgcm |
\texttt{tile}???\texttt{.mitgrid} where ??? is \texttt{001} through |
67 |
executable is run. These files are provided in the example |
\texttt{006} (one per subdomain), must be in the working directory |
68 |
experiments for cubed sphere configurations with 32$\times$32 cube |
when the MITgcm executable is run. These files are provided in the |
69 |
sides and are non-trivial to generate -- please contact MITgcm |
example experiments for cubed sphere configurations with |
70 |
support if you want to generate files for other configurations. |
32$\times$32 cube sides and are non-trivial to generate -- please |
71 |
This is lame. ?? \\ |
contact MITgcm support if you want to generate files for other |
72 |
|
configurations. \\ |
73 |
|
|
74 |
|
$\bullet$ As always when compiling MITgcm, the file \texttt{SIZE.h} |
75 |
|
must be placed where \texttt{genmake2} will find it. In particular |
76 |
|
for the exch2, the domain decompositin specified in \texttt{SIZE.h} |
77 |
|
must correspond with the particular configuration's topology |
78 |
|
specified in \texttt{W2\_EXCH2\_TOPOLOGY.h} and |
79 |
|
\texttt{w2\_e2setup.F}. Domain decomposition issues particular to |
80 |
|
exch2 are addressed in Section \ref{sec:topogen}: ``Generating |
81 |
|
Topology Files for exch2''; a more general background on the subject |
82 |
|
relvant to MITgcm is presented in Section |
83 |
|
\ref{sect:specifying_a_decomposition}: ``Specifying a |
84 |
|
decomposition''.\\ |
85 |
|
|
86 |
As of the time of writing the following examples use exch2 and may be |
As of the time of writing the following examples use exch2 and may be |
87 |
used for guidance: |
used for guidance: |
97 |
|
|
98 |
|
|
99 |
|
|
100 |
\subsection{Generating Topology Files} |
\subsection{Generating Topology Files for exch2} |
101 |
|
\label{sec:topogen} |
102 |
|
|
103 |
|
Alternate cubed sphere topologies may be created using the Matlab |
104 |
|
scripts in \texttt{utils/exch2/matlab-topology-generator}. Running the |
105 |
|
m-file \texttt{driver} from the Matlab prompt (without passing any |
106 |
|
function parameters) generates exch2 topology files |
107 |
|
\texttt{W2\_EXCH2\_TOPOLOGY.h} and \texttt{w2\_e2setup.F} in the |
108 |
|
working directory and displays via Matlab a figure of the topology. |
109 |
|
The other m-files in the directory are subroutines of \texttt{driver} |
110 |
|
and should not be run except for development purposes. \\ |
111 |
|
|
112 |
|
The parameters that determine the dimensions and topology of the |
113 |
|
generated configuration are nr, nb, ng, tnx and tny, and all are |
114 |
|
assigned early in the script. |
115 |
|
|
116 |
|
The first three determine the size of the subdomains (cube faces) and |
117 |
|
hence the size of the overall domain. Each one determines the number |
118 |
|
of grid points, and therefore the resolution, along the subdomain |
119 |
|
sides in a ``great circle'' around each axis of the cube. At the time |
120 |
|
of this writing MITgcm requires these three parameters to be equal, |
121 |
|
but they provide for future releases of MITgcm to accomodate different |
122 |
|
resolutions around the axes to allow (for example) greater resolution |
123 |
|
around the equator.\\ |
124 |
|
|
125 |
|
The parameters tnx and tny determine the dimensions of the tiles into |
126 |
|
which the subdomains are decomposed, and must evenly divide the |
127 |
|
integer assigned to nr, nb and ng. The result is a rectangular tiling |
128 |
|
of the subdomain. Figure \ref{fig:24tile} shows one possible topology |
129 |
|
for a twenty-four tile cube, and figure \ref{fig:12tile} shows one for |
130 |
|
twelve tiles. \\ |
131 |
|
|
132 |
|
\begin{figure} |
133 |
|
\begin{center} |
134 |
|
\resizebox{4in}{!}{ |
135 |
|
\includegraphics{part6/s24t_16x16.ps} |
136 |
|
} |
137 |
|
\end{center} |
138 |
|
\caption{Plot of cubed sphere topology with a 32$\times$32 grid and |
139 |
|
twenty-four tiles (tnx=16, tny=16) |
140 |
|
} \label{fig:24tile} |
141 |
|
\end{figure} |
142 |
|
|
143 |
|
\begin{figure} |
144 |
|
\begin{center} |
145 |
|
\resizebox{4in}{!}{ |
146 |
|
\includegraphics{part6/s12t_16x32.ps} |
147 |
|
} |
148 |
|
\end{center} |
149 |
|
\caption{Plot of cubed sphere topology with a 32$\times$32 grid and |
150 |
|
twelve tiles (tnx=16, tny=32) |
151 |
|
} \label{fig:12tile} |
152 |
|
\end{figure} |
153 |
|
|
154 |
|
Tiles can be selected from the topology to be omitted from being |
155 |
|
allocated memory and processors. This kind of tuning is useful in |
156 |
|
ocean modeling for omitting tiles that fall entirely on land. The |
157 |
|
tiles omitted are specified in the file \texttt{blanklist.txt} by |
158 |
|
their tile number in the topology, separated by a newline. \\ |
159 |
|
|
160 |
|
|
161 |
|
|
162 |
|
|
163 |
|
|
164 |
|
|
165 |
\subsection{Key Variables} |
\subsection{Key Variables} |
166 |
|
|