s_examples/held_suarez_cs/inp_data.templ


%\subsubsection{File {\it input/data}}
%\label{www:tutorials}

This file, reproduced completely below, specifies the main parameters 
for the experiment. 
The parameters that are significant for this configuration are:

\begin{itemize}

\item Lines PUT_LINE_NB:tRef=,
\begin{verbatim} 
 tRef=295.2, 295.5, 295.9, 296.3, 296.7, 297.1, 297.6, 298.1, 298.7, 299.3,
\end{verbatim} 
$\cdots$ \\
set reference values for potential temperature (in Kelvin units) 
at each model level.
The entries are ordered like model level, from surface up to the top.
Density is calculated from anomalies at each level evaluated
with respect to the reference values set here.
\\ \fbox{
\begin{minipage}{5.0in}
{\it S/R INI\_THETA}({\it ini\_theta.F})
\end{minipage}
}


\item Line PUT_LINE_NB:no_slip_sides=,
\begin{verbatim}
 no_slip_sides=.FALSE.,
\end{verbatim}
this line selects a free-slip lateral boundary condition for
the horizontal Laplacian friction operator 
e.g. $\frac{\partial u}{\partial y}$=0 along boundaries in $y$ and
$\frac{\partial v}{\partial x}$=0 along boundaries in $x$.

\item Lines PUT_LINE_NB:no_slip_bottom=,
\begin{verbatim}
 no_slip_bottom=.FALSE.,
\end{verbatim}
this line selects a free-slip boundary condition at the top,
in the vertical Laplacian friction operator 
e.g. $\frac{\partial u}{\partial p} = \frac{\partial v}{\partial p} = 0$

\item Line PUT_LINE_NB:buoyancyRelation=,
\begin{verbatim}
 buoyancyRelation='ATMOSPHERIC',
\end{verbatim}
this line sets the type of fluid and the type of vertical coordinate to use,
which, in this case, is air with a pressure like coordinate ($p$ or $p^*$).

\item Line PUT_LINE_NB:eosType=,
\begin{verbatim}
 eosType='IDEALG',
\end{verbatim}
Selects the Ideal gas equation of state.
%\\ \fbox{
%\begin{minipage}{5.0in}
%{\it S/R FIND\_RHO}~({\it find\_rho.F})\\
%{\it S/R FIND\_ALPHA}~({\it find\_alpha.F})
%\end{minipage}
%}

\item Line PUT_LINE_NB:implicitFreeSurface=,
\begin{verbatim}
 implicitFreeSurface=.TRUE.,
\end{verbatim}
Selects the way the barotropic equation is solved, using here the implicit 
free-surface formulation.
\\ \fbox{
\begin{minipage}{5.0in}
{\it S/R SOLVE\_FOR\_PRESSURE}~({\it solve\_for\_pressure.F})
\end{minipage}
}

\item Line PUT_LINE_NB:exactConserv=,
\begin{verbatim}
 exactConserv=.TRUE.,
\end{verbatim}
Explicitly calculate again the surface pressure changes from 
the divergence of the vertically integrated horizontal flow,
after the implicit free surface solver and filters are applied.
\\ \fbox{
\begin{minipage}{5.0in}
{\it S/R INTEGR\_CONTINUITY}~({\it integr\_continuity.F})
\end{minipage}
}

\item Line PUT_LINE_NB:nonlinFreeSurf=
 and Line PUT_LINE_NB:select_rStar=,
\begin{verbatim}
 nonlinFreeSurf=4,
 select_rStar=2,
\end{verbatim}
Select the Non-Linear free surface formulation, using $r^*$ vertical coordinate
(here $p^*$).
Note that, except for the default ($= 0$), other values of those 2 parameters
are only permitted for testing/debuging purpose.
\\ \fbox{
\begin{minipage}{5.0in}
{\it S/R CALC\_R\_STAR}~({\it calc\_r\_star.F})\\
{\it S/R UPDATE\_R\_STAR}~({\it update\_r\_star.F})
\end{minipage}
}

\item Line PUT_LINE_NB:uniformLin_PhiSurf=
\begin{verbatim}
 uniformLin_PhiSurf=.FALSE.,
\end{verbatim}
Select the linear relation between surface geopotential anomaly 
and surface pressure anomaly to be evaluated from 
$\frac{\partial \Phi_s}{\partial p_s} = 1/\rho(\theta_{Ref})$
(see section \ref{sec:phi-freesurface}).
Note that using the default (=TRUE), the constant $1/\rho_0$ is
used instead, and is not necessary consistent with other 
parts of the geopotential that relies on $\theta_{Ref}$.
\\ \fbox{
\begin{minipage}{5.0in}
{\it S/R INI\_LINEAR\_PHISURF}~({\it ini\_linear\_phisurf.F})
\end{minipage}
}

\item Line PUT_LINE_NB:saltStepping= and Line PUT_LINE_NB:momViscosity=
\begin{verbatim}
 saltStepping=.FALSE.,
 momViscosity=.FALSE.,
\end{verbatim}
Do not step forward Water vapour and do not compute viscous terms.
This allow to save some computer time.

\item Line PUT_LINE_NB:vectorInvariantMomentum=
\begin{verbatim}
 vectorInvariantMomentum=.TRUE.,
\end{verbatim}
Select the vector-invariant form to solve the momentum equation.
\\ \fbox{
\begin{minipage}{5.0in}
{\it S/R MOM\_VECINV}~({\it mom\_vecinv.F})
\end{minipage}
}

\item Line PUT_LINE_NB:staggerTimeStep=
\begin{verbatim}
 staggerTimeStep=.TRUE.,
\end{verbatim}
Select the staggered time-stepping (rather than syncronous time stepping).

\item Line PUT_LINE_NB:readBinaryPrec= and PUT_LINE_NB:writeBinaryPrec=
\begin{verbatim}
 readBinaryPrec=64,
 writeBinaryPrec=64,
\end{verbatim}
Sets format for reading binary input datasets and writing output fields to
use 64-bit representation for floating-point numbers.
\\ \fbox{
\begin{minipage}{5.0in}
{\it S/R READ\_WRITE\_FLD}~({\it read\_write\_fld.F})\\
{\it S/R READ\_WRITE\_REC}~({\it read\_write\_rec.F})
\end{minipage}
}

\item Line PUT_LINE_NB:cg2dMaxIters=,
\begin{verbatim}
 cg2dMaxIters=200,
\end{verbatim}
Sets maximum number of iterations the two-dimensional, conjugate
gradient solver will use, {\bf irrespective of convergence 
criteria being met}.
\\ \fbox{
\begin{minipage}{5.0in}
{\it S/R CG2D}~({\it cg2d.F})
\end{minipage}
}

\item Line PUT_LINE_NB:cg2dTargetResWunit=,
\begin{verbatim}
 cg2dTargetResWunit=1.E-17,
\end{verbatim}
Sets the tolerance (in units of $\omega$) which the two-dimensional, 
conjugate gradient solver will use to test for convergence in equation 
%- note: Description of Conjugate gradient method (& related params) is missing
%  in the mean time, substitute this eq ref:
\ref{eq:elliptic-backward-free-surface} %\ref{eq:eg-hs-congrad_2d_resid} 
to $1 \times 10^{-17} Pa/s$.
Solver will iterate until tolerance falls below this value or until the 
maximum number of solver iterations is reached.
\\ \fbox{
\begin{minipage}{5.0in}
{\it S/R CG2D}~({\it cg2d.F})
\end{minipage}
}

\item Line PUT_LINE_NB:deltaT=,
\begin{verbatim}
 deltaT=450.,
\end{verbatim}
Sets the timestep $\Delta t$ used in the model to
$450~{\rm s}$ ($= 1/8 {\rm h}$).
\\ \fbox{
\begin{minipage}{5.0in}
{\it S/R TIMESTEP}({\it timestep.F})\\
{\it S/R TIMESTEP\_TRACER}({\it timestep\_tracer.F})
\end{minipage}
}

\item Line PUT_LINE_NB:startTime=,
\begin{verbatim}
 startTime=124416000.,
\end{verbatim}
Sets the starting time, in seconds, for the model time counter.
A non-zero starting time requires to read the initial state
from a pickup file. By default the pickup file is named according 
to the integer number ({\it nIter0}) of time steps 
in the {\it startTime} value ($ nIter0 = startTime / deltaT $).

\item Line PUT_LINE_NB:#nTimeSteps=,
\begin{verbatim}
#nTimeSteps=69120,
\end{verbatim}
A commented out setting for the length of the simulation 
(in number of time-step) that corresponds to 1 year simulation.

\item Line PUT_LINE_NB:nTimeSteps= and PUT_LINE_NB:monitorFreq=,
\begin{verbatim}
 nTimeSteps=16,
 monitorFreq=1.,
\end{verbatim}
Sets the length of the simulation (in number of time-step)
and the frequency (in seconds) for "monitor" output.
to 16 iterations and 1 seconds respectively. This choice
corresponds to a short simulation test.

\item Line PUT_LINE_NB:pChkptFreq=,
\begin{verbatim}
 pChkptFreq=31104000.,
\end{verbatim}
Sets the time interval, in seconds, bewteen 2 consecutive 
"permanent" pickups ("permanent checkpoint frequency")
that are used to restart the simuilation, to 1 year.

\item Line PUT_LINE_NB:chkptFreq=,
\begin{verbatim}
 chkptFreq=2592000.,
\end{verbatim}
Sets the time interval, in seconds, bewteen 2 consecutive 
"temporary" pickups ("checkpoint frequency") to 1 month. 
The "temporary" pickup file name is alternatively "ckptA" 
and "ckptB" ; thoses pickup (as opposed to the permanent ones)
are designed to be over-written by the model as the simulation
progresses.

\item Line PUT_LINE_NB:dumpFreq=,
\begin{verbatim}
 dumpFreq=2592000.,
\end{verbatim}
Set the frequencies (in seconds) for the snap-shot output
to 1 month.

\item Line PUT_LINE_NB:#monitorFreq=,
\begin{verbatim}
#monitorFreq=43200.,
\end{verbatim}
A commented out line setting the frequency (in seconds) for the 
"monitor" output to 12.h. This frequency fits 
better the longer simulation of 1 year.

\item Line PUT_LINE_NB:usingCurvilinearGrid=,
\begin{verbatim}
 usingCurvilinearGrid=.TRUE.,
\end{verbatim}
Set the horizontal type of grid to Curvilinear-Grid.

\item Line PUT_LINE_NB:horizGridFile=,
\begin{verbatim}
 horizGridFile='grid_cs32',
\end{verbatim}
Set the root for the grid file name to "{\it grid\_cs32}".
The grid-file names are derived from the root, adding a 
suffix with the face number (e.g.: {\it .face001.bin}, 
{\it .face002.bin} $\cdots$ )
\\ \fbox{
\begin{minipage}{5.0in}
{\it S/R INI\_CURVILINEAR\_GRID}~({\it ini\_curvilinear\_grid.F})
\end{minipage}
}

\item Lines PUT_LINE_NB:delR= and PUT_LINE_NB:Ro_SeaLevel=,
\begin{verbatim}
 delR=20*50.E2,
 Ro_SeaLevel=1.E5,
\end{verbatim}
Those 2 lines define the vertical discretization, in pressure units.
The $1^{rst}$ one sets the increments in pressure units (Pa),
to 20 equally thick levels of $50 \times 10^2 {\rm Pa}$ each.
The $2^{nd}$ one sets the reference pressure at the sea-level, 
to $10^5 {\rm Pa}$. This define the origin (interface $k=1$) 
of the vertical pressure axis, with decreasing pressure 
as the level index $k$ increases.
\\ \fbox{
\begin{minipage}{5.0in}
{\it S/R INI\_VERTICAL\_GRID}~({\it ini\_vertical\_grid.F})
\end{minipage}
}

\item Line PUT_LINE_NB:#topoFile=,
\begin{verbatim}
#topoFile='topo.cs.bin'
\end{verbatim}
This commented out line would allow to set the file name 
of a 2-D orography file, in meters units, to '{\it topo.cs.bin}'.
\\ \fbox{
\begin{minipage}{5.0in}
{\it S/R INI\_DEPTH}~({\it ini\_depth.F})
\end{minipage}
}

\end{itemize}

\noindent other lines in the file {\it input/data} are standard values
that are described in the MITgcm Getting Started and MITgcm Parameters
notes.

\begin{small}
\input{s_examples/held_suarez_cs/input/data}
\end{small}
1
2	%\subsubsection{File {\it input/data}}
3	%\label{www:tutorials}
4
5	This file, reproduced completely below, specifies the main parameters
6	for the experiment.
7	The parameters that are significant for this configuration are:
8
9	\begin{itemize}
10
11	\item Lines PUT_LINE_NB:tRef=,
12	\begin{verbatim}
13	tRef=295.2, 295.5, 295.9, 296.3, 296.7, 297.1, 297.6, 298.1, 298.7, 299.3,
14	\end{verbatim}
15	$\cdots$ \\
16	set reference values for potential temperature (in Kelvin units)
17	at each model level.
18	The entries are ordered like model level, from surface up to the top.
19	Density is calculated from anomalies at each level evaluated
20	with respect to the reference values set here.
21	\\ \fbox{
22	\begin{minipage}{5.0in}
23	{\it S/R INI\_THETA}({\it ini\_theta.F})
24	\end{minipage}
25	}
26
27
28	\item Line PUT_LINE_NB:no_slip_sides=,
29	\begin{verbatim}
30	no_slip_sides=.FALSE.,
31	\end{verbatim}
32	this line selects a free-slip lateral boundary condition for
33	the horizontal Laplacian friction operator
34	e.g. $\frac{\partial u}{\partial y}$=0 along boundaries in $y$ and
35	$\frac{\partial v}{\partial x}$=0 along boundaries in $x$.
36
37	\item Lines PUT_LINE_NB:no_slip_bottom=,
38	\begin{verbatim}
39	no_slip_bottom=.FALSE.,
40	\end{verbatim}
41	this line selects a free-slip boundary condition at the top,
42	in the vertical Laplacian friction operator
43	e.g. $\frac{\partial u}{\partial p} = \frac{\partial v}{\partial p} = 0$
44
45	\item Line PUT_LINE_NB:buoyancyRelation=,
46	\begin{verbatim}
47	buoyancyRelation='ATMOSPHERIC',
48	\end{verbatim}
49	this line sets the type of fluid and the type of vertical coordinate to use,
50	which, in this case, is air with a pressure like coordinate ($p$ or $p^*$).
51
52	\item Line PUT_LINE_NB:eosType=,
53	\begin{verbatim}
54	eosType='IDEALG',
55	\end{verbatim}
56	Selects the Ideal gas equation of state.
57	%\\ \fbox{
58	%\begin{minipage}{5.0in}
59	%{\it S/R FIND\_RHO}~({\it find\_rho.F})\\
60	%{\it S/R FIND\_ALPHA}~({\it find\_alpha.F})
61	%\end{minipage}
62	%}
63
64	\item Line PUT_LINE_NB:implicitFreeSurface=,
65	\begin{verbatim}
66	implicitFreeSurface=.TRUE.,
67	\end{verbatim}
68	Selects the way the barotropic equation is solved, using here the implicit
69	free-surface formulation.
70	\\ \fbox{
71	\begin{minipage}{5.0in}
72	{\it S/R SOLVE\_FOR\_PRESSURE}~({\it solve\_for\_pressure.F})
73	\end{minipage}
74	}
75
76	\item Line PUT_LINE_NB:exactConserv=,
77	\begin{verbatim}
78	exactConserv=.TRUE.,
79	\end{verbatim}
80	Explicitly calculate again the surface pressure changes from
81	the divergence of the vertically integrated horizontal flow,
82	after the implicit free surface solver and filters are applied.
83	\\ \fbox{
84	\begin{minipage}{5.0in}
85	{\it S/R INTEGR\_CONTINUITY}~({\it integr\_continuity.F})
86	\end{minipage}
87	}
88
89	\item Line PUT_LINE_NB:nonlinFreeSurf=
90	and Line PUT_LINE_NB:select_rStar=,
91	\begin{verbatim}
92	nonlinFreeSurf=4,
93	select_rStar=2,
94	\end{verbatim}
95	Select the Non-Linear free surface formulation, using $r^*$ vertical coordinate
96	(here $p^*$).
97	Note that, except for the default ($= 0$), other values of those 2 parameters
98	are only permitted for testing/debuging purpose.
99	\\ \fbox{
100	\begin{minipage}{5.0in}
101	{\it S/R CALC\_R\_STAR}~({\it calc\_r\_star.F})\\
102	{\it S/R UPDATE\_R\_STAR}~({\it update\_r\_star.F})
103	\end{minipage}
104	}
105
106	\item Line PUT_LINE_NB:uniformLin_PhiSurf=
107	\begin{verbatim}
108	uniformLin_PhiSurf=.FALSE.,
109	\end{verbatim}
110	Select the linear relation between surface geopotential anomaly
111	and surface pressure anomaly to be evaluated from
112	$\frac{\partial \Phi_s}{\partial p_s} = 1/\rho(\theta_{Ref})$
113	(see section \ref{sec:phi-freesurface}).
114	Note that using the default (=TRUE), the constant $1/\rho_0$ is
115	used instead, and is not necessary consistent with other
116	parts of the geopotential that relies on $\theta_{Ref}$.
117	\\ \fbox{
118	\begin{minipage}{5.0in}
119	{\it S/R INI\_LINEAR\_PHISURF}~({\it ini\_linear\_phisurf.F})
120	\end{minipage}
121	}
122
123	\item Line PUT_LINE_NB:saltStepping= and Line PUT_LINE_NB:momViscosity=
124	\begin{verbatim}
125	saltStepping=.FALSE.,
126	momViscosity=.FALSE.,
127	\end{verbatim}
128	Do not step forward Water vapour and do not compute viscous terms.
129	This allow to save some computer time.
130
131	\item Line PUT_LINE_NB:vectorInvariantMomentum=
132	\begin{verbatim}
133	vectorInvariantMomentum=.TRUE.,
134	\end{verbatim}
135	Select the vector-invariant form to solve the momentum equation.
136	\\ \fbox{
137	\begin{minipage}{5.0in}
138	{\it S/R MOM\_VECINV}~({\it mom\_vecinv.F})
139	\end{minipage}
140	}
141
142	\item Line PUT_LINE_NB:staggerTimeStep=
143	\begin{verbatim}
144	staggerTimeStep=.TRUE.,
145	\end{verbatim}
146	Select the staggered time-stepping (rather than syncronous time stepping).
147
148	\item Line PUT_LINE_NB:readBinaryPrec= and PUT_LINE_NB:writeBinaryPrec=
149	\begin{verbatim}
150	readBinaryPrec=64,
151	writeBinaryPrec=64,
152	\end{verbatim}
153	Sets format for reading binary input datasets and writing output fields to
154	use 64-bit representation for floating-point numbers.
155	\\ \fbox{
156	\begin{minipage}{5.0in}
157	{\it S/R READ\_WRITE\_FLD}~({\it read\_write\_fld.F})\\
158	{\it S/R READ\_WRITE\_REC}~({\it read\_write\_rec.F})
159	\end{minipage}
160	}
161
162	\item Line PUT_LINE_NB:cg2dMaxIters=,
163	\begin{verbatim}
164	cg2dMaxIters=200,
165	\end{verbatim}
166	Sets maximum number of iterations the two-dimensional, conjugate
167	gradient solver will use, {\bf irrespective of convergence
168	criteria being met}.
169	\\ \fbox{
170	\begin{minipage}{5.0in}
171	{\it S/R CG2D}~({\it cg2d.F})
172	\end{minipage}
173	}
174
175	\item Line PUT_LINE_NB:cg2dTargetResWunit=,
176	\begin{verbatim}
177	cg2dTargetResWunit=1.E-17,
178	\end{verbatim}
179	Sets the tolerance (in units of $\omega$) which the two-dimensional,
180	conjugate gradient solver will use to test for convergence in equation
181	%- note: Description of Conjugate gradient method (& related params) is missing
182	% in the mean time, substitute this eq ref:
183	\ref{eq:elliptic-backward-free-surface} %\ref{eq:eg-hs-congrad_2d_resid}
184	to $1 \times 10^{-17} Pa/s$.
185	Solver will iterate until tolerance falls below this value or until the
186	maximum number of solver iterations is reached.
187	\\ \fbox{
188	\begin{minipage}{5.0in}
189	{\it S/R CG2D}~({\it cg2d.F})
190	\end{minipage}
191	}
192
193	\item Line PUT_LINE_NB:deltaT=,
194	\begin{verbatim}
195	deltaT=450.,
196	\end{verbatim}
197	Sets the timestep $\Delta t$ used in the model to
198	$450~{\rm s}$ ($= 1/8 {\rm h}$).
199	\\ \fbox{
200	\begin{minipage}{5.0in}
201	{\it S/R TIMESTEP}({\it timestep.F})\\
202	{\it S/R TIMESTEP\_TRACER}({\it timestep\_tracer.F})
203	\end{minipage}
204	}
205
206	\item Line PUT_LINE_NB:startTime=,
207	\begin{verbatim}
208	startTime=124416000.,
209	\end{verbatim}
210	Sets the starting time, in seconds, for the model time counter.
211	A non-zero starting time requires to read the initial state
212	from a pickup file. By default the pickup file is named according
213	to the integer number ({\it nIter0}) of time steps
214	in the {\it startTime} value ($ nIter0 = startTime / deltaT $).
215
216	\item Line PUT_LINE_NB:#nTimeSteps=,
217	\begin{verbatim}
218	#nTimeSteps=69120,
219	\end{verbatim}
220	A commented out setting for the length of the simulation
221	(in number of time-step) that corresponds to 1 year simulation.
222
223	\item Line PUT_LINE_NB:nTimeSteps= and PUT_LINE_NB:monitorFreq=,
224	\begin{verbatim}
225	nTimeSteps=16,
226	monitorFreq=1.,
227	\end{verbatim}
228	Sets the length of the simulation (in number of time-step)
229	and the frequency (in seconds) for "monitor" output.
230	to 16 iterations and 1 seconds respectively. This choice
231	corresponds to a short simulation test.
232
233	\item Line PUT_LINE_NB:pChkptFreq=,
234	\begin{verbatim}
235	pChkptFreq=31104000.,
236	\end{verbatim}
237	Sets the time interval, in seconds, bewteen 2 consecutive
238	"permanent" pickups ("permanent checkpoint frequency")
239	that are used to restart the simuilation, to 1 year.
240
241	\item Line PUT_LINE_NB:chkptFreq=,
242	\begin{verbatim}
243	chkptFreq=2592000.,
244	\end{verbatim}
245	Sets the time interval, in seconds, bewteen 2 consecutive
246	"temporary" pickups ("checkpoint frequency") to 1 month.
247	The "temporary" pickup file name is alternatively "ckptA"
248	and "ckptB" ; thoses pickup (as opposed to the permanent ones)
249	are designed to be over-written by the model as the simulation
250	progresses.
251
252	\item Line PUT_LINE_NB:dumpFreq=,
253	\begin{verbatim}
254	dumpFreq=2592000.,
255	\end{verbatim}
256	Set the frequencies (in seconds) for the snap-shot output
257	to 1 month.
258
259	\item Line PUT_LINE_NB:#monitorFreq=,
260	\begin{verbatim}
261	#monitorFreq=43200.,
262	\end{verbatim}
263	A commented out line setting the frequency (in seconds) for the
264	"monitor" output to 12.h. This frequency fits
265	better the longer simulation of 1 year.
266
267	\item Line PUT_LINE_NB:usingCurvilinearGrid=,
268	\begin{verbatim}
269	usingCurvilinearGrid=.TRUE.,
270	\end{verbatim}
271	Set the horizontal type of grid to Curvilinear-Grid.
272
273	\item Line PUT_LINE_NB:horizGridFile=,
274	\begin{verbatim}
275	horizGridFile='grid_cs32',
276	\end{verbatim}
277	Set the root for the grid file name to "{\it grid\_cs32}".
278	The grid-file names are derived from the root, adding a
279	suffix with the face number (e.g.: {\it .face001.bin},
280	{\it .face002.bin} $\cdots$ )
281	\\ \fbox{
282	\begin{minipage}{5.0in}
283	{\it S/R INI\_CURVILINEAR\_GRID}~({\it ini\_curvilinear\_grid.F})
284	\end{minipage}
285	}
286
287	\item Lines PUT_LINE_NB:delR= and PUT_LINE_NB:Ro_SeaLevel=,
288	\begin{verbatim}
289	delR=20*50.E2,
290	Ro_SeaLevel=1.E5,
291	\end{verbatim}
292	Those 2 lines define the vertical discretization, in pressure units.
293	The $1^{rst}$ one sets the increments in pressure units (Pa),
294	to 20 equally thick levels of $50 \times 10^2 {\rm Pa}$ each.
295	The $2^{nd}$ one sets the reference pressure at the sea-level,
296	to $10^5 {\rm Pa}$. This define the origin (interface $k=1$)
297	of the vertical pressure axis, with decreasing pressure
298	as the level index $k$ increases.
299	\\ \fbox{
300	\begin{minipage}{5.0in}
301	{\it S/R INI\_VERTICAL\_GRID}~({\it ini\_vertical\_grid.F})
302	\end{minipage}
303	}
304
305	\item Line PUT_LINE_NB:#topoFile=,
306	\begin{verbatim}
307	#topoFile='topo.cs.bin'
308	\end{verbatim}
309	This commented out line would allow to set the file name
310	of a 2-D orography file, in meters units, to '{\it topo.cs.bin}'.
311	\\ \fbox{
312	\begin{minipage}{5.0in}
313	{\it S/R INI\_DEPTH}~({\it ini\_depth.F})
314	\end{minipage}
315	}
316
317	\end{itemize}
318
319	\noindent other lines in the file {\it input/data} are standard values
320	that are described in the MITgcm Getting Started and MITgcm Parameters
321	notes.
322
323	\begin{small}
324	\input{s_examples/held_suarez_cs/input/data}
325	\end{small}