Changeset 3764 for branches/2012/dev_MERGE_2012/DOC
- Timestamp:
- 2013-01-23T15:33:04+01:00 (12 years ago)
- Location:
- branches/2012/dev_MERGE_2012/DOC/TexFiles
- Files:
-
- 9 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_MERGE_2012/DOC/TexFiles/Biblio/Biblio.bib
r3680 r3764 1412 1412 } 1413 1413 1414 @ ARTICLE{Hunke2008,1414 @TECHREPORT{Hunke2008, 1415 1415 author = {E.C. Hunke and W.H. Lipscomb}, 1416 1416 title = {CICE: the Los Alamos sea ice model documentation and software user's manual, 1417 1417 Version 4.0}, 1418 institution = { Los Alamos National Laboratory, N.M.}, 1418 1419 publisher = {LA-CC-06-012, Los Alamos National Laboratory, N.M.}, 1419 1420 year = {2008} -
branches/2012/dev_MERGE_2012/DOC/TexFiles/Chapters/Chap_CFG.tex
r3294 r3764 31 31 32 32 % ================================================================ 33 % 1D model functionality33 % 1D model configuration 34 34 % ================================================================ 35 35 \section{Water column model: 1D model (C1D) (\key{c1d})} … … 48 48 49 49 The methodology is based on the use of the zoom functionality over the smallest possible 50 domain : a 3 x3 domain centred on the grid point of interest (see \S\ref{MISC_zoom}),50 domain : a 3x3 domain centred on the grid point of interest (see \S\ref{MISC_zoom}), 51 51 with some extra routines. There is no need to define a new mesh, bathymetry, 52 52 initial state or forcing, since the 1D model will use those of the configuration it is a zoom of. 53 The chosen grid point is set in par\_oce.F90module by setting the \jp{jpizoom} and \jp{jpjzoom}53 The chosen grid point is set in \mdl{par\_oce} module by setting the \jp{jpizoom} and \jp{jpjzoom} 54 54 parameters to the indices of the location of the chosen grid point. 55 55 56 The 1D model has some specifies. First, all the horizontal derivatives are assumed to be zero. 57 Therefore a simplified \rou{step} routine is used (\rou{step\_c1d}) in which both lateral tendancy 58 terms and lateral physics are not called, and the vertical velocity is zero (so far, no attempt at 59 introducing a Ekman pumping velocity has been made). 60 Second, the two components of the velocity are moved on a $T$-point. 61 This requires a specific treatment of the Coriolis term (see \rou{dyncor\_c1d}) and of the 62 dynamic time stepping (\rou{dynnxt\_c1d}). 63 All the relevant modules can be found in the NEMOGCM/NEMO/OPA\_SRC/C1D directory of 56 The 1D model has some specifies. First, all the horizontal derivatives are assumed to be zero, and 57 second, the two components of the velocity are moved on a $T$-point. 58 Therefore, defining \key{c1d} changes five main things in the code behaviour: 59 \begin{description} 60 \item[(1)] the lateral boundary condition routine (\rou{lbc\_lnk}) set the value of the central column 61 of the 3x3 domain is imposed over the whole domain ; 62 \item[(3)] a call to \rou{lbc\_lnk} is systematically done when reading input data ($i.e.$ in \mdl{iom}) ; 63 \item[(3)] a simplified \rou{stp} routine is used (\rou{stp\_c1d}, see \mdl{step\_c1d} module) in which 64 both lateral tendancy terms and lateral physics are not called ; 65 \item[(4)] the vertical velocity is zero (so far, no attempt at introducing a Ekman pumping velocity 66 has been made) ; 67 \item[(5)] a simplified treatment of the Coriolis term is performed as $U$- and $V$-points are the same 68 (see \mdl{dyncor\_c1d}). 69 \end{description} 70 All the relevant \textit{\_c1d} modules can be found in the NEMOGCM/NEMO/OPA\_SRC/C1D directory of 64 71 the \NEMO distribution. 65 72 … … 206 213 % ------------------------------------------------------------------------------------------------------------- 207 214 \section{GYRE family: double gyre basin (\key{gyre})} 208 \label{ MISC_config_gyre}215 \label{CFG_gyre} 209 216 210 217 The GYRE configuration \citep{Levy_al_OM10} have been built to simulated -
branches/2012/dev_MERGE_2012/DOC/TexFiles/Chapters/Chap_DIA.tex
r3680 r3764 1018 1018 In addition, a series of diagnostics has been added in the \mdl{diaar5}. 1019 1019 They corresponds to outputs that are required for AR5 simulations 1020 (see Section \ref{ MISC_steric} below for one of them).1020 (see Section \ref{DIA_steric} below for one of them). 1021 1021 Activating those outputs requires to define the \key{diaar5} CPP key. 1022 1022 \\ -
branches/2012/dev_MERGE_2012/DOC/TexFiles/Chapters/Chap_DOM.tex
r3680 r3764 499 499 Hybridation of the three main coordinates are available: $s-z$ or $s-zps$ coordinate 500 500 (Fig.~\ref{Fig_z_zps_s_sps}d and \ref{Fig_z_zps_s_sps}e). When using the variable 501 volume option \key{vvl} )($i.e.$ non-linear free surface), the coordinate follow the501 volume option \key{vvl} ($i.e.$ non-linear free surface), the coordinate follow the 502 502 time-variation of the free surface so that the transformation is time dependent: 503 503 $z(i,j,k,t)$ (Fig.~\ref{Fig_z_zps_s_sps}f). This option can be used with full step -
branches/2012/dev_MERGE_2012/DOC/TexFiles/Chapters/Chap_DYN.tex
r3294 r3764 127 127 This is of paramount importance. Replacing $T$ by the number $1$ in the tracer equation and summing 128 128 over the water column must lead to the sea surface height equation otherwise tracer content 129 will not be conserved \ ref{Griffies_al_MWR01, LeclairMadec2009}.129 will not be conserved \citep{Griffies_al_MWR01, Leclair_Madec_OM09}. 130 130 131 131 The vertical velocity is computed by an upward integration of the horizontal … … 189 189 the relative vorticity term and horizontal kinetic energy for the planetary vorticity 190 190 term (MIX scheme) ; or conserving both the potential enstrophy of horizontally non-divergent 191 flow and horizontal kinetic energy (EEN scheme) (see Appendix~\ref{Apdx_C_vor _zad}). In the191 flow and horizontal kinetic energy (EEN scheme) (see Appendix~\ref{Apdx_C_vorEEN}). In the 192 192 case of ENS, ENE or MIX schemes the land sea mask may be slightly modified to ensure the 193 193 consistency of vorticity term with analytical equations (\textit{ln\_dynvor\_con}=true). … … 331 331 This EEN scheme in fact combines the conservation properties of the ENS and ENE schemes. 332 332 It conserves both total energy and potential enstrophy in the limit of horizontally 333 nondivergent flow ($i.e.$ $\chi$=$0$) (see Appendix~\ref{Apdx_C_vor _zad}).333 nondivergent flow ($i.e.$ $\chi$=$0$) (see Appendix~\ref{Apdx_C_vorEEN}). 334 334 Applied to a realistic ocean configuration, it has been shown that it leads to a significant 335 335 reduction of the noise in the vertical velocity field \citep{Le_Sommer_al_OM09}. … … 938 938 is the \textit{before} velocity in time, except for the pure vertical component 939 939 that appears when a tensor of rotation is used. This latter term is solved 940 implicitly together with the vertical diffusion term (see \S\ref{ DOM_nxt})940 implicitly together with the vertical diffusion term (see \S\ref{STP}) 941 941 942 942 At the lateral boundaries either free slip, no slip or partial slip boundary … … 1066 1066 scheme (\np{ln\_zdfexp}=true) using a time splitting technique 1067 1067 (\np{nn\_zdfexp} $>$ 1) or $(b)$ a backward (or implicit) time differencing scheme 1068 (\np{ln\_zdfexp}=false) (see \S\ref{ DOM_nxt}). Note that namelist variables1068 (\np{ln\_zdfexp}=false) (see \S\ref{STP}). Note that namelist variables 1069 1069 \np{ln\_zdfexp} and \np{nn\_zdfexp} apply to both tracers and dynamics. 1070 1070 -
branches/2012/dev_MERGE_2012/DOC/TexFiles/Chapters/Chap_TRA.tex
r3308 r3764 264 264 transport) rather than TVD. The TVD scheme is implemented in the \mdl{traadv\_tvd} module. 265 265 266 For stability reasons (see \S\ref{ DOM_nxt}),266 For stability reasons (see \S\ref{STP}), 267 267 $\tau _u^{cen2}$ is evaluated in (\ref{Eq_tra_adv_tvd}) using the \textit{now} tracer while $\tau _u^{ups}$ 268 268 is evaluated using the \textit{before} tracer. In other words, the advective part of … … 337 337 \np{ln\_traadv\_ubs}=true. 338 338 339 For stability reasons (see \S\ref{ DOM_nxt}),339 For stability reasons (see \S\ref{STP}), 340 340 the first term in \eqref{Eq_tra_adv_ubs} (which corresponds to a second order centred scheme) 341 341 is evaluated using the \textit{now} tracer (centred in time) while the … … 451 451 except for the pure vertical component that appears when a rotation tensor 452 452 is used. This latter term is solved implicitly together with the 453 vertical diffusion term (see \S\ref{ DOM_nxt}).453 vertical diffusion term (see \S\ref{STP}). 454 454 455 455 % ------------------------------------------------------------------------------------------------------------- -
branches/2012/dev_MERGE_2012/DOC/TexFiles/Chapters/Chap_ZDF.tex
r3294 r3764 120 120 \end{equation} 121 121 122 is computed from the wind stress vector $|\tau|$ and the reference den dity $ \rho_o$.122 is computed from the wind stress vector $|\tau|$ and the reference density $ \rho_o$. 123 123 The final $h_{e}$ is further constrained by the adjustable bounds \np{rn\_mldmin} and \np{rn\_mldmax}. 124 124 Once $h_{e}$ is computed, the vertical eddy coefficients within $h_{e}$ are set to … … 1188 1188 \includegraphics[width=0.90\textwidth]{./TexFiles/Figures/Fig_ZDF_M2_K1_tmx.pdf} 1189 1189 \caption{ \label{Fig_ZDF_M2_K1_tmx} 1190 (a) M2 and (b) K 2internal wave drag energy from \citet{Carrere_Lyard_GRL03} ($W/m^2$). }1190 (a) M2 and (b) K1 internal wave drag energy from \citet{Carrere_Lyard_GRL03} ($W/m^2$). } 1191 1191 \end{center} \end{figure} 1192 1192 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 1205 1205 1206 1206 When \np{ln\_tmx\_itf}=true, the two key parameters $q$ and $F(z)$ are adjusted following 1207 the parameterisation developed by \ ref{Koch-Larrouy_al_GRL07}:1207 the parameterisation developed by \citet{Koch-Larrouy_al_GRL07}: 1208 1208 1209 1209 First, the Indonesian archipelago is a complex geographic region … … 1219 1219 Second, the vertical structure function, $F(z)$, is no more associated 1220 1220 with a bottom intensification of the mixing, but with a maximum of 1221 energy available within the thermocline. \ ref{Koch-Larrouy_al_GRL07}1221 energy available within the thermocline. \citet{Koch-Larrouy_al_GRL07} 1222 1222 have suggested that the vertical distribution of the energy dissipation 1223 1223 proportional to $N^2$ below the core of the thermocline and to $N$ above. … … 1236 1236 and vertical distributions of the mixing are adequately prescribed 1237 1237 \citep{Koch-Larrouy_al_GRL07, Koch-Larrouy_al_OD08a, Koch-Larrouy_al_OD08b}. 1238 Note also that such a parameterisation has a s ugnificant impact on the behaviour1238 Note also that such a parameterisation has a significant impact on the behaviour 1239 1239 of global coupled GCMs \citep{Koch-Larrouy_al_CD10}. 1240 1240 -
branches/2012/dev_MERGE_2012/DOC/TexFiles/Namelist/namasm
r3294 r3764 3 3 !----------------------------------------------------------------------- 4 4 ln_bkgwri = .false. ! Logical switch for writing out background state 5 ln_trjwri = .false. ! Logical switch for writing out state trajectory6 5 ln_trainc = .false. ! Logical switch for applying tracer increments 7 6 ln_dyninc = .false. ! Logical switch for applying velocity increments … … 14 13 nitiaufin = 15 ! Timestep of end of IAU interval in [0,nitend-nit000-1] 15 14 niaufn = 0 ! Type of IAU weighting function 16 nittrjfrq = 0 ! Frequency of trajectory output for 4D-VAR17 15 ln_salfix = .false. ! Logical switch for ensuring that the sa > salfixmin 18 16 salfixmin = -9999 ! Minimum salinity after applying the increments -
branches/2012/dev_MERGE_2012/DOC/TexFiles/Namelist/namdyn_vor
r3294 r3764 2 2 &namdyn_vor ! option of physics/algorithm (not control by CPP keys) 3 3 !----------------------------------------------------------------------- 4 ln_dynvor_ene = .false. ! en strophyconserving scheme5 ln_dynvor_ens = .false. ! en ergy conserving scheme4 ln_dynvor_ene = .false. ! energy conserving scheme 5 ln_dynvor_ens = .false. ! enstrophy conserving scheme 6 6 ln_dynvor_mix = .false. ! mixed scheme 7 7 ln_dynvor_een = .true. ! energy & enstrophy scheme
Note: See TracChangeset
for help on using the changeset viewer.