\hypertarget{namespacemom__fixed__initialization}{}\section{mom\+\_\+fixed\+\_\+initialization Module Reference}
\label{namespacemom__fixed__initialization}\index{mom\+\_\+fixed\+\_\+initialization@{mom\+\_\+fixed\+\_\+initialization}}


\subsection{Detailed Description}
Initializes fixed aspects of the model, such as horizontal grid metrics, topography and Coriolis. \subsection*{Functions/\+Subroutines}
\begin{DoxyCompactItemize}
\item 
subroutine, public \mbox{\hyperlink{namespacemom__fixed__initialization_a79fad1bd49f85eaa0579def2d9d4ae53}{mom\+\_\+initialize\+\_\+fixed}} (G, US, O\+BC, PF, write\+\_\+geom, output\+\_\+dir)
\begin{DoxyCompactList}\small\item\em M\+O\+M\+\_\+initialize\+\_\+fixed sets up time-\/invariant quantities related to M\+O\+M6\textquotesingle{}s horizontal grid, bathymetry, and the Coriolis parameter. \end{DoxyCompactList}\item 
subroutine, public \mbox{\hyperlink{namespacemom__fixed__initialization_ae980c841ada6296d07df08442398e87b}{mom\+\_\+initialize\+\_\+topography}} (D, max\+\_\+depth, G, PF, US)
\begin{DoxyCompactList}\small\item\em M\+O\+M\+\_\+initialize\+\_\+topography makes the appropriate call to set up the bathymetry. At this point the topography is in units of \mbox{[}m\mbox{]}, but this can be changed later. \end{DoxyCompactList}\end{DoxyCompactItemize}


\subsection{Function/\+Subroutine Documentation}
\mbox{\Hypertarget{namespacemom__fixed__initialization_a79fad1bd49f85eaa0579def2d9d4ae53}\label{namespacemom__fixed__initialization_a79fad1bd49f85eaa0579def2d9d4ae53}} 
\index{mom\+\_\+fixed\+\_\+initialization@{mom\+\_\+fixed\+\_\+initialization}!mom\+\_\+initialize\+\_\+fixed@{mom\+\_\+initialize\+\_\+fixed}}
\index{mom\+\_\+initialize\+\_\+fixed@{mom\+\_\+initialize\+\_\+fixed}!mom\+\_\+fixed\+\_\+initialization@{mom\+\_\+fixed\+\_\+initialization}}
\subsubsection{\texorpdfstring{mom\+\_\+initialize\+\_\+fixed()}{mom\_initialize\_fixed()}}
{\footnotesize\ttfamily subroutine, public mom\+\_\+fixed\+\_\+initialization\+::mom\+\_\+initialize\+\_\+fixed (\begin{DoxyParamCaption}\item[{type(dyn\+\_\+horgrid\+\_\+type), intent(inout)}]{G,  }\item[{type(unit\+\_\+scale\+\_\+type), intent(in)}]{US,  }\item[{type(ocean\+\_\+obc\+\_\+type), pointer}]{O\+BC,  }\item[{type(param\+\_\+file\+\_\+type), intent(in)}]{PF,  }\item[{logical, intent(in)}]{write\+\_\+geom,  }\item[{character(len=$\ast$), intent(in)}]{output\+\_\+dir }\end{DoxyParamCaption})}



M\+O\+M\+\_\+initialize\+\_\+fixed sets up time-\/invariant quantities related to M\+O\+M6\textquotesingle{}s horizontal grid, bathymetry, and the Coriolis parameter. 


\begin{DoxyParams}[1]{Parameters}
\mbox{\tt in,out}  & {\em g} & The ocean\textquotesingle{}s grid structure.\\
\hline
\mbox{\tt in}  & {\em us} & A dimensional unit scaling type\\
\hline
 & {\em obc} & Open boundary structure.\\
\hline
\mbox{\tt in}  & {\em pf} & A structure indicating the open file to parse for model parameter values.\\
\hline
\mbox{\tt in}  & {\em write\+\_\+geom} & If true, write grid geometry files.\\
\hline
\mbox{\tt in}  & {\em output\+\_\+dir} & The directory into which to write files. \\
\hline
\end{DoxyParams}


Definition at line 57 of file M\+O\+M\+\_\+fixed\+\_\+initialization.\+F90.


\begin{DoxyCode}
57   \textcolor{keywordtype}{type}(dyn\_horgrid\_type),  \textcolor{keywordtype}{intent(inout)} :: G\textcolor{comment}{    !< The ocean's grid structure.}
58   \textcolor{keywordtype}{type}(unit\_scale\_type),   \textcolor{keywordtype}{intent(in)}    :: US\textcolor{comment}{   !< A dimensional unit scaling type}
59   \textcolor{keywordtype}{type}(ocean\_OBC\_type),    \textcolor{keywordtype}{pointer}       :: OBC\textcolor{comment}{  !< Open boundary structure.}
60   \textcolor{keywordtype}{type}(param\_file\_type),   \textcolor{keywordtype}{intent(in)}    :: PF\textcolor{comment}{   !< A structure indicating the open file}
61 \textcolor{comment}{                                                 !! to parse for model parameter values.}
62   \textcolor{keywordtype}{logical},                 \textcolor{keywordtype}{intent(in)}    :: write\_geom\textcolor{comment}{ !< If true, write grid geometry files.}
63   \textcolor{keywordtype}{character(len=*)},        \textcolor{keywordtype}{intent(in)}    :: output\_dir\textcolor{comment}{ !< The directory into which to write files.}
64 
65   \textcolor{comment}{! Local}
66   \textcolor{keywordtype}{character(len=200)} :: inputdir   \textcolor{comment}{! The directory where NetCDF input files are.}
67   \textcolor{keywordtype}{character(len=200)} :: config
68   \textcolor{keywordtype}{character(len=40)}  :: mdl = \textcolor{stringliteral}{"MOM\_fixed\_initialization"} \textcolor{comment}{! This module's name.}
69   \textcolor{keywordtype}{logical} :: debug
70 \textcolor{comment}{! This include declares and sets the variable "version".}
71 \textcolor{preprocessor}{#include "version\_variable.h"}
72 \textcolor{preprocessor}{}
73   \textcolor{keyword}{call }calltree\_enter(\textcolor{stringliteral}{"MOM\_initialize\_fixed(), MOM\_fixed\_initialization.F90"})
74   \textcolor{keyword}{call }log\_version(pf, mdl, version, \textcolor{stringliteral}{""})
75   \textcolor{keyword}{call }get\_param(pf, mdl, \textcolor{stringliteral}{"DEBUG"}, debug, default=.false.)
76 
77   \textcolor{keyword}{call }get\_param(pf, mdl, \textcolor{stringliteral}{"INPUTDIR"}, inputdir, &
78          \textcolor{stringliteral}{"The directory in which input files are found."}, default=\textcolor{stringliteral}{"."})
79   inputdir = slasher(inputdir)
80 
81   \textcolor{comment}{! Set up the parameters of the physical domain (i.e. the grid), G}
82   \textcolor{keyword}{call }set\_grid\_metrics(g, pf, us)
83 
84   \textcolor{comment}{! Set up the bottom depth, G%bathyT either analytically or from file}
85   \textcolor{comment}{! This also sets G%max\_depth based on the input parameter MAXIMUM\_DEPTH,}
86   \textcolor{comment}{! or, if absent, is diagnosed as G%max\_depth = max( G%D(:,:) )}
87   \textcolor{keyword}{call }mom\_initialize\_topography(g%bathyT, g%max\_depth, g, pf, us)
88 \textcolor{comment}{!  call rescale\_dyn\_horgrid\_bathymetry(G, US%Z\_to\_m)}
89 
90   \textcolor{comment}{! To initialize masks, the bathymetry in halo regions must be filled in}
91   \textcolor{keyword}{call }pass\_var(g%bathyT, g%Domain)
92 
93   \textcolor{comment}{! Determine the position of any open boundaries}
94   \textcolor{keyword}{call }open\_boundary\_config(g, us, pf, obc)
95 
96   \textcolor{comment}{! Make bathymetry consistent with open boundaries}
97   \textcolor{keyword}{call }open\_boundary\_impose\_normal\_slope(obc, g, g%bathyT)
98 
99   \textcolor{comment}{! This call sets masks that prohibit flow over any point interpreted as land}
100   \textcolor{keyword}{call }initialize\_masks(g, pf, us)
101 
102   \textcolor{comment}{! Make OBC mask consistent with land mask}
103   \textcolor{keyword}{call }open\_boundary\_impose\_land\_mask(obc, g, g%areaCu, g%areaCv, us)
104 
105   \textcolor{keywordflow}{if} (debug) \textcolor{keywordflow}{then}
106     \textcolor{keyword}{call }hchksum(g%bathyT, \textcolor{stringliteral}{'MOM\_initialize\_fixed: depth '}, g%HI, haloshift=1, scale=us%Z\_to\_m)
107     \textcolor{keyword}{call }hchksum(g%mask2dT, \textcolor{stringliteral}{'MOM\_initialize\_fixed: mask2dT '}, g%HI)
108     \textcolor{keyword}{call }uvchksum(\textcolor{stringliteral}{'MOM\_initialize\_fixed: mask2dC[uv]'}, g%mask2dCu, &
109                   g%mask2dCv, g%HI)
110     \textcolor{keyword}{call }qchksum(g%mask2dBu, \textcolor{stringliteral}{'MOM\_initialize\_fixed: mask2dBu '}, g%HI)
111 \textcolor{keywordflow}{  endif}
112 
113   \textcolor{comment}{! Modulate geometric scales according to geography.}
114   \textcolor{keyword}{call }get\_param(pf, mdl, \textcolor{stringliteral}{"CHANNEL\_CONFIG"}, config, &
115                  \textcolor{stringliteral}{"A parameter that determines which set of channels are \(\backslash\)n"}//&
116                  \textcolor{stringliteral}{"restricted to specific  widths.  Options are:\(\backslash\)n"}//&
117                  \textcolor{stringliteral}{" \(\backslash\)t none - All channels have the grid width.\(\backslash\)n"}//&
118                  \textcolor{stringliteral}{" \(\backslash\)t global\_1deg - Sets 16 specific channels appropriate \(\backslash\)n"}//&
119                  \textcolor{stringliteral}{" \(\backslash\)t\(\backslash\)t for a 1-degree model, as used in CM2G.\(\backslash\)n"}//&
120                  \textcolor{stringliteral}{" \(\backslash\)t list - Read the channel locations and widths from a \(\backslash\)n"}//&
121                  \textcolor{stringliteral}{" \(\backslash\)t\(\backslash\)t text file, like MOM\_channel\_list in the MOM\_SIS \(\backslash\)n"}//&
122                  \textcolor{stringliteral}{" \(\backslash\)t\(\backslash\)t test case.\(\backslash\)n"}//&
123                  \textcolor{stringliteral}{" \(\backslash\)t file - Read open face widths everywhere from a \(\backslash\)n"}//&
124                  \textcolor{stringliteral}{" \(\backslash\)t\(\backslash\)t NetCDF file on the model grid."}, &
125                  default=\textcolor{stringliteral}{"none"})
126   \textcolor{keywordflow}{select case} ( trim(config) )
127     \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"none"})
128     \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"list"}) ; \textcolor{keyword}{call }reset\_face\_lengths\_list(g, pf, us)
129     \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"file"}) ; \textcolor{keyword}{call }reset\_face\_lengths\_file(g, pf, us)
130     \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"global\_1deg"}) ; \textcolor{keyword}{call }reset\_face\_lengths\_named(g, pf, trim(config), us)
131 \textcolor{keywordflow}{    case default} ; \textcolor{keyword}{call }mom\_error(fatal, \textcolor{stringliteral}{"MOM\_initialize\_fixed: "}// &
132       \textcolor{stringliteral}{"Unrecognized channel configuration "}//trim(config))
133 \textcolor{keywordflow}{  end select}
134 
135   \textcolor{comment}{!   This call sets the topography at velocity points.}
136   \textcolor{keywordflow}{if} (g%bathymetry\_at\_vel) \textcolor{keywordflow}{then}
137     \textcolor{keyword}{call }get\_param(pf, mdl, \textcolor{stringliteral}{"VELOCITY\_DEPTH\_CONFIG"}, config, &
138                    \textcolor{stringliteral}{"A string that determines how the topography is set at "}//&
139                    \textcolor{stringliteral}{"velocity points. This may be 'min' or 'max'."}, &
140                    default=\textcolor{stringliteral}{"max"})
141     \textcolor{keywordflow}{select case} ( trim(config) )
142       \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"max"}) ; \textcolor{keyword}{call }set\_velocity\_depth\_max(g)
143       \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"min"}) ; \textcolor{keyword}{call }set\_velocity\_depth\_min(g)
144 \textcolor{keywordflow}{      case default} ; \textcolor{keyword}{call }mom\_error(fatal, \textcolor{stringliteral}{"MOM\_initialize\_fixed: "}// &
145         \textcolor{stringliteral}{"Unrecognized velocity depth configuration "}//trim(config))
146 \textcolor{keywordflow}{    end select}
147 \textcolor{keywordflow}{  endif}
148 
149 \textcolor{comment}{!    Calculate the value of the Coriolis parameter at the latitude   !}
150 \textcolor{comment}{!  of the q grid points [s-1].}
151   \textcolor{keyword}{call }mom\_initialize\_rotation(g%CoriolisBu, g, pf, us=us)
152 \textcolor{comment}{!   Calculate the components of grad f (beta)}
153   \textcolor{keyword}{call }mom\_calculate\_grad\_coriolis(g%dF\_dx, g%dF\_dy, g, us=us)
154   \textcolor{keywordflow}{if} (debug) \textcolor{keywordflow}{then}
155     \textcolor{keyword}{call }qchksum(g%CoriolisBu, \textcolor{stringliteral}{"MOM\_initialize\_fixed: f "}, g%HI, scale=us%s\_to\_T)
156     \textcolor{keyword}{call }hchksum(g%dF\_dx, \textcolor{stringliteral}{"MOM\_initialize\_fixed: dF\_dx "}, g%HI, scale=us%m\_to\_L*us%s\_to\_T)
157     \textcolor{keyword}{call }hchksum(g%dF\_dy, \textcolor{stringliteral}{"MOM\_initialize\_fixed: dF\_dy "}, g%HI, scale=us%m\_to\_L*us%s\_to\_T)
158 \textcolor{keywordflow}{  endif}
159 
160   \textcolor{keyword}{call }initialize\_grid\_rotation\_angle(g, pf)
161 
162 \textcolor{comment}{! Compute global integrals of grid values for later use in scalar diagnostics !}
163   \textcolor{keyword}{call }compute\_global\_grid\_integrals(g, us=us)
164 
165 \textcolor{comment}{! Write out all of the grid data used by this run.}
166   \textcolor{keywordflow}{if} (write\_geom) \textcolor{keyword}{call }write\_ocean\_geometry\_file(g, pf, output\_dir, us=us)
167 
168   \textcolor{keyword}{call }calltree\_leave(\textcolor{stringliteral}{'MOM\_initialize\_fixed()'})
169 
\end{DoxyCode}
\mbox{\Hypertarget{namespacemom__fixed__initialization_ae980c841ada6296d07df08442398e87b}\label{namespacemom__fixed__initialization_ae980c841ada6296d07df08442398e87b}} 
\index{mom\+\_\+fixed\+\_\+initialization@{mom\+\_\+fixed\+\_\+initialization}!mom\+\_\+initialize\+\_\+topography@{mom\+\_\+initialize\+\_\+topography}}
\index{mom\+\_\+initialize\+\_\+topography@{mom\+\_\+initialize\+\_\+topography}!mom\+\_\+fixed\+\_\+initialization@{mom\+\_\+fixed\+\_\+initialization}}
\subsubsection{\texorpdfstring{mom\+\_\+initialize\+\_\+topography()}{mom\_initialize\_topography()}}
{\footnotesize\ttfamily subroutine, public mom\+\_\+fixed\+\_\+initialization\+::mom\+\_\+initialize\+\_\+topography (\begin{DoxyParamCaption}\item[{real, dimension(g\%isd\+:g\%ied,g\%jsd\+:g\%jed), intent(out)}]{D,  }\item[{real, intent(out)}]{max\+\_\+depth,  }\item[{type(dyn\+\_\+horgrid\+\_\+type), intent(in)}]{G,  }\item[{type(param\+\_\+file\+\_\+type), intent(in)}]{PF,  }\item[{type(unit\+\_\+scale\+\_\+type), intent(in), optional}]{US }\end{DoxyParamCaption})}



M\+O\+M\+\_\+initialize\+\_\+topography makes the appropriate call to set up the bathymetry. At this point the topography is in units of \mbox{[}m\mbox{]}, but this can be changed later. 


\begin{DoxyParams}[1]{Parameters}
\mbox{\tt in}  & {\em g} & The dynamic horizontal grid type\\
\hline
\mbox{\tt out}  & {\em d} & Ocean bottom depth \mbox{[}m\mbox{]}\\
\hline
\mbox{\tt in}  & {\em pf} & Parameter file structure\\
\hline
\mbox{\tt out}  & {\em max\+\_\+depth} & Maximum depth of model \mbox{[}m\mbox{]}\\
\hline
\mbox{\tt in}  & {\em us} & A dimensional unit scaling type \\
\hline
\end{DoxyParams}


Definition at line 175 of file M\+O\+M\+\_\+fixed\+\_\+initialization.\+F90.


\begin{DoxyCode}
175   \textcolor{keywordtype}{type}(dyn\_horgrid\_type),           \textcolor{keywordtype}{intent(in)}  :: G\textcolor{comment}{  !< The dynamic horizontal grid type}
176   \textcolor{keywordtype}{real}, \textcolor{keywordtype}{dimension(G%isd:G%ied,G%jsd:G%jed)}, &
177                                     \textcolor{keywordtype}{intent(out)} :: D\textcolor{comment}{  !< Ocean bottom depth [m]}
178   \textcolor{keywordtype}{type}(param\_file\_type),            \textcolor{keywordtype}{intent(in)}  :: PF\textcolor{comment}{ !< Parameter file structure}
179   \textcolor{keywordtype}{real},                             \textcolor{keywordtype}{intent(out)} :: max\_depth\textcolor{comment}{ !< Maximum depth of model [m]}
180   \textcolor{keywordtype}{type}(unit\_scale\_type),  \textcolor{keywordtype}{optional}, \textcolor{keywordtype}{intent(in)}  :: US\textcolor{comment}{ !< A dimensional unit scaling type}
181 
182   \textcolor{comment}{! This subroutine makes the appropriate call to set up the bottom depth.}
183   \textcolor{comment}{! This is a separate subroutine so that it can be made public and shared with}
184   \textcolor{comment}{! the ice-sheet code or other components.}
185 
186   \textcolor{comment}{! Local variables}
187   \textcolor{keywordtype}{real} :: m\_to\_Z, Z\_to\_m  \textcolor{comment}{! Dimensional rescaling factors}
188   \textcolor{keywordtype}{character(len=40)}  :: mdl = \textcolor{stringliteral}{"MOM\_initialize\_topography"} \textcolor{comment}{! This subroutine's name.}
189   \textcolor{keywordtype}{character(len=200)} :: config
190 
191   m\_to\_z = 1.0 ; \textcolor{keywordflow}{if} (\textcolor{keyword}{present}(us)) m\_to\_z = us%m\_to\_Z
192   z\_to\_m = 1.0 ; \textcolor{keywordflow}{if} (\textcolor{keyword}{present}(us)) z\_to\_m = us%Z\_to\_m
193 
194   \textcolor{keyword}{call }get\_param(pf, mdl, \textcolor{stringliteral}{"TOPO\_CONFIG"}, config, &
195                  \textcolor{stringliteral}{"This specifies how bathymetry is specified: \(\backslash\)n"}//&
196                  \textcolor{stringliteral}{" \(\backslash\)t file - read bathymetric information from the file \(\backslash\)n"}//&
197                  \textcolor{stringliteral}{" \(\backslash\)t\(\backslash\)t specified by (TOPO\_FILE).\(\backslash\)n"}//&
198                  \textcolor{stringliteral}{" \(\backslash\)t flat - flat bottom set to MAXIMUM\_DEPTH. \(\backslash\)n"}//&
199                  \textcolor{stringliteral}{" \(\backslash\)t bowl - an analytically specified bowl-shaped basin \(\backslash\)n"}//&
200                  \textcolor{stringliteral}{" \(\backslash\)t\(\backslash\)t ranging between MAXIMUM\_DEPTH and MINIMUM\_DEPTH. \(\backslash\)n"}//&
201                  \textcolor{stringliteral}{" \(\backslash\)t spoon - a similar shape to 'bowl', but with an vertical \(\backslash\)n"}//&
202                  \textcolor{stringliteral}{" \(\backslash\)t\(\backslash\)t wall at the southern face. \(\backslash\)n"}//&
203                  \textcolor{stringliteral}{" \(\backslash\)t halfpipe - a zonally uniform channel with a half-sine \(\backslash\)n"}//&
204                  \textcolor{stringliteral}{" \(\backslash\)t\(\backslash\)t profile in the meridional direction. \(\backslash\)n"}//&
205                  \textcolor{stringliteral}{" \(\backslash\)t bbuilder - build topography from list of functions. \(\backslash\)n"}//&
206                  \textcolor{stringliteral}{" \(\backslash\)t benchmark - use the benchmark test case topography. \(\backslash\)n"}//&
207                  \textcolor{stringliteral}{" \(\backslash\)t Neverworld - use the Neverworld test case topography. \(\backslash\)n"}//&
208                  \textcolor{stringliteral}{" \(\backslash\)t DOME - use a slope and channel configuration for the \(\backslash\)n"}//&
209                  \textcolor{stringliteral}{" \(\backslash\)t\(\backslash\)t DOME sill-overflow test case. \(\backslash\)n"}//&
210                  \textcolor{stringliteral}{" \(\backslash\)t ISOMIP - use a slope and channel configuration for the \(\backslash\)n"}//&
211                  \textcolor{stringliteral}{" \(\backslash\)t\(\backslash\)t ISOMIP test case. \(\backslash\)n"}//&
212                  \textcolor{stringliteral}{" \(\backslash\)t DOME2D - use a shelf and slope configuration for the \(\backslash\)n"}//&
213                  \textcolor{stringliteral}{" \(\backslash\)t\(\backslash\)t DOME2D gravity current/overflow test case. \(\backslash\)n"}//&
214                  \textcolor{stringliteral}{" \(\backslash\)t Kelvin - flat but with rotated land mask.\(\backslash\)n"}//&
215                  \textcolor{stringliteral}{" \(\backslash\)t seamount - Gaussian bump for spontaneous motion test case.\(\backslash\)n"}//&
216                  \textcolor{stringliteral}{" \(\backslash\)t dumbbell - Sloshing channel with reservoirs on both ends.\(\backslash\)n"}//&
217                  \textcolor{stringliteral}{" \(\backslash\)t shelfwave - exponential slope for shelfwave test case.\(\backslash\)n"}//&
218                  \textcolor{stringliteral}{" \(\backslash\)t Phillips - ACC-like idealized topography used in the Phillips config.\(\backslash\)n"}//&
219                  \textcolor{stringliteral}{" \(\backslash\)t dense - Denmark Strait-like dense water formation and overflow.\(\backslash\)n"}//&
220                  \textcolor{stringliteral}{" \(\backslash\)t USER - call a user modified routine."}, &
221                  fail\_if\_missing=.true.)
222   max\_depth = -1.e9*m\_to\_z ; \textcolor{keyword}{call }read\_param(pf, \textcolor{stringliteral}{"MAXIMUM\_DEPTH"}, max\_depth, scale=m\_to\_z)
223   \textcolor{keywordflow}{select case} ( trim(config) )
224     \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"file"});      \textcolor{keyword}{call }initialize\_topography\_from\_file(d, g, pf, us)
225     \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"flat"});      \textcolor{keyword}{call }initialize\_topography\_named(d, g, pf, config, max\_depth, us)
226     \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"spoon"});     \textcolor{keyword}{call }initialize\_topography\_named(d, g, pf, config, max\_depth, us)
227     \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"bowl"});      \textcolor{keyword}{call }initialize\_topography\_named(d, g, pf, config, max\_depth, us)
228     \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"halfpipe"});  \textcolor{keyword}{call }initialize\_topography\_named(d, g, pf, config, max\_depth, us)
229     \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"DOME"});      \textcolor{keyword}{call }dome\_initialize\_topography(d, g, pf, max\_depth, us)
230     \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"ISOMIP"});    \textcolor{keyword}{call }isomip\_initialize\_topography(d, g, pf, max\_depth, us)
231     \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"bbuilder"});  \textcolor{keyword}{call }basin\_builder\_topography(d, g, pf, max\_depth)
232     \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"benchmark"}); \textcolor{keyword}{call }benchmark\_initialize\_topography(d, g, pf, max\_depth, us)
233     \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"Neverworld"},\textcolor{stringliteral}{"Neverland"}); \textcolor{keyword}{call }neverworld\_initialize\_topography(d, g, pf, max\_depth)
234     \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"DOME2D"});    \textcolor{keyword}{call }dome2d\_initialize\_topography(d, g, pf, max\_depth)
235     \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"Kelvin"});    \textcolor{keyword}{call }kelvin\_initialize\_topography(d, g, pf, max\_depth, us)
236     \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"sloshing"});  \textcolor{keyword}{call }sloshing\_initialize\_topography(d, g, pf, max\_depth)
237     \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"seamount"});  \textcolor{keyword}{call }seamount\_initialize\_topography(d, g, pf, max\_depth)
238     \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"dumbbell"});  \textcolor{keyword}{call }dumbbell\_initialize\_topography(d, g, pf, max\_depth)
239     \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"shelfwave"}); \textcolor{keyword}{call }shelfwave\_initialize\_topography(d, g, pf, max\_depth, us)
240     \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"Phillips"});  \textcolor{keyword}{call }phillips\_initialize\_topography(d, g, pf, max\_depth, us)
241     \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"dense"});     \textcolor{keyword}{call }dense\_water\_initialize\_topography(d, g, pf, max\_depth)
242     \textcolor{keywordflow}{case} (\textcolor{stringliteral}{"USER"});      \textcolor{keyword}{call }user\_initialize\_topography(d, g, pf, max\_depth, us)
243 \textcolor{keywordflow}{    case default} ;      \textcolor{keyword}{call }mom\_error(fatal,\textcolor{stringliteral}{"MOM\_initialize\_topography: "}// &
244       \textcolor{stringliteral}{"Unrecognized topography setup '"}//trim(config)//\textcolor{stringliteral}{"'"})
245 \textcolor{keywordflow}{  end select}
246   \textcolor{keywordflow}{if} (max\_depth>0.) \textcolor{keywordflow}{then}
247     \textcolor{keyword}{call }log\_param(pf, mdl, \textcolor{stringliteral}{"MAXIMUM\_DEPTH"}, max\_depth*z\_to\_m, &
248                    \textcolor{stringliteral}{"The maximum depth of the ocean."}, units=\textcolor{stringliteral}{"m"})
249   \textcolor{keywordflow}{else}
250     max\_depth = diagnosemaximumdepth(d,g)
251     \textcolor{keyword}{call }log\_param(pf, mdl, \textcolor{stringliteral}{"!MAXIMUM\_DEPTH"}, max\_depth*z\_to\_m, &
252                    \textcolor{stringliteral}{"The (diagnosed) maximum depth of the ocean."}, units=\textcolor{stringliteral}{"m"}, like\_default=.true.)
253 \textcolor{keywordflow}{  endif}
254   \textcolor{keywordflow}{if} (trim(config) /= \textcolor{stringliteral}{"DOME"}) \textcolor{keywordflow}{then}
255     \textcolor{keyword}{call }limit\_topography(d, g, pf, max\_depth, us)
256 \textcolor{keywordflow}{  endif}
257 
\end{DoxyCode}
