--- manual/s_examples/deep_convection/convection.tex 2002/05/16 15:54:37 1.3
+++ manual/s_examples/deep_convection/convection.tex 2010/08/27 13:25:31 1.9
@@ -1,6 +1,12 @@
\section{Surface Driven Convection}
\label{www:tutorials}
\label{sect:eg-bconv}
+\begin{rawhtml}
+
+\end{rawhtml}
+\begin{center}
+(in directory: {\it verification/tutorial\_deep\_convection/})
+\end{center}
\bodytext{bgcolor="#FFFFFFFF"}
@@ -17,7 +23,7 @@
\begin{center}
\resizebox{7.5cm}{5.5cm}{
\includegraphics*[0.2in,0.7in][10.5in,10.5in]
- {part3/case_studies/doubly_periodic_convection/simulation_config.eps} }
+ {s_examples/deep_convection/simulation_config.eps} }
\end{center}
\caption{Schematic of simulation domain
for the surface driven convection experiment. The domain is doubly periodic
@@ -26,9 +32,11 @@
\label{FIG:eg-bconv-simulation_config}
\end{figure}
-This experiment, figure \ref{FIG:eg-bconv-simulation_config}, showcasing MITgcm's non-hydrostatic capability, was designed to explore
+This experiment, figure \ref{FIG:eg-bconv-simulation_config}, showcasing MITgcm's non-hydrostatic
+capability, was designed to explore
the temporal and spatial characteristics of convection plumes as they might exist during a
-period of oceanic deep convection. It is
+period of oceanic deep convection. The files for this experiment can be found in the verification
+directory under tutorial\_deep\_convection. It is
\begin{itemize}
\item non-hydrostatic
@@ -82,7 +90,7 @@
\begin{center}
\resizebox{15cm}{10cm}{
\includegraphics*[0.2in,0.7in][10.5in,10.5in]
- {part3/case_studies/doubly_periodic_convection/verticalsection.ps} }
+ {s_examples/deep_convection/verticalsection.ps} }
\end{center}
\caption{
}
@@ -96,7 +104,7 @@
\begin{center}
\resizebox{10cm}{10cm}{
\includegraphics*[0.2in,0.7in][10.5in,10.5in]
- {part3/case_studies/doubly_periodic_convection/surfacesection.ps} }
+ {s_examples/deep_convection/surfacesection.ps} }
\end{center}
\caption{
}
@@ -265,7 +273,7 @@
\begin{rawhtml}\end{rawhtml}
\begin{small}
-\input{part3/case_studies/doubly_periodic_convection/code/SIZE.h}
+\input{s_examples/deep_convection/code/SIZE.h}
\end{small}
\begin{rawhtml}\end{rawhtml}
@@ -280,29 +288,29 @@
\item Line 4,
\begin{verbatim}
- 4 tRef=20*20.0,
+ 4 tRef=20*20.0,
\end{verbatim}
this line sets
the initial and reference values of potential temperature at each model
-level in units of $^{\circ}$C. Here the value is arbitrary since, in this case, the
+level in units of $^{\circ}\mathrm{C}$. Here the value is arbitrary since, in this case, the
flow evolves independently of the absolute magnitude of the reference temperature.
For each depth level the initial and reference profiles will be uniform in
$x$ and $y$. The values specified are read into the
variable
{\bf
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
tRef
\begin{rawhtml} \end{rawhtml}
}
in the model code, by procedure
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R INI\_PARMS ({\it ini\_parms.F})
\begin{rawhtml} \end{rawhtml}.
}
The temperature field is initialised, by procedure
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R INI\_THETA ({\it ini\_theta.F})
\begin{rawhtml} \end{rawhtml}.
}
@@ -310,7 +318,7 @@
\item Line 5,
\begin{verbatim}
- 5 sRef=20*35.0,
+ 5 sRef=20*35.0,
\end{verbatim}
this line sets the initial and reference values of salinity at each model
level in units of ppt. In this case salinity is set to an (arbitrary) uniform value of
@@ -320,19 +328,19 @@
$x$ and $y$. The values specified are read into the
variable
{\bf
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
sRef
\begin{rawhtml} \end{rawhtml}
}
in the model code, by procedure
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R INI\_PARMS ({\it ini\_parms.F})
}
\begin{rawhtml} \end{rawhtml}.
The salinity field is initialised, by procedure
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R INI\_SALT ({\it ini\_salt.F})
\begin{rawhtml} \end{rawhtml}.
}
@@ -340,30 +348,30 @@
\item Line 6,
\begin{verbatim}
- 6 viscAh=0.1,
+ 6 viscAh=0.1,
\end{verbatim}
this line sets the horizontal laplacian dissipation coefficient to
0.1 ${\rm m^{2}s^{-1}}$. Boundary conditions
for this operator are specified later.
The variable
{\bf
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
viscAh
\begin{rawhtml} \end{rawhtml}
}
is read in the routine
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R INI\_PARMS ({\it ini\_params.F})
\begin{rawhtml} \end{rawhtml}
} and applied in routines
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R CALC\_MOM\_RHS ({\it calc\_mom\_rhs.F})
\begin{rawhtml} \end{rawhtml}
} and
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R CALC\_GW ({\it calc\_gw.F})
\begin{rawhtml} \end{rawhtml}
}.
@@ -371,32 +379,32 @@
\item Line 7,
\begin{verbatim}
- 7 viscAz=0.1,
+ 7 viscAz=0.1,
\end{verbatim}
this line sets the vertical laplacian frictional dissipation coefficient to
0.1 ${\rm m^{2}s^{-1}}$. Boundary conditions
for this operator are specified later.
The variable
{\bf
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
viscAz
\begin{rawhtml} \end{rawhtml}
}
is read in the routine
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R INI\_PARMS ({\it ini\_parms.F})
\begin{rawhtml} \end{rawhtml}
}
and is copied into model general vertical coordinate variable
{\bf
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
viscAr
\begin{rawhtml} \end{rawhtml}
}. At each time step, the viscous term contribution to the momentum equations
is calculated in routine
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R CALC\_DIFFUSIVITY ({\it calc\_diffusivity.F})
\begin{rawhtml} \end{rawhtml}
}.
@@ -412,18 +420,18 @@
$\frac{\partial v}{\partial x}$=0 along boundaries in $x$.
The variable
{\bf
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
no\_slip\_sides
\begin{rawhtml} \end{rawhtml}
}
is read in the routine
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R INI\_PARMS ({\it ini\_parms.F})
\begin{rawhtml} \end{rawhtml}
} and the boundary condition is evaluated in routine
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R CALC\_MOM\_RHS ({\it calc\_mom\_rhs.F})
\begin{rawhtml} \end{rawhtml}
}.
@@ -438,18 +446,18 @@
e.g. $u=v=0$ at $z=-H$, where $H$ is the local depth of the domain.
The variable
{\bf
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
no\_slip\_bottom
\begin{rawhtml} \end{rawhtml}
}
is read in the routine
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R INI\_PARMS ({\it ini\_parms.F})
\begin{rawhtml} \end{rawhtml}
} and is applied in the routine
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R CALC\_MOM\_RHS ({\it calc\_mom\_rhs.F})
\begin{rawhtml} \end{rawhtml}
}.
@@ -464,18 +472,18 @@
all boundaries.
The variable
{\bf
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
diffKhT
\begin{rawhtml} \end{rawhtml}
}
is read in the routine
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R INI\_PARMS ({\it ini\_parms.F})
\begin{rawhtml} \end{rawhtml}
} and used in routine
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R CALC\_GT ({\it calc\_gt.F})
\begin{rawhtml} \end{rawhtml}
}.
@@ -489,24 +497,24 @@
operator is $\frac{\partial}{\partial z}$ = 0 on all boundaries.
The variable
{\bf
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
diffKzT
\begin{rawhtml} \end{rawhtml}
}
is read in the routine
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R INI\_PARMS ({\it ini\_parms.F})
\begin{rawhtml} \end{rawhtml}
}.
It is copied into model general vertical coordinate variable
{\bf
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
diffKrT
\begin{rawhtml} \end{rawhtml}
} which is used in routine
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R CALC\_DIFFUSIVITY ({\it calc\_diffusivity.F})
\begin{rawhtml} \end{rawhtml}
}.
@@ -522,18 +530,18 @@
all boundaries.
The variable
{\bf
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
diffKsT
\begin{rawhtml} \end{rawhtml}
}
is read in the routine
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R INI\_PARMS ({\it ini\_parms.F})
\begin{rawhtml} \end{rawhtml}
} and used in routine
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R CALC\_GS ({\it calc\_gs.F})
\begin{rawhtml} \end{rawhtml}
}.
@@ -548,24 +556,24 @@
operator is $\frac{\partial}{\partial z}$ = 0 on all boundaries.
The variable
{\bf
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
diffKzS
\begin{rawhtml} \end{rawhtml}
}
is read in the routine
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R INI\_PARMS ({\it ini\_parms.F})
\begin{rawhtml} \end{rawhtml}
}.
It is copied into model general vertical coordinate variable
{\bf
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
diffKrS
\begin{rawhtml} \end{rawhtml}
} which is used in routine
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R CALC\_DIFFUSIVITY ({\it calc\_diffusivity.F})
\begin{rawhtml} \end{rawhtml}
}.
@@ -594,19 +602,19 @@
to $2 \times 10^{-4}$ $^o$ C$^{-1}$.
The variable
{\bf
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
tAlpha
\begin{rawhtml} \end{rawhtml}
}
is read in the routine
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R INI\_PARMS ({\it ini\_parms.F})
\begin{rawhtml} \end{rawhtml}
}.
The routine
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R FIND\_RHO ({\it find\_rho.F})
\begin{rawhtml} \end{rawhtml}
} makes use of {\bf tAlpha}.
@@ -672,7 +680,7 @@
gradient solver will use to 40, {\bf irrespective of the convergence
criteria being met}. Used in routine
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R CG3D ({\it cg3d.F})
\begin{rawhtml} \end{rawhtml}
}.
@@ -690,7 +698,7 @@
tolerance falls below this value or until the maximum number of
solver iterations is reached. Used in routine
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R CG3D ({\it cg3d.F})
\begin{rawhtml} \end{rawhtml}
}.
@@ -759,19 +767,19 @@
surface heat flux file used in the example.
The variable
{\bf
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
Qsurf
\begin{rawhtml} \end{rawhtml}
}
is read in the routine
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R INI\_PARMS ({\it ini\_parms.F})
\begin{rawhtml} \end{rawhtml}
}
and applied in
{\it
-\begin{rawhtml} \end{rawhtml}
+\begin{rawhtml} \end{rawhtml}
S/R EXTERNAL\_FORCING\_SURF ({\it external\_forcing\_surf.F})
\begin{rawhtml} \end{rawhtml}
} where the flux is converted to a temperature tendency.
@@ -782,7 +790,7 @@
\begin{rawhtml}\end{rawhtml}
\begin{small}
-\input{part3/case_studies/doubly_periodic_convection/input/data}
+\input{s_examples/deep_convection/input/data}
\end{small}
\begin{rawhtml}\end{rawhtml}
@@ -816,7 +824,7 @@
\begin{center}
% \resizebox{15cm}{10cm}{
% \includegraphics*[0.2in,0.7in][10.5in,10.5in]
-% {part3/case_studies/doubly_periodic_convection/Qsurf.ps} }
+% {s_examples/deep_convection/Qsurf.ps} }
\end{center}
\caption{
}