New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 14403 for NEMO/branches/2021/dev_r14312_MPI_Interface/doc/latex/NEMO/subfiles – NEMO

Ignore:
Timestamp:
2021-02-05T13:13:14+01:00 (3 years ago)
Author:
smasson
Message:

dev_r14312_MPI_Interface: merge with trunk@14402, #2598

Location:
NEMO/branches/2021/dev_r14312_MPI_Interface/doc/latex/NEMO/subfiles
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2021/dev_r14312_MPI_Interface/doc/latex/NEMO/subfiles

    • Property svn:ignore
      •  

        old new  
         1*.aux 
         2*.bbl 
         3*.blg 
         4*.fdb* 
         5*.fls 
         6*.idx 
         7*.ilg 
        18*.ind 
        2 *.ilg 
         9*.lo* 
         10*.out 
         11*.pdf 
         12*.pyg 
         13*.tdo 
         14*.toc 
         15*.xdv 
         16cache* 
  • NEMO/branches/2021/dev_r14312_MPI_Interface/doc/latex/NEMO/subfiles/apdx_DOMAINcfg.tex

    r14303 r14403  
    66\label{apdx:DOMCFG} 
    77 
    8 %    {\em 4.0} & {\em Andrew Coward} & {\em Created at v4.0 from materials removed from chap\_DOM that are still relevant to the \forcode{DOMAINcfg} tool and which illustrate and explain the choices to be made by the user when setting up new domains }  \\ 
    9  
    108\chaptertoc 
    119 
     
    1412{\footnotesize 
    1513  \begin{tabularx}{\textwidth}{l||X|X} 
    16     Release & Author(s) & Modifications \\ 
    17     \hline 
    18     {\em   next}& {\em Pierre Mathiot} & {\em add ice shelf and closed sea option description } \\ 
    19     {\em   4.0} & {\em Andrew Coward}  & {\em Created at v4.0 from materials removed from chap\_DOM that are still relevant to the \forcode{DOMAINcfg} tool and which illustrate and explain the choices to be made by the user when setting up new domains }  \\ 
    20     {\em   3.6} & {\em ...} & {\em ...} \\ 
    21     {\em   3.4} & {\em ...} & {\em ...} \\ 
    22     {\em <=3.4} & {\em ...} & {\em ...} 
     14    Release     & Author(s)            & Modifications                                                 \\ 
     15    \hline 
     16    {\em  next} & {\em Pierre Mathiot} & {\em Add ice shelf and closed sea option description        } \\ 
     17    {\em   4.0} & {\em  Andrew Coward} & {\em Creation from materials removed from \autoref{chap:DOM} 
     18                                              that are still relevant to the DOMAINcfg tool 
     19                                              when setting up new domains                            } 
    2320  \end{tabularx} 
    2421} 
     
    592589  the \textit{isfdraft\_meter} file (Netcdf format). This file need to include the \textit{isf\_draft} variable.  
    593590  A positive value will mean ice shelf/ocean or ice shelf bedrock interface below the reference 0m ssh.  
    594   The exact shape of the ice shelf cavity (grounding line position and minimum thickness of the water column under an ice shelf, ...) can be specify in \nam{zgr_isf}{zgr_isf}. 
     591  The exact shape of the ice shelf cavity (grounding line position and minimum thickness of the water column under an ice shelf, ...) can be specify in \nam{zgr_isf}{zgr\_isf}. 
    595592 
    596593\begin{listing} 
     
    616613\end{listing} 
    617614 
    618    The options available to define the shape of the under ice shelf cavities are listed in \nam{zgr_isf}{zgr_isf} (\texttt{DOMAINcfg} only, \autoref{lst:namzgr_isf}). 
    619  
    620    \subsection{Model ice shelf draft definition} 
    621    \label{subsec:zgrisf_isfd} 
    622  
    623    First of all, the tool make sure, the ice shelf draft ($h_{isf}$) is sensible and compatible with the bathymetry. 
    624    There are 3 compulsory steps to achieve this: 
    625  
    626    \begin{description} 
    627    \item{\np{rn_isfdep_min}{rn\_isfdep\_min}:} this is the minimum ice shelf draft. This is to make sure there is no ridiculous thin ice shelf. If \np{rn_isfdep_min}{rn\_isfdep\_min} is smaller than the surface level, \np{rn_isfdep_min}{rn\_isfdep\_min} is set to $e3t\_1d(1)$.  
    628    Where $h_{isf} < MAX(e3t\_1d(1),\np{rn_isfdep_min}{rn\_isfdep\_min}$), $h_{isf}$ is set to \np{rn_isfdep_min}{rn\_isfdep\_min}. 
    629  
    630    \item{\np{rn_glhw_min}{rn\_glhw\_min}:} This parameter is used to define the grounding line position. 
    631    Where the difference between the bathymetry and the ice shelf draft is smaller than \np{rn_glhw_min}{rn\_glhw\_min}, the cell are grounded (ie masked).  
    632    This step is needed to take into account possible small mismatch between ice shelf draft value and bathymetry value (sources are coming from different grid, different data processes, rounding error, ...). 
    633  
    634    \item{\np{rn_isfhw_min}{rn\_isfhw\_min}:} This parameter is the minimum water column thickness in the cavity.  
    635    Where the water column thickness is lower than \np{rn_isfhw_min}{rn\_isfhw\_min}, the ice shelf draft is adjusted to match this criterion.  
    636    If for any reason, this adjustement break the minimum ice shelf draft allowed (\np{rn_isfdep_min}{rn\_isfdep\_min}), the cell is masked. 
    637    \end{description} 
    638  
    639    Once all these adjustements are made, if the water column thickness contains one cell wide channels, these channels can be closed using \np{ln_isfchannel}{ln\_isfchannel}.   
     615   The options available to define the shape of the under ice shelf cavities are listed in \nam{zgr_isf}{zgr\_isf} (\texttt{DOMAINcfg} only, \autoref{lst:namzgr_isf}). 
     616 
     617\subsection{Model ice shelf draft definition} 
     618\label{subsec:zgrisf_isfd} 
     619 
     620First of all, the tool make sure, the ice shelf draft ($h_{isf}$) is sensible and compatible with the bathymetry. 
     621There are 3 compulsory steps to achieve this: 
     622 
     623\begin{description} 
     624\item{\np{rn_isfdep_min}{rn\_isfdep\_min}:} this is the minimum ice shelf draft. This is to make sure there is no ridiculous thin ice shelf. If \np{rn_isfdep_min}{rn\_isfdep\_min} is smaller than the surface level, \np{rn_isfdep_min}{rn\_isfdep\_min} is set to $e3t\_1d(1)$.  
     625  Where $h_{isf} < MAX(e3t\_1d(1),rn\_isfdep\_min)$, $h_{isf}$ is set to \np{rn_isfdep_min}{rn\_isfdep\_min}. 
     626 
     627\item{\np{rn_glhw_min}{rn\_glhw\_min}:} This parameter is used to define the grounding line position. 
     628  Where the difference between the bathymetry and the ice shelf draft is smaller than \np{rn_glhw_min}{rn\_glhw\_min}, the cell are grounded (ie masked).  
     629  This step is needed to take into account possible small mismatch between ice shelf draft value and bathymetry value (sources are coming from different grid, different data processes, rounding error, ...). 
     630 
     631\item{\np{rn_isfhw_min}{rn\_isfhw\_min}:} This parameter is the minimum water column thickness in the cavity.  
     632  Where the water column thickness is lower than \np{rn_isfhw_min}{rn\_isfhw\_min}, the ice shelf draft is adjusted to match this criterion.  
     633  If for any reason, this adjustement break the minimum ice shelf draft allowed (\np{rn_isfdep_min}{rn\_isfdep\_min}), the cell is masked. 
     634\end{description} 
     635 
     636Once all these adjustements are made, if the water column thickness contains one cell wide channels, these channels can be closed using \np{ln_isfchannel}{ln\_isfchannel}.   
    640637  
    641    \subsection{Model top level definition} 
    642    After the definition of the ice shelf draft, the tool defines the top level.  
    643    The compulsory criterion is that the water column needs at least 2 wet cells in the water column at U- and V-points. 
    644    To do so, if there one cell wide water column, the tools adjust the ice shelf draft to fillful the requierement.\\ 
    645  
    646    The process is the following: 
    647    \begin{description} 
    648    \item{step 1:} The top level is defined in the same way as the bottom level is defined. 
    649    \item{step 2:} The isolated grid point in the bathymetry are filled (as it is done in a domain without ice shelf) 
    650    \item{step 3:} The tools make sure, the top level is above or equal to the bottom level 
    651    \item{step 4:} If the water column at a U- or V- point is one wet cell wide, the ice shelf draft is adjusted. So the actual top cell become fully open and the new 
    652    top cell thickness is set to the minimum cell thickness allowed (following the same logic as for the bottom partial cell). This step is iterated 4 times to ensure the condition is fullfill along the 4 sides of the cell. 
    653    \end{description} 
    654  
    655    In case of steep slope and shallow water column, it likely that 2 cells are disconnected (bathymetry above its neigbourging ice shelf draft).  
    656    The option \np{ln_isfconnect}{ln\_isfconnect} allow the tool to force the connection between these 2 cells. 
    657    Some limiters in meter or levels on the digging allowed by the tool are available (respectively, \np{rn_zisfmax}{rn\_zisfmax} or \np{rn_kisfmax}{rn\_kisfmax}). 
    658    This will prevent the formation of subglacial lakes at the expense of long vertical pipe to connect cells at very different levels. 
    659  
    660    \subsection{Subglacial lakes} 
    661    Despite careful setting of your ice shelf draft and bathymetry input file as well as setting described in \autoref{subsec:zgrisf_isfd}, some situation are unavoidable. 
    662    For exemple if you setup your ice shelf draft and bathymetry to do ocean/ice sheet coupling,  
    663    you may decide to fill the whole antarctic with a bathymetry and an ice shelf draft value (ice/bedrock interface depth when grounded).  
    664    If you do so, the subglacial lakes will show up (Vostock for example). An other possibility is with coarse vertical resolution, some ice shelves could be cut in 2 parts:  
    665    one connected to the main ocean and an other one closed which can be considered as a subglacial sea be the model.\\ 
    666  
    667    The namelist option \np{ln_isfsubgl}{ln\_isfsubgl} allow you to remove theses subglacial lakes. 
    668    This may be useful for esthetical reason or for stability reasons: 
    669  
    670    \begin{description} 
    671    \item $\bullet$ In a subglacial lakes, in case of very weak circulation (often the case), the only heat flux is the conductive heat flux through the ice sheet.  
    672          This will lead to constant freezing until water reaches -20C.  
    673          This is one of the defitiency of the 3 equation melt formulation (for details on this formulation, see: \autoref{sec:isf}). 
    674    \item $\bullet$ In case of coupling with an ice sheet model,  
    675          the ssh in the subglacial lakes and the main ocean could be very different (ssh initial adjustement for example),  
    676          and so if for any reason both a connected at some point, the model is likely to fall over.\\ 
    677    \end{description} 
     638\subsection{Model top level definition} 
     639After the definition of the ice shelf draft, the tool defines the top level.  
     640The compulsory criterion is that the water column needs at least 2 wet cells in the water column at U- and V-points. 
     641To do so, if there one cell wide water column, the tools adjust the ice shelf draft to fillful the requierement.\\ 
     642 
     643The process is the following: 
     644\begin{description} 
     645\item{step 1:} The top level is defined in the same way as the bottom level is defined. 
     646\item{step 2:} The isolated grid point in the bathymetry are filled (as it is done in a domain without ice shelf) 
     647\item{step 3:} The tools make sure, the top level is above or equal to the bottom level 
     648\item{step 4:} If the water column at a U- or V- point is one wet cell wide, the ice shelf draft is adjusted. So the actual top cell become fully open and the new 
     649  top cell thickness is set to the minimum cell thickness allowed (following the same logic as for the bottom partial cell). This step is iterated 4 times to ensure the condition is fullfill along the 4 sides of the cell. 
     650\end{description} 
     651 
     652In case of steep slope and shallow water column, it likely that 2 cells are disconnected (bathymetry above its neigbourging ice shelf draft).  
     653The option \np{ln_isfconnect}{ln\_isfconnect} allow the tool to force the connection between these 2 cells. 
     654Some limiters in meter or levels on the digging allowed by the tool are available (respectively, \np{rn_zisfmax}{rn\_zisfmax} or \np{rn_kisfmax}{rn\_kisfmax}). 
     655This will prevent the formation of subglacial lakes at the expense of long vertical pipe to connect cells at very different levels. 
     656 
     657\subsection{Subglacial lakes} 
     658Despite careful setting of your ice shelf draft and bathymetry input file as well as setting described in \autoref{subsec:zgrisf_isfd}, some situation are unavoidable. 
     659For exemple if you setup your ice shelf draft and bathymetry to do ocean/ice sheet coupling,  
     660you may decide to fill the whole antarctic with a bathymetry and an ice shelf draft value (ice/bedrock interface depth when grounded).  
     661If you do so, the subglacial lakes will show up (Vostock for example). An other possibility is with coarse vertical resolution, some ice shelves could be cut in 2 parts:  
     662one connected to the main ocean and an other one closed which can be considered as a subglacial sea be the model.\\ 
     663 
     664The namelist option \np{ln_isfsubgl}{ln\_isfsubgl} allow you to remove theses subglacial lakes. 
     665This may be useful for esthetical reason or for stability reasons: 
     666 
     667\begin{description} 
     668\item $\bullet$ In a subglacial lakes, in case of very weak circulation (often the case), the only heat flux is the conductive heat flux through the ice sheet.  
     669  This will lead to constant freezing until water reaches -20C.  
     670  This is one of the defitiency of the 3 equation melt formulation (for details on this formulation, see: \autoref{sec:isf}). 
     671\item $\bullet$ In case of coupling with an ice sheet model,  
     672  the ssh in the subglacial lakes and the main ocean could be very different (ssh initial adjustement for example),  
     673  and so if for any reason both a connected at some point, the model is likely to fall over.\\ 
     674\end{description} 
    678675 
    679676\section{Closed sea definition} 
     
    707704\end{listing} 
    708705 
    709 The options available to define the closed seas and how closed sea net fresh water input will be redistributed by NEMO are listed in \nam{clo}{dom_clo} (\texttt{DOMAINcfg} only). 
     706The options available to define the closed seas and how closed sea net fresh water input will be redistributed by NEMO are listed in \nam{dom_clo}{dom\_clo} (\texttt{DOMAINcfg} only). 
    710707The individual definition of each closed sea is managed by \np{sn_lake}{sn\_lake}. In this fields the user needs to define:\\ 
    711708   \begin{description} 
  • NEMO/branches/2021/dev_r14312_MPI_Interface/doc/latex/NEMO/subfiles/chap_SBC.tex

    r14303 r14403  
    11791179%% ================================================================================================= 
    11801180\section[Ice Shelf (ISF)]{Interaction with ice shelves (ISF)} 
    1181 \label{sec:isf} 
     1181\label{sec:SBC_isf} 
    11821182 
    11831183\begin{listing} 
  • NEMO/branches/2021/dev_r14312_MPI_Interface/doc/latex/NEMO/subfiles/chap_TRA.tex

    r14303 r14403  
    733733  (see \autoref{sec:SBC_rnf} for further detail of how it acts on temperature and salinity tendencies) 
    734734\item [\textit{fwfisf}] The mass flux associated with ice shelf melt, 
    735   (see \autoref{sec:isf} for further details on how the ice shelf melt is computed and applied). 
     735  (see \autoref{sec:SBC_isf} for further details on how the ice shelf melt is computed and applied). 
    736736\end{labeling} 
    737737 
Note: See TracChangeset for help on using the changeset viewer.