- Timestamp:
- 2016-07-19T10:38:35+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5549_BDY_ZEROGRAD/DOC/TexFiles/Chapters/Chap_DIA.tex
- Property svn:keywords set to Id
r5515 r6808 2 2 % Chapter I/O & Diagnostics 3 3 % ================================================================ 4 \chapter{Ou put and Diagnostics (IOM, DIA, TRD, FLO)}4 \chapter{Output and Diagnostics (IOM, DIA, TRD, FLO)} 5 5 \label{DIA} 6 6 \minitoc 7 7 8 8 \newpage 9 $\ $\newline % force a new li gne9 $\ $\newline % force a new line 10 10 11 11 % ================================================================ … … 33 33 "\textit{grep -i numout}" in the source code directory. 34 34 35 By default, diagnostic output files are written in NetCDF format but an IEEE binary output format, called DIMG, can be cho osen by defining \key{dimgout}.35 By default, diagnostic output files are written in NetCDF format but an IEEE binary output format, called DIMG, can be chosen by defining \key{dimgout}. 36 36 37 37 Since version 3.2, when defining \key{iomput}, an I/O server has been added which provides more flexibility in the choice of the fields to be written as well as how the writing work is distributed over the processors in massively parallel computing. The complete description of the use of this I/O server is presented in the next section. … … 57 57 \item The choice of output frequencies that can be different for each file (including real months and years). 58 58 \item The choice of file contents; includes complete flexibility over which data are written in which files (the same data can be written in different files). 59 \item The possibility to split output files at a cho osen frequency.59 \item The possibility to split output files at a chosen frequency. 60 60 \item The possibility to extract a vertical or an horizontal subdomain. 61 61 \item The choice of the temporal operation to perform, e.g.: average, accumulate, instantaneous, min, max and once. … … 1319 1319 1320 1320 % ------------------------------------------------------------------------------------------------------------- 1321 % 25 hour mean and hourly Surface, Mid and Bed 1322 % ------------------------------------------------------------------------------------------------------------- 1323 \section{25 hour mean output for tidal models } 1324 1325 A module is available to compute a crudely detided M2 signal by obtaining a 25 hour mean. 1326 The 25 hour mean is available for daily runs by summing up the 25 hourly instantananeous hourly values from 1327 midnight at the start of the day to midight at the day end. 1328 This diagnostic is actived with the logical $ln\_dia25h$ 1329 1330 %------------------------------------------nam_dia25h------------------------------------------------------ 1331 \namdisplay{nam_dia25h} 1332 %---------------------------------------------------------------------------------------------------------- 1333 1334 \section{Top Middle and Bed hourly output } 1335 1336 A module is available to output the surface (top), mid water and bed diagnostics of a set of standard variables. 1337 This can be a useful diagnostic when hourly or sub-hourly output is required in high resolution tidal outputs. 1338 The tidal signal is retained but the overall data usage is cut to just three vertical levels. Also the bottom level 1339 is calculated for each cell. 1340 This diagnostic is actived with the logical $ln\_diatmb$ 1341 1342 %------------------------------------------nam_diatmb----------------------------------------------------- 1343 \namdisplay{nam_diatmb} 1344 %---------------------------------------------------------------------------------------------------------- 1345 1346 % ------------------------------------------------------------------------------------------------------------- 1321 1347 % Sections transports 1322 1348 % ------------------------------------------------------------------------------------------------------------- … … 1522 1548 \\ 1523 1549 1550 \section{Courant numbers} 1551 Courant numbers provide a theoretical indication of the model's numerical stability. The advective Courant numbers can be calculated according to 1552 \begin{equation} 1553 \label{eq:CFL} 1554 C_u = |u|\frac{\rdt}{e_{1u}}, \quad C_v = |v|\frac{\rdt}{e_{2v}}, \quad C_w = |w|\frac{\rdt}{e_{3w}} 1555 \end{equation} 1556 in the zonal, meridional and vertical directions respectively. The vertical component is included although it is not strictly valid as the vertical velocity is calculated from the continuity equation rather than as a prognostic variable. Physically this represents the rate at which information is propogated across a grid cell. Values greater than 1 indicate that information is propagated across more than one grid cell in a single time step. 1557 1558 The variables can be activated by setting the \np{nn\_diacfl} namelist parameter to 1 in the \ngn{namctl} namelist. The diagnostics will be written out to an ascii file named cfl\_diagnostics.ascii. In this file the maximum value of $C_u$, $C_v$, and $C_w$ are printed at each timestep along with the coordinates of where the maximum value occurs. At the end of the model run the maximum value of $C_u$, $C_v$, and $C_w$ for the whole model run is printed along with the coordinates of each. The maximum values from the run are also copied to the ocean.output file. 1524 1559 1525 1560
Note: See TracChangeset
for help on using the changeset viewer.