Changeset 9407 for branches/2017/dev_merge_2017/DOC/tex_sub/chap_SBC.tex
- Timestamp:
- 2018-03-15T17:40:35+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/DOC/tex_sub/chap_SBC.tex
r9394 r9407 5 5 % ================================================================ 6 6 \chapter{Surface Boundary Condition (SBC, ISF, ICB) } 7 \label{ SBC}7 \label{chap:SBC} 8 8 \minitoc 9 9 … … 40 40 need not be supplied on the model grid. Instead a file of coordinates and weights can 41 41 be 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}). 43 43 If the Interpolation on the Fly option is used, input data belonging to land points (in the native grid), 44 44 can be masked to avoid spurious results in proximity of the coasts as large sea-land gradients characterize … … 65 65 Next the scheme for interpolation on the fly is described. 66 66 Finally, 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}),67 One of these is modification by icebergs (see \autoref{sec: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 \autoref{sec:SBC_isf}), 69 69 which provides additional sources of fresh water. 70 70 … … 74 74 % ================================================================ 75 75 \section{Surface boundary condition for the ocean} 76 \label{ SBC_general}76 \label{sec:SBC_general} 77 77 78 78 The surface ocean stress is the stress exerted by the wind and the sea-ice 79 79 on 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}).80 computation of the momentum vertical mixing trend (see \autoref{eq:dynzdf_sbc} in \autoref{sec:DYN_zdf}). 81 81 As such, it has to be provided as a 2D vector interpolated 82 82 onto the horizontal velocity ocean mesh, $i.e.$ resolved onto the model … … 88 88 plus the heat content of the mass exchange with the atmosphere and sea-ice). 89 89 It 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}).90 the first level temperature time evolution equation (see \autoref{eq:tra_sbc} 91 and \autoref{eq:tra_sbc_lin} in \autoref{subsec:TRA_sbc}). 92 92 The latter is the penetrative part of the heat flux. It is applied as a 3D 93 93 trends of the temperature equation (\mdl{traqsr} module) when \np{ln\_traqsr}\forcode{ = .true.}. 94 94 The way the light penetrates inside the water column is generally a sum of decreasing 95 exponentials (see \ S\ref{TRA_qsr}).95 exponentials (see \autoref{subsec:TRA_qsr}). 96 96 97 97 The surface freshwater budget is provided by the \textit{emp} field. … … 130 130 The ocean model provides, at each time step, to the surface module (\mdl{sbcmod}) 131 131 the surface currents, temperature and salinity. 132 These variables are averaged over \np{nn\_fsbc} time-step (\ ref{Tab_ssm}),132 These variables are averaged over \np{nn\_fsbc} time-step (\autoref{tab:ssm}), 133 133 and it is these averaged fields which are used to computes the surface fluxes 134 134 at a frequency of \np{nn\_fsbc} time-step. … … 144 144 Sea surface salinty & sss\_m & $psu$ & T \\ \hline 145 145 \end{tabular} 146 \caption{ \protect\label{ Tab_ssm}146 \caption{ \protect\label{tab:ssm} 147 147 Ocean variables provided by the ocean to the surface module (SBC). 148 148 The variable are averaged over nn{\_}fsbc time step, … … 158 158 % ================================================================ 159 159 \section{Input data generic interface} 160 \label{ SBC_input}160 \label{sec:SBC_input} 161 161 162 162 A generic interface has been introduced to manage the way input data (2D or 3D fields, … … 181 181 182 182 The 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-fly184 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 184 interpolation is used, a file of weights must be supplied (see \autoref{subsec:SBC_iof}). 185 185 186 186 Note that when an input data is archived on a disc which is accessible directly … … 193 193 % ------------------------------------------------------------------------------------------------------------- 194 194 \subsection{Input data specification (\protect\mdl{fldread})} 195 \label{ SBC_fldread}195 \label{subsec:SBC_fldread} 196 196 197 197 The structure associated with an input variable contains the following information: … … 205 205 This stem will be completed automatically by the model, with the addition of a '.nc' at its end 206 206 and 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 whether207 Tab.\autoref{tab:fldread} provides the resulting file name in all possible cases according to whether 208 208 it is a climatological file or not, and to the open/close frequency (see below for definition). 209 209 … … 218 218 \end{tabular} 219 219 \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, 221 221 as a function of the Open/close frequency. The stem name is assumed to be 'fn'. 222 222 For 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 … … 259 259 260 260 \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}.261 which is described in \autoref{subsec:SBC_iof}. 262 262 263 263 \end{description} … … 301 301 % ------------------------------------------------------------------------------------------------------------- 302 302 \subsection{Interpolation on-the-fly} 303 \label{ SBC_iof}303 \label{subsec:SBC_iof} 304 304 305 305 Interpolation on the Fly allows the user to supply input files required … … 325 325 326 326 \subsubsection{Bilinear interpolation} 327 \label{ SBC_iof_bilinear}327 \label{subsec:SBC_iof_bilinear} 328 328 329 329 The input weights file in this case has two sets of variables: src01, src02, … … 347 347 348 348 \subsubsection{Bicubic interpolation} 349 \label{ SBC_iof_bicubic}349 \label{subsec:SBC_iof_bicubic} 350 350 351 351 Again there are two sets of variables: "src" and "wgt". … … 363 363 364 364 \subsubsection{Implementation} 365 \label{ SBC_iof_imp}365 \label{subsec:SBC_iof_imp} 366 366 367 367 To activate this option, a non-empty string should be supplied in the weights filename column … … 398 398 399 399 \subsubsection{Limitations} 400 \label{ SBC_iof_lim}400 \label{subsec:SBC_iof_lim} 401 401 402 402 \begin{enumerate} … … 412 412 413 413 \subsubsection{Utilities} 414 \label{ SBC_iof_util}414 \label{subsec:SBC_iof_util} 415 415 416 416 % to be completed … … 422 422 % ------------------------------------------------------------------------------------------------------------- 423 423 \subsection{Standalone surface boundary condition scheme} 424 \label{ SAS_iof}424 \label{subsec:SAS_iof} 425 425 426 426 %---------------------------------------namsbc_ana-------------------------------------------------- … … 482 482 % ================================================================ 483 483 \section{Analytical formulation (\protect\mdl{sbcana})} 484 \label{ SBC_ana}484 \label{sec:SBC_ana} 485 485 486 486 %---------------------------------------namsbc_ana-------------------------------------------------- … … 506 506 % ================================================================ 507 507 \section{Flux formulation (\protect\mdl{sbcflx})} 508 \label{ SBC_flx}508 \label{sec:SBC_flx} 509 509 %------------------------------------------namsbc_flx---------------------------------------------------- 510 510 \forfile{../namelists/namsbc_flx} … … 516 516 read in the file, the time frequency at which it is given (in hours), and a logical 517 517 setting 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.518 for this field. See \autoref{subsec:SBC_fldread} for a more detailed description of the parameters. 519 519 520 520 Note 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 its521 restoring term to observed SST and/or SSS. See \autoref{subsec:SBC_ssr} for its 522 522 specification. 523 523 … … 528 528 \section[Bulk formulation {(\textit{sbcblk\{\_core,\_clio,\_mfs\}.F90})}] 529 529 {Bulk formulation {(\protect\mdl{sbcblk\_core}, \protect\mdl{sbcblk\_clio}, \protect\mdl{sbcblk\_mfs})}} 530 \label{ SBC_blk}530 \label{sec:SBC_blk} 531 531 532 532 In the bulk formulation, the surface boundary condition fields are computed … … 545 545 % ------------------------------------------------------------------------------------------------------------- 546 546 \subsection{CORE formulea (\protect\mdl{sbcblk\_core}, \protect\np{ln\_core}\forcode{ = .true.})} 547 \label{ SBC_blk_core}547 \label{subsec:SBC_blk_core} 548 548 %------------------------------------------namsbc_core---------------------------------------------------- 549 549 %\forfile{../namelists/namsbc_core} … … 566 566 567 567 %--------------------------------------------------TABLE-------------------------------------------------- 568 \begin{table}[htbp] \label{ Tab_CORE}568 \begin{table}[htbp] \label{tab:CORE} 569 569 \begin{center} 570 570 \begin{tabular}{|l|c|c|c|} … … 609 609 % ------------------------------------------------------------------------------------------------------------- 610 610 \subsection{CLIO formulea (\protect\mdl{sbcblk\_clio}, \protect\np{ln\_clio}\forcode{ = .true.})} 611 \label{ SBC_blk_clio}611 \label{subsec:SBC_blk_clio} 612 612 %------------------------------------------namsbc_clio---------------------------------------------------- 613 613 %\forfile{../namelists/namsbc_clio} … … 623 623 624 624 %--------------------------------------------------TABLE-------------------------------------------------- 625 \begin{table}[htbp] \label{ Tab_CLIO}625 \begin{table}[htbp] \label{tab:CLIO} 626 626 \begin{center} 627 627 \begin{tabular}{|l|l|l|l|} … … 643 643 As for the flux formulation, information about the input data required by the 644 644 model is provided in the namsbc\_blk\_core or namsbc\_blk\_clio 645 namelist (see \ S\ref{SBC_fldread}).645 namelist (see \autoref{subsec:SBC_fldread}). 646 646 647 647 % ------------------------------------------------------------------------------------------------------------- … … 649 649 % ------------------------------------------------------------------------------------------------------------- 650 650 \subsection{MFS formulea (\protect\mdl{sbcblk\_mfs}, \protect\np{ln\_mfs}\forcode{ = .true.})} 651 \label{ SBC_blk_mfs}651 \label{subsec:SBC_blk_mfs} 652 652 %------------------------------------------namsbc_mfs---------------------------------------------------- 653 653 %\forfile{../namelists/namsbc_mfs} … … 687 687 % ================================================================ 688 688 \section{Coupled formulation (\protect\mdl{sbccpl})} 689 \label{ SBC_cpl}689 \label{sec:SBC_cpl} 690 690 %------------------------------------------namsbc_cpl---------------------------------------------------- 691 691 \forfile{../namelists/namsbc_cpl} … … 725 725 % ================================================================ 726 726 \section{Atmospheric pressure (\protect\mdl{sbcapr})} 727 \label{ SBC_apr}727 \label{sec:SBC_apr} 728 728 %------------------------------------------namsbc_apr---------------------------------------------------- 729 729 \forfile{../namelists/namsbc_apr} … … 737 737 pressure is further transformed into an equivalent inverse barometer sea surface height, 738 738 $\eta_{ib}$, using: 739 \begin{equation} \label{ SBC_ssh_ib}739 \begin{equation} \label{eq:SBC_ssh_ib} 740 740 \eta_{ib} = - \frac{1}{g\,\rho_o} \left( P_{atm} - P_o \right) 741 741 \end{equation} … … 759 759 % ================================================================ 760 760 \section{Tidal potential (\protect\mdl{sbctide})} 761 \label{ SBC_tide}761 \label{sec:SBC_tide} 762 762 763 763 %------------------------------------------nam_tide--------------------------------------- … … 814 814 % ================================================================ 815 815 \section{River runoffs (\protect\mdl{sbcrnf})} 816 \label{ SBC_rnf}816 \label{sec:SBC_rnf} 817 817 %------------------------------------------namsbc_rnf---------------------------------------------------- 818 818 \forfile{../namelists/namsbc_rnf} … … 826 826 %coastal modelling and becomes more and more often open ocean and climate modelling 827 827 %\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}.}. 829 829 830 830 … … 847 847 more common in open ocean and climate modelling 848 848 \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}.}.849 required to properly represent the diurnal cycle \citep{Bernie_al_JC05}. see also \autoref{fig:SBC_dcy}.}. 850 850 851 851 As such from V~3.3 onwards it is possible to add river runoff through a non-zero depth, and for the … … 929 929 % ================================================================ 930 930 \section{Ice shelf melting (\protect\mdl{sbcisf})} 931 \label{ SBC_isf}931 \label{sec:SBC_isf} 932 932 %------------------------------------------namsbc_isf---------------------------------------------------- 933 933 \forfile{../namelists/namsbc_isf} … … 1006 1006 The fw addition due to the ice shelf melting is, at each relevant depth level, added to the horizontal divergence 1007 1007 (\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.1008 See the runoff section \autoref{sec:SBC_rnf} for all the details about the divergence correction. 1009 1009 1010 1010 1011 1011 \section{Ice sheet coupling} 1012 \label{ SBC_iscpl}1012 \label{sec:SBC_iscpl} 1013 1013 %------------------------------------------namsbc_iscpl---------------------------------------------------- 1014 1014 \forfile{../namelists/namsbc_iscpl} … … 1048 1048 % ================================================================ 1049 1049 \section{Handling of icebergs (ICB)} 1050 \label{ ICB_icebergs}1050 \label{sec:ICB_icebergs} 1051 1051 %------------------------------------------namberg---------------------------------------------------- 1052 1052 \forfile{../namelists/namberg} … … 1113 1113 % ================================================================ 1114 1114 \section{Miscellaneous options} 1115 \label{ SBC_misc}1115 \label{sec:SBC_misc} 1116 1116 1117 1117 % ------------------------------------------------------------------------------------------------------------- … … 1119 1119 % ------------------------------------------------------------------------------------------------------------- 1120 1120 \subsection{Diurnal cycle (\protect\mdl{sbcdcy})} 1121 \label{ SBC_dcy}1121 \label{subsec:SBC_dcy} 1122 1122 %------------------------------------------namsbc_rnf---------------------------------------------------- 1123 1123 %\forfile{../namelists/namsbc} … … 1127 1127 \begin{figure}[!t] \begin{center} 1128 1128 \includegraphics[width=0.8\textwidth]{Fig_SBC_diurnal} 1129 \caption{ \protect\label{ Fig_SBC_diurnal}1129 \caption{ \protect\label{fig:SBC_diurnal} 1130 1130 Example of recontruction of the diurnal cycle variation of short wave flux 1131 1131 from daily mean values. The reconstructed diurnal cycle (black line) is chosen … … 1149 1149 can be found in the appendix~A of \cite{Bernie_al_CD07}. The algorithm preserve the daily 1150 1150 mean 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}).1151 of the analytical cycle over this time step (\autoref{fig:SBC_diurnal}). 1152 1152 The use of diurnal cycle reconstruction requires the input SWF to be daily 1153 1153 ($i.e.$ a frequency of 24 and a time interpolation set to true in \np{sn\_qsr} namelist parameter). 1154 1154 Furthermore, it is recommended to have a least 8 surface module time step per day, 1155 1155 that 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~hours1156 is given in \autoref{fig:SBC_dcy} for a 12 reconstructed diurnal cycle, one every 2~hours 1157 1157 (from 1am to 11pm). 1158 1158 … … 1160 1160 \begin{figure}[!t] \begin{center} 1161 1161 \includegraphics[width=0.7\textwidth]{Fig_SBC_dcy} 1162 \caption{ \protect\label{ Fig_SBC_dcy}1162 \caption{ \protect\label{fig:SBC_dcy} 1163 1163 Example of recontruction of the diurnal cycle variation of short wave flux 1164 1164 from daily mean values on an ORCA2 grid with a time sampling of 2~hours (from 1am to 11pm). … … 1176 1176 % ------------------------------------------------------------------------------------------------------------- 1177 1177 \subsection{Rotation of vector pairs onto the model grid directions} 1178 \label{ SBC_rotation}1178 \label{subsec:SBC_rotation} 1179 1179 1180 1180 When using a flux (\np{ln\_flx}\forcode{ = .true.}) or bulk (\np{ln\_clio}\forcode{ = .true.} or \np{ln\_core}\forcode{ = .true.}) formulation, … … 1195 1195 % ------------------------------------------------------------------------------------------------------------- 1196 1196 \subsection{Surface restoring to observed SST and/or SSS (\protect\mdl{sbcssr})} 1197 \label{ SBC_ssr}1197 \label{subsec:SBC_ssr} 1198 1198 %------------------------------------------namsbc_ssr---------------------------------------------------- 1199 1199 \forfile{../namelists/namsbc_ssr} … … 1203 1203 n forced mode using a flux formulation (\np{ln\_flx}\forcode{ = .true.}), a 1204 1204 feedback 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} 1206 1206 Q_{ns} = Q_{ns}^o + \frac{dQ}{dT} \left( \left. T \right|_{k=1} - SST_{Obs} \right) 1207 1207 \end{equation} … … 1216 1216 equivalent freshwater flux, it takes the following expression : 1217 1217 1218 \begin{equation} \label{ Eq_sbc_dmp_emp}1218 \begin{equation} \label{eq:sbc_dmp_emp} 1219 1219 \textit{emp} = \textit{emp}_o + \gamma_s^{-1} e_{3t} \frac{ \left(\left.S\right|_{k=1}-SSS_{Obs}\right)} 1220 1220 {\left.S\right|_{k=1}} … … 1226 1226 $\left.S\right|_{k=1}$ is the model surface layer salinity and $\gamma_s$ is a negative 1227 1227 feedback 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 atmosphere1228 physical justification for the feedback term in \autoref{eq:sbc_dmp_emp} as the atmosphere 1229 1229 does not care about ocean surface salinity \citep{Madec1997}. The SSS restoring 1230 1230 term should be viewed as a flux correction on freshwater fluxes to reduce the … … 1235 1235 % ------------------------------------------------------------------------------------------------------------- 1236 1236 \subsection{Handling of ice-covered area (\textit{sbcice\_...})} 1237 \label{ SBC_ice-cover}1237 \label{subsec:SBC_ice-cover} 1238 1238 1239 1239 The presence at the sea surface of an ice covered area modifies all the fluxes … … 1264 1264 1265 1265 \subsection{Interface to CICE (\protect\mdl{sbcice\_cice})} 1266 \label{ SBC_cice}1266 \label{subsec:SBC_cice} 1267 1267 1268 1268 It is now possible to couple a regional or global NEMO configuration (without AGRIF) to the CICE sea-ice … … 1291 1291 % ------------------------------------------------------------------------------------------------------------- 1292 1292 \subsection{Freshwater budget control (\protect\mdl{sbcfwb})} 1293 \label{ SBC_fwb}1293 \label{subsec:SBC_fwb} 1294 1294 1295 1295 For global ocean simulation it can be useful to introduce a control of the mean sea … … 1313 1313 \subsection[Neutral drag coeff. from external wave model (\protect\mdl{sbcwave})] 1314 1314 {Neutral drag coefficient from external wave model (\protect\mdl{sbcwave})} 1315 \label{ SBC_wave}1315 \label{subsec:SBC_wave} 1316 1316 %------------------------------------------namwave---------------------------------------------------- 1317 1317 \forfile{../namelists/namsbc_wave} … … 1322 1322 The \mdl{sbcwave} module containing the routine \np{sbc\_wave} reads the 1323 1323 namelist \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})1324 the miscellanous options allowed by Input Data generic Interface see \autoref{sec:SBC_input}) 1325 1325 and a 2D field of neutral drag coefficient. 1326 1326 Then 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.