/[MITgcm]/manual/s_phys_pkgs/text/shelfice.tex

Diff of /manual/s_phys_pkgs/text/shelfice.tex

Parent Directory | Revision Log | View Revision Graph Revision Graph | View Patch Patch

-revision 1.2 by mlosch,
Thu May 14 15:35:48 2009 UTC
+revision 1.5 by mlosch,
Thu Dec  1 11:59:46 2011 UTC
 Line 11 
 Authors: Martin Losch, Jean-Michel Campi
  %----------------------------------------------------------------------
  \subsubsection{Introduction
- \label{sec:pkg:exf:intro}}
+ \label{sec:pkg:shelfice:intro}}
  Package ``shelfice'' provides a thermodynamic model for basal melting
 Line 26 
 A description of key subroutines is give
  \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}.
  %----------------------------------------------------------------------
 Line 56 
 schemes, be aware, that non-local parame
  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
  \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
 Line 94 
 and \code{data.shelfice} (read in \code{
  ~ \\
  %
  A package is switched on/off at run-time by setting
- (e.g. for SHELFICE) \code{useSHELFICE = .TRUE.} in \code{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
 Line 129 
 Table~\ref{tab:pkg:shelfice:runtimeparms
        &  %---ref---
        \\
        SHELFICEtopoFile         & UNSET
-       &   under-ice topography of ice shelfes
+       &   under-ice topography of ice shelves
        &  %---ref---
        \\
        SHELFICElatentHeat       &  334.0E+03
 Line 137 
 Table~\ref{tab:pkg:shelfice:runtimeparms
        &  %---ref---
        \\
        SHELFICEHeatCapacity\_Cp & 2000.0E+00
-       &   latent heat of fusion ($c_{p,I}$)
+       &   specific heat capacity of ice ($c_{p,I}$)
        &  %---ref---
        \\
        rhoShelfIce              &  917.0E+00
 Line 199 
 Table~\ref{tab:pkg:shelfice:runtimeparms
  }
  \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
-Line 390 
 $T_{f}$ which is computed from a linear
+Line 407 
 $T_{f}$ which is computed from a linear
  \begin{equation}
    \label{eq:helmerfreeze}
      T_{f} = (0.0901 - 0.0575\ S_{b})^{\circ}
-     - 7.61 \times 10^{-4}\frac{^{\circ}}{\text{dBar}}\ p_{b}
+     - 7.61 \times 10^{-4}\frac{\text{K}}{\text{dBar}}\ p_{b}
  \end{equation}
  with the salinity $S_{b}$ and the pressure $p_{b}$ (in dBar) in the
  cell at the ice-water interface. From the salt budget, the salt flux
-Line 505 
 C ...
+Line 522 
 C ...
  C | | |-EXTERNAL_FORCING_T :: Problem specific forcing for temperature.
  C | | |-SHELFICE_FORCING_T :: apply heat fluxes from ice shelf model
  C ...
- C | | |-EXTERNAL_FORCING_S :: Problem specific forcing for temperature.
+ C | | |-EXTERNAL_FORCING_S :: Problem specific forcing for salinity.
  C | | |-SHELFICE_FORCING_S :: apply fresh water fluxes from ice shelf model
  C ...
  C | |-DYNAMICS             :: Momentum equations driver.
-Line 535 
 Diagnostics output is available via the
+Line 552 
 Diagnostics output is available via the
  Available output fields are summarized in
  Table \ref{tab:pkg:shelfice:diagnostics}.
- \begin{table}[h!]
+ \begin{table}[!ht]
  \centering
  \label{tab:pkg:shelfice:diagnostics}
  {\footnotesize
-Line 569 
 Table \ref{tab:pkg:shelfice:diagnostics}
+Line 586 
 Table \ref{tab:pkg:shelfice:diagnostics}
  %%% Local Variables:
  %%% mode: latex
- %%% TeX-master: "../manual"
+ %%% TeX-master: "../../manual"
  %%% End:

 Legend:



Removed from v.1.2
 


changed lines


 
Added in v.1.5
 Legend:



Removed from v.1.2
 


changed lines


 
Added in v.1.5
-Removed from v.1.2
+Added in v.1.5

	ViewVC Help
Powered by ViewVC 1.1.22