--- manual/s_examples/deep_convection/convection.tex 2002/05/16 15:54:37 1.3 +++ manual/s_examples/deep_convection/convection.tex 2003/07/30 13:42:52 1.4 @@ -290,19 +290,19 @@ $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}. } @@ -320,19 +320,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}. } @@ -347,23 +347,23 @@ 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} }. @@ -378,25 +378,25 @@ 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 +412,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 +438,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 +464,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 +489,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 +522,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 +548,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 +594,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 +672,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 +690,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 +759,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.