Changeset 2376 for branches/nemo_v3_3_beta/DOC/TexFiles/Chapters
- Timestamp:
- 2010-11-11T18:01:29+01:00 (14 years ago)
- Location:
- branches/nemo_v3_3_beta/DOC/TexFiles/Chapters
- Files:
-
- 1 added
- 15 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/nemo_v3_3_beta/DOC/TexFiles/Chapters/Annex_E.tex
r2282 r2376 300 300 \begin{center} 301 301 \includegraphics[width=0.70\textwidth]{./TexFiles/Figures/Fig_ISO_triad.pdf} 302 \caption{Triads used in the Griffies's like iso-neutral diffision scheme for 302 \caption{ \label{Fig_ISO_triad} 303 Triads used in the Griffies's like iso-neutral diffision scheme for 303 304 $u$-component (upper panel) and $w$-component (lower panel).} 304 305 \end{center} -
branches/nemo_v3_3_beta/DOC/TexFiles/Chapters/Annex_ISO.tex
r2285 r2376 100 100 \begin{figure}[h] \begin{center} 101 101 \includegraphics[width=0.90\textwidth]{./TexFiles/Figures/Fig_triad_fluxes} 102 \caption{(a) Arrangement of triads $S_i$ and tracer gradients to 103 give lateral tracer flux from box $i,k$ to $i+1,k$ (b) Triads 104 $S'_i$ and tracer gradients to give vertical tracer flux from 105 box $i,k$ to $i,k+1$.} 102 \caption{ \label{Fig_ISO_triad} 103 (a) Arrangement of triads $S_i$ and tracer gradients to 104 give lateral tracer flux from box $i,k$ to $i+1,k$ 105 (b) Triads $S'_i$ and tracer gradients to give vertical tracer flux from 106 box $i,k$ to $i,k+1$.} 106 107 \label{Fig_triad} 107 108 \end{center} \end{figure} … … 168 169 \begin{figure}[h] \begin{center} 169 170 \includegraphics[width=0.60\textwidth]{./TexFiles/Figures/Fig_qcells} 170 \caption{Triad notation for quarter cells.T-cells are inside 171 \caption{ \label{Fig_ISO_triad_notation} 172 Triad notation for quarter cells.T-cells are inside 171 173 boxes, while the $i+\half,k$ u-cell is shaded in green and the 172 174 $i,k+\half$ w-cell is shaded in pink.} -
branches/nemo_v3_3_beta/DOC/TexFiles/Chapters/Chap_DOM.tex
r2349 r2376 1 2 1 % ================================================================ 3 2 % Chapter 2 Ñ Space and Time Domain (DOM) … … 40 39 41 40 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 42 \begin{figure}[!tb] \label{Fig_cell}\begin{center}41 \begin{figure}[!tb] \begin{center} 43 42 \includegraphics[width=0.90\textwidth]{./TexFiles/Figures/Fig_cell.pdf} 44 \caption{Arrangement of variables. $t$ indicates scalar points where temperature, 43 \caption{ \label{Fig_cell} 44 Arrangement of variables. $t$ indicates scalar points where temperature, 45 45 salinity, density, pressure and horizontal divergence are defined. ($u$,$v$,$w$) 46 46 indicates vector points, and $f$ indicates vorticity points where both relative and … … 80 80 as the sum of the relevant scale factors (see \eqref{DOM_bar}) in the next section). 81 81 82 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 82 83 \begin{table}[!tb] \label{Tab_cell} 83 84 \begin{center} \begin{tabular}{|p{46pt}|p{56pt}|p{56pt}|p{56pt}|} … … 92 93 fw & $i+1/2$ & $j+1/2$ & $k+1/2$ \\ \hline 93 94 \end{tabular} 94 \caption{Location of grid-points as a function of integer or integer and a half value 95 of the column, line or level. This indexing is only used for the writing of the semi- 96 discrete equation. In the code, the indexing uses integer values only and has a 97 reverse direction in the vertical (see \S\ref{DOM_Num_Index})} 95 \caption{ \label{Tab_cell} 96 Location of grid-points as a function of integer or integer and a half value of the column, 97 line or level. This indexing is only used for the writing of the semi-discrete equation. 98 In the code, the indexing uses integer values only and has a reverse direction 99 in the vertical (see \S\ref{DOM_Num_Index})} 98 100 \end{center} 99 101 \end{table} 102 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 100 103 101 104 % ------------------------------------------------------------------------------------------------------------- … … 206 209 207 210 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 208 \begin{figure}[!tb] \label{Fig_index_hor}\begin{center}211 \begin{figure}[!tb] \begin{center} 209 212 \includegraphics[width=0.90\textwidth]{./TexFiles/Figures/Fig_index_hor.pdf} 210 \caption{ Horizontal integer indexing used in the \textsc{Fortran} code. The dashed211 area indicates the cell in which variables contained in arrays have the same212 $i$- and $j$-indices}213 \caption{ \label{Fig_index_hor} 214 Horizontal integer indexing used in the \textsc{Fortran} code. The dashed area indicates 215 the cell in which variables contained in arrays have the same $i$- and $j$-indices} 213 216 \end{center} \end{figure} 214 217 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 256 259 257 260 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 258 \begin{figure}[!pt] \label{Fig_index_vert}\begin{center}261 \begin{figure}[!pt] \begin{center} 259 262 \includegraphics[width=.90\textwidth]{./TexFiles/Figures/Fig_index_vert.pdf} 260 \caption{Vertical integer indexing used in the \textsc{Fortran } code. Note that 263 \caption{ \label{Fig_index_vert} 264 Vertical integer indexing used in the \textsc{Fortran } code. Note that 261 265 the $k$-axis is orientated downward. The dashed area indicates the cell in 262 266 which variables contained in arrays have the same $k$-index.} … … 364 368 Fig.~\ref{Fig_zgr_e3}. 365 369 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 366 \begin{figure}[!t] \label{Fig_zgr_e3}\begin{center}370 \begin{figure}[!t] \begin{center} 367 371 \includegraphics[width=0.90\textwidth]{./TexFiles/Figures/Fig_zgr_e3.pdf} 368 \caption{Comparison of (a) traditional definitions of grid-point position and grid-size 369 in the vertical, and (b) analytically derived grid-point position and scale factors. For 370 both grids here, the same $w$-point depth has been chosen but in (a) the 372 \caption{ \label{Fig_zgr_e3} 373 Comparison of (a) traditional definitions of grid-point position and grid-size in the vertical, 374 and (b) analytically derived grid-point position and scale factors. 375 For both grids here, the same $w$-point depth has been chosen but in (a) the 371 376 $t$-points are set half way between $w$-points while in (b) they are defined from 372 377 an analytical function: $z(k)=5\,(i-1/2)^3 - 45\,(i-1/2)^2 + 140\,(i-1/2) - 150$. … … 471 476 472 477 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 473 \begin{figure}[!tb] \label{Fig_z_zps_s_sps}\begin{center}478 \begin{figure}[!tb] \begin{center} 474 479 \includegraphics[width=1.0\textwidth]{./TexFiles/Figures/Fig_z_zps_s_sps.pdf} 475 \caption{The ocean bottom as seen by the model: 480 \caption{ \label{Fig_z_zps_s_sps} 481 The ocean bottom as seen by the model: 476 482 (a) $z$-coordinate with full step, 477 483 (b) $z$-coordinate with partial step, … … 575 581 576 582 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 577 \begin{figure}[!tb] \label{Fig_zgr}\begin{center}583 \begin{figure}[!tb] \begin{center} 578 584 \includegraphics[width=0.90\textwidth]{./TexFiles/Figures/Fig_zgr.pdf} 579 \caption{Default vertical mesh for ORCA2: 30 ocean levels (L30). Vertical level 580 functions for (a) T-point depth and (b) the associated scale factor as computed 585 \caption{ \label{Fig_zgr} 586 Default vertical mesh for ORCA2: 30 ocean levels (L30). Vertical level functions for 587 (a) T-point depth and (b) the associated scale factor as computed 581 588 from \eqref{DOM_zgr_ana} using \eqref{DOM_zgr_coef} in $z$-coordinate.} 582 589 \end{center} \end{figure} … … 651 658 652 659 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 653 \begin{table} \label{Tab_orca_zgr} 654 \begin{center} \begin{tabular}{c||r|r|r|r} 660 \begin{table} \begin{center} \begin{tabular}{c||r|r|r|r} 655 661 \hline 656 662 \textbf{LEVEL}& \textbf{gdept}& \textbf{gdepw}& \textbf{e3t }& \textbf{e3w } \\ \hline … … 687 693 31 & \textbf{5250.23}& 5000.00 & \textbf{500.56} & 500.33 \\ \hline 688 694 \end{tabular} \end{center} 689 \caption{ Default vertical mesh in $z$-coordinate for 30 layers ORCA2 configuration690 as computed from \eqref{DOM_zgr_ana} using the coefficients given in691 \eqref{DOM_zgr_coef}}695 \caption{ \label{Tab_orca_zgr} 696 Default vertical mesh in $z$-coordinate for 30 layers ORCA2 configuration as computed 697 from \eqref{DOM_zgr_ana} using the coefficients given in \eqref{DOM_zgr_coef}} 692 698 \end{table} 693 699 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 766 772 767 773 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 768 \begin{figure}[!tb] \label{Fig_sco_function}\begin{center}774 \begin{figure}[!tb] \begin{center} 769 775 \includegraphics[width=1.0\textwidth]{./TexFiles/Figures/Fig_sco_function.pdf} 770 \caption{Examples of the stretching function applied to a sea mont; from left to right: 776 \caption{ \label{Fig_sco_function} 777 Examples of the stretching function applied to a sea mont; from left to right: 771 778 surface, surface and bottom, and bottom intensified resolutions} 772 779 \end{center} \end{figure} -
branches/nemo_v3_3_beta/DOC/TexFiles/Chapters/Chap_DYN.tex
r2349 r2376 291 291 292 292 %>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 293 \begin{figure}[!ht] \label{Fig_DYN_een_triad} 294 \begin{center} 293 \begin{figure}[!ht] \begin{center} 295 294 \includegraphics[width=0.70\textwidth]{./TexFiles/Figures/Fig_DYN_een_triad.pdf} 296 \caption{Triads used in the energy and enstrophy conserving scheme (een) for 295 \caption{ \label{Fig_DYN_een_triad} 296 Triads used in the energy and enstrophy conserving scheme (een) for 297 297 $u$-component (upper panel) and $v$-component (lower panel).} 298 \end{center} 299 \end{figure} 298 \end{center} \end{figure} 300 299 %>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 301 300 … … 769 768 770 769 %> > > > > > > > > > > > > > > > > > > > > > > > > > > > 771 \begin{figure}[!t] \label{Fig_DYN_dynspg_ts} 772 \begin{center} 770 \begin{figure}[!t] \begin{center} 773 771 \includegraphics[width=0.90\textwidth]{./TexFiles/Figures/Fig_DYN_dynspg_ts.pdf} 774 \caption{Schematic of the split-explicit time stepping scheme for the external 772 \caption{ \label{Fig_DYN_dynspg_ts} 773 Schematic of the split-explicit time stepping scheme for the external 775 774 and internal modes. Time increases to the right. 776 775 Internal mode time steps (which are also the model time steps) are denoted … … 790 789 velocity. The model time stepping scheme can then be achieved by a baroclinic 791 790 leap-frog time step that carries the surface height from $t-\rdt$ to $t+\rdt$. } 792 \end{center} 793 \end{figure} 791 \end{center} \end{figure} 794 792 %> > > > > > > > > > > > > > > > > > > > > > > > > > > > 795 793 -
branches/nemo_v3_3_beta/DOC/TexFiles/Chapters/Chap_LBC.tex
r2349 r2376 51 51 52 52 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 53 \begin{figure}[!t] \label{Fig_LBC_uv}\begin{center}53 \begin{figure}[!t] \begin{center} 54 54 \includegraphics[width=0.90\textwidth]{./TexFiles/Figures/Fig_LBC_uv.pdf} 55 \caption {Lateral boundary (thick line) at T-level. The velocity normal to the56 55 \caption{ \label{Fig_LBC_uv} 56 Lateral boundary (thick line) at T-level. The velocity normal to the boundary is set to zero.} 57 57 \end{center} \end{figure} 58 58 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 74 74 75 75 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 76 \begin{figure}[!p] \ label{Fig_LBC_shlat} \begin{center}76 \begin{figure}[!p] \begin{center} 77 77 \includegraphics[width=0.90\textwidth]{./TexFiles/Figures/Fig_LBC_shlat.pdf} 78 \caption {lateral boundary condition (a) free-slip ($rn\_shlat=0$) ; (b) no-slip ($rn\_shlat=2$) 78 \caption{ \label{Fig_LBC_shlat} 79 lateral boundary condition (a) free-slip ($rn\_shlat=0$) ; (b) no-slip ($rn\_shlat=2$) 79 80 ; (c) "partial" free-slip ($0<rn\_shlat<2$) and (d) "strong" no-slip ($2<rn\_shlat$). 80 81 Implied "ghost" velocity inside land area is display in grey. } 81 \end{center} \end{figure}82 \end{center} \end{figure} 82 83 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 83 84 … … 192 193 193 194 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 194 \begin{figure}[!t] \label{Fig_LBC_jperio}\begin{center}195 \begin{figure}[!t] \begin{center} 195 196 \includegraphics[width=1.0\textwidth]{./TexFiles/Figures/Fig_LBC_jperio.pdf} 196 \caption {setting of (a) east-west cyclic (b) symmetric across the equator boundary conditions.} 197 \caption{ \label{Fig_LBC_jperio} 198 setting of (a) east-west cyclic (b) symmetric across the equator boundary conditions.} 197 199 \end{center} \end{figure} 198 200 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 209 211 210 212 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 211 \begin{figure}[!t] \label{Fig_North_Fold_T}\begin{center}213 \begin{figure}[!t] \begin{center} 212 214 \includegraphics[width=0.90\textwidth]{./TexFiles/Figures/Fig_North_Fold_T.pdf} 213 \caption {North fold boundary with a $T$-point pivot and cyclic east-west boundary condition 214 ($jperio=4$), as used in ORCA 2, 1/4, and 1/12. Pink shaded area corresponds to the inner 215 domain mask (see text). } 215 \caption{ \label{Fig_North_Fold_T} 216 North fold boundary with a $T$-point pivot and cyclic east-west boundary condition 217 ($jperio=4$), as used in ORCA 2, 1/4, and 1/12. Pink shaded area corresponds 218 to the inner domain mask (see text). } 216 219 \end{center} \end{figure} 217 220 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 274 277 275 278 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 276 \begin{figure}[!t] \label{Fig_mpp}\begin{center}279 \begin{figure}[!t] \begin{center} 277 280 \includegraphics[width=0.90\textwidth]{./TexFiles/Figures/Fig_mpp.pdf} 278 \caption {Positioning of a sub-domain when massively parallel processing is used. } 281 \caption{ \label{Fig_mpp} 282 Positioning of a sub-domain when massively parallel processing is used. } 279 283 \end{center} \end{figure} 280 284 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 363 367 364 368 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 365 \begin{figure}[!ht] \label{Fig_mppini2}\begin{center}369 \begin{figure}[!ht] \begin{center} 366 370 \includegraphics[width=0.90\textwidth]{./TexFiles/Figures/Fig_mppini2.pdf} 367 \caption {Example of Atlantic domain defined for the CLIPPER projet. Initial grid is 368 composed of 773 x 1236 horizontal points. (a) the domain is split onto 9 \time 20 369 subdomains (jpni=9, jpnj=20). 52 subdomains are land areas. (b) 52 subdomains 370 are eliminated (white rectangles) and the resulting number of processors really 371 used during the computation is jpnij=128.} 371 \caption { \label{Fig_mppini2} 372 Example of Atlantic domain defined for the CLIPPER projet. Initial grid is 373 composed of 773 x 1236 horizontal points. 374 (a) the domain is split onto 9 \time 20 subdomains (jpni=9, jpnj=20). 375 52 subdomains are land areas. 376 (b) 52 subdomains are eliminated (white rectangles) and the resulting number 377 of processors really used during the computation is jpnij=128.} 372 378 \end{center} \end{figure} 373 379 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 444 450 445 451 %--------------------------------------------------TABLE-------------------------------------------------- 446 447 \begin{table}[htbp] \label{Tab_obc_param} 448 \begin{center} 449 \begin{tabular}{|l|c|c|c|} 452 \begin{table}[htbp] \begin{center} \begin{tabular}{|l|c|c|c|} 450 453 \hline 451 454 Boundary and & Constant index & Starting index (d\'{e}but) & Ending index (fin) \\ … … 464 467 lp\_obc\_north & $j$-index of a $v$ point & $i$ of a $T$ point & $i$ of a $T$ point \\ 465 468 \hline 466 \end{tabular} 467 \ end{center}468 \caption{Names of different indices relating to the open boundaries. In the case469 \end{tabular} \end{center} 470 \caption{ \label{Tab_obc_param} 471 Names of different indices relating to the open boundaries. In the case 469 472 of a completely open ocean domain with four ocean boundaries, the parameters 470 473 take exactly the values indicated.} 471 474 \end{table} 475 %------------------------------------------------------------------------------------------------------------ 472 476 473 477 The open boundaries must be along coordinate lines. On the C-grid, the boundary … … 496 500 497 501 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 498 \begin{figure}[!t] \label{Fig_obc_north}\begin{center}502 \begin{figure}[!t] \begin{center} 499 503 \includegraphics[width=0.70\textwidth]{./TexFiles/Figures/Fig_obc_north.pdf} 500 \caption {Localization of the North open boundary points.}501 \end{center} 502 \end{ figure}504 \caption{ \label{Fig_obc_north} 505 Localization of the North open boundary points.} 506 \end{center} \end{figure} 503 507 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 504 508 … … 559 563 560 564 %--------------------------------------------------TABLE-------------------------------------------------- 561 562 \begin{table}[htbp] \label{Tab_obc_ind} 563 \begin{center} 564 \begin{tabular}{|l|c|c|c|c|c|} 565 \begin{table}[htbp] \begin{center} \begin{tabular}{|l|c|c|c|c|c|} 565 566 \hline 566 567 OBC & Variable & file name & Index & Start & end \\ … … 581 582 & V & obcnorth\_V.nc & $je$-2 & $ib$+1 & $ie-1$ \\ 582 583 \hline 583 \end{tabular} 584 \ end{center}585 \caption{Requirements for creating open boundary files from a global configuration,584 \end{tabular} \end{center} 585 \caption{ \label{Tab_obc_ind} 586 Requirements for creating open boundary files from a global configuration, 586 587 appropriate for the subdomain of indices $ib:ie$, $jb:je$. ``Index'' designates the 587 588 $i$ or $j$ index along which the $u$ of $v$ boundary point is situated in the global … … 590 591 $-F$ $-d\;y,je-2$ $-d\;x,ib+1,ie-1$ } 591 592 \end{table} 593 %----------------------------------------------------------------------------------------------------------- 592 594 593 595 It is assumed that the open boundary files contain the variables for the period of … … 878 880 879 881 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 880 \begin{figure}[!t] \label{Fig_LBC_bdy_geom}\begin{center}882 \begin{figure}[!t] \begin{center} 881 883 \includegraphics[width=1.0\textwidth]{./TexFiles/Figures/Fig_LBC_bdy_geom.pdf} 882 \caption {Example of geometry of unstructured open boundary} 884 \caption { \label{Fig_LBC_bdy_geom} 885 Example of geometry of unstructured open boundary} 883 886 \end{center} \end{figure} 884 887 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 923 926 924 927 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 925 \begin{figure}[!t] \label{Fig_LBC_nc_header}\begin{center}928 \begin{figure}[!t] \begin{center} 926 929 \includegraphics[width=1.0\textwidth]{./TexFiles/Figures/Fig_LBC_nc_header.pdf} 927 \caption {Example of header of netcdf input data file for BDY} 930 \caption { \label{Fig_LBC_nc_header} 931 Example of header of netcdf input data file for BDY} 928 932 \end{center} \end{figure} 929 933 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> -
branches/nemo_v3_3_beta/DOC/TexFiles/Chapters/Chap_LDF.tex
r2349 r2376 353 353 354 354 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 355 \begin{figure}[!ht] \label{Fig_LDF_ZDF1}\begin{center}355 \begin{figure}[!ht] \begin{center} 356 356 \includegraphics[width=0.70\textwidth]{./TexFiles/Figures/Fig_LDF_ZDF1.pdf} 357 \caption {averaging procedure for isopycnal slope computation.} 358 \end{center} \end{figure} 357 \caption { \label{Fig_LDF_ZDF1} 358 averaging procedure for isopycnal slope computation.} 359 \end{center} \end{figure} 359 360 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 360 361 … … 380 381 381 382 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 382 \begin{figure}[!ht] \label{Fig_eiv_slp}\begin{center}383 \begin{figure}[!ht] \begin{center} 383 384 \includegraphics[width=0.70\textwidth]{./TexFiles/Figures/Fig_eiv_slp.pdf} 384 \caption {Vertical profile of the slope used for lateral mixing in the mixed layer : 385 \caption { \label{Fig_eiv_slp} 386 Vertical profile of the slope used for lateral mixing in the mixed layer : 385 387 \textit{(a)} in the real ocean the slope is the iso-neutral slope in the ocean interior, 386 388 which has to be adjusted at the surface boundary (i.e. it must tend to zero at the -
branches/nemo_v3_3_beta/DOC/TexFiles/Chapters/Chap_MISC.tex
r2364 r2376 58 58 59 59 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 60 \begin{figure}[!tbp] \label{Fig_MISC_strait_hand}\begin{center}60 \begin{figure}[!tbp] \begin{center} 61 61 \includegraphics[width=0.80\textwidth]{./TexFiles/Figures/Fig_Gibraltar.pdf} 62 62 \includegraphics[width=0.80\textwidth]{./TexFiles/Figures/Fig_Gibraltar2.pdf} 63 \caption {Example of the Gibraltar strait defined in a $1\deg \times 1\deg$ mesh. 63 \caption{ \label{Fig_MISC_strait_hand} 64 Example of the Gibraltar strait defined in a $1\deg \times 1\deg$ mesh. 64 65 \textit{Top}: using partially open cells. The meridional scale factor at $v$-point 65 66 is reduced on both sides of the strait to account for the real width of the strait … … 134 135 135 136 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 136 \begin{figure}[!ht] \label{Fig_LBC_zoom}\begin{center}137 \begin{figure}[!ht] \begin{center} 137 138 \includegraphics[width=0.90\textwidth]{./TexFiles/Figures/Fig_LBC_zoom.pdf} 138 \caption {Position of a model domain compared to the data input domain when the zoom functionality is used.} 139 \caption{ \label{Fig_LBC_zoom} 140 Position of a model domain compared to the data input domain when the zoom functionality is used.} 139 141 \end{center} \end{figure} 140 142 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 141 143 142 143 % ================================================================144 % 1D model functionality145 % ================================================================146 \section{Water column model: 1D model (\key{c1d})}147 \label{MISC_1D}148 149 The 1D model option simulates a stand alone water column within the 3D \NEMO system.150 It can be applied to the ocean alone or to the ocean-ice system and can include passive tracers151 or a biogeochemical model. It is set up by defining the \key{c1d} CPP key.152 The 1D model is a very useful tool153 \textit{(a)} to learn about the physics and numerical treatment of vertical mixing processes ;154 \textit{(b)} to investigate suitable parameterisations of unresolved turbulence (wind steering,155 langmuir circulation, skin layers) ;156 \textit{(c)} to compare the behaviour of different vertical mixing schemes ;157 \textit{(d)} to perform sensitivity studies on the vertical diffusion at a particular point of an ocean domain ;158 \textit{(d)} to produce extra diagnostics, without the large memory requirement of the full 3D model.159 160 The methodology is based on the use of the zoom functionality over the smallest possible161 domain : a 3 x 3 domain centred on the grid point of interest (see \S\ref{MISC_zoom}),162 with some extra routines. There is no need to define a new mesh, bathymetry,163 initial state or forcing, since the 1D model will use those of the configuration it is a zoom of.164 The chosen grid point is set in par\_oce.F90 module by setting the jpizoom and jpjzoom165 parameters to the indices of the location of the chosen grid point.166 144 167 145 % ================================================================ … … 260 238 The "bit comparison" option has been introduced in order to be able to check that 261 239 mono-processor and multi-processor runs give exactly the same results. 240 %THIS is to be updated with the mpp_sum_glo introduced in v3.3 241 % nn_bit_cmp today only check that the nn_cla = 0 (no cross land advection) 262 242 263 243 $\bullet$ Benchmark (\np{nn\_bench}). This option defines a benchmark run based on 264 a GYRE configuration in which the resolution remains the same whatever the domain265 size. This allows a very large model domain to be used, just by changing the domain266 size (\jp{jpiglo}, \jp{jpjglo}) and without adjusting either the time-step or the physical267 parameterisations.244 a GYRE configuration (see \S\ref{CFG_gyre}) in which the resolution remains the same 245 whatever the domain size. This allows a very large model domain to be used, just by 246 changing the domain size (\jp{jpiglo}, \jp{jpjglo}) and without adjusting either the time-step 247 or the physical parameterisations. 268 248 269 249 … … 607 587 volume ratio of each processing region. 608 588 609 \begin{table} 610 \begin{tab ular}{lrrr}589 %------------------------------------------TABLE---------------------------------------------------- 590 \begin{table} \begin{tabular}{lrrr} 611 591 Filename & NetCDF3 & NetCDF4 & Reduction\\ 612 592 &filesize & filesize & \% \\ … … 638 618 ORCA2\_2d\_grid\_W\_0007.nc & 4416 & 1368 & 70\%\\ 639 619 \end{tabular} 640 \caption{ \label{Tab_NC4} Filesize comparison between NetCDF3 and NetCDF4641 with chunking and compression}620 \caption{ \label{Tab_NC4} 621 Filesize comparison between NetCDF3 and NetCDF4 with chunking and compression} 642 622 \end{table} 623 %---------------------------------------------------------------------------------------------------- 643 624 644 625 Since version 3.2, an I/O server has been added which provides more … … 758 739 %------------------------------------------------------------------------------------------------------------- 759 740 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 760 \begin{figure}[!t] \label{Fig_mask_subasins}\begin{center}741 \begin{figure}[!t] \begin{center} 761 742 \includegraphics[width=1.0\textwidth]{./TexFiles/Figures/Fig_mask_subasins.pdf} 762 \caption {Decomposition of the World Ocean (here ORCA2) into sub-basin used in to compute 743 \caption{ \label{Fig_mask_subasins} 744 Decomposition of the World Ocean (here ORCA2) into sub-basin used in to compute 763 745 the heat and salt transports as well as the meridional stream-function: Atlantic basin (red), 764 746 Pacific basin (green), Indian basin (bleue), Indo-Pacific basin (bleue+green). … … 931 913 932 914 % ================================================================ 933 % predefined configurations934 % ================================================================935 \section{predefined configurations}936 \label{MISC_config}937 938 There is several predefined ocean configuration which use is controlled by a specific CPP key.939 940 The key set the domain sizes (\jp{jpiglo}, \jp{jpjglo}, \jp{jpk}), the mesh and the bathymetry,941 and, in some cases, add to the model physics some specific treatments.942 943 % -------------------------------------------------------------------------------------------------------------944 % ORCA family configurations945 % -------------------------------------------------------------------------------------------------------------946 \subsection{ORCA family: global ocean with tripolar grid}947 \label{MISC_config_orca}948 949 The NEMO system is provided with four built-in ORCA configurations which differ in the950 horizontal resolution used:951 \begin{description}952 \item[\key{orca\_r4}] \jp{cp\_cfg}~=~orca ; \jp{jp\_cfg}~=~4953 \item[\key{orca\_r2}] \jp{cp\_cfg}~=~orca ; \jp{jp\_cfg}~=~2954 \item[\key{orca\_r1}] \jp{cp\_cfg}~=~orca ; \jp{jp\_cfg}~=~1955 \item[\key{orca\_r05}] \jp{cp\_cfg}~=~orca ; \jp{jp\_cfg}~=~05956 \item[\key{orca\_r025}] \jp{cp\_cfg}~=~orca ; \jp{jp\_cfg}~=~025957 \end{description}958 959 \subsubsection{ORCA mesh}960 961 %>>>>>>>>>>>>>>>>>>>>>>>>>>>>962 \begin{figure}[!t] \label{Fig_MISC_ORCA_msh} \begin{center}963 \includegraphics[width=0.98\textwidth]{./TexFiles/Figures/Fig_ORCA_NH_mesh.pdf}964 \caption {ORCA mesh conception. The departure from an isotropic Mercator grid start poleward of 20\deg N.965 The two "north pole" are the foci of a series of embedded ellipses (blue curves)966 which are determined analytically and form the i-lines of the ORCA mesh (pseudo latitudes).967 Then, following \citet{Madec_Imbard_CD96}, the normal to the series of ellipses (red curves) is computed968 which provide the j-lines of the mesh (pseudo longitudes).969 }970 \end{center} \end{figure}971 %>>>>>>>>>>>>>>>>>>>>>>>>>>>>972 973 974 %>>>>>>>>>>>>>>>>>>>>>>>>>>>>975 \begin{figure}[!tbp] \label{Fig_MISC_ORCA_e1e2} \begin{center}976 \includegraphics[width=1.0\textwidth]{./TexFiles/Figures/Fig_ORCA_NH_msh05_e1_e2.pdf}977 \includegraphics[width=0.80\textwidth]{./TexFiles/Figures/Fig_ORCA_aniso.pdf}978 \caption {\textit{Top}: Horizontal scale factors ($e_1$, $e_2$) and979 \textit{Bottom}: ratio of anisotropy ($e_1 / e_2$)980 for ORCA 0.5\deg ~mesh. South of 20\deg N a Mercator grid is used ($e_1 = e_2$)981 so that the anisotropy ratio is 1. Poleward of 20\deg N, the two "north pole"982 introduce a weak anisotropy over the ocean areas ($< 1.2$) except in vicinity of Victoria Island983 (Canadian Arctic Archipelago). }984 \end{center} \end{figure}985 %>>>>>>>>>>>>>>>>>>>>>>>>>>>>986 987 %--------------------------------------------------TABLE--------------------------------------------------988 \begin{table}[htbp] \label{Tab_ORCA}989 \begin{center}990 \begin{tabular}{ccccc}991 key & \jp{jp\_cfg} & \jp{jpiglo} & \jp{jpiglo} & \\992 \hline \hline993 \key{orca\_r4} & 4 & 92 & 76 & \\994 \key{orca\_r2} & 2 & 182 & 149 & \\995 %\key{orca\_r1} & 1 & 362 & 511 & \\996 \key{orca\_r05} & 05 & 722 & 261 & \\997 \key{orca\_r025} & 025 & 1442 & 1021 & \\998 %\key{orca\_r8} & 8 & 2882 & 2042 & \\999 %\key{orca\_r12} & 12 & 4322 & 3062 & \\1000 \hline1001 \hline1002 \end{tabular}1003 \caption {Set of predefined ORCA parameters. }1004 \end{center}1005 \end{table}1006 %--------------------------------------------------------------------------------------------------------------1007 1008 The tripolar grid used in ORCA configuration ....1009 1010 NB: the two north poles position has been chosen to minimise the anisotropy ratio in1011 the Gulf Stream and kuroshio areas, two highly turbulent regions.1012 1013 ORCA~2 : a $2\deg$ zonal resolution, and a meridional resolution varying from $0.5\deg$ at the1014 equator to $2\deg cos\phi$ south of $20\deg$S (Fig. 1). The grid features two points of convergence in the1015 Northern Hemisphere, both situated on continents. Minimum resolution in high latitudes is about1016 65~km in the Arctic and 50~km in the Antarctic. Local mesh refinements are applied to the1017 Mediterranean, Red, Black and Caspian Seas. None of them appears to be of particular1018 importance for the study of high latitude climate, but the fine resolution is needed in order to have1019 their local circulation and their role in the World Ocean's circulation considered correctly.1020 1021 1022 1023 ORCA2-LIM (global ocean sea-ice configuration \citep{Timmermann_al_OM05}.1024 The horizontal mesh is based on a $2\deg \times 2\deg$ Mercator grid ($i.e.$ same zonal and1025 meridional grid spacing) which has been modified poleward1026 of $20\deg$N in order to include two numerical inland poles \citep{Murray_JCP96}.1027 This modification is semi-analytical \citep{Madec_Imbard_CD96}1028 and based on a series of embedded ellipses. It insures that the mesh remains1029 close to isotropy and that the smallest grid cell is along Antarctica.1030 In order to refine the meridional resolution up to $0.5\deg$ at the equator,1031 additional local transformations were applied with in the Tropics.1032 Local mesh refinements are also applied to the Mediterranean, Red, Black1033 and Caspian Seas so that the resolution is $1\deg \time 1\deg$ there.1034 There are 31 levels in the vertical, with the highest resolution (10m)1035 in the upper 150m. The bottom topography and the coastlines are derived1036 from the global atlas of Smith and Sandwell (1997).1037 1038 \key{orca\_lev10} 10 time more vertical levels1039 1040 \key{agrif} : ORCA2-LIM plus an AGRIF zoom over the Agulhas current area1041 1042 \key{arctic}, \key{antarctic} (not used in ORCA\_R4)1043 1044 1045 We thus only provide a brief introduction in this chapter.1046 The global coupled ocean-ice configuration is very similar to that used as part of the climate1047 model developed at GFDL for the 4th IPCC assessment of climate change (Griffies et al., 2005;1048 Gnanadesikan et al., 2006).1049 The ORCA2-LIM configuration is also the basis for the \NEMO contribution to the1050 Coordinate Ocean-ice Reference Experiments (COREs) documented in \citet{Griffies_al_OM09}.1051 These experiments employ the boundary forcing from \citet{Large_Yeager_Rep04} (see \S\ref{SBC_blk_core}),1052 which was developed for the purpose of running global coupled ocean-ice simulations without an1053 interactive atmosphere. This \citet{Large_Yeager_Rep04} dataset is available through the GFDL web1054 site \footnote{http://nomads.gfdl.noaa.gov/nomads/forms/mom4/CORE.html}.1055 The "normal year" of \citet{Large_Yeager_Rep04} has been chosen of the \NEMO distribution1056 since release v3.3.1057 1058 % -------------------------------------------------------------------------------------------------------------1059 % GYRE family configuration1060 % -------------------------------------------------------------------------------------------------------------1061 \subsection{GYRE family: double gyre basin (\key{gyre})}1062 \label{MISC_config_gyre}1063 1064 The GYRE configuration \citep{Levy_al_OM10} have been built to simulated1065 the seasonal cycle of a double-gyre box model. It consist in an idealized domain1066 similar to that used in the studies of \citet{Drijfhout_JPO94} and \citet{Hazeleger_Drijfhout_JPO98,1067 Hazeleger_Drijfhout_JPO99, Hazeleger_Drijfhout_JGR00, Hazeleger_Drijfhout_JPO00},1068 over which an analytical seasonal forcing is applied. This allows to investigate the1069 spontaneous generation of a large number of interacting, transient mesoscale eddies1070 and their contribution to the large scale circulation.1071 1072 The domain geometry is a closed rectangular basin on the $\beta$-plane centred1073 at $\sim 30\deg$N and rotated by 45\deg, 3180~km long, 2120~km wide1074 and 4~km deep (Fig.~\ref{Fig_MISC_strait_hand}).1075 The domain is bounded by vertical walls and by a ßat bottom. The configuration is1076 meant to represent an idealized North Atlantic or North Pacific basin.1077 The circulation is forced by analytical profiles of wind and buoyancy ßuxes.1078 The applied forcings vary seasonally in a sinusoidal manner between winter1079 and summer extrema \citep{Levy_al_OM10}.1080 The wind stress is zonal and its curl changes sign at 22\deg N and 36\deg N.1081 It forces a subpolar gyre in the north, a subtropical gyre in the wider part of the domain1082 and a small recirculation gyre in the southern corner.1083 The net heat ßux takes the form of a restoring toward a zonal apparent air1084 temperature profile. A portion of the net heat ßux which comes from the solar radiation1085 is allowed to penetrate within the water column.1086 The fresh water ßux is also prescribed and varies zonally.1087 It is determined such as, at each time step, the basin-integrated ßux is zero.1088 The basin is initialised at rest with vertical profiles of temperature and salinity1089 uniformly applied to the whole domain.1090 1091 The GYRE configuration is set through the \key{gyre} CPP key. Its horizontal resolution1092 (and thus the size of the domain) is determined by setting \jp{jp\_cfg} in \hf{par\_GYRE} file: \\1093 \jp{jpiglo} $= 30 \times$ \jp{jp\_cfg} + 2 \\1094 \jp{jpjglo} $= 20 \times$ \jp{jp\_cfg} + 2 \\1095 Obviously, the namelist parameters have to be adjusted to the chosen resolution.1096 In the vertical, GYRE uses the default 30 ocean levels (\jp{jpk}=31) (Fig.~\ref{Fig_zgr}).1097 1098 The GYRE configuration is also used in benchmark test as it is very simple to increase1099 its resolution and as it does not requires any input file. For example, keeping a same model size1100 on each processor while increasing the number of processor used is very easy, even though the1101 physical integrity of the solution can be compromised.1102 1103 %>>>>>>>>>>>>>>>>>>>>>>>>>>>>1104 \begin{figure}[!t] \label{Fig_GYRE} \begin{center}1105 \includegraphics[width=1.0\textwidth]{./TexFiles/Figures/Fig_GYRE.pdf}1106 \caption {Snapshot of relative vorticity at the surface of the model domain1107 in GYRE R9, R27 and R54. From \citet{Levy_al_OM10}.}1108 \end{center} \end{figure}1109 %>>>>>>>>>>>>>>>>>>>>>>>>>>>>1110 1111 % -------------------------------------------------------------------------------------------------------------1112 % EEL family configuration1113 % -------------------------------------------------------------------------------------------------------------1114 \subsection{EEL family: periodic channel}1115 \label{MISC_config_EEL}1116 1117 \begin{description}1118 \item[\key{eel\_r2}]1119 \item[\key{eel\_r5}]1120 \item[\key{eel\_r6}]1121 \end{description}1122 1123 % -------------------------------------------------------------------------------------------------------------1124 % POMME configuration1125 % -------------------------------------------------------------------------------------------------------------1126 \subsection{POMME: mid-latitude sub-domain}1127 \label{MISC_config_POMME}1128 1129 1130 \key{pomme\_r025}1131 1132 1133 -
branches/nemo_v3_3_beta/DOC/TexFiles/Chapters/Chap_Model_Basics.tex
r2349 r2376 6 6 \label{PE} 7 7 \minitoc 8 9 8 10 9 \newpage … … 114 113 115 114 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 116 \begin{figure}[!ht] \label{Fig_ocean_bc}\begin{center}115 \begin{figure}[!ht] \begin{center} 117 116 \includegraphics[width=0.90\textwidth]{./TexFiles/Figures/Fig_I_ocean_bc.pdf} 118 \caption{The ocean is bounded by two surfaces, $z=-H(i,j)$ and $z=\eta(i,j,k,t)$, where $H$ 119 is the depth of the sea floor and $\eta$ the height of the sea surface. Both $H$ and $\eta $ 120 are referenced to $z=0$.} 117 \caption{ \label{Fig_ocean_bc} 118 The ocean is bounded by two surfaces, $z=-H(i,j)$ and $z=\eta(i,j,t)$, where $H$ 119 is the depth of the sea floor and $\eta$ the height of the sea surface. 120 Both $H$ and $\eta$ are referenced to $z=0$.} 121 121 \end{center} \end{figure} 122 122 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 167 167 168 168 169 \newpage170 $\ $\newline % force a new ligne169 %\newpage 170 %$\ $\newline % force a new ligne 171 171 172 172 % ================================================================ … … 371 371 372 372 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 373 \begin{figure}[!tb] \label{Fig_referential}\begin{center}373 \begin{figure}[!tb] \begin{center} 374 374 \includegraphics[width=0.60\textwidth]{./TexFiles/Figures/Fig_I_earth_referential.pdf} 375 \caption{the geographical coordinate system $(\lambda,\varphi,z)$ and the curvilinear 375 \caption{ \label{Fig_referential} 376 the geographical coordinate system $(\lambda,\varphi,z)$ and the curvilinear 376 377 coordinate system (\textbf{i},\textbf{j},\textbf{k}). } 377 378 \end{center} \end{figure} … … 703 704 \label{PE_gco} 704 705 705 %\gmcomment{706 706 The ocean domain presents a huge diversity of situation in the vertical. First the ocean surface is a time dependent surface (moving surface). Second the ocean floor depends on the geographical position, varying from more than 6,000 meters in abyssal trenches to zero at the coast. Last but not least, the ocean stratification exerts a strong barrier to vertical motions and mixing. 707 707 Therefore, in order to represent the ocean with respect to the first point a space and time dependent vertical coordinate that follows the variation of the sea surface height $e.g.$ an $z$*-coordinate; for the second point, a space variation to fit the change of bottom topography $e.g.$ a terrain-following or $\sigma$-coordinate; and for the third point, one will be tempted to use a space and time dependent coordinate that follows the isopycnal surfaces, $e.g.$ an isopycnic coordinate. … … 717 717 The coordinate is also sometime referenced as an adaptive coordinate \citep{Hofmeister_al_OM09}, since the coordinate system is adapted in the course of the simulation. Its most often used implementation is via an ALE algorithm, in which a pure lagrangian step is followed by regridding and remapping steps, the later step implicitly embedding the vertical advection \citep{Hirt_al_JCP74, Chassignet_al_JPO03, White_al_JCP09}. Here we follow the \citep{Kasahara_MWR74} strategy : a regridding step (an update of the vertical coordinate) followed by an eulerian step with an explicit computation of vertical advection relative to the moving s-surfaces. 718 718 719 A key point here is that the $s$-coordinate depends on $(i,j)$ ==> horizontal pressure gradient... 720 721 the generalized vertical coordinates used in ocean modelling are not orthogonal, which contrasts with many other applications in mathematical physics. Hence, it is useful to keep in mind the following properties that may seem odd on initial encounter. 722 723 the horizontal velocity in ocean models measures motions in the horizontal plane, perpendicular to the local gravitational field. That is, horizontal velocity is mathematically the same regardless the vertical coordinate, be it geopotential, isopycnal, pressure, or terrain following. The key motivation for maintaining the same horizontal velocity component is that the hydrostatic and geostrophic balances are dominant in the large-scale ocean. Use of an alternative quasi-horizontal velocity, for example one oriented parallel to the generalized surface, would lead to unacceptable numerical errors. Correspondingly, the vertical direction is anti-parallel to the gravitational force in all of the coordinate systems. We do not choose the alternative of a quasi-vertical direction oriented normal to the surface of a constant generalized vertical coordinate. 724 725 It is the method used to measure transport across the generalized vertical coordinate surfaces which differs between the vertical coordinate choices. That is, computation of the dia-surface velocity component represents the fundamental distinction between the various coordinates. In some models, such as geopotential, pressure, 726 and terrain following, this transport is typically diagnosed from volume or mass conservation. In other models, such as isopycnal layered models, this transport is prescribed based on assumptions about the physical processes producing a flux across the layer interfaces. 727 728 729 In this section we first establish the PE in the generalised vertical $s$-coordinate, then we discuss the particular cases available in \NEMO, namely $z$, $z$*, $s$, and $\tilde z$. 719 %\gmcomment{ 720 721 %A key point here is that the $s$-coordinate depends on $(i,j)$ ==> horizontal pressure gradient... 722 723 the generalized vertical coordinates used in ocean modelling are not orthogonal, 724 which contrasts with many other applications in mathematical physics. 725 Hence, it is useful to keep in mind the following properties that may seem 726 odd on initial encounter. 727 728 The horizontal velocity in ocean models measures motions in the horizontal plane, 729 perpendicular to the local gravitational field. That is, horizontal velocity is mathematically 730 the same regardless the vertical coordinate, be it geopotential, isopycnal, pressure, 731 or terrain following. The key motivation for maintaining the same horizontal velocity 732 component is that the hydrostatic and geostrophic balances are dominant in the large-scale ocean. 733 Use of an alternative quasi-horizontal velocity, for example one oriented parallel 734 to the generalized surface, would lead to unacceptable numerical errors. 735 Correspondingly, the vertical direction is anti-parallel to the gravitational force in all 736 of the coordinate systems. We do not choose the alternative of a quasi-vertical 737 direction oriented normal to the surface of a constant generalized vertical coordinate. 738 739 It is the method used to measure transport across the generalized vertical coordinate 740 surfaces which differs between the vertical coordinate choices. That is, computation 741 of the dia-surface velocity component represents the fundamental distinction between 742 the various coordinates. In some models, such as geopotential, pressure, and 743 terrain following, this transport is typically diagnosed from volume or mass conservation. 744 In other models, such as isopycnal layered models, this transport is prescribed based 745 on assumptions about the physical processes producing a flux across the layer interfaces. 746 747 748 In this section we first establish the PE in the generalised vertical $s$-coordinate, 749 then we discuss the particular cases available in \NEMO, namely $z$, $z$*, $s$, and $\tilde z$. 730 750 %} 731 751 … … 821 841 822 842 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 823 \begin{figure}[!b] \label{Fig_z_zstar}\begin{center}843 \begin{figure}[!b] \begin{center} 824 844 \includegraphics[width=1.0\textwidth]{./TexFiles/Figures/Fig_z_zstar.pdf} 825 \caption{(a) $z$-coordinate in linear free-surface case ; (b) $z-$coordinate in non-linear 826 free surface case (c) re-scaled height coordinate (become popular as the \textit{z*-}coordinate 845 \caption{ \label{Fig_z_zstar} 846 (a) $z$-coordinate in linear free-surface case ; 847 (b) $z-$coordinate in non-linear free surface case ; 848 (c) re-scaled height coordinate (become popular as the \textit{z*-}coordinate 827 849 \citep{Adcroft_Campin_OM04} ).} 828 850 \end{center} \end{figure} -
branches/nemo_v3_3_beta/DOC/TexFiles/Chapters/Chap_Model_Basics_zstar.tex
r996 r2376 117 117 118 118 %> > > > > > > > > > > > > > > > > > > > > > > > > > > > 119 \begin{figure}[!t] \label{Fig_DYN_dynspg_ts} 120 \begin{center} 119 \begin{figure}[!t] \begin{center} 121 120 \includegraphics[width=0.90\textwidth]{./Figures/Fig_DYN_dynspg_ts.pdf} 122 \caption{Schematic of the split-explicit time stepping scheme for the barotropic and baroclinic modes, after \citet{Griffies2004}. Time increases to the right. Baroclinic time steps are denoted by $t-\Delta t$, $t, t+\Delta t$, and $t+2\Delta t$. The curved line represents a leap-frog time step, and the smaller barotropic time steps $N \Delta t=2\Delta t$ are denoted by the zig-zag line. The vertically integrated forcing \textbf{M}(t) computed at baroclinic time step t represents the interaction between the barotropic and baroclinic motions. While keeping the total depth, tracer, and freshwater forcing fields fixed, a leap-frog integration carries the surface height and vertically integrated velocity from t to $t+2 \Delta t$ using N barotropic time steps of length $\Delta t$. Time averaging the barotropic fields over the N+1 time steps (endpoints included) centers the vertically integrated velocity at the baroclinic timestep $t+\Delta t$. A baroclinic leap-frog time step carries the surface height to $t+\Delta t$ using the convergence of the time averaged vertically integrated velocity taken from baroclinic time step t. } 121 \caption{ \label{Fig_DYN_dynspg_ts} 122 Schematic of the split-explicit time stepping scheme for the barotropic and baroclinic modes, 123 after \citet{Griffies2004}. Time increases to the right. Baroclinic time steps are denoted by 124 $t-\Delta t$, $t, t+\Delta t$, and $t+2\Delta t$. The curved line represents a leap-frog time step, 125 and the smaller barotropic time steps $N \Delta t=2\Delta t$ are denoted by the zig-zag line. 126 The vertically integrated forcing \textbf{M}(t) computed at baroclinic time step t represents 127 the interaction between the barotropic and baroclinic motions. While keeping the total depth, 128 tracer, and freshwater forcing fields fixed, a leap-frog integration carries the surface height 129 and vertically integrated velocity from t to $t+2 \Delta t$ using N barotropic time steps of length 130 $\Delta t$. Time averaging the barotropic fields over the N+1 time steps (endpoints included) 131 centers the vertically integrated velocity at the baroclinic timestep $t+\Delta t$. 132 A baroclinic leap-frog time step carries the surface height to $t+\Delta t$ using the convergence 133 of the time averaged vertically integrated velocity taken from baroclinic time step t. } 123 134 \end{center} 124 135 \end{figure} -
branches/nemo_v3_3_beta/DOC/TexFiles/Chapters/Chap_OBS.tex
r2349 r2376 678 678 \subsubsection{Geographical distribution of observations among processors} 679 679 680 \begin{figure} 681 \begin{ center}680 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 681 \begin{figure} \begin{center} 682 682 \includegraphics[width=10cm,height=12cm,angle=-90.]{./TexFiles/Figures/Fig_ASM_obsdist_local} 683 \end{center} 684 \caption{Example of the distribution of observations with the geographical 685 distribution of observational data.} 686 \label{fig:obslocal} 687 \end{figure} 683 \caption{ \label{fig:obslocal} 684 Example of the distribution of observations with the geographical distribution of observational data.} 685 \end{center} \end{figure} 686 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 688 687 689 688 This is the simplest option in which the observations are distributed according … … 706 705 \subsubsection{Round-robin distribution of observations among processors} 707 706 708 \begin{figure} 709 \begin{ center}707 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 708 \begin{figure} \begin{center} 710 709 \includegraphics[width=10cm,height=12cm,angle=-90.]{./TexFiles/Figures/Fig_ASM_obsdist_global} 711 \end{center} 712 \caption{Example of the distribution of observations with the round-robin 713 distribution of observational data.} 714 \label{fig:obsglobal} 715 \end{figure} 710 \caption{ \label{fig:obsglobal} 711 Example of the distribution of observations with the round-robin distribution of observational data.} 712 \end{center} \end{figure} 713 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 716 714 717 715 An alternative approach is to distribute the observations equally -
branches/nemo_v3_3_beta/DOC/TexFiles/Chapters/Chap_SBC.tex
r2366 r2376 29 29 (\np{ln\_core}~=~true) or CLIO (\np{ln\_clio}~=~true) bulk formulae) and a coupled 30 30 formulation (exchanges with a atmospheric model via the OASIS coupler) 31 (\np{ln\_cpl}~=~true). The optional atmospheric pressure can be used either 32 to force ocean and ice dynamics (\np{ln\_apr\_dyn}~=~true), or in the bulk 33 formulae computation (\np{ln\_apr\_dyn}~=~true) 34 \footnote{None of the two current bulk formulea (CLIO and CORE) uses the 35 atmospheric pressure field.}. 31 (\np{ln\_cpl}~=~true). When used, the atmospheric pressure forces both 32 ocean and ice dynamics (\np{ln\_apr\_dyn}~=~true) 33 \footnote{The surface pressure field could be use in bulk formulae, nevertheless 34 none of the current bulk formulea (CLIO and CORE) uses the it.}. 36 35 The frequency at which the six or seven fields have to be updated is the \np{nn\_fsbc} 37 36 namelist parameter. … … 39 38 need not be supplied on the model grid. Instead a file of coordinates and weights can 40 39 be supplied which maps the data from the supplied grid to the model points 41 (so called "Interpolation on the Fly" ).40 (so called "Interpolation on the Fly", see \S\ref{SBC_iof}). 42 41 In addition, the resulting fields can be further modified using several namelist options. 43 42 These options control the rotation of vector components supplied relative to an east-north … … 52 51 53 52 In this chapter, we first discuss where the surface boundary condition appears in the 54 model equations. Then we present the four ways of providing the surface boundary condition. 53 model equations. Then we present the four ways of providing the surface boundary condition, 54 followed by the description of the atmospheric pressure and the river runoff. 55 55 Next the scheme for interpolation on the fly is described. 56 56 Finally, the different options that further modify the fluxes applied to the ocean are discussed. … … 157 157 158 158 %-------------------------------------------------TABLE--------------------------------------------------- 159 \begin{table}[tb] \label{Tab_ssm} 160 \begin{center} 161 \begin{tabular}{|l|l|l|l|} 159 \begin{table}[tb] \begin{center} \begin{tabular}{|l|l|l|l|} 162 160 \hline 163 161 Variable description & Model variable & Units & point \\ \hline … … 167 165 Sea surface salinty & sss\_m & $psu$ & T \\ \hline 168 166 \end{tabular} 169 \caption{Ocean variables provided by the ocean to the surface module (SBC). 167 \caption{ \label{Tab_ssm} 168 Ocean variables provided by the ocean to the surface module (SBC). 170 169 The variable are averaged over nf{\_}sbc time step, $i.e.$ the frequency of 171 170 computation of surface fluxes.} 172 \end{center} 173 \end{table} 171 \end{center} \end{table} 174 172 %-------------------------------------------------------------------------------------------------------------- 175 173 176 177 178 %\colorbox{yellow}{Penser a} mettre dans le restant l'info nf{\_}sbc ET nf{\_}sbc*rdt de sorte de reinitialiser la moyenne si on change la frequence ou le pdt 174 %\colorbox{yellow}{Penser a} mettre dans le restant l'info nn{\_}fsbc ET nn{\_}fsbc*rdt de sorte de reinitialiser la moyenne si on change la frequence ou le pdt 179 175 180 176 … … 385 381 %------------------------------------------------------------------------------------------------------------- 386 382 387 The optional atmospheric pressure can be used either to force ocean and ice dynamics 388 (\np{ln\_apr\_dyn}~=~true), or in the bulk formulae computation (\np{ln\_apr\_dyn}~=~true). 389 The input atmospheric forcing is interpolated in time to the model time step, and optionally 390 in space when interpolation on-the-fly is used. When used to force the dynamics, it is further 391 transformed into an equivalent inverse barometer sea surface height, $\eta_{ib}$, using: 383 The optional atmospheric pressure can be used to force ocean and ice dynamics 384 (\np{ln\_apr\_dyn}~=~true, \textit{namsbc} namelist ). 385 The input atmospheric forcing defined via \np{sn\_apr} structure (\textit{namsbc\_apr} namelist) 386 can be interpolated in time to the model time step, and even in space when the 387 interpolation on-the-fly is used. When used to force the dynamics, the atmospheric 388 pressure is further transformed into an equivalent inverse barometer sea surface height, 389 $\eta_{ib}$, using: 392 390 \begin{equation} \label{SBC_ssh_ib} 393 391 \eta_{ib} = - \frac{1}{g\,\rho_o} \left( P_{atm} - P_o \right) … … 398 396 $\eta_{ib}$ is kept to zero at all time step. 399 397 400 Agradient of $\eta_{ib}$ is added to the RHS of the ocean momentum equation398 The gradient of $\eta_{ib}$ is added to the RHS of the ocean momentum equation 401 399 (see \mdl{dynspg} for the ocean). For sea-ice, the sea surface height, $\eta_m$, 402 400 which is provided to the sea ice model is set to $\eta - \eta_{ib}$ (see \mdl{sbcssr} module). 403 Furthermore, $\eta_{ib}$ can be set in the output. This simplifies the altirmetry data 404 and model comparison as inverse barometer sea surface height is usually removed 405 from thise date prior to their distribution. 401 $\eta_{ib}$ can be set in the output. This can simplify the altirmetry data and model comparison 402 as inverse barometer sea surface height is usually removed from these date prior to their distribution. 406 403 407 404 % ================================================================ … … 522 519 523 520 521 % ================================================================ 522 % Interpolation on the Fly 523 % ================================================================ 524 525 \section [Interpolation on the Fly] {Interpolation on the Fly} 526 \label{SBC_iof} 527 528 Interpolation on the Fly allows the user to supply input files required 529 for the surface forcing on grids other than the model grid. 530 To do this he or she must supply, in addition to the source data file, 531 a file of weights to be used to interpolate from the data grid to the model 532 grid. 533 The original development of this code used the SCRIP package (freely available 534 under a copyright agreement from http://climate.lanl.gov/Software/SCRIP). 535 In principle, any package can be used to generate the weights, but the 536 variables in the input weights file must have the same names and meanings as 537 assumed by the model. 538 Two methods are currently available: bilinear and bicubic interpolation. 539 540 \subsection{Bilinear Interpolation} 541 \label{SBC_iof_bilinear} 542 543 The input weights file in this case has two sets of variables: src01, src02, 544 src03, src04 and wgt01, wgt02, wgt03, wgt04. 545 The "src" variables correspond to the point in the input grid to which the weight 546 "wgt" is to be applied. Each src value is an integer corresponding to the index of a 547 point in the input grid when written as a one dimensional array. For example, for an input grid 548 of size 5x10, point (3,2) is referenced as point 8, since (2-1)*5+3=8. 549 There are four of each variable because bilinear interpolation uses the four points defining 550 the grid box containing the point to be interpolated. 551 All of these arrays are on the model grid, so that values src01(i,j) and 552 wgt01(i,j) are used to generate a value for point (i,j) in the model. 553 554 Symbolically, the algorithm used is: 555 556 \begin{equation} 557 f_{m}(i,j) = f_{m}(i,j) + \sum_{k=1}^{4} {wgt(k)f(idx(src(k)))} 558 \end{equation} 559 where function idx() transforms a one dimensional index src(k) into a two dimensional index, 560 and wgt(1) corresponds to variable "wgt01" for example. 561 562 \subsection{Bicubic Interpolation} 563 \label{SBC_iof_bicubic} 564 565 Again there are two sets of variables: "src" and "wgt". 566 But in this case there are 16 of each. 567 The symbolic algorithm used to calculate values on the model grid is now: 568 569 \begin{equation*} \begin{split} 570 f_{m}(i,j) = f_{m}(i,j) +& \sum_{k=1}^{4} {wgt(k)f(idx(src(k)))} 571 + \sum_{k=5}^{8} {wgt(k)\left.\frac{\partial f}{\partial i}\right| _{idx(src(k))} } \\ 572 +& \sum_{k=9}^{12} {wgt(k)\left.\frac{\partial f}{\partial j}\right| _{idx(src(k))} } 573 + \sum_{k=13}^{16} {wgt(k)\left.\frac{\partial ^2 f}{\partial i \partial j}\right| _{idx(src(k))} } 574 \end{split} 575 \end{equation*} 576 The gradients here are taken with respect to the horizontal indices and not distances since the spatial dependency has been absorbed into the weights. 577 578 \subsection{Implementation} 579 \label{SBC_iof_imp} 580 581 To activate this option, a non-empty string should be supplied in the weights filename column 582 of the relevant namelist; if this is left as an empty string no action is taken. 583 In the model, weights files are read in and stored in a structured type (WGT) in the fldread 584 module, as and when they are first required. 585 This initialisation procedure determines whether the input data grid should be treated 586 as cyclical or not by inspecting a global attribute stored in the weights input file. 587 This attribute must be called "ew\_wrap" and be of integer type. 588 If it is negative, the input non-model grid is assumed not to be cyclic. 589 If zero or greater, then the value represents the number of columns that overlap. 590 $E.g.$ if the input grid has columns at longitudes 0, 1, 2, .... , 359, then ew\_wrap should be set to 0; 591 if longitudes are 0.5, 2.5, .... , 358.5, 360.5, 362.5, ew\_wrap should be 2. 592 If the model does not find attribute ew\_wrap, then a value of -999 is assumed. 593 In this case the \rou{fld\_read} routine defaults ew\_wrap to value 0 and therefore the grid 594 is assumed to be cyclic with no overlapping columns. 595 (In fact this only matters when bicubic interpolation is required.) 596 Note that no testing is done to check the validity in the model, since there is no way 597 of knowing the name used for the longitude variable, 598 so it is up to the user to make sure his or her data is correctly represented. 599 600 Next the routine reads in the weights. 601 Bicubic interpolation is assumed if it finds a variable with name "src05", otherwise 602 bilinear interpolation is used. The WGT structure includes dynamic arrays both for 603 the storage of the weights (on the model grid), and when required, for reading in 604 the variable to be interpolated (on the input data grid). 605 The size of the input data array is determined by examining the values in the "src" 606 arrays to find the minimum and maximum i and j values required. 607 Since bicubic interpolation requires the calculation of gradients at each point on the grid, 608 the corresponding arrays are dimensioned with a halo of width one grid point all the way around. 609 When the array of points from the data file is adjacent to an edge of the data grid, 610 the halo is either a copy of the row/column next to it (non-cyclical case), or is a copy 611 of one from the first few columns on the opposite side of the grid (cyclical case). 612 613 \subsection{Limitations} 614 \label{SBC_iof_lim} 615 616 \begin{description} 617 \item 618 The case where input data grids are not logically rectangular has not been tested. 619 \item 620 This code is not guaranteed to produce positive definite answers from positive definite inputs. 621 \item 622 The cyclic condition is only applied on left and right columns, and not to top and bottom rows. 623 \item 624 The gradients across the ends of a cyclical grid assume that the grid spacing between the two columns involved are consistent with the weights used. 625 \item 626 Neither interpolation scheme is conservative. 627 (There is a conservative scheme available in SCRIP, but this has not been implemented.) 628 \end{description} 629 630 \subsection{Utilities} 631 \label{SBC_iof_util} 632 633 % to be completed 634 A set of utilities to create a weights file for a rectilinear input grid is available 635 (see the directory NEMOGCM/TOOLS/WEIGHTS). 636 637 % ================================================================ 638 % Miscellanea options 639 % ================================================================ 640 \section{Miscellaneous options} 641 \label{SBC_misc} 642 643 % ------------------------------------------------------------------------------------------------------------- 644 % Diurnal cycle 645 % ------------------------------------------------------------------------------------------------------------- 646 \subsection [Diurnal cycle (\textit{sbcdcy})] 647 {Diurnal cycle (\mdl{sbcdcy})} 648 \label{SBC_dcy} 649 %------------------------------------------namsbc_rnf---------------------------------------------------- 650 %\namdisplay{namsbc} 651 %------------------------------------------------------------------------------------------------------------- 652 524 653 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 525 \begin{figure}[!t] \label{Fig_SBC_diurnal}\begin{center}654 \begin{figure}[!t] \begin{center} 526 655 \includegraphics[width=0.8\textwidth]{./TexFiles/Figures/Fig_SBC_diurnal.pdf} 527 \caption{Example of recontruction of the diurnal cycle variation of short wave flux 656 \caption{ \label{Fig_SBC_diurnal} 657 Example of recontruction of the diurnal cycle variation of short wave flux 528 658 from daily mean values. The reconstructed diurnal cycle (black line) is chosen 529 659 as the mean value of the analytical cycle (blue line) over a time step, not … … 531 661 \end{center} \end{figure} 532 662 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 533 534 % ================================================================535 % Diurnal cycle536 % ================================================================537 \section [Diurnal cycle (\textit{sbcdcy})]538 {Diurnal cycle (\mdl{sbcdcy})}539 \label{SBC_dcy}540 %------------------------------------------namsbc_rnf----------------------------------------------------541 %\namdisplay{namsbc}542 %-------------------------------------------------------------------------------------------------------------543 663 544 664 \cite{Bernie_al_JC05} have shown that to capture 90$\%$ of the diurnal variability of … … 565 685 566 686 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 567 \begin{figure}[!t] \label{Fig_SBC_dcy}\begin{center}687 \begin{figure}[!t] \begin{center} 568 688 \includegraphics[width=0.7\textwidth]{./TexFiles/Figures/Fig_SBC_dcy.pdf} 569 \caption{Example of recontruction of the diurnal cycle variation of short wave flux 689 \caption{ \label{Fig_SBC_dcy} 690 Example of recontruction of the diurnal cycle variation of short wave flux 570 691 from daily mean values on an ORCA2 grid with a time sampling of 2~hours (from 1am to 11pm). 571 692 The display is on (i,j) plane. } … … 577 698 appear due to an inconsistency between the scale of the vertical resolution 578 699 and the forcing acting on that scale. 579 580 % ================================================================581 % Interpolation on the Fly582 % ================================================================583 584 \section [Interpolation on the Fly] {Interpolation on the Fly}585 \label{SBC_iof}586 587 Interpolation on the Fly allows the user to supply input files required588 for the surface forcing on grids other than the model grid.589 To do this he or she must supply, in addition to the source data file,590 a file of weights to be used to interpolate from the data grid to the model591 grid.592 The original development of this code used the SCRIP package (freely available593 under a copyright agreement from http://climate.lanl.gov/Software/SCRIP).594 In principle, any package can be used to generate the weights, but the595 variables in the input weights file must have the same names and meanings as596 assumed by the model.597 Two methods are currently available: bilinear and bicubic interpolation.598 599 \subsection{Bilinear Interpolation}600 \label{SBC_iof_bilinear}601 602 The input weights file in this case has two sets of variables: src01, src02,603 src03, src04 and wgt01, wgt02, wgt03, wgt04.604 The "src" variables correspond to the point in the input grid to which the weight605 "wgt" is to be applied. Each src value is an integer corresponding to the index of a606 point in the input grid when written as a one dimensional array. For example, for an input grid607 of size 5x10, point (3,2) is referenced as point 8, since (2-1)*5+3=8.608 There are four of each variable because bilinear interpolation uses the four points defining609 the grid box containing the point to be interpolated.610 All of these arrays are on the model grid, so that values src01(i,j) and611 wgt01(i,j) are used to generate a value for point (i,j) in the model.612 613 Symbolically, the algorithm used is:614 615 \begin{equation}616 f_{m}(i,j) = f_{m}(i,j) + \sum_{k=1}^{4} {wgt(k)f(idx(src(k)))}617 \end{equation}618 where function idx() transforms a one dimensional index src(k) into a two dimensional index,619 and wgt(1) corresponds to variable "wgt01" for example.620 621 \subsection{Bicubic Interpolation}622 \label{SBC_iof_bicubic}623 624 Again there are two sets of variables: "src" and "wgt".625 But in this case there are 16 of each.626 The symbolic algorithm used to calculate values on the model grid is now:627 628 \begin{equation*} \begin{split}629 f_{m}(i,j) = f_{m}(i,j) +& \sum_{k=1}^{4} {wgt(k)f(idx(src(k)))}630 + \sum_{k=5}^{8} {wgt(k)\left.\frac{\partial f}{\partial i}\right| _{idx(src(k))} } \\631 +& \sum_{k=9}^{12} {wgt(k)\left.\frac{\partial f}{\partial j}\right| _{idx(src(k))} }632 + \sum_{k=13}^{16} {wgt(k)\left.\frac{\partial ^2 f}{\partial i \partial j}\right| _{idx(src(k))} }633 \end{split}634 \end{equation*}635 The gradients here are taken with respect to the horizontal indices and not distances since the spatial dependency has been absorbed into the weights.636 637 \subsection{Implementation}638 \label{SBC_iof_imp}639 640 To activate this option, a non-empty string should be supplied in the weights filename column641 of the relevant namelist; if this is left as an empty string no action is taken.642 In the model, weights files are read in and stored in a structured type (WGT) in the fldread643 module, as and when they are first required.644 This initialisation procedure determines whether the input data grid should be treated645 as cyclical or not by inspecting a global attribute stored in the weights input file.646 This attribute must be called "ew\_wrap" and be of integer type.647 If it is negative, the input non-model grid is assumed not to be cyclic.648 If zero or greater, then the value represents the number of columns that overlap.649 $E.g.$ if the input grid has columns at longitudes 0, 1, 2, .... , 359, then ew\_wrap should be set to 0;650 if longitudes are 0.5, 2.5, .... , 358.5, 360.5, 362.5, ew\_wrap should be 2.651 If the model does not find attribute ew\_wrap, then a value of -999 is assumed.652 In this case the \rou{fld\_read} routine defaults ew\_wrap to value 0 and therefore the grid653 is assumed to be cyclic with no overlapping columns.654 (In fact this only matters when bicubic interpolation is required.)655 Note that no testing is done to check the validity in the model, since there is no way656 of knowing the name used for the longitude variable,657 so it is up to the user to make sure his or her data is correctly represented.658 659 Next the routine reads in the weights.660 Bicubic interpolation is assumed if it finds a variable with name "src05", otherwise661 bilinear interpolation is used. The WGT structure includes dynamic arrays both for662 the storage of the weights (on the model grid), and when required, for reading in663 the variable to be interpolated (on the input data grid).664 The size of the input data array is determined by examining the values in the "src"665 arrays to find the minimum and maximum i and j values required.666 Since bicubic interpolation requires the calculation of gradients at each point on the grid,667 the corresponding arrays are dimensioned with a halo of width one grid point all the way around.668 When the array of points from the data file is adjacent to an edge of the data grid,669 the halo is either a copy of the row/column next to it (non-cyclical case), or is a copy670 of one from the first few columns on the opposite side of the grid (cyclical case).671 672 \subsection{Limitations}673 \label{SBC_iof_lim}674 675 \begin{description}676 \item677 The case where input data grids are not logically rectangular has not been tested.678 \item679 This code is not guaranteed to produce positive definite answers from positive definite inputs.680 \item681 The cyclic condition is only applied on left and right columns, and not to top and bottom rows.682 \item683 The gradients across the ends of a cyclical grid assume that the grid spacing between the two columns involved are consistent with the weights used.684 \item685 Neither interpolation scheme is conservative.686 (There is a conservative scheme available in SCRIP, but this has not been implemented.)687 \end{description}688 689 \subsection{Utilities}690 \label{SBC_iof_util}691 692 % to be completed693 A set of utilities to create a weights file for a rectilinear input grid is available.694 695 % ================================================================696 % Miscellanea options697 % ================================================================698 \section{Miscellaneous options}699 \label{SBC_misc}700 700 701 701 % ------------------------------------------------------------------------------------------------------------- -
branches/nemo_v3_3_beta/DOC/TexFiles/Chapters/Chap_STP.tex
r2282 r2376 195 195 %\gmcomment{ 196 196 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 197 \begin{figure}[!t] \label{Fig_TimeStep_flowchart}\begin{center}197 \begin{figure}[!t] \begin{center} 198 198 \includegraphics[width=0.7\textwidth]{./TexFiles/Figures/Fig_TimeStepping_flowchart.pdf} 199 \caption{Sketch of the leapfrog time stepping sequence in \NEMO from \citet{Leclair_Madec_OM09}. 199 \caption{ \label{Fig_TimeStep_flowchart} 200 Sketch of the leapfrog time stepping sequence in \NEMO from \citet{Leclair_Madec_OM09}. 200 201 The use of a semi-implicit computation of the hydrostatic pressure gradient requires 201 202 the tracer equation to be stepped forward prior to the momentum equation. … … 286 287 287 288 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 288 \begin{figure}[!t] \label{Fig_MLF_forcing}\begin{center}289 \begin{figure}[!t] \begin{center} 289 290 \includegraphics[width=0.90\textwidth]{./TexFiles/Figures/Fig_MLF_forcing.pdf} 290 \caption{Illustration of forcing integration methods. 291 \caption{ \label{Fig_MLF_forcing} 292 Illustration of forcing integration methods. 291 293 (top) ''Traditional'' formulation : the forcing is defined at the same time as the variable 292 294 to which it is applied (integer value of the time step index) and it is applied over a $2\rdt$ period. -
branches/nemo_v3_3_beta/DOC/TexFiles/Chapters/Chap_TRA.tex
r2349 r2376 90 90 the continuity equation which is used to calculate the vertical velocity. 91 91 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 92 \begin{figure}[!t] \label{Fig_adv_scheme}\begin{center}92 \begin{figure}[!t] \begin{center} 93 93 \includegraphics[width=0.9\textwidth]{./TexFiles/Figures/Fig_adv_scheme.pdf} 94 \caption{Schematic representation of some ways used to evaluate the tracer value 94 \caption{ \label{Fig_adv_scheme} 95 Schematic representation of some ways used to evaluate the tracer value 95 96 at $u$-point and the amount of tracer exchanged between two neighbouring grid 96 97 points. Upsteam biased scheme (ups): the upstream value is used and the black … … 836 837 837 838 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 838 \begin{figure}[!t] \label{Fig_traqsr_irradiance}\begin{center}839 \begin{figure}[!t] \begin{center} 839 840 \includegraphics[width=1.0\textwidth]{./TexFiles/Figures/Fig_TRA_Irradiance.pdf} 840 \caption{Penetration profile of the downward solar irradiance 841 calculated by four models. Two waveband chlorophyll-independent formulation (blue), 842 a chlorophyll-dependent monochromatic formulation (green), 4 waveband RGB formulation (red), 841 \caption{ \label{Fig_traqsr_irradiance} 842 Penetration profile of the downward solar irradiance calculated by four models. 843 Two waveband chlorophyll-independent formulation (blue), a chlorophyll-dependent 844 monochromatic formulation (green), 4 waveband RGB formulation (red), 843 845 61 waveband Morel (1988) formulation (black) for a chlorophyll concentration of 844 846 (a) Chl=0.05 mg/m$^3$ and (b) Chl=0.5 mg/m$^3$. From \citet{Lengaigne_al_CD07}.} … … 856 858 %-------------------------------------------------------------------------------------------------------------- 857 859 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 858 \begin{figure}[!t] \label{Fig_geothermal}\begin{center}860 \begin{figure}[!t] \begin{center} 859 861 \includegraphics[width=1.0\textwidth]{./TexFiles/Figures/Fig_TRA_geoth.pdf} 860 \caption{Geothermal Heat flux (in $mW.m^{-2}$) used by \cite{Emile-Geay_Madec_OS09}. 862 \caption{ \label{Fig_geothermal} 863 Geothermal Heat flux (in $mW.m^{-2}$) used by \cite{Emile-Geay_Madec_OS09}. 861 864 It is inferred from the age of the sea floor and the formulae of \citet{Stein_Stein_Nat92}.} 862 865 \end{center} \end{figure} … … 963 966 964 967 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 965 \begin{figure}[!t] \label{Fig_bbl}\begin{center}968 \begin{figure}[!t] \begin{center} 966 969 \includegraphics[width=0.7\textwidth]{./TexFiles/Figures/Fig_BBL_adv.pdf} 967 \caption{Advective/diffusive Bottom Boundary Layer. The BBL parameterisation is 970 \caption{ \label{Fig_bbl} 971 Advective/diffusive Bottom Boundary Layer. The BBL parameterisation is 968 972 activated when $\rho^i_{kup}$ is larger than $\rho^{i+1}_{kdnw}$. 969 973 Red arrows indicate the additional overturning circulation due to the advective BBL. … … 1316 1320 1317 1321 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 1318 \begin{figure}[!p] \label{Fig_Partial_step_scheme}\begin{center}1322 \begin{figure}[!p] \begin{center} 1319 1323 \includegraphics[width=0.9\textwidth]{./TexFiles/Figures/Partial_step_scheme.pdf} 1320 \caption{ Discretisation of the horizontal difference and average of tracers in the $z$-partial step coordinate (\np{ln\_zps}=true) in the case $( e3w_k^{i+1} - e3w_k^i )>0$. A linear interpolation is used to estimate $\widetilde{T}_k^{i+1}$, the tracer value at the depth of the shallower tracer point of the two adjacent bottom $T$-points. The horizontal difference is then given by: $\delta _{i+1/2} T_k= \widetilde{T}_k^{\,i+1} -T_k^{\,i}$ and the average by: $\overline{T}_k^{\,i+1/2}= ( \widetilde{T}_k^{\,i+1/2} - T_k^{\,i} ) / 2$. } 1324 \caption{ \label{Fig_Partial_step_scheme} 1325 Discretisation of the horizontal difference and average of tracers in the $z$-partial 1326 step coordinate (\np{ln\_zps}=true) in the case $( e3w_k^{i+1} - e3w_k^i )>0$. 1327 A linear interpolation is used to estimate $\widetilde{T}_k^{i+1}$, the tracer value 1328 at the depth of the shallower tracer point of the two adjacent bottom $T$-points. 1329 The horizontal difference is then given by: $\delta _{i+1/2} T_k= \widetilde{T}_k^{\,i+1} -T_k^{\,i}$ 1330 and the average by: $\overline{T}_k^{\,i+1/2}= ( \widetilde{T}_k^{\,i+1/2} - T_k^{\,i} ) / 2$. } 1321 1331 \end{center} \end{figure} 1322 1332 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> -
branches/nemo_v3_3_beta/DOC/TexFiles/Chapters/Chap_ZDF.tex
r2349 r2376 154 154 The choice of $P_{rt}$ is controlled by the \np{nn\_pdl} namelist parameter. 155 155 156 At the sea surface, the value of $\bar{e}$ is prescribed from the wind 157 stress field as $\bar{e}_o = e_{bb} |\tau| / \rho_o$, with $e_{bb}$ the \np{rn\_ebb} 158 namelist parameter. The default value of $e_{bb}$ is 3.75. \citep{Gaspar1990}), 159 however a much larger value can be used when taking into account the 160 surface wave breaking (see below Eq. \eqref{ZDF_Esbc}). 161 The bottom value of TKE is assumed to be equal to the value of the level just above. 162 The time integration of the $\bar{e}$ equation may formally lead to negative values 163 because the numerical scheme does not ensure its positivity. To overcome this 164 problem, a cut-off in the minimum value of $\bar{e}$ is used (\np{rn\_emin} 165 namelist parameter). Following \citet{Gaspar1990}, the cut-off value is set 166 to $\sqrt{2}/2~10^{-6}~m^2.s^{-2}$. This allows the subsequent formulations 167 to match that of \citet{Gargett1984} for the diffusion in the thermocline and 168 deep ocean : $K_\rho = 10^{-3} / N$. 169 In addition, a cut-off is applied on $K_m$ and $K_\rho$ to avoid numerical 170 instabilities associated with too weak vertical diffusion. They must be 171 specified at least larger than the molecular values, and are set through 172 \np{rn\_avm0} and \np{rn\_avt0} (namzdf namelist, see \S\ref{ZDF_cst}). 173 174 \subsubsection{Turbulent length scale} 156 175 For computational efficiency, the original formulation of the turbulent length 157 176 scales proposed by \citet{Gaspar1990} has been simplified. Four formulations … … 187 206 mixing length scales as (and note that here we use numerical indexing): 188 207 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 189 \begin{figure}[!t] \ label{Fig_mixing_length} \begin{center}208 \begin{figure}[!t] \begin{center} 190 209 \includegraphics[width=1.00\textwidth]{./TexFiles/Figures/Fig_mixing_length.pdf} 191 \caption {Illustration of the mixing length computation. } 210 \caption{ \label{Fig_mixing_length} 211 Illustration of the mixing length computation. } 192 212 \end{center} 193 213 \end{figure} … … 204 224 $i.e.$ $l^{(k)} = \sqrt {2 {\bar e}^{(k)} / {N^2}^{(k)} }$. 205 225 206 In the \np{nn\_mxl} =2 case, the dissipation and mixing length scales take the same226 In the \np{nn\_mxl}~=~2 case, the dissipation and mixing length scales take the same 207 227 value: $ l_k= l_\epsilon = \min \left(\ l_{up} \;,\; l_{dwn}\ \right)$, while in the 208 \np{nn\_mxl} =2case, the dissipation and mixing turbulent length scales are give228 \np{nn\_mxl}~=~3 case, the dissipation and mixing turbulent length scales are give 209 229 as in \citet{Gaspar1990}: 210 230 \begin{equation} \label{Eq_tke_mxl_gaspar} … … 215 235 \end{equation} 216 236 217 At the sea surface the value of $\bar{e}$ is prescribed from the wind 218 stress field: $\bar{e}=rn\_ebb\;\left| \tau \right|$ (\np{rn\_ebb}=60 by default) 219 with a minimal threshold of \np{rn\_emin0}$=10^{-4}~m^2.s^{-2}$ (namelist 220 parameters). Its value at the bottom of the ocean is assumed to be 221 equal to the value of the level just above. The time integration of the 222 $\bar{e}$ equation may formally lead to negative values because the 223 numerical scheme does not ensure its positivity. To overcome this 224 problem, a cut-off in the minimum value of $\bar{e}$ is used (\np{rn\_emin} 225 namelist parameter). Following \citet{Gaspar1990}, the cut-off value is set 226 to $\sqrt{2}/2~10^{-6}~m^2.s^{-2}$. This allows the subsequent formulations 227 to match that of \citet{Gargett1984} for the diffusion in the thermocline and 228 deep ocean : $K_\rho = 10^{-3} / N$. 229 In addition, a cut-off is applied on $K_m$ and $K_\rho$ to avoid numerical 230 instabilities associated with too weak vertical diffusion. They must be 231 specified at least larger than the molecular values, and are set through 232 \np{rn\_avm0} and \np{rn\_avt0} (namzdf namelist, see \S\ref{ZDF_cst}). 233 234 % ------------------------------------------------------------------------------------------------------------- 235 % TKE Turbulent Closure Scheme : new organization to energetic considerations 237 At the ocean surface, a non zero length scale is set through the \np{rn\_lmin0} namelist 238 parameter. Usually the surface scale is given by $l_o = \kappa \,z_o$ 239 where $\kappa = 0.4$ is von Karman's constant and $z_o$ the roughness 240 parameter of the surface. Assuming $z_o=0.1$~m \citep{Craig_Banner_JPO94} 241 leads to a 0.04~m, the default value of \np{rn\_lsurf}. In the ocean interior 242 a minimum length scale is set to recover the molecular viscosity when $\bar{e}$ 243 reach its minimum value ($1.10^{-6}= C_k\, l_{min} \,\sqrt{\bar{e}_{min}}$ ). 244 245 246 \subsubsection{Surface wave breaking parameterization} 247 %-----------------------------------------------------------------------% 248 249 Following \citet{Mellor_Blumberg_JPO04}, the TKE turbulence closure model has been modified 250 to include the effect of surface wave breaking energetics. This results in a reduction of summertime 251 surface temperature when the mixed layer is relatively shallow. The \citet{Mellor_Blumberg_JPO04} 252 modifications acts on surface length scale and TKE values and air-sea drag coefficient. 253 The latter concerns the bulk formulea and is not discussed here. 254 255 Following \citet{Craig_Banner_JPO94}, the boundary condition on surface TKE value is : 256 \begin{equation} \label{ZDF_Esbc} 257 \bar{e}_o = \frac{1}{2}\,\left( 15.8\,\alpha_{CB} \right)^{2/3} \,\frac{|\tau|}{\rho_o} 258 \end{equation} 259 where $\alpha_{CB}$ is the \citet{Craig_Banner_JPO94} constant of proportionality 260 which depends on the ''wave age'', ranging from 57 for mature waves to 146 for 261 younger waves \citep{Mellor_Blumberg_JPO04}. 262 The boundary condition on the turbulent length scale follows the Charnock's relation: 263 \begin{equation} \label{ZDF_Lsbc} 264 l_o = \kappa \beta \,\frac{|\tau|}{g\,\rho_o} 265 \end{equation} 266 where $\kappa=0.40$ is the von Karman constant, and $\beta$ is the Charnock's constant. 267 \citet{Mellor_Blumberg_JPO04} suggest $\beta = 2.10^{5}$ the value chosen by \citet{Stacey_JPO99} 268 citing observation evidence, and $\alpha_{CB} = 100$ the Craig and Banner's value. 269 As the surface boundary condition on TKE is prescribed through $\bar{e}_o = e_{bb} |\tau| / \rho_o$, 270 with $e_{bb}$ the \np{rn\_ebb} namelist parameter, setting \np{rn\_ebb}~=~67.83 corresponds 271 to $\alpha_{CB} = 100$. further setting \np{ln\_lsurf} to true applies \eqref{ZDF_Lsbc} 272 as surface boundary condition on length scale, with $\beta$ hard coded to the Stacet's value. 273 Note that a minimal threshold of \np{rn\_emin0}$=10^{-4}~m^2.s^{-2}$ (namelist parameters) 274 is applied on surface $\bar{e}$ value. 275 276 277 \subsubsection{Langmuir cells} 278 %--------------------------------------% 279 Langmuir circulations (LC) can be described as ordered large-scale vertical motions 280 in the surface layer of the oceans. Although LC have nothing to do with convection, 281 the circulation pattern is rather similar to so-called convective rolls in the atmospheric 282 boundary layer. The detailed physics behind LC is described in, for example, 283 \citet{Craik_Leibovich_JFM76}. The prevailing explanation is that LC arise from 284 a nonlinear interaction between the Stokes drift and wind drift currents. 285 286 Here we introduced in the TKE turbulent closure the simple parameterization of 287 Langmuir circulations proposed by \citep{Axell_JGR02} for a $k-\epsilon$ turbulent closure. 288 The parameterization, tuned against large-eddy simulation, includes the whole effect 289 of LC in an extra source terms of TKE, $P_{LC}$. 290 The presence of $P_{LC}$ in \eqref{Eq_zdftke_e}, the TKE equation, is controlled 291 by setting \np{ln\_lc} to \textit{true} in the namtke namelist. 292 293 By making an analogy with the characteristic convective velocity scale 294 ($e.g.$, \citet{D'Alessio_al_JPO98}), $P_{LC}$ is assumed to be : 295 \begin{equation} 296 P_{LC}(z) = \frac{w_{LC}^3(z)}{H_{LC}} 297 \end{equation} 298 where $w_{LC}(z)$ is the vertical velocity profile of LC, and $H_{LC}$ is the LC depth. 299 With no information about the wave field, $w_{LC}$ is assumed to be proportional to 300 the Stokes drift $u_s = 0.377\,\,|\tau|^{1/2}$, where $|\tau|$ is the surface wind stress module 301 \footnote{Following \citet{Li_Garrett_JMR93}, the surface Stoke drift velocity 302 may be expressed as $u_s = 0.016 \,|U_{10m}|$. Assuming an air density of 303 $\rho_a=1.22 \,Kg/m^3$ and a drag coefficient of $1.5~10^{-3}$ give the expression 304 used of $u_s$ as a function of the module of surface stress}. 305 For the vertical variation, $w_{LC}$ is assumed to be zero at the surface as well as 306 at a finite depth $H_{LC}$ (which is often close to the mixed layer depth), and simply 307 varies as a sine function in between (a first-order profile for the Langmuir cell structures). 308 The resulting expression for $w_{LC}$ is : 309 \begin{equation} 310 w_{LC} = \begin{cases} 311 c_{LC} \,u_s \,\sin(- \pi\,z / H_{LC} ) & \text{if $-z \leq H_{LC}$} \\ 312 0 & \text{otherwise} 313 \end{cases} 314 \end{equation} 315 where $c_{LC} = 0.15$ has been chosen by \citep{Axell_JGR02} as a good compromise 316 to fit LES data. The chosen value yields maximum vertical velocities $w_{LC}$ of the order 317 of a few centimeters per second. The value of $c_{LC}$ is set through the \np{rn\_lc} 318 namelist parameter, having in mind that it should stay between 0.15 and 0.54 \citep{Axell_JGR02}. 319 320 The $H_{LC}$ is estimated in a similar way as the turbulent length scale of TKE equations: 321 $H_{LC}$ is depth to which a water parcel with kinetic energy due to Stoke drift 322 can reach on its own by converting its kinetic energy to potential energy, according to 323 \begin{equation} 324 - \int_{-H_{LC}}^0 { N^2\;z \;dz} = \frac{1}{2} u_s^2 325 \end{equation} 326 327 328 %\subsubsection{Mixing just below the mixed layer} 329 %---------------------------------------------------------------% 330 331 % add here a description of "penetration of TKE" and the associated namelist parameters 332 333 % ------------------------------------------------------------------------------------------------------------- 334 % TKE discretization considerations 236 335 % ------------------------------------------------------------------------------------------------------------- 237 336 \subsection{TKE discretization considerations (\key{zdftke})} … … 239 338 240 339 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 241 \begin{figure}[!t] \label{Fig_TKE_time_scheme}\begin{center}340 \begin{figure}[!t] \begin{center} 242 341 \includegraphics[width=1.00\textwidth]{./TexFiles/Figures/Fig_ZDF_TKE_time_scheme.pdf} 243 \caption {Illustration of the TKE time integration and its links to the momentum and tracer time integration. } 342 \caption{ \label{Fig_TKE_time_scheme} 343 Illustration of the TKE time integration and its links to the momentum and tracer time integration. } 244 344 \end{center} 245 345 \end{figure} … … 389 489 390 490 %--------------------------------------------------TABLE-------------------------------------------------- 391 \begin{table}[htbp] \label{Tab_GLS} 392 \begin{center} 491 \begin{table}[htbp] \begin{center} 393 492 %\begin{tabular}{cp{70pt}cp{70pt}cp{70pt}cp{70pt}cp{70pt}cp{70pt}c} 394 493 \begin{tabular}{ccccc} … … 408 507 \hline 409 508 \end{tabular} 410 \caption {Set of predefined GLS parameters, or equivalently predefined turbulence models available with \key{zdfgls} and controlled by the \np{nn\_clos} namelist parameter.} 411 \end{center} 412 \end{table} 509 \caption{ \label{Tab_GLS} 510 Set of predefined GLS parameters, or equivalently predefined turbulence models available 511 with \key{zdfgls} and controlled by the \np{nn\_clos} namelist parameter.} 512 \end{center} \end{table} 413 513 %-------------------------------------------------------------------------------------------------------------- 414 514 … … 417 517 value near physical boundaries (logarithmic boundary layer law). $C_{\mu}$ and $C_{\mu'}$ 418 518 are calculated from stability function proposed by \citet{Galperin_al_JAS88}, or by \citet{Kantha_Clayson_1994} 419 or one of the two functions suggested by \citet{Canuto_2001} (\np{nn\_stab\_func} = 0, 1, 2 or 3, resp.}). The value of $C_{0\mu}$ depends of the choice of the stability function. 519 or one of the two functions suggested by \citet{Canuto_2001} (\np{nn\_stab\_func} = 0, 1, 2 or 3, resp.}). 520 The value of $C_{0\mu}$ depends of the choice of the stability function. 420 521 421 522 The surface and bottom boundary condition on both $\bar{e}$ and $\psi$ can be calculated 422 523 thanks to Dirichlet or Neumann condition through \np{nn\_tkebc\_surf} and \np{nn\_tkebc\_bot}, resp. 423 The wave effect on the mixing could be also being considered \citep{Craig_Banner_1994}. 524 As for TKE closure , the wave effect on the mixing is considered when \np{ln\_crban}~=~true 525 \citep{Craig_Banner_JPO94, Mellor_Blumberg_JPO04}. The \np{rn\_crban} namelist parameter 526 is $\alpha_{CB}$ in \eqref{ZDF_Esbc} and \np{rn\_charn} provides the value of $\beta$ in \eqref{ZDF_Lsbc}. 424 527 425 528 The $\psi$ equation is known to fail in stably stratified flows, and for this reason … … 433 536 if \np{ln\_length\_lim}=true, and the $c_{lim}$ is set to the \np{rn\_clim\_galp} value. 434 537 538 The time and space discretization of the GLS equations follows the same energetic 539 consideration as for the TKE case described in \S\ref{ZDF_tke_ene} \citep{Burchard_OM02}. 540 Examples of performance of the 4 turbulent closure scheme can be found in \citet{Warner_al_OM05}. 541 435 542 % ------------------------------------------------------------------------------------------------------------- 436 543 % K Profile Parametrisation (KPP) … … 479 586 480 587 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 481 \begin{figure}[!htb] \label{Fig_npc}\begin{center}588 \begin{figure}[!htb] \begin{center} 482 589 \includegraphics[width=0.90\textwidth]{./TexFiles/Figures/Fig_npc.pdf} 483 \caption {Example of an unstable density profile treated by the non penetrative 590 \caption{ \label{Fig_npc} 591 Example of an unstable density profile treated by the non penetrative 484 592 convective adjustment algorithm. $1^{st}$ step: the initial profile is checked from 485 593 the surface to the bottom. It is found to be unstable between levels 3 and 4. … … 641 749 642 750 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 643 \begin{figure}[!t] \label{Fig_zdfddm}\begin{center}751 \begin{figure}[!t] \begin{center} 644 752 \includegraphics[width=0.99\textwidth]{./TexFiles/Figures/Fig_zdfddm.pdf} 645 \caption {From \citet{Merryfield1999} : (a) Diapycnal diffusivities $A_f^{vT}$ 753 \caption{ \label{Fig_zdfddm} 754 From \citet{Merryfield1999} : (a) Diapycnal diffusivities $A_f^{vT}$ 646 755 and $A_f^{vS}$ for temperature and salt in regions of salt fingering. Heavy 647 756 curves denote $A^{\ast v} = 10^{-3}~m^2.s^{-1}$ and thin curves … … 986 1095 987 1096 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 988 \begin{figure}[!t] \label{Fig_ZDF_M2_K1_tmx}\begin{center}1097 \begin{figure}[!t] \begin{center} 989 1098 \includegraphics[width=0.90\textwidth]{./TexFiles/Figures/Fig_ZDF_M2_K1_tmx.pdf} 990 \caption {(a) M2 and (b) K2 internal wave drag energy from \citet{Carrere_Lyard_GRL03} ($W/m^2$). }991 \end{center} 992 \end{ figure}1099 \caption{ \label{Fig_ZDF_M2_K1_tmx} 1100 (a) M2 and (b) K2 internal wave drag energy from \citet{Carrere_Lyard_GRL03} ($W/m^2$). } 1101 \end{center} \end{figure} 993 1102 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 994 1103 -
branches/nemo_v3_3_beta/DOC/TexFiles/Chapters/Introduction.tex
r2349 r2376 103 103 is a module related to the TRAcers equation, computing the Lateral DiFfussion. 104 104 The complete list of module names is presented in Appendix~\ref{Apdx_D}. 105 Furthermore, modules are organized in a few directories 106 that correspond to their category, as indicated by the first three letters of their name.105 Furthermore, modules are organized in a few directories that correspond to their category, 106 as indicated by the first three letters of their name (Tab.~\ref{Tab_chap}). 107 107 108 108 The manual mirrors the organization of the model. 109 109 After the presentation of the continuous equations (Chapter \ref{PE}), the following chapters 110 refer to specific terms of the equations each associated with a group of modules .110 refer to specific terms of the equations each associated with a group of modules (Tab.~\ref{Tab_chap}). 111 111 112 112 113 \begin{table}[htbp] \label{tab1} 113 %--------------------------------------------------TABLE-------------------------------------------------- 114 \begin{table}[!t] 114 115 %\begin{center} \begin{tabular}{|p{143pt}|l|l|} \hline 115 \begin{center} \begin{tabular}{|l|l|l|} \hline116 \begin{center} \begin{tabular}{|l|l|l|} \hline 116 117 Chapter \ref{STP} & - & model time STePping environment \\ \hline 117 118 Chapter \ref{DOM} & DOM & model DOMain \\ \hline … … 121 122 Chapter \ref{LBC} & LBC & Lateral Boundary Conditions (also OBC and BDY) \\ \hline 122 123 Chapter \ref{LDF} & LDF & Lateral DiFfusion (parameterisations) \\ \hline 123 Chapter \ref{ZDF} & ZDF & vertical (Z) DiFfusion \\ \hline124 Chapter \ref{ZDF} & ZDF & vertical (Z) DiFfusion (parameterisations) \\ \hline 124 125 Chapter \ref{OBS} & OBS & OBServation and model comparison \\ \hline 125 126 Chapter \ref{ASM} & ASM & ASsimilation increment \\ \hline 126 Chapter \ref{MISC} & ... & Miscellaneous topics (DIA, DTA, IOM, SOL, TRD, FLO...) \\ \hline 127 \end{tabular} \end{center} 128 \end{table} 127 Chapter \ref{MISC} & ... & Miscellaneous topics (DIA, DTA, IOM, \\ 128 & & SOL, TRD, FLO...) \\ \hline 129 Chapter \ref{CFG} & - & predefined configurations \\ \hline 130 \end{tabular} 131 \caption{ \label{Tab_chap} 132 Organization of Chapters which miminc the one of the model directories. } 133 \end{center} \end{table} 134 %-------------------------------------------------------------------------------------------------------------- 129 135 130 \vspace{1cm} Nota Bene : \vspace{0.25cm}131 136 132 137 \subsubsection{Changes between releases} 133 138 NEMO/OPA, like all research tools, is in perpetual evolution. The present document describes 134 139 the OPA version include in the release 3.3 of NEMO. This release differs significantly 135 from version 8, documented in \citet{Madec1998}. 140 from version 8, documented in \citet{Madec1998}.\\ 136 141 137 142 $\bullet$ The main modifications from OPA v8 and NEMO/OPA v3.2 are :\\ 143 \\ 138 144 (1) transition to full native \textsc{Fortran} 90, deep code restructuring and drastic 139 145 reduction of CPP keys; \\ … … 161 167 \vspace{1cm} 162 168 $\bullet$ The main modifications from NEMO/OPA v3.2 and v3.2 are :\\ 169 \\ 163 170 (1) introduction of a modified leapfrog-Asselin filter time stepping scheme \citep{Leclair_Madec_OM09}; \\ 164 171 (2) additional scheme for iso-neutral mixing \citep{Griffies_al_JPO98}, although it is still a "work in progress"; \\
Note: See TracChangeset
for help on using the changeset viewer.