diff --git a/latex/NEMO/subfiles/chap_misc.tex b/latex/NEMO/subfiles/chap_misc.tex index 71abb9585b85b749501ee90f72fe67d89d4b6f06..9a3cd47de79acff2cc99eb239d09990727e84bdf 100644 --- a/latex/NEMO/subfiles/chap_misc.tex +++ b/latex/NEMO/subfiles/chap_misc.tex @@ -175,70 +175,6 @@ These points can be masked using the namelist option \np{ln_mask_csundef}{ln\_ma The masks needed for the closed sea can be created using the DOMAINcfg tool in the utils/tools/DOMAINcfg directory. See \autoref{sec:clocfg} for details on the usage of definition of the closed sea masks. -%% ================================================================================================= -\section{Sub-domain functionality} -\label{sec:MISC_zoom} - -%% ================================================================================================= -\subsection{Simple subsetting of input files via NetCDF attributes} - -The extended grids for use with the under-shelf ice cavities will result in redundant rows -around Antarctica if the ice cavities are not active. A simple mechanism for subsetting -input files associated with the extended domains has been implemented to avoid the need to -maintain different sets of input fields for use with or without active ice cavities. This -subsetting operates for the j-direction only and works by optionally looking for and using -a global file attribute (named: \np{open_ocean_jstart}{open\_ocean\_jstart}) to determine the starting j-row -for input. The use of this option is best explained with an example: -\medskip - -\noindent Consider an ORCA1 -configuration using the extended grid domain configuration file: \textit{eORCA1\_domcfg.nc} -This file define a horizontal domain of 362x332. The first row with -open ocean wet points in the non-isf bathymetry for this set is row 42 (\fortran\ indexing) -then the formally correct setting for \np{open_ocean_jstart}{open\_ocean\_jstart} is 41. Using this value as -the first row to be read will result in a 362x292 domain which is the same size as the -original ORCA1 domain. Thus the extended domain configuration file can be used with all -the original input files for ORCA1 if the ice cavities are not active (\np[=.false.]{ln_isfcav}{ln\_isfcav}). Full instructions for achieving this are: - -\begin{itemize} -\item Add the new attribute to any input files requiring a j-row offset, i.e: -\begin{cmds} -ncatted -a open_ocean_jstart,global,a,d,41 eORCA1_domcfg.nc -\end{cmds} - -\item Add the logical switch \np{ln_use_jattr}{ln\_use\_jattr} to \nam{cfg}{cfg} in the configuration -namelist (if it is not already there) and set \forcode{.true.} -\end{itemize} - -\noindent Note that with this option, the j-size of the global domain is (extended -j-size minus \np{open_ocean_jstart}{open\_ocean\_jstart} + 1 ) and this must match the \texttt{jpjglo} value -for the configuration. This means an alternative version of \textit{eORCA1\_domcfg.nc} must -be created for when \np{ln_use_jattr}{ln\_use\_jattr} is active. The \texttt{ncap2} tool provides a -convenient way of achieving this: - -\begin{cmds} -ncap2 -s 'jpjglo=292' eORCA1_domcfg.nc nORCA1_domcfg.nc -\end{cmds} - -The domain configuration file is unique in this respect since it also contains the value of \texttt{jpjglo} -that is read and used by the model. -Any other global, 2D and 3D, netcdf, input field can be prepared for use in a reduced domain by adding the -\texttt{open\_ocean\_jstart} attribute to the file's global attributes. -In particular this is true for any field that is read by \NEMO\ using the following optional argument to -the appropriate call to \np{iom_get}{iom\_get}. - -\begin{forlines} -lrowattr=ln_use_jattr -\end{forlines} - -Currently, only the domain configuration variables make use of this optional argument so -this facility is of little practical use except for tests where no other external input -files are needed or you wish to use an extended domain configuration with inputs from -earlier, non-extended configurations. Alternatively, it should be possible to exclude -empty rows for extended domain, forced ocean runs using interpolation on the fly, by -adding the optional argument to \texttt{iom\_get} calls for the weights and initial -conditions. Experimenting with this remains an exercise for the user. - %% ================================================================================================= \section[Accuracy and reproducibility (\textit{lib\_fortran.F90})]{Accuracy and reproducibility (\protect\mdl{lib\_fortran})} \label{sec:MISC_fortran}