- Timestamp:
- 2016-02-01T03:35:04+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/nemo_v3_6_STABLE/DOC/TexFiles/Chapters/Chap_DIA.tex
r5515 r6275 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 % ================================================================ … … 48 48 49 49 50 Since version 3.2, iomput is the NEMO output interface of choice. It has been designed to be simple to use, flexible and efficient. The two main purposes of iomput are: 50 Since version 3.2, iomput is the NEMO output interface of choice. 51 It has been designed to be simple to use, flexible and efficient. 52 The two main purposes of iomput are: 51 53 \begin{enumerate} 52 54 \item The complete and flexible control of the output files through external XML files adapted by the user from standard templates. … … 1116 1118 % ------------------------------------------------------------------------------------------------------------- 1117 1119 \section[Tracer/Dynamics Trends (TRD)] 1118 {Tracer/Dynamics Trends (\key{trdtra}, \key{trddyn}, \\ 1119 \key{trddvor}, \key{trdmld})} 1120 {Tracer/Dynamics Trends (\ngn{namtrd})} 1120 1121 \label{DIA_trd} 1121 1122 … … 1124 1125 %------------------------------------------------------------------------------------------------------------- 1125 1126 1126 When \key{trddyn} and/or \key{trddyn} CPP variables are defined, each 1127 trend of the dynamics and/or temperature and salinity time evolution equations 1128 is stored in three-dimensional arrays just after their computation ($i.e.$ at the end 1129 of each $dyn\cdots.F90$ and/or $tra\cdots.F90$ routines). Options are defined by 1130 \ngn{namtrd} namelist variables. These trends are then 1131 used in \mdl{trdmod} (see TRD directory) every \textit{nn\_trd } time-steps. 1132 1133 What is done depends on the CPP keys defined: 1127 Each trend of the dynamics and/or temperature and salinity time evolution equations 1128 can be send to \mdl{trddyn} and/or \mdl{trdtra} modules (see TRD directory) just after their computation 1129 ($i.e.$ at the end of each $dyn\cdots.F90$ and/or $tra\cdots.F90$ routines). 1130 This capability is controlled by options offered in \ngn{namtrd} namelist. 1131 Note that the output are done with xIOS, and therefore the \key{IOM} is required. 1132 1133 What is done depends on the \ngn{namtrd} logical set to \textit{true}: 1134 1134 \begin{description} 1135 \item[\key{trddyn}, \key{trdtra}] : a check of the basin averaged properties of the momentum 1136 and/or tracer equations is performed ; 1137 \item[\key{trdvor}] : a vertical summation of the moment tendencies is performed, 1138 then the curl is computed to obtain the barotropic vorticity tendencies which are output ; 1139 \item[\key{trdmld}] : output of the tracer tendencies averaged vertically 1140 either over the mixed layer (\np{nn\_ctls}=0), 1141 or over a fixed number of model levels (\np{nn\_ctls}$>$1 provides the number of level), 1142 or over a spatially varying but temporally fixed number of levels (typically the base 1143 of the winter mixed layer) read in \ifile{ctlsurf\_idx} (\np{nn\_ctls}=1) ; 1135 \item[\np{ln\_glo\_trd}] : at each \np{nn\_trd} time-step a check of the basin averaged properties 1136 of the momentum and tracer equations is performed. This also includes a check of $T^2$, $S^2$, 1137 $\tfrac{1}{2} (u^2+v2)$, and potential energy time evolution equations properties ; 1138 \item[\np{ln\_dyn\_trd}] : each 3D trend of the evolution of the two momentum components is output ; 1139 \item[\np{ln\_dyn\_mxl}] : each 3D trend of the evolution of the two momentum components averaged 1140 over the mixed layer is output ; 1141 \item[\np{ln\_vor\_trd}] : a vertical summation of the moment tendencies is performed, 1142 then the curl is computed to obtain the barotropic vorticity tendencies which are output ; 1143 \item[\np{ln\_KE\_trd}] : each 3D trend of the Kinetic Energy equation is output ; 1144 \item[\np{ln\_tra\_trd}] : each 3D trend of the evolution of temperature and salinity is output ; 1145 \item[\np{ln\_tra\_mxl}] : each 2D trend of the evolution of temperature and salinity averaged 1146 over the mixed layer is output ; 1144 1147 \end{description} 1145 1146 The units in the output file can be changed using the \np{nn\_ucf} namelist parameter.1147 For example, in case of salinity tendency the units are given by PSU/s/\np{nn\_ucf}.1148 Setting \np{nn\_ucf}=86400 ($i.e.$ the number of second in a day) provides the tendencies in PSU/d.1149 1150 When \key{trdmld} is defined, two time averaging procedure are proposed.1151 Setting \np{ln\_trdmld\_instant} to \textit{true}, a simple time averaging is performed,1152 so that the resulting tendency is the contribution to the change of a quantity between1153 the two instantaneous values taken at the extremities of the time averaging period.1154 Setting \np{ln\_trdmld\_instant} to \textit{false}, a double time averaging is performed,1155 so that the resulting tendency is the contribution to the change of a quantity between1156 two \textit{time mean} values. The later option requires the use of an extra file, \ifile{restart\_mld}1157 (\np{ln\_trdmld\_restart}=true), to restart a run.1158 1159 1148 1160 1149 Note that the mixed layer tendency diagnostic can also be used on biogeochemical models 1161 1150 via the \key{trdtrc} and \key{trdmld\_trc} CPP keys. 1151 1152 \textbf{Note that} in the current version (v3.6), many changes has been introduced but not fully tested. 1153 In particular, options associated with \np{ln\_dyn\_mxl}, \np{ln\_vor\_trd}, and \np{ln\_tra\_mxl} 1154 are not working, and none of the option have been tested with variable volume ($i.e.$ \key{vvl} defined). 1155 1162 1156 1163 1157 % ------------------------------------------------------------------------------------------------------------- … … 1280 1274 \label{DIA_diag_harm} 1281 1275 1282 A module is available to compute the amplitude and phase for tidal waves.1283 This diagnostic is actived with \key{diaharm}.1284 1285 1276 %------------------------------------------namdia_harm---------------------------------------------------- 1286 1277 \namdisplay{namdia_harm} 1287 1278 %---------------------------------------------------------------------------------------------------------- 1288 1279 1289 Concerning the on-line Harmonic analysis, some parameters are available in namelist 1290 \ngn{namdia\_harm} : 1291 1292 - \texttt{nit000\_han} is the first time step used for harmonic analysis 1293 1294 - \texttt{nitend\_han} is the last time step used for harmonic analysis 1295 1296 - \texttt{nstep\_han} is the time step frequency for harmonic analysis 1297 1298 - \texttt{nb\_ana} is the number of harmonics to analyse 1299 1300 - \texttt{tname} is an array with names of tidal constituents to analyse 1301 1302 \texttt{nit000\_han} and \texttt{nitend\_han} must be between \texttt{nit000} and \texttt{nitend} of the simulation. 1280 A module is available to compute the amplitude and phase of tidal waves. 1281 This on-line Harmonic analysis is actived with \key{diaharm}. 1282 Some parameters are available in namelist \ngn{namdia\_harm} : 1283 1284 - \np{nit000\_han} is the first time step used for harmonic analysis 1285 1286 - \np{nitend\_han} is the last time step used for harmonic analysis 1287 1288 - \np{nstep\_han} is the time step frequency for harmonic analysis 1289 1290 - \np{nb\_ana} is the number of harmonics to analyse 1291 1292 - \np{tname} is an array with names of tidal constituents to analyse 1293 1294 \np{nit000\_han} and \np{nitend\_han} must be between \np{nit000} and \np{nitend} of the simulation. 1303 1295 The restart capability is not implemented. 1304 1296 1305 The Harmonic analysis solve th isequation:1297 The Harmonic analysis solve the following equation: 1306 1298 \begin{equation} 1307 1299 h_{i} - A_{0} + \sum^{nb\_ana}_{j=1}[A_{j}cos(\nu_{j}t_{j}-\phi_{j})] = e_{i} … … 1324 1316 \label{DIA_diag_dct} 1325 1317 1326 A module is available to compute the transport of volume, heat and salt through sections. This diagnostic1327 is actived with \key{diadct}.1318 A module is available to compute the transport of volume, heat and salt through sections. 1319 This diagnostic is actived with \key{diadct}. 1328 1320 1329 1321 Each section is defined by the coordinates of its 2 extremities. The pathways between them are contructed … … 1347 1339 %------------------------------------------------------------------------------------------------------------- 1348 1340 1349 \ texttt{nn\_dct}: frequency of instantaneous transports computing1350 1351 \ texttt{nn\_dctwri}: frequency of writing ( mean of instantaneous transports )1352 1353 \ texttt{nn\_debug}: debugging of the section1341 \np{nn\_dct}: frequency of instantaneous transports computing 1342 1343 \np{nn\_dctwri}: frequency of writing ( mean of instantaneous transports ) 1344 1345 \np{nn\_debug}: debugging of the section 1354 1346 1355 1347 \subsubsection{ To create a binary file containing the pathway of each section } … … 1482 1474 the \key{diahth} CPP key: 1483 1475 1484 - the mixed layer depth (based on a density criterion , \citet{de_Boyer_Montegut_al_JGR04}) (\mdl{diahth})1476 - the mixed layer depth (based on a density criterion \citep{de_Boyer_Montegut_al_JGR04}) (\mdl{diahth}) 1485 1477 1486 1478 - the turbocline depth (based on a turbulent mixing coefficient criterion) (\mdl{diahth})
Note: See TracChangeset
for help on using the changeset viewer.