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 9407 for branches/2017/dev_merge_2017/DOC/tex_sub/chap_SBC.tex – NEMO

Ignore:
Timestamp:
2018-03-15T17:40:35+01:00 (6 years ago)
Author:
nicolasmartin
Message:

Complete refactoring of cross-referencing

  • Use of \autoref instead of simple \ref for contextual text depending on target type
  • creation of few prefixes for marker to identify the type reference: apdx|chap|eq|fig|sec|subsec|tab
File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_merge_2017/DOC/tex_sub/chap_SBC.tex

    r9394 r9407  
    55% ================================================================ 
    66\chapter{Surface Boundary Condition (SBC, ISF, ICB) } 
    7 \label{SBC} 
     7\label{chap:SBC} 
    88\minitoc 
    99 
     
    4040need not be supplied on the model grid. Instead a file of coordinates and weights can  
    4141be supplied which maps the data from the supplied grid to the model points  
    42 (so called "Interpolation on the Fly", see \S\ref{SBC_iof}). 
     42(so called "Interpolation on the Fly", see \autoref{subsec:SBC_iof}). 
    4343If the Interpolation on the Fly option is used, input data belonging to land points (in the native grid), 
    4444can be masked to avoid spurious results in proximity of the coasts  as large sea-land gradients characterize 
     
    6565Next the scheme for interpolation on the fly is described. 
    6666Finally, the different options that further modify the fluxes applied to the ocean are discussed. 
    67 One of these is modification by icebergs (see \S\ref{ICB_icebergs}), which act as drifting sources of fresh water. 
    68 Another example of modification is that due to the ice shelf melting/freezing (see \S\ref{SBC_isf}),  
     67One of these is modification by icebergs (see \autoref{sec:ICB_icebergs}), which act as drifting sources of fresh water. 
     68Another example of modification is that due to the ice shelf melting/freezing (see \autoref{sec:SBC_isf}),  
    6969which provides additional sources of fresh water. 
    7070 
     
    7474% ================================================================ 
    7575\section{Surface boundary condition for the ocean} 
    76 \label{SBC_general} 
     76\label{sec:SBC_general} 
    7777 
    7878The surface ocean stress is the stress exerted by the wind and the sea-ice  
    7979on the ocean. It is applied in \mdl{dynzdf} module as a surface boundary condition of the  
    80 computation of the momentum vertical mixing trend (see \eqref{Eq_dynzdf_sbc} in \S\ref{DYN_zdf}). 
     80computation of the momentum vertical mixing trend (see \autoref{eq:dynzdf_sbc} in \autoref{sec:DYN_zdf}). 
    8181As such, it has to be provided as a 2D vector interpolated  
    8282onto the horizontal velocity ocean mesh, $i.e.$ resolved onto the model  
     
    8888plus the heat content of the mass exchange with the atmosphere and sea-ice).  
    8989It is applied in \mdl{trasbc} module as a surface boundary condition trend of  
    90 the first level temperature time evolution equation (see \eqref{Eq_tra_sbc}  
    91 and \eqref{Eq_tra_sbc_lin} in \S\ref{TRA_sbc}).  
     90the first level temperature time evolution equation (see \autoref{eq:tra_sbc}  
     91and \autoref{eq:tra_sbc_lin} in \autoref{subsec:TRA_sbc}).  
    9292The latter is the penetrative part of the heat flux. It is applied as a 3D  
    9393trends of the temperature equation (\mdl{traqsr} module) when \np{ln\_traqsr}\forcode{ = .true.}. 
    9494The way the light penetrates inside the water column is generally a sum of decreasing  
    95 exponentials (see \S\ref{TRA_qsr}).  
     95exponentials (see \autoref{subsec:TRA_qsr}).  
    9696 
    9797The surface freshwater budget is provided by the \textit{emp} field. 
     
    130130The ocean model provides, at each time step, to the surface module (\mdl{sbcmod})  
    131131the surface currents, temperature and salinity.   
    132 These variables are averaged over \np{nn\_fsbc} time-step (\ref{Tab_ssm}),  
     132These variables are averaged over \np{nn\_fsbc} time-step (\autoref{tab:ssm}),  
    133133and it is these averaged fields which are used to computes the surface fluxes  
    134134at a frequency of \np{nn\_fsbc} time-step. 
     
    144144Sea surface salinty              & sss\_m & $psu$        & T \\   \hline 
    145145\end{tabular} 
    146 \caption{  \protect\label{Tab_ssm}    
     146\caption{  \protect\label{tab:ssm}    
    147147Ocean variables provided by the ocean to the surface module (SBC).  
    148148The variable are averaged over nn{\_}fsbc time step,  
     
    158158% ================================================================ 
    159159\section{Input data generic interface} 
    160 \label{SBC_input} 
     160\label{sec:SBC_input} 
    161161 
    162162A generic interface has been introduced to manage the way input data (2D or 3D fields,  
     
    181181 
    182182The only constraints are that the input file is a NetCDF file, the file name follows a nomenclature  
    183 (see \S\ref{SBC_fldread}), the period it cover is one year, month, week or day, and, if on-the-fly  
    184 interpolation is used, a file of weights must be supplied (see \S\ref{SBC_iof}). 
     183(see \autoref{subsec:SBC_fldread}), the period it cover is one year, month, week or day, and, if on-the-fly  
     184interpolation is used, a file of weights must be supplied (see \autoref{subsec:SBC_iof}). 
    185185 
    186186Note that when an input data is archived on a disc which is accessible directly  
     
    193193% ------------------------------------------------------------------------------------------------------------- 
    194194\subsection{Input data specification (\protect\mdl{fldread})} 
    195 \label{SBC_fldread} 
     195\label{subsec:SBC_fldread} 
    196196 
    197197The structure associated with an input variable contains the following information: 
     
    205205This stem will be completed automatically by the model, with the addition of a '.nc' at its end  
    206206and by date information and possibly a prefix (when using AGRIF).  
    207 Tab.\ref{Tab_fldread} provides the resulting file name in all possible cases according to whether  
     207Tab.\autoref{tab:fldread} provides the resulting file name in all possible cases according to whether  
    208208it is a climatological file or not, and to the open/close frequency (see below for definition).  
    209209 
     
    218218\end{tabular} 
    219219\end{center} 
    220 \caption{ \protect\label{Tab_fldread}   naming nomenclature for climatological or interannual input file,  
     220\caption{ \protect\label{tab:fldread}   naming nomenclature for climatological or interannual input file,  
    221221as a function of the Open/close frequency. The stem name is assumed to be 'fn'.  
    222222For weekly files, the 'LLL' corresponds to the first three letters of the first day of the week ($i.e.$ 'sun','sat','fri','thu','wed','tue','mon'). The 'YYYY', 'MM' and 'DD' should be replaced by the  
     
    259259 
    260260\item[Others]: 'weights filename', 'pairing rotation' and 'land/sea mask' are associted with on-the-fly interpolation  
    261 which is described in \S\ref{SBC_iof}. 
     261which is described in \autoref{subsec:SBC_iof}. 
    262262 
    263263\end{description} 
     
    301301% ------------------------------------------------------------------------------------------------------------- 
    302302\subsection{Interpolation on-the-fly} 
    303 \label{SBC_iof} 
     303\label{subsec:SBC_iof} 
    304304 
    305305Interpolation on the Fly allows the user to supply input files required 
     
    325325 
    326326\subsubsection{Bilinear interpolation} 
    327 \label{SBC_iof_bilinear} 
     327\label{subsec:SBC_iof_bilinear} 
    328328 
    329329The input weights file in this case has two sets of variables: src01, src02, 
     
    347347 
    348348\subsubsection{Bicubic interpolation} 
    349 \label{SBC_iof_bicubic} 
     349\label{subsec:SBC_iof_bicubic} 
    350350 
    351351Again there are two sets of variables: "src" and "wgt". 
     
    363363 
    364364\subsubsection{Implementation} 
    365 \label{SBC_iof_imp} 
     365\label{subsec:SBC_iof_imp} 
    366366 
    367367To activate this option, a non-empty string should be supplied in the weights filename column  
     
    398398 
    399399\subsubsection{Limitations} 
    400 \label{SBC_iof_lim} 
     400\label{subsec:SBC_iof_lim} 
    401401 
    402402\begin{enumerate}   
     
    412412 
    413413\subsubsection{Utilities} 
    414 \label{SBC_iof_util} 
     414\label{subsec:SBC_iof_util} 
    415415 
    416416% to be completed 
     
    422422% ------------------------------------------------------------------------------------------------------------- 
    423423\subsection{Standalone surface boundary condition scheme} 
    424 \label{SAS_iof} 
     424\label{subsec:SAS_iof} 
    425425 
    426426%---------------------------------------namsbc_ana-------------------------------------------------- 
     
    482482% ================================================================ 
    483483\section{Analytical formulation (\protect\mdl{sbcana})} 
    484 \label{SBC_ana} 
     484\label{sec:SBC_ana} 
    485485 
    486486%---------------------------------------namsbc_ana-------------------------------------------------- 
     
    506506% ================================================================ 
    507507\section{Flux formulation (\protect\mdl{sbcflx})} 
    508 \label{SBC_flx} 
     508\label{sec:SBC_flx} 
    509509%------------------------------------------namsbc_flx---------------------------------------------------- 
    510510\forfile{../namelists/namsbc_flx}  
     
    516516read in the file, the time frequency at which it is given (in hours), and a logical  
    517517setting whether a time interpolation to the model time step is required  
    518 for this field. See \S\ref{SBC_fldread} for a more detailed description of the parameters. 
     518for this field. See \autoref{subsec:SBC_fldread} for a more detailed description of the parameters. 
    519519 
    520520Note that in general, a flux formulation is used in associated with a  
    521 restoring term to observed SST and/or SSS. See \S\ref{SBC_ssr} for its  
     521restoring term to observed SST and/or SSS. See \autoref{subsec:SBC_ssr} for its  
    522522specification. 
    523523 
     
    528528\section[Bulk formulation {(\textit{sbcblk\{\_core,\_clio,\_mfs\}.F90})}] 
    529529         {Bulk formulation {(\protect\mdl{sbcblk\_core}, \protect\mdl{sbcblk\_clio}, \protect\mdl{sbcblk\_mfs})}} 
    530 \label{SBC_blk} 
     530\label{sec:SBC_blk} 
    531531 
    532532In the bulk formulation, the surface boundary condition fields are computed  
     
    545545% ------------------------------------------------------------------------------------------------------------- 
    546546\subsection{CORE formulea (\protect\mdl{sbcblk\_core}, \protect\np{ln\_core}\forcode{ = .true.})} 
    547 \label{SBC_blk_core} 
     547\label{subsec:SBC_blk_core} 
    548548%------------------------------------------namsbc_core---------------------------------------------------- 
    549549%\forfile{../namelists/namsbc_core} 
     
    566566 
    567567%--------------------------------------------------TABLE-------------------------------------------------- 
    568 \begin{table}[htbp]   \label{Tab_CORE} 
     568\begin{table}[htbp]   \label{tab:CORE} 
    569569\begin{center} 
    570570\begin{tabular}{|l|c|c|c|} 
     
    609609% ------------------------------------------------------------------------------------------------------------- 
    610610\subsection{CLIO formulea (\protect\mdl{sbcblk\_clio}, \protect\np{ln\_clio}\forcode{ = .true.})} 
    611 \label{SBC_blk_clio} 
     611\label{subsec:SBC_blk_clio} 
    612612%------------------------------------------namsbc_clio---------------------------------------------------- 
    613613%\forfile{../namelists/namsbc_clio} 
     
    623623 
    624624%--------------------------------------------------TABLE-------------------------------------------------- 
    625 \begin{table}[htbp]   \label{Tab_CLIO} 
     625\begin{table}[htbp]   \label{tab:CLIO} 
    626626\begin{center} 
    627627\begin{tabular}{|l|l|l|l|} 
     
    643643As for the flux formulation, information about the input data required by the  
    644644model is provided in the namsbc\_blk\_core or namsbc\_blk\_clio  
    645 namelist (see \S\ref{SBC_fldread}).  
     645namelist (see \autoref{subsec:SBC_fldread}).  
    646646 
    647647% ------------------------------------------------------------------------------------------------------------- 
     
    649649% ------------------------------------------------------------------------------------------------------------- 
    650650\subsection{MFS formulea (\protect\mdl{sbcblk\_mfs}, \protect\np{ln\_mfs}\forcode{ = .true.})} 
    651 \label{SBC_blk_mfs} 
     651\label{subsec:SBC_blk_mfs} 
    652652%------------------------------------------namsbc_mfs---------------------------------------------------- 
    653653%\forfile{../namelists/namsbc_mfs} 
     
    687687% ================================================================ 
    688688\section{Coupled formulation (\protect\mdl{sbccpl})} 
    689 \label{SBC_cpl} 
     689\label{sec:SBC_cpl} 
    690690%------------------------------------------namsbc_cpl---------------------------------------------------- 
    691691\forfile{../namelists/namsbc_cpl}  
     
    725725% ================================================================ 
    726726\section{Atmospheric pressure (\protect\mdl{sbcapr})} 
    727 \label{SBC_apr} 
     727\label{sec:SBC_apr} 
    728728%------------------------------------------namsbc_apr---------------------------------------------------- 
    729729\forfile{../namelists/namsbc_apr}  
     
    737737pressure is further transformed into an equivalent inverse barometer sea surface height,  
    738738$\eta_{ib}$, using: 
    739 \begin{equation} \label{SBC_ssh_ib} 
     739\begin{equation} \label{eq:SBC_ssh_ib} 
    740740   \eta_{ib} = -  \frac{1}{g\,\rho_o}  \left( P_{atm} - P_o \right)  
    741741\end{equation} 
     
    759759% ================================================================ 
    760760\section{Tidal potential (\protect\mdl{sbctide})} 
    761 \label{SBC_tide} 
     761\label{sec:SBC_tide} 
    762762 
    763763%------------------------------------------nam_tide--------------------------------------- 
     
    814814% ================================================================ 
    815815\section{River runoffs (\protect\mdl{sbcrnf})} 
    816 \label{SBC_rnf} 
     816\label{sec:SBC_rnf} 
    817817%------------------------------------------namsbc_rnf---------------------------------------------------- 
    818818\forfile{../namelists/namsbc_rnf}  
     
    826826%coastal modelling and becomes more and more often open ocean and climate modelling  
    827827%\footnote{At least a top cells thickness of 1~meter and a 3 hours forcing frequency are 
    828 %required to properly represent the diurnal cycle \citep{Bernie_al_JC05}. see also \S\ref{SBC_dcy}.}. 
     828%required to properly represent the diurnal cycle \citep{Bernie_al_JC05}. see also \autoref{fig:SBC_dcy}.}. 
    829829 
    830830 
     
    847847more common in open ocean and climate modelling  
    848848\footnote{At least a top cells thickness of 1~meter and a 3 hours forcing frequency are 
    849 required to properly represent the diurnal cycle \citep{Bernie_al_JC05}. see also \S\ref{SBC_dcy}.}. 
     849required to properly represent the diurnal cycle \citep{Bernie_al_JC05}. see also \autoref{fig:SBC_dcy}.}. 
    850850 
    851851As such from V~3.3 onwards it is possible to add river runoff through a non-zero depth, and for the  
     
    929929% ================================================================ 
    930930\section{Ice shelf melting (\protect\mdl{sbcisf})} 
    931 \label{SBC_isf} 
     931\label{sec:SBC_isf} 
    932932%------------------------------------------namsbc_isf---------------------------------------------------- 
    933933\forfile{../namelists/namsbc_isf} 
     
    10061006The fw addition due to the ice shelf melting is, at each relevant depth level, added to the horizontal divergence  
    10071007(\textit{hdivn}) in the subroutine \rou{sbc\_isf\_div}, called from \mdl{divcur}.  
    1008 See the runoff section \ref{SBC_rnf} for all the details about the divergence correction.  
     1008See the runoff section \autoref{sec:SBC_rnf} for all the details about the divergence correction.  
    10091009 
    10101010 
    10111011\section{Ice sheet coupling} 
    1012 \label{SBC_iscpl} 
     1012\label{sec:SBC_iscpl} 
    10131013%------------------------------------------namsbc_iscpl---------------------------------------------------- 
    10141014\forfile{../namelists/namsbc_iscpl} 
     
    10481048% ================================================================ 
    10491049\section{Handling of icebergs (ICB)} 
    1050 \label{ICB_icebergs} 
     1050\label{sec:ICB_icebergs} 
    10511051%------------------------------------------namberg---------------------------------------------------- 
    10521052\forfile{../namelists/namberg} 
     
    11131113% ================================================================ 
    11141114\section{Miscellaneous options} 
    1115 \label{SBC_misc} 
     1115\label{sec:SBC_misc} 
    11161116 
    11171117% ------------------------------------------------------------------------------------------------------------- 
     
    11191119% ------------------------------------------------------------------------------------------------------------- 
    11201120\subsection{Diurnal cycle (\protect\mdl{sbcdcy})} 
    1121 \label{SBC_dcy} 
     1121\label{subsec:SBC_dcy} 
    11221122%------------------------------------------namsbc_rnf---------------------------------------------------- 
    11231123%\forfile{../namelists/namsbc}  
     
    11271127\begin{figure}[!t]    \begin{center} 
    11281128\includegraphics[width=0.8\textwidth]{Fig_SBC_diurnal} 
    1129 \caption{ \protect\label{Fig_SBC_diurnal}     
     1129\caption{ \protect\label{fig:SBC_diurnal}     
    11301130Example of recontruction of the diurnal cycle variation of short wave flux   
    11311131from daily mean values. The reconstructed diurnal cycle (black line) is chosen  
     
    11491149can be found in the appendix~A of \cite{Bernie_al_CD07}. The algorithm preserve the daily  
    11501150mean incomming SWF as the reconstructed SWF at a given time step is the mean value  
    1151 of the analytical cycle over this time step (Fig.\ref{Fig_SBC_diurnal}).  
     1151of the analytical cycle over this time step (\autoref{fig:SBC_diurnal}).  
    11521152The use of diurnal cycle reconstruction requires the input SWF to be daily  
    11531153($i.e.$ a frequency of 24 and a time interpolation set to true in \np{sn\_qsr} namelist parameter). 
    11541154Furthermore, it is recommended to have a least 8 surface module time step per day, 
    11551155that is  $\rdt \ nn\_fsbc < 10,800~s = 3~h$. An example of recontructed SWF  
    1156 is given in Fig.\ref{Fig_SBC_dcy} for a 12 reconstructed diurnal cycle, one every 2~hours  
     1156is given in \autoref{fig:SBC_dcy} for a 12 reconstructed diurnal cycle, one every 2~hours  
    11571157(from 1am to 11pm). 
    11581158 
     
    11601160\begin{figure}[!t]  \begin{center} 
    11611161\includegraphics[width=0.7\textwidth]{Fig_SBC_dcy} 
    1162 \caption{ \protect\label{Fig_SBC_dcy}    
     1162\caption{ \protect\label{fig:SBC_dcy}    
    11631163Example of recontruction of the diurnal cycle variation of short wave flux   
    11641164from daily mean values on an ORCA2 grid with a time sampling of 2~hours (from 1am to 11pm).  
     
    11761176% ------------------------------------------------------------------------------------------------------------- 
    11771177\subsection{Rotation of vector pairs onto the model grid directions} 
    1178 \label{SBC_rotation} 
     1178\label{subsec:SBC_rotation} 
    11791179 
    11801180When using a flux (\np{ln\_flx}\forcode{ = .true.}) or bulk (\np{ln\_clio}\forcode{ = .true.} or \np{ln\_core}\forcode{ = .true.}) formulation,  
     
    11951195% ------------------------------------------------------------------------------------------------------------- 
    11961196\subsection{Surface restoring to observed SST and/or SSS (\protect\mdl{sbcssr})} 
    1197 \label{SBC_ssr} 
     1197\label{subsec:SBC_ssr} 
    11981198%------------------------------------------namsbc_ssr---------------------------------------------------- 
    11991199\forfile{../namelists/namsbc_ssr}  
     
    12031203n forced mode using a flux formulation (\np{ln\_flx}\forcode{ = .true.}), a  
    12041204feedback term \emph{must} be added to the surface heat flux $Q_{ns}^o$: 
    1205 \begin{equation} \label{Eq_sbc_dmp_q} 
     1205\begin{equation} \label{eq:sbc_dmp_q} 
    12061206Q_{ns} = Q_{ns}^o + \frac{dQ}{dT} \left( \left. T \right|_{k=1} - SST_{Obs} \right) 
    12071207\end{equation} 
     
    12161216equivalent freshwater flux, it takes the following expression : 
    12171217 
    1218 \begin{equation} \label{Eq_sbc_dmp_emp} 
     1218\begin{equation} \label{eq:sbc_dmp_emp} 
    12191219\textit{emp} = \textit{emp}_o + \gamma_s^{-1} e_{3t}  \frac{  \left(\left.S\right|_{k=1}-SSS_{Obs}\right)} 
    12201220                                             {\left.S\right|_{k=1}} 
     
    12261226$\left.S\right|_{k=1}$ is the model surface layer salinity and $\gamma_s$ is a negative  
    12271227feedback coefficient which is provided as a namelist parameter. Unlike heat flux, there is no  
    1228 physical justification for the feedback term in \ref{Eq_sbc_dmp_emp} as the atmosphere  
     1228physical justification for the feedback term in \autoref{eq:sbc_dmp_emp} as the atmosphere  
    12291229does not care about ocean surface salinity \citep{Madec1997}. The SSS restoring  
    12301230term should be viewed as a flux correction on freshwater fluxes to reduce the  
     
    12351235% ------------------------------------------------------------------------------------------------------------- 
    12361236\subsection{Handling of ice-covered area  (\textit{sbcice\_...})} 
    1237 \label{SBC_ice-cover} 
     1237\label{subsec:SBC_ice-cover} 
    12381238 
    12391239The presence at the sea surface of an ice covered area modifies all the fluxes  
     
    12641264 
    12651265\subsection{Interface to CICE (\protect\mdl{sbcice\_cice})} 
    1266 \label{SBC_cice} 
     1266\label{subsec:SBC_cice} 
    12671267 
    12681268It is now possible to couple a regional or global NEMO configuration (without AGRIF) to the CICE sea-ice 
     
    12911291% ------------------------------------------------------------------------------------------------------------- 
    12921292\subsection{Freshwater budget control (\protect\mdl{sbcfwb})} 
    1293 \label{SBC_fwb} 
     1293\label{subsec:SBC_fwb} 
    12941294 
    12951295For global ocean simulation it can be useful to introduce a control of the mean sea  
     
    13131313\subsection[Neutral drag coeff. from external wave model (\protect\mdl{sbcwave})] 
    13141314            {Neutral drag coefficient from external wave model (\protect\mdl{sbcwave})} 
    1315 \label{SBC_wave} 
     1315\label{subsec:SBC_wave} 
    13161316%------------------------------------------namwave---------------------------------------------------- 
    13171317\forfile{../namelists/namsbc_wave} 
     
    13221322The \mdl{sbcwave} module containing the routine \np{sbc\_wave} reads the 
    13231323namelist \ngn{namsbc\_wave} (for external data names, locations, frequency, interpolation and all  
    1324 the miscellanous options allowed by Input Data generic Interface see \S\ref{SBC_input})  
     1324the miscellanous options allowed by Input Data generic Interface see \autoref{sec:SBC_input})  
    13251325and a 2D field of neutral drag coefficient.  
    13261326Then using the routine TURB\_CORE\_1Z or TURB\_CORE\_2Z, and starting from the neutral drag coefficent provided,  
Note: See TracChangeset for help on using the changeset viewer.