Skip to content
Snippets Groups Projects
Commit e370ce8d authored by Andrew Coward's avatar Andrew Coward
Browse files

Removed Sub-domain functionality section from chap_misc.tex since this option...

Removed Sub-domain functionality section from chap_misc.tex since this option has been removed (practically at v4.2 and completely at v5.0)
parent a7def739
1 merge request!16Resolve "Removal of ln_use_jattr option"
Pipeline #2857 failed with stage
in 6 minutes and 21 seconds
......@@ -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}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment