1 |
\section{Surface Driven Convection} |
\section{Surface Driven Convection} |
2 |
|
\label{www:tutorials} |
3 |
\label{sect:eg-bconv} |
\label{sect:eg-bconv} |
4 |
|
|
5 |
\bodytext{bgcolor="#FFFFFFFF"} |
\bodytext{bgcolor="#FFFFFFFF"} |
40 |
\end{itemize} |
\end{itemize} |
41 |
|
|
42 |
\subsection{Overview} |
\subsection{Overview} |
43 |
|
\label{www:tutorials} |
44 |
|
|
45 |
The model domain consists of an approximately 3 |
The model domain consists of an approximately 3 |
46 |
km square by 1 km deep box of initially |
km square by 1 km deep box of initially |
110 |
in a binary data file generated using the Matlab script {\it input/gendata.m}. |
in a binary data file generated using the Matlab script {\it input/gendata.m}. |
111 |
|
|
112 |
\subsection{Equations solved} |
\subsection{Equations solved} |
113 |
|
\label{www:tutorials} |
114 |
|
|
115 |
The model is configured in nonhydrostatic form, that is, all terms in the Navier |
The model is configured in nonhydrostatic form, that is, all terms in the Navier |
116 |
Stokes equations are retained and the pressure field is found, subject to appropriate |
Stokes equations are retained and the pressure field is found, subject to appropriate |
180 |
\\ |
\\ |
181 |
|
|
182 |
\subsection{Discrete numerical configuration} |
\subsection{Discrete numerical configuration} |
183 |
|
\label{www:tutorials} |
184 |
|
|
185 |
The domain is discretised with a uniform grid spacing in each direction. There are 64 |
The domain is discretised with a uniform grid spacing in each direction. There are 64 |
186 |
grid cells in directions $x$ and $y$ and 20 vertical levels thus the domain |
grid cells in directions $x$ and $y$ and 20 vertical levels thus the domain |
187 |
comprises a total of just over 80 000 gridpoints. |
comprises a total of just over 80 000 gridpoints. |
188 |
|
|
189 |
\subsection{Numerical stability criteria and other considerations} |
\subsection{Numerical stability criteria and other considerations} |
190 |
|
\label{www:tutorials} |
191 |
|
|
192 |
For a heat flux of 800 Wm$^2$ and a rotation rate of $10^{-4}$ s$^{-1}$ the |
For a heat flux of 800 Wm$^2$ and a rotation rate of $10^{-4}$ s$^{-1}$ the |
193 |
plume-scale can be expected to be a few hundred meters guiding our choice of grid |
plume-scale can be expected to be a few hundred meters guiding our choice of grid |
213 |
correlated over 50 m. |
correlated over 50 m. |
214 |
|
|
215 |
\subsection{Experiment configuration} |
\subsection{Experiment configuration} |
216 |
|
\label{www:tutorials} |
217 |
|
|
218 |
The model configuration for this experiment resides under the directory |
The model configuration for this experiment resides under the directory |
219 |
{\it verification/convection/}. The experiment files |
{\it verification/convection/}. The experiment files |
230 |
experiment. Below we describe these experiment-specific customisations. |
experiment. Below we describe these experiment-specific customisations. |
231 |
|
|
232 |
\subsubsection{File {\it code/CPP\_EEOPTIONS.h}} |
\subsubsection{File {\it code/CPP\_EEOPTIONS.h}} |
233 |
|
\label{www:tutorials} |
234 |
|
|
235 |
This file uses standard default values and does not contain |
This file uses standard default values and does not contain |
236 |
customisations for this experiment. |
customisations for this experiment. |
237 |
|
|
238 |
\subsubsection{File {\it code/CPP\_OPTIONS.h}} |
\subsubsection{File {\it code/CPP\_OPTIONS.h}} |
239 |
|
\label{www:tutorials} |
240 |
|
|
241 |
This file uses standard default values and does not contain |
This file uses standard default values and does not contain |
242 |
customisations for this experiment. |
customisations for this experiment. |
243 |
|
|
244 |
\subsubsection{File {\it code/SIZE.h}} |
\subsubsection{File {\it code/SIZE.h}} |
245 |
|
\label{www:tutorials} |
246 |
|
|
247 |
Three lines are customized in this file. These prescribe the domain grid dimensions. |
Three lines are customized in this file. These prescribe the domain grid dimensions. |
248 |
\begin{itemize} |
\begin{itemize} |
270 |
\begin{rawhtml}</PRE>\end{rawhtml} |
\begin{rawhtml}</PRE>\end{rawhtml} |
271 |
|
|
272 |
\subsubsection{File {\it input/data}} |
\subsubsection{File {\it input/data}} |
273 |
|
\label{www:tutorials} |
274 |
|
|
275 |
This file, reproduced completely below, specifies the main parameters |
This file, reproduced completely below, specifies the main parameters |
276 |
for the experiment. The parameters that are significant for this configuration |
for the experiment. The parameters that are significant for this configuration |
290 |
$x$ and $y$. The values specified are read into the |
$x$ and $y$. The values specified are read into the |
291 |
variable |
variable |
292 |
{\bf |
{\bf |
293 |
\begin{rawhtml} <A href=../../../code_reference/vdb/names/OK.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/names/OK.htm> \end{rawhtml} |
294 |
tRef |
tRef |
295 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
296 |
} |
} |
297 |
in the model code, by procedure |
in the model code, by procedure |
298 |
{\it |
{\it |
299 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
300 |
S/R INI\_PARMS ({\it ini\_parms.F}) |
S/R INI\_PARMS ({\it ini\_parms.F}) |
301 |
\begin{rawhtml} </A>\end{rawhtml}. |
\begin{rawhtml} </A>\end{rawhtml}. |
302 |
} |
} |
303 |
The temperature field is initialised, by procedure |
The temperature field is initialised, by procedure |
304 |
{\it |
{\it |
305 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/OK.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/OK.htm> \end{rawhtml} |
306 |
S/R INI\_THETA ({\it ini\_theta.F}) |
S/R INI\_THETA ({\it ini\_theta.F}) |
307 |
\begin{rawhtml} </A>\end{rawhtml}. |
\begin{rawhtml} </A>\end{rawhtml}. |
308 |
} |
} |
320 |
$x$ and $y$. The values specified are read into the |
$x$ and $y$. The values specified are read into the |
321 |
variable |
variable |
322 |
{\bf |
{\bf |
323 |
\begin{rawhtml} <A href=../../../code_reference/vdb/names/OK.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/names/OK.htm> \end{rawhtml} |
324 |
sRef |
sRef |
325 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
326 |
} |
} |
327 |
in the model code, by procedure |
in the model code, by procedure |
328 |
{\it |
{\it |
329 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
330 |
S/R INI\_PARMS ({\it ini\_parms.F}) |
S/R INI\_PARMS ({\it ini\_parms.F}) |
331 |
} |
} |
332 |
\begin{rawhtml} </A>\end{rawhtml}. |
\begin{rawhtml} </A>\end{rawhtml}. |
333 |
The salinity field is initialised, by procedure |
The salinity field is initialised, by procedure |
334 |
{\it |
{\it |
335 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/OK.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/OK.htm> \end{rawhtml} |
336 |
S/R INI\_SALT ({\it ini\_salt.F}) |
S/R INI\_SALT ({\it ini\_salt.F}) |
337 |
\begin{rawhtml} </A>\end{rawhtml}. |
\begin{rawhtml} </A>\end{rawhtml}. |
338 |
} |
} |
347 |
for this operator are specified later. |
for this operator are specified later. |
348 |
The variable |
The variable |
349 |
{\bf |
{\bf |
350 |
\begin{rawhtml} <A href=../../../code_reference/vdb/names/SI.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/names/SI.htm> \end{rawhtml} |
351 |
viscAh |
viscAh |
352 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
353 |
} |
} |
354 |
is read in the routine |
is read in the routine |
355 |
{\it |
{\it |
356 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
357 |
S/R INI\_PARMS ({\it ini\_params.F}) |
S/R INI\_PARMS ({\it ini\_params.F}) |
358 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
359 |
} and applied in routines |
} and applied in routines |
360 |
{\it |
{\it |
361 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
362 |
S/R CALC\_MOM\_RHS ({\it calc\_mom\_rhs.F}) |
S/R CALC\_MOM\_RHS ({\it calc\_mom\_rhs.F}) |
363 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
364 |
} and |
} and |
365 |
{\it |
{\it |
366 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
367 |
S/R CALC\_GW ({\it calc\_gw.F}) |
S/R CALC\_GW ({\it calc\_gw.F}) |
368 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
369 |
}. |
}. |
378 |
for this operator are specified later. |
for this operator are specified later. |
379 |
The variable |
The variable |
380 |
{\bf |
{\bf |
381 |
\begin{rawhtml} <A href=../../../code_reference/vdb/names/ZQ.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/names/ZQ.htm> \end{rawhtml} |
382 |
viscAz |
viscAz |
383 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
384 |
} |
} |
385 |
is read in the routine |
is read in the routine |
386 |
{\it |
{\it |
387 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
388 |
S/R INI\_PARMS ({\it ini\_parms.F}) |
S/R INI\_PARMS ({\it ini\_parms.F}) |
389 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
390 |
} |
} |
391 |
and is copied into model general vertical coordinate variable |
and is copied into model general vertical coordinate variable |
392 |
{\bf |
{\bf |
393 |
\begin{rawhtml} <A href=../../../code_reference/vdb/names/PF.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/names/PF.htm> \end{rawhtml} |
394 |
viscAr |
viscAr |
395 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
396 |
}. At each time step, the viscous term contribution to the momentum equations |
}. At each time step, the viscous term contribution to the momentum equations |
397 |
is calculated in routine |
is calculated in routine |
398 |
{\it |
{\it |
399 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
400 |
S/R CALC\_DIFFUSIVITY ({\it calc\_diffusivity.F}) |
S/R CALC\_DIFFUSIVITY ({\it calc\_diffusivity.F}) |
401 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
402 |
}. |
}. |
412 |
$\frac{\partial v}{\partial x}$=0 along boundaries in $x$. |
$\frac{\partial v}{\partial x}$=0 along boundaries in $x$. |
413 |
The variable |
The variable |
414 |
{\bf |
{\bf |
415 |
\begin{rawhtml} <A href=../../../code_reference/vdb/names/UT.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/names/UT.htm> \end{rawhtml} |
416 |
no\_slip\_sides |
no\_slip\_sides |
417 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
418 |
} |
} |
419 |
is read in the routine |
is read in the routine |
420 |
{\it |
{\it |
421 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
422 |
S/R INI\_PARMS ({\it ini\_parms.F}) |
S/R INI\_PARMS ({\it ini\_parms.F}) |
423 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
424 |
} and the boundary condition is evaluated in routine |
} and the boundary condition is evaluated in routine |
425 |
{\it |
{\it |
426 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
427 |
S/R CALC\_MOM\_RHS ({\it calc\_mom\_rhs.F}) |
S/R CALC\_MOM\_RHS ({\it calc\_mom\_rhs.F}) |
428 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
429 |
}. |
}. |
438 |
e.g. $u=v=0$ at $z=-H$, where $H$ is the local depth of the domain. |
e.g. $u=v=0$ at $z=-H$, where $H$ is the local depth of the domain. |
439 |
The variable |
The variable |
440 |
{\bf |
{\bf |
441 |
\begin{rawhtml} <A href=../../../code_reference/vdb/names/UK.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/names/UK.htm> \end{rawhtml} |
442 |
no\_slip\_bottom |
no\_slip\_bottom |
443 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
444 |
} |
} |
445 |
is read in the routine |
is read in the routine |
446 |
{\it |
{\it |
447 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
448 |
S/R INI\_PARMS ({\it ini\_parms.F}) |
S/R INI\_PARMS ({\it ini\_parms.F}) |
449 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
450 |
} and is applied in the routine |
} and is applied in the routine |
451 |
{\it |
{\it |
452 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
453 |
S/R CALC\_MOM\_RHS ({\it calc\_mom\_rhs.F}) |
S/R CALC\_MOM\_RHS ({\it calc\_mom\_rhs.F}) |
454 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
455 |
}. |
}. |
464 |
all boundaries. |
all boundaries. |
465 |
The variable |
The variable |
466 |
{\bf |
{\bf |
467 |
\begin{rawhtml} <A href=../../../code_reference/vdb/names/RC.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/names/RC.htm> \end{rawhtml} |
468 |
diffKhT |
diffKhT |
469 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
470 |
} |
} |
471 |
is read in the routine |
is read in the routine |
472 |
{\it |
{\it |
473 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
474 |
S/R INI\_PARMS ({\it ini\_parms.F}) |
S/R INI\_PARMS ({\it ini\_parms.F}) |
475 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
476 |
} and used in routine |
} and used in routine |
477 |
{\it |
{\it |
478 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
479 |
S/R CALC\_GT ({\it calc\_gt.F}) |
S/R CALC\_GT ({\it calc\_gt.F}) |
480 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
481 |
}. |
}. |
489 |
operator is $\frac{\partial}{\partial z}$ = 0 on all boundaries. |
operator is $\frac{\partial}{\partial z}$ = 0 on all boundaries. |
490 |
The variable |
The variable |
491 |
{\bf |
{\bf |
492 |
\begin{rawhtml} <A href=../../../code_reference/vdb/names/ZT.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/names/ZT.htm> \end{rawhtml} |
493 |
diffKzT |
diffKzT |
494 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
495 |
} |
} |
496 |
is read in the routine |
is read in the routine |
497 |
{\it |
{\it |
498 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
499 |
S/R INI\_PARMS ({\it ini\_parms.F}) |
S/R INI\_PARMS ({\it ini\_parms.F}) |
500 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
501 |
}. |
}. |
502 |
It is copied into model general vertical coordinate variable |
It is copied into model general vertical coordinate variable |
503 |
{\bf |
{\bf |
504 |
\begin{rawhtml} <A href=../../../code_reference/vdb/names/PD.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/names/PD.htm> \end{rawhtml} |
505 |
diffKrT |
diffKrT |
506 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
507 |
} which is used in routine |
} which is used in routine |
508 |
{\it |
{\it |
509 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
510 |
S/R CALC\_DIFFUSIVITY ({\it calc\_diffusivity.F}) |
S/R CALC\_DIFFUSIVITY ({\it calc\_diffusivity.F}) |
511 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
512 |
}. |
}. |
522 |
all boundaries. |
all boundaries. |
523 |
The variable |
The variable |
524 |
{\bf |
{\bf |
525 |
\begin{rawhtml} <A href=../../../code_reference/vdb/names/RC.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/names/RC.htm> \end{rawhtml} |
526 |
diffKsT |
diffKsT |
527 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
528 |
} |
} |
529 |
is read in the routine |
is read in the routine |
530 |
{\it |
{\it |
531 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
532 |
S/R INI\_PARMS ({\it ini\_parms.F}) |
S/R INI\_PARMS ({\it ini\_parms.F}) |
533 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
534 |
} and used in routine |
} and used in routine |
535 |
{\it |
{\it |
536 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
537 |
S/R CALC\_GS ({\it calc\_gs.F}) |
S/R CALC\_GS ({\it calc\_gs.F}) |
538 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
539 |
}. |
}. |
548 |
operator is $\frac{\partial}{\partial z}$ = 0 on all boundaries. |
operator is $\frac{\partial}{\partial z}$ = 0 on all boundaries. |
549 |
The variable |
The variable |
550 |
{\bf |
{\bf |
551 |
\begin{rawhtml} <A href=../../../code_reference/vdb/names/ZT.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/names/ZT.htm> \end{rawhtml} |
552 |
diffKzS |
diffKzS |
553 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
554 |
} |
} |
555 |
is read in the routine |
is read in the routine |
556 |
{\it |
{\it |
557 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
558 |
S/R INI\_PARMS ({\it ini\_parms.F}) |
S/R INI\_PARMS ({\it ini\_parms.F}) |
559 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
560 |
}. |
}. |
561 |
It is copied into model general vertical coordinate variable |
It is copied into model general vertical coordinate variable |
562 |
{\bf |
{\bf |
563 |
\begin{rawhtml} <A href=../../../code_reference/vdb/names/PD.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/names/PD.htm> \end{rawhtml} |
564 |
diffKrS |
diffKrS |
565 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
566 |
} which is used in routine |
} which is used in routine |
567 |
{\it |
{\it |
568 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
569 |
S/R CALC\_DIFFUSIVITY ({\it calc\_diffusivity.F}) |
S/R CALC\_DIFFUSIVITY ({\it calc\_diffusivity.F}) |
570 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
571 |
}. |
}. |
594 |
to $2 \times 10^{-4}$ $^o$ C$^{-1}$. |
to $2 \times 10^{-4}$ $^o$ C$^{-1}$. |
595 |
The variable |
The variable |
596 |
{\bf |
{\bf |
597 |
\begin{rawhtml} <A href=../../../code_reference/vdb/names/ZV.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/names/ZV.htm> \end{rawhtml} |
598 |
tAlpha |
tAlpha |
599 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
600 |
} |
} |
601 |
is read in the routine |
is read in the routine |
602 |
{\it |
{\it |
603 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
604 |
S/R INI\_PARMS ({\it ini\_parms.F}) |
S/R INI\_PARMS ({\it ini\_parms.F}) |
605 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
606 |
}. |
}. |
607 |
The routine |
The routine |
608 |
{\it |
{\it |
609 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
610 |
S/R FIND\_RHO ({\it find\_rho.F}) |
S/R FIND\_RHO ({\it find\_rho.F}) |
611 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
612 |
} makes use of {\bf tAlpha}. |
} makes use of {\bf tAlpha}. |
672 |
gradient solver will use to 40, {\bf irrespective of the convergence |
gradient solver will use to 40, {\bf irrespective of the convergence |
673 |
criteria being met}. Used in routine |
criteria being met}. Used in routine |
674 |
{\it |
{\it |
675 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
676 |
S/R CG3D ({\it cg3d.F}) |
S/R CG3D ({\it cg3d.F}) |
677 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
678 |
}. |
}. |
690 |
tolerance falls below this value or until the maximum number of |
tolerance falls below this value or until the maximum number of |
691 |
solver iterations is reached. Used in routine |
solver iterations is reached. Used in routine |
692 |
{\it |
{\it |
693 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
694 |
S/R CG3D ({\it cg3d.F}) |
S/R CG3D ({\it cg3d.F}) |
695 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
696 |
}. |
}. |
759 |
surface heat flux file used in the example. |
surface heat flux file used in the example. |
760 |
The variable |
The variable |
761 |
{\bf |
{\bf |
762 |
\begin{rawhtml} <A href=../../../code_reference/vdb/names/179.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/names/179.htm> \end{rawhtml} |
763 |
Qsurf |
Qsurf |
764 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
765 |
} |
} |
766 |
is read in the routine |
is read in the routine |
767 |
{\it |
{\it |
768 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
769 |
S/R INI\_PARMS ({\it ini\_parms.F}) |
S/R INI\_PARMS ({\it ini\_parms.F}) |
770 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
771 |
} |
} |
772 |
and applied in |
and applied in |
773 |
{\it |
{\it |
774 |
\begin{rawhtml} <A href=../../../code_reference/vdb/code/94.htm> \end{rawhtml} |
\begin{rawhtml} <A href=../code_reference/vdb/code/94.htm> \end{rawhtml} |
775 |
S/R EXTERNAL\_FORCING\_SURF ({\it external\_forcing\_surf.F}) |
S/R EXTERNAL\_FORCING\_SURF ({\it external\_forcing\_surf.F}) |
776 |
\begin{rawhtml} </A>\end{rawhtml} |
\begin{rawhtml} </A>\end{rawhtml} |
777 |
} where the flux is converted to a temperature tendency. |
} where the flux is converted to a temperature tendency. |
788 |
|
|
789 |
|
|
790 |
\subsubsection{File {\it input/data.pkg}} |
\subsubsection{File {\it input/data.pkg}} |
791 |
|
\label{www:tutorials} |
792 |
|
|
793 |
This file uses standard default values and does not contain |
This file uses standard default values and does not contain |
794 |
customisations for this experiment. |
customisations for this experiment. |
795 |
|
|
796 |
\subsubsection{File {\it input/eedata}} |
\subsubsection{File {\it input/eedata}} |
797 |
|
\label{www:tutorials} |
798 |
|
|
799 |
This file uses standard default values and does not contain |
This file uses standard default values and does not contain |
800 |
customisations for this experiment. |
customisations for this experiment. |
801 |
|
|
802 |
|
|
803 |
\subsubsection{File {\it input/Qsurf.bin}} |
\subsubsection{File {\it input/Qsurf.bin}} |
804 |
|
\label{www:tutorials} |
805 |
|
|
806 |
The file {\it input/Qsurf.bin} specifies a two-dimensional ($x,y$) |
The file {\it input/Qsurf.bin} specifies a two-dimensional ($x,y$) |
807 |
map of heat flux values where |
map of heat flux values where |
824 |
\end{figure} |
\end{figure} |
825 |
|
|
826 |
\subsection{Running the example} |
\subsection{Running the example} |
827 |
|
\label{www:tutorials} |
828 |
|
|
829 |
\subsubsection{Code download} |
\subsubsection{Code download} |
830 |
|
\label{www:tutorials} |
831 |
|
|
832 |
In order to run the examples you must first download the code distribution. |
In order to run the examples you must first download the code distribution. |
833 |
Instructions for downloading the code can be found in \ref{sect:obtainingCode}. |
Instructions for downloading the code can be found in \ref{sect:obtainingCode}. |
834 |
|
|
835 |
\subsubsection{Experiment Location} |
\subsubsection{Experiment Location} |
836 |
|
\label{www:tutorials} |
837 |
|
|
838 |
This example experiments is located under the release sub-directory |
This example experiments is located under the release sub-directory |
839 |
|
|
841 |
{\it verification/convection/ } |
{\it verification/convection/ } |
842 |
|
|
843 |
\subsubsection{Running the Experiment} |
\subsubsection{Running the Experiment} |
844 |
|
\label{www:tutorials} |
845 |
|
|
846 |
To run the experiment |
To run the experiment |
847 |
|
|