--- manual/s_phys_pkgs/text/shelfice.tex 2009/05/14 10:19:12 1.1 +++ manual/s_phys_pkgs/text/shelfice.tex 2010/08/30 23:09:22 1.4 @@ -1,4 +1,4 @@ -% $Header: /home/ubuntu/mnt/e9_copy/manual/s_phys_pkgs/text/shelfice.tex,v 1.1 2009/05/14 10:19:12 mlosch Exp $ +% $Header: /home/ubuntu/mnt/e9_copy/manual/s_phys_pkgs/text/shelfice.tex,v 1.4 2010/08/30 23:09:22 jmc Exp $ % $Name: $ \subsection{SHELFICE Package} @@ -11,7 +11,7 @@ %---------------------------------------------------------------------- \subsubsection{Introduction -\label{sec:pkg:exf:intro}} +\label{sec:pkg:shelfice:intro}} Package ``shelfice'' provides a thermodynamic model for basal melting @@ -20,13 +20,13 @@ CPP options enable or disable different aspects of the package (Section \ref{sec:pkg:shelfice:config}). Run-Time options, flags, filenames and field-related dates/times are -set in \texttt{data.shelfice} +set in \code{data.shelfice} (Section \ref{sec:pkg:shelfice:runtime}). A description of key subroutines is given in Section \ref{sec:pkg:shelfice:subroutines}. Input fields, units and sign conventions are summarized in Section \ref{sec:pkg:shelfice:fields_units}, and available diagnostics -output is listed in Section \ref{sec:pkg:shelfice:fields_diagnostics}. +output is listed in Section \ref{sec:pkg:shelfice:diagnostics}. %---------------------------------------------------------------------- @@ -41,11 +41,11 @@ \begin{itemize} % \item -using the \texttt{packages.conf} file by adding \texttt{shelfice} to it, +using the \code{packages.conf} file by adding \code{shelfice} to it, % \item -or using \texttt{genmake2} adding -\texttt{-enable=shelfice} or \texttt{-disable=shelfice} switches +or using \code{genmake2} adding +\code{-enable=shelfice} or \code{-disable=shelfice} switches % \item \textit{required packages and CPP options}: \\ @@ -56,14 +56,14 @@ off, e.g.\ for KPP (\ref{sec:pkg:kpp}). % \end{itemize} -(see Section \ref{sect:buildingCode}). +(see Section \ref{sec:buildingCode}). Parts of the SHELFICE code can be enabled or disabled at compile time via CPP preprocessor flags. These options are set -\texttt{SHELFICE\_OPTIONS.h}. +\code{SHELFICE\_OPTIONS.h}. Table \ref{tab:pkg:shelfice:cpp} summarizes these options. -\begin{table}[h!] +\begin{table}[!ht] \centering \label{tab:pkg:shelfice:cpp} {\footnotesize @@ -71,14 +71,14 @@ \hline \textbf{CPP option} & \textbf{Description} \\ \hline \hline - \texttt{ALLOW\_SHELFICE\_DEBUG} & + \code{ALLOW\_SHELFICE\_DEBUG} & Include code for enhanged diagnosis \\ - \texttt{ALLOW\_ISOMIP\_TD} & + \code{ALLOW\_ISOMIP\_TD} & Include code for simplifed ISOMIP thermodynamics \\ \hline \end{tabular} } - \caption{Available CPP-flags to be set in \texttt{SHELFICE\_OPTIONS.h}} + \caption{Available CPP-flags to be set in \code{SHELFICE\_OPTIONS.h}} \end{table} %---------------------------------------------------------------------- @@ -87,20 +87,20 @@ \label{sec:pkg:shelfice:runtime}} Run-time parameters are set in files -\texttt{data.pkg} (read in \texttt{packages\_readparms.F}), -and \texttt{data.shelfice} (read in \texttt{shelfice\_readparms.F}). +\code{data.pkg} (read in \code{packages\_readparms.F}), +and \code{data.shelfice} (read in \code{shelfice\_readparms.F}). \paragraph{Enabling the package} ~ \\ % A package is switched on/off at run-time by setting -(e.g. for SHELFICE) \texttt{useSHELFICE = .TRUE.} in \texttt{data.pkg}. +(e.g. for SHELFICE) \code{useSHELFICE = .TRUE.} in \code{data.pkg}. \paragraph{General flags and parameters} ~ \\ % Table~\ref{tab:pkg:shelfice:runtimeparms} lists all run-time parameters. -\begin{table}[h!] +\begin{table}[!ht] \caption{Run-time parameters and default values \label{tab:pkg:shelfice:runtimeparms}} {\footnotesize @@ -129,7 +129,7 @@ & %---ref--- \\ SHELFICEtopoFile & UNSET - & under-ice topography of ice shelfes + & under-ice topography of ice shelves & %---ref--- \\ SHELFICElatentHeat & 334.0E+03 @@ -199,7 +199,24 @@ } \end{table} +\paragraph{Input fields and units\label{sec:pkg:shelfice:fields_units}} +\begin{description} +\item[\code{SHEFLICEtopoFile}:] under-ice topography of ice shelves in + meters; upwards is positive, that as for the bathymetry files, + negative values are required for topography below the sea-level; +\item[\code{SHEFLICEloadAnomalyFile}:] pressure load anomaly at the bottom of + the ice shelves in pressure units (Pa); this field is absolutely + required to avoid large excursions of the free surface during + initial adjustment processes; obtained by integrating an approximate + density from the surface at $z=0$ down to the bottom of the last + fully dry cell within the ice shelf, see + Eq.~(\ref{eq:surfacepressure}); however, the file + \code{SHEFLICEloadAnomalyFile} must not be $p_{top}$, but + $p_{top}-g\sum_{k'=1}^{n-1}\rho_{0}\Delta{z}_{k'}$, with + $\rho_{0}=$~\code{rhoConst}, so that in the absenses of a $\rho^{*}$ + that is different from $\rho_{0}$, the anomaly is zero. +\end{description} %---------------------------------------------------------------------- \subsubsection{Description @@ -228,7 +245,7 @@ of the water column is the atmospheric pressure $p_{a}$ plus the weight of the ice-shelf. It is this weight of the ice-shelf that has to be provided as a boundary condition at the top of the water column -(in run-time parameter \texttt{SHELFICEloadAnomalyFile}). +(in run-time parameter \code{SHELFICEloadAnomalyFile}). The weight is conveniently computed by integrating a density profile $\rho^*$, that is constant in time and corresponds to the sea-water replaced by ice, from $z=0$ to a ``reference'' ice-shelf draft at @@ -277,7 +294,7 @@ \end{equation} where $H(k'-k)=1$ for $k'