2020-04-17T17:06:11+02:00 (8 months ago)

#2444: update of the documentation (prior to changes suggested by Dave)

1 edited


  • NEMO/branches/2020/ticket_2444/doc/latex/NEMO/subfiles/apdx_DOMAINcfg.tex

    r11693 r12769  
    319319  Defining the bathymetry also defines the coastline: where the bathymetry is zero, 
    320320  no wet levels are defined (all levels are masked). 
    322   The \ifile{isfdraft\_meter} file (Netcdf format) provides the ice shelf draft (positive, in meters) at 
    323   each grid point of the model grid. 
    324   This file is only needed if \np[=.true.]{ln_isfcav}{ln\_isfcav}. 
    325   Defining the ice shelf draft will also define the ice shelf edge and the grounding line position. 
    530525This option is described in the Report by Levier \textit{et al.} (2007), available on the \NEMO\ web site. 
     527\section{Ice shelf cavity definition} 
     530  If the under ice shelf seas are opened (\np{ln_isfcav}{ln\_isfcav}), the depth of the ice shelf/ocean interface has to be include in  
     531  the \ifile{isfdraft\_meter} file (Netcdf format). This file need to include the \ifile{isf\_draft} variable.  
     532  A positive value will me an ice shelf/ocean or ice shelf bedrock interface below the reference 0m ssh.  
     533  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{lst:namzgr_isf}. 
     537  \nlst{namzgr_isf_domcfg} 
     538  \caption{\forcode{&namzgr_isf}} 
     539  \label{lst:namzgr_isf} 
     542   The options available to define the shape of the under ice shelf cavities are listed in \nam{namzgr_isf}{namzgr\_isf} (\texttt{DOMAINcfg} only, \autoref{lst:namzgr_isf}). 
     544   \subsection{Model ice shelf draft definition} 
     545   \label{subsec:zgrisf_isfd} 
     547   First of all, the tool make sure, the ice shelf draft ($h_{isf}$) is sensible and compatible with the bathymetry. 
     548   There are 3 compulsory steps to achieve this: 
     550   \begin{description} 
     551   \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)$.  
     552   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}. 
     554   \item{\np{rn_glhw_min}{rn\_glhw\_min}:} This parameter is used to define the grounding line position. 
     555   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).  
     556   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, ...). 
     558   \item{\np{rn_isfhw_min}{rn\_isfhw\_min}:} This parameter is minimum water column thickness in the cavity.  
     559   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.  
     560   If for any reason, this adjustement break the minimum ice shelf draft allowed (\np{rn_isfdep_min}{rn\_isfdep\_min}), the cell is masked. 
     561   \end{description} 
     563   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}.   
     565   \subsection{Model top level definition} 
     566   After the definition of the ice shelf draft, the tool defines the top level.  
     567   The compulsory criterion is that the water column needs at least 2 wet cells in the water column at U- and V-points. 
     568   To do so, if there one cell wide water column, the tools adjust the ice shelf draft to fillful the requierement.\\ 
     570   The process is the following: 
     571   \begin{description} 
     572   \item{step 1:} The top level is defined in the same way as the bottom level is defined. 
     573   \item{step 2:} The isolated grid point in the bathymetry are filled (as it is done in a domain without ice shelf) 
     574   \item{step 3:} The tools make sure, the top level is above or equal to the bottom level 
     575   \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 
     576   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. 
     577   \end{description} 
     579   In case of steep slope and shallow water column, it likely that 2 cells are disconnected (bathymetry above its neigbourg ice shelf draft).  
     580   The option \np{ln_isfconnect}{ln\_isfconnect} allow the tool to force the connection between these 2 cells. 
     581   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}). 
     582   This will prevent the formation of subglacial lakes at the expense of long vertical pipe to connect cells at very different levels. 
     584   \subsection{Subglacial lakes} 
     585   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. 
     586   For exemple if you setup your ice shelf draft and bathymetry to do ocean/ice sheet coupling,  
     587   you may decide to fill the whole antarctic with a bathymetry and an ice shelf draft value (ice/bedrock interface depth when grounded).  
     588   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:  
     589   one connected to the main ocean and an other one closed which can be considered as a subglacial sea be the model.\\ 
     591   The namelist option \np{ln_isfsubgl}{ln\_isfsubgl} allow you to remove theses subglacial lakes. 
     592   This may be useful for esthetical reason or for stability reasons: 
     594   \begin{description} 
     595   \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.  
     596         This will lead to constant freezing until water reaches -20C.  
     597         This is one of the defitiency of the 3 equation melt formulation (for details on this formulation, see: \autoref{sec:isf}). 
     598   \item $\bullet$ In case of coupling with an ice sheet model,  
     599         the ssh in the subglacial lakes and the main ocean could be very different (ssh initial adjustement for example),  
     600         and so if for any reason both a connected at some point, the model is likely to fall over.\\ 
     601   \end{description} 
     603\section{Closed sea definition} 
     607  \nlst{namclo_domcfg} 
     608  \caption{\forcode{&namclo}} 
     609  \label{lst:namclo} 
     612The 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} (\texttt{DOMAINcfg} only, \autoref{lst:namclo}). 
     613The individual definition of each closed sea is managed by \np{sn_lake}{sn\_lake}. In this fields the user needs to defined:\\ 
     614   \begin{description} 
     615   \item $\bullet$    the name of the closed sea (print output purposes). 
     616   \item $\bullet$    the seed location to define the area of the closed sea (if seed on land because not present in this configuration, this closed sea will be ignored).\\ 
     617   \item $\bullet$    the seed location for the target area. 
     618   \item $\bullet$    the type of target area ('local','coast' or 'global'). See point 6 for definition of these cases. 
     619   \item $\bullet$    the type of redistribution scheme for the net fresh water flux over the closed sea (as a runoff in a target area, as emp in a target area, as emp globally). For the runoff case, if the net fwf is negative, it will be redistribut globally. 
     620   \item $\bullet$    the radius of the target area (not used for the 'global' case). So the target defined by a 'local' target area of a radius of 100km, for example, correspond to all the wet points within this radius. The coastal case will return only the coastal point within the specifid radius. 
     621   \item $\bullet$    the target id. This target id is used to group multiple lakes into the same river ouflow (Great Lakes for example). 
     622   \end{description} 
     624The closed sea module defines a number of masks in the \ifile{domain\_cfg} output: 
     625   \begin{description} 
     626   \item[\textit{mask\_opensea}:] a mask of the main ocean without all the closed seas closed. This mask is defined by a flood filling algorithm with an initial seed (localisation defined by \np{rn_lon_opnsea}{rn\_lon\_opnsea} and \np{rn_lat_opnsea}{rn\_lat\_opnsea}). 
     627   \item[\textit{mask\_csglo}, \textit{mask\_csrnf}, \textit{mask\_csemp}:] a mask of all the closed seas defined in the namelist by \np{sn_lake}{sn\_lake} for each redistribution scheme. The total number of defined closed seas has to be defined in \np{nn_closea}{nn\_closea}. 
     628   \item[\textit{mask\_csgrpglo}, \textit{mask\_csgrprnf}, \textit{mask\_csgrpemp}:] a mask of all the closed seas and targets grouped by target id for each type of redistribution scheme. 
     629   \item[\textit{mask\_csundef}:] a mask of all the closed sea not defined in \np{sn_lake}{sn\_lake}. This will allows NEMO to mask them if needed or to inform the user of potential minor issues in its bathymetry. 
     630   \end{description} 
Note: See TracChangeset for help on using the changeset viewer.