Changeset 3989 for branches/2013/dev_r3853_CNRS9_ConfSetting/DOC
- Timestamp:
- 2013-07-24T11:48:35+02:00 (11 years ago)
- Location:
- branches/2013/dev_r3853_CNRS9_ConfSetting/DOC
- Files:
-
- 55 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/NEMO_book.tex
r3294 r3989 177 177 \newcommand{\rou} [1] {\textit{#1}\index{Routines!#1}} %module (routine) 178 178 \newcommand{\hf} [1] {\textit{#1.h90}\index{h90 file!#1}} %module (h90 files) 179 \newcommand{\np} [1] {\textit{#1}\index{Namelist parameters!#1}} %namelist parameter (nampar) 179 \newcommand{\ngn} [1] {\textit{#1}\index{Namelist Group Name!#1}} %namelist name (nampar) 180 \newcommand{\np} [1] {\textit{#1}\index{Namelist variables!#1}} %namelist variable 180 181 \newcommand{\jp} [1] {\textit{#1}\index{Model parameters!#1}} %model parameter (jp) 181 182 \newcommand{\pp} [1] {\textit{#1}\index{Model parameters!#1}} %namelist parameter (pp) … … 245 246 % ================================================================ 246 247 % ================================================================ 248 247 249 \begin{document} 248 250 … … 291 293 \include{./TexFiles/Chapters/Chap_ZDF} % Vertical diffusion 292 294 293 \include{./TexFiles/Chapters/Chap_DIA} % Miscellaneous topics295 \include{./TexFiles/Chapters/Chap_DIA} % Outputs and Diagnostics 294 296 295 297 \include{./TexFiles/Chapters/Chap_OBS} % Observation operator -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Annex_ISO.tex
r3297 r3989 7 7 \minitoc 8 8 \pagebreak 9 \section{Choice of namelist parameters}9 \section{Choice of \ngn{namtra\_ldf} namelist parameters} 10 10 %-----------------------------------------nam_traldf------------------------------------------------------ 11 11 \namdisplay{namtra_ldf} … … 26 26 \np{rn\_aeiv\_0}. If 2D-varying coefficients are set with 27 27 \key{traldf\_c2d} then $A_l$ is reduced in proportion with horizontal 28 scale factor according to \eqref{Eq_title} \footnote{Except in global 29 $0.5^{\circ}$ runs (\key{orca\_r05})with \key{traldf\_eiv}, where28 scale factor according to \eqref{Eq_title} \footnote{Except in global ORCA 29 $0.5^{\circ}$ runs with \key{traldf\_eiv}, where 30 30 $A_l$ is set like $A_e$ but with a minimum vale of 31 31 $100\;\mathrm{m}^2\;\mathrm{s}^{-1}$}. In idealised setups with 32 32 \key{traldf\_c2d}, $A_e$ is reduced similarly, but if \key{traldf\_eiv} 33 is set in the global configurations \key{orca\_r2}, \key{orca\_r1} or 34 \key{orca\_r05} with \key{traldf\_c2d}, a horizontally varying $A_e$ is 33 is set in the global configurations with \key{traldf\_c2d}, a horizontally varying $A_e$ is 35 34 instead set from the Held-Larichev parameterisation\footnote{In this 36 35 case, $A_e$ at low latitudes $|\theta|<20^{\circ}$ is further -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Chap_ASM.tex
r3294 r3989 18 18 assimilation code. The code can also output model background fields which are used 19 19 as an input to data assimilation code. This is all controlled by the namelist 20 \textit{ nam\_asminc}. There is a brief description of all the namelist options20 \textit{\ngn{nam\_asminc} }. There is a brief description of all the namelist options 21 21 provided. To build the ASM code \key{asminc} must be set. 22 22 … … 125 125 \label{ASM_details} 126 126 127 Here we show an example namelist and the header of an example assimilation127 Here we show an example \ngn{namasm} namelist and the header of an example assimilation 128 128 increments file on the ORCA2 grid. 129 129 -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Chap_CFG.tex
r3764 r3989 1 1 % ================================================================ 2 % Chapter ÑConfigurations2 % Chapter � Configurations 3 3 % ================================================================ 4 4 \chapter{Configurations} … … 16 16 17 17 18 The purpose of this part of the manual is to introduce the \NEMO predefined configuration.18 The purpose of this part of the manual is to introduce the \NEMO reference configurations. 19 19 These configurations are offered as means to explore various numerical and physical options, 20 20 thus allowing the user to verify that the code is performing in a manner consistent with that 21 21 we are running. This form of verification is critical as one adopts the code for his or her particular 22 22 research purposes. The test cases also provide a sense for some of the options available 23 in the code, though by no means are all options exercised in the predefined configurations. 24 25 26 %There is several predefined ocean configuration which use is controlled by a specific CPP key. 27 28 %The key set the domain sizes (\jp{jpiglo}, \jp{jpjglo}, \jp{jpk}), the mesh and the bathymetry, 29 %and, in some cases, add to the model physics some specific treatments. 30 23 in the code, though by no means are all options exercised in the reference configurations. 24 25 Configuration is defined mainly through the \ngn{namcfg} namelist variables: 26 %------------------------------------------namcfg---------------------------------------------------- 27 \namdisplay{namcfg} 28 %------------------------------------------------------------------------------------------------------------- 31 29 32 30 % ================================================================ 33 31 % 1D model configuration 34 32 % ================================================================ 35 \section{Water column model: 1D model (C1D) (\key{c1d}) }33 \section{Water column model: 1D model (C1D) (\key{c1d}) } 36 34 \label{CFG_c1d} 37 35 38 36 The 1D model option simulates a stand alone water column within the 3D \NEMO system. 39 37 It can be applied to the ocean alone or to the ocean-ice system and can include passive tracers 40 or a biogeochemical model. It is set up by defining the \key{c1d} CPP key. 38 or a biogeochemical model. It is set up by defining the position of the 1D water column in the grid 39 (see \textit{CONFIG/SHARED/namelist\_ref} ). 41 40 The 1D model is a very useful tool 42 41 \textit{(a)} to learn about the physics and numerical treatment of vertical mixing processes ; … … 48 47 49 48 The methodology is based on the use of the zoom functionality over the smallest possible 50 domain : a 3x3 domain cent red on the grid point of interest (see \S\ref{MISC_zoom}),49 domain : a 3x3 domain centered on the grid point of interest, 51 50 with some extra routines. There is no need to define a new mesh, bathymetry, 52 51 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 \ mdl{par\_oce} module by setting the \jp{jpizoom} and \jp{jpjzoom}52 The chosen grid point is set in \textit{\ngn{namcfg}} namelist by setting the \np{jpizoom} and \np{jpjzoom} 54 53 parameters to the indices of the location of the chosen grid point. 55 54 … … 76 75 % ORCA family configurations 77 76 % ================================================================ 78 \section{ORCA family: global ocean with tripolar grid (\key{orca\_rX})}77 \section{ORCA family: global ocean with tripolar grid } 79 78 \label{CFG_orca} 80 79 … … 82 81 the LIM sea-ice model (ORCA-LIM) and possibly with PISCES biogeochemical model 83 82 (ORCA-LIM-PISCES), using various resolutions. 83 The appropriate \textit{\&namcfg} namelist is available in \textit{CONFIG/ORCA2\_LIM/EXP00/namelist\_cfg} 84 for ORCA2 and in \textit{CONFIG/SHARED/README\_other\_configurations\_namelist\_namcfg} 85 for other resolutions 84 86 85 87 … … 147 149 The NEMO system is provided with five built-in ORCA configurations which differ in the 148 150 horizontal resolution. The value of the resolution is given by the resolution at the Equator 149 expressed in degrees. Each of configuration is set through a CPP key, \key{orca\_rX}150 (with X being an indicator of the resolution), which setthe grid size and configuration151 name parameters (Tab. ~\ref{Tab_ORCA}).151 expressed in degrees. Each of configuration is set through the \textit{\ngn{namcfg}} namelist, 152 which sets the grid size and configuration 153 name parameters (Tab. \ref{Tab_ORCA}). 152 154 . 153 155 … … 155 157 \begin{table}[!t] \begin{center} 156 158 \begin{tabular}{p{4cm} c c c c} 157 CPP key & \jp{jp\_cfg} & \jp{jpiglo} & \jp{jpiglo} & \\159 Horizontal Grid & \np{jp\_cfg} & \np{jpiglo} & \np{jpjglo} & \\ 158 160 \hline \hline 159 \ key{orca\_r4}& 4 & 92 & 76 & \\160 \ key{orca\_r2}& 2 & 182 & 149 & \\161 \ key{orca\_r1}& 1 & 362 & 292 & \\162 \ key{orca\_r05}& 05 & 722 & 511 & \\163 \ key{orca\_r025}& 025 & 1442 & 1021 & \\161 \~4\deg & 4 & 92 & 76 & \\ 162 \~2\deg & 2 & 182 & 149 & \\ 163 \~1\deg & 1 & 362 & 292 & \\ 164 \~0.5\deg & 05 & 722 & 511 & \\ 165 \~0.25\deg & 025 & 1442 & 1021 & \\ 164 166 %\key{orca\_r8} & 8 & 2882 & 2042 & \\ 165 167 %\key{orca\_r12} & 12 & 4322 & 3062 & \\ … … 168 170 \caption{ \label{Tab_ORCA} 169 171 Set of predefined parameters for ORCA family configurations. 170 In all cases, the name of the configuration is set to "orca" ($i.e.$ \ jp{cp\_cfg}~=~orca). }172 In all cases, the name of the configuration is set to "orca" ($i.e.$ \np{cp\_cfg}~=~orca). } 171 173 \end{center} 172 174 \end{table} … … 197 199 in the upper 150m (see Tab.~\ref{Tab_orca_zgr} and Fig.~\ref{Fig_zgr}). 198 200 The bottom topography and the coastlines are derived from the global atlas of Smith and Sandwell (1997). 199 The default forcing employthe boundary forcing from \citet{Large_Yeager_Rep04} (see \S\ref{SBC_blk_core}),201 The default forcing uses the boundary forcing from \citet{Large_Yeager_Rep04} (see \S\ref{SBC_blk_core}), 200 202 which was developed for the purpose of running global coupled ocean-ice simulations 201 203 without an interactive atmosphere. This \citet{Large_Yeager_Rep04} dataset is available … … 205 207 206 208 ORCA\_R2 pre-defined configuration can also be run with an AGRIF zoom over the Agulhas 207 current area ( \key{agrif} defined) and, by setting the key \key{arctic} or \key{antarctic}, 209 current area ( \key{agrif} defined) and, by setting the appropriate variables in 210 \textit{\&namcfg}, see \textit{CONFIG/SHARED/namelist\_ref} 208 211 a regional Arctic or peri-Antarctic configuration is extracted from an ORCA\_R2 or R05 configurations 209 212 using sponge layers at open boundaries. … … 212 215 % GYRE family: double gyre basin 213 216 % ------------------------------------------------------------------------------------------------------------- 214 \section{GYRE family: double gyre basin (\key{gyre})}217 \section{GYRE family: double gyre basin } 215 218 \label{CFG_gyre} 216 219 217 The GYRE configuration \citep{Levy_al_OM10} ha ve been built to simulated218 the seasonal cycle of a double-gyre box model. It consist in an idealized domain220 The GYRE configuration \citep{Levy_al_OM10} has been built to simulate 221 the seasonal cycle of a double-gyre box model. It consists in an idealized domain 219 222 similar to that used in the studies of \citet{Drijfhout_JPO94} and \citet{Hazeleger_Drijfhout_JPO98, 220 223 Hazeleger_Drijfhout_JPO99, Hazeleger_Drijfhout_JGR00, Hazeleger_Drijfhout_JPO00}, … … 242 245 uniformly applied to the whole domain. 243 246 244 The GYRE configuration is set through the \key{gyre} CPP key. Its horizontal resolution 245 (and thus the size of the domain) is determined by setting \jp{jp\_cfg} in \hf{par\_GYRE} file: \\ 246 \jp{jpiglo} $= 30 \times$ \jp{jp\_cfg} + 2 \\ 247 \jp{jpjglo} $= 20 \times$ \jp{jp\_cfg} + 2 \\ 248 Obviously, the namelist parameters have to be adjusted to the chosen resolution. 249 In the vertical, GYRE uses the default 30 ocean levels (\jp{jpk}=31) (Fig.~\ref{Fig_zgr}). 247 The GYRE configuration is set through the \textit{\&namcfg} namelist defined in the reference 248 configuration \textit{CONFIG/GYRE/EXP00/namelist\_cfg}. Its horizontal resolution 249 (and thus the size of the domain) is determined by setting \np{jp\_cfg} : \\ 250 \np{jpiglo} $= 30 \times$ \np{jp\_cfg} + 2 \\ 251 \np{jpjglo} $= 20 \times$ \np{jp\_cfg} + 2 \\ 252 Obviously, the namelist parameters have to be adjusted to the chosen resolution, see the Configurations 253 pages on the NEMO web site (Using NEMO\/Configurations) . 254 In the vertical, GYRE uses the default 30 ocean levels (\pp{jpk}=31) (Fig.~\ref{Fig_zgr}). 250 255 251 256 The GYRE configuration is also used in benchmark test as it is very simple to increase … … 270 275 271 276 \begin{description} 272 \item[ \key{eel\_r2}] to be described....273 \item[ \key{eel\_r5}]274 \item[ \key{eel\_r6}]277 \item[eel\_r2] to be described.... 278 \item[eel\_r5] 279 \item[eel\_r6] 275 280 \end{description} 276 281 The appropriate \textit{\&namcfg} namelists are available in 282 \textit{CONFIG/SHARED/README\_other\_configurations\_namelist\_namcfg} 277 283 % ------------------------------------------------------------------------------------------------------------- 278 284 % AMM configuration 279 285 % ------------------------------------------------------------------------------------------------------------- 280 \section{AMM: atlantic margin configuration (\key{amm\_12km})}286 \section{AMM: atlantic margin configuration } 281 287 \label{MISC_config_AMM} 282 288 283 289 The AMM, Atlantic Margins Model, is a regional model covering the 284 290 Northwest European Shelf domain on a regular lat-lon grid at 285 approximately 12km horizontal resolution. The key \key{amm\_12km} 286 is used to create the correct dimensions of the AMM domain. 291 approximately 12km horizontal resolution. The appropriate 292 \textit{\&namcfg} namelist is available in \textit{CONFIG/AMM12/EXP00/namelist\_cfg}. 293 It is used to build the correct dimensions of the AMM domain. 287 294 288 295 This configuration tests several features of NEMO functionality specific -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Chap_DIA.tex
r3764 r3989 1 1 % ================================================================ 2 % Chapter ÑI/O & Diagnostics2 % Chapter � I/O & Diagnostics 3 3 % ================================================================ 4 4 \chapter{Ouput and Diagnostics (IOM, DIA, TRD, FLO)} … … 609 609 set via an equivalent and identically named namelist to \textit{namnc4} 610 610 in \np{xmlio\_server.def}. Typically this namelist serves the mean files 611 whilst the \n p{ namnc4} in the main namelist file continues to serve the611 whilst the \ngn{ namnc4} in the main namelist file continues to serve the 612 612 restart files. This duplication is unfortunate but appropriate since, if 613 613 using io\_servers, the domain sizes of the individual files produced by the … … 631 631 trend of the dynamics and/or temperature and salinity time evolution equations 632 632 is stored in three-dimensional arrays just after their computation ($i.e.$ at the end 633 of each $dyn\cdots.F90$ and/or $tra\cdots.F90$ routines). These trends are then 633 of each $dyn\cdots.F90$ and/or $tra\cdots.F90$ routines). Options are defined by 634 \ngn{namtrd} namelist variables. These trends are then 634 635 used in \mdl{trdmod} (see TRD directory) every \textit{nn\_trd } time-steps. 635 636 … … 675 676 The on-line computation of floats advected either by the three dimensional velocity 676 677 field or constraint to remain at a given depth ($w = 0$ in the computation) have been 677 introduced in the system during the CLIPPER project. The algorithm used is based 678 introduced in the system during the CLIPPER project. Options are defined by \ngn{namflo} 679 namelis variables. The algorithm used is based 678 680 either on the work of \cite{Blanke_Raynaud_JPO97} (default option), or on a $4^th$ 679 681 Runge-Hutta algorithm (\np{ln\_flork4}=true). Note that the \cite{Blanke_Raynaud_JPO97} … … 687 689 688 690 689 In case of Ariane convention, input filename is \np{ "init\_float\_ariane"}. Its format is:691 In case of Ariane convention, input filename is \np{init\_float\_ariane}. Its format is: 690 692 691 693 \texttt{ I J K nisobfl itrash itrash } … … 709 711 710 712 711 In the other case ( longitude and latitude ), input filename is \np{"init\_float"}. Its format is:713 In the other case ( longitude and latitude ), input filename is init\_float. Its format is: 712 714 713 715 \texttt{ Long Lat depth nisobfl ngrpfl itrash} … … 732 734 733 735 \np{jpnfl} is the total number of floats during the run. 734 When initial positions are read in a restart file ( \np{ln\_rstflo = .TRUE.}), \np{jpnflnewflo}736 When initial positions are read in a restart file ( \np{ln\_rstflo}= .TRUE. ), \np{jpnflnewflo} 735 737 can be added in the initialization file. 736 738 … … 740 742 is the frequency of creation of the float restart file. 741 743 742 Output data can be written in ascii files (\np{ln\_flo\_ascii = .TRUE.}). In that case,743 output filename is \np{is trajec\_float}.744 745 Another possiblity of writing format is Netcdf (\np{ln\_flo\_ascii = .FALSE.}). There are 2 possibilities:746 747 - if (\key{iomput}) is used, outputs are selected in \np{iodef.xml}. Here it is an example of specification744 Output data can be written in ascii files (\np{ln\_flo\_ascii} = .TRUE. ). In that case, 745 output filename is trajec\_float. 746 747 Another possiblity of writing format is Netcdf (\np{ln\_flo\_ascii} = .FALSE. ). There are 2 possibilities: 748 749 - if (\key{iomput}) is used, outputs are selected in iodef.xml. Here it is an example of specification 748 750 to put in files description section: 749 751 … … 768 770 769 771 770 - if (\key{iomput}) is not used, a file called \np{trajec\_float.nc}will be created by IOIPSL library.772 - if (\key{iomput}) is not used, a file called trajec\_float.nc will be created by IOIPSL library. 771 773 772 774 … … 789 791 %---------------------------------------------------------------------------------------------------------- 790 792 791 Concerning the on-line Harmonic analysis, some parameters are available in namelist: 793 Concerning the on-line Harmonic analysis, some parameters are available in namelist 794 \ngn{namdia\_harm} : 792 795 793 796 - \texttt{nit000\_han} is the first time step used for harmonic analysis … … 841 844 842 845 843 Namelist parameters control how frequently the flows are summed and the time scales over which844 they are averaged, as well as the level of output for debugging:846 Namelist variables in \ngn{namdct} control how frequently the flows are summed 847 and the time scales over which they are averaged, as well as the level of output for debugging: 845 848 846 849 %------------------------------------------namdct---------------------------------------------------- … … 992 995 993 996 The poleward heat and salt transports, their advective and diffusive component, and 994 the meriodional stream function can be computed on-line in \mdl{diaptr} by setting995 \np{ln\_diaptr} to true (see the \textit{ namptr} namelist below).997 the meriodional stream function can be computed on-line in \mdl{diaptr} 998 \np{ln\_diaptr} to true (see the \textit{\ngn{namptr} } namelist below). 996 999 When \np{ln\_subbas}~=~true, transports and stream function are computed 997 1000 for the Atlantic, Indian, Pacific and Indo-Pacific Oceans (defined north of 30\deg S) -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Chap_DOM.tex
r3764 r3989 1 1 % ================================================================ 2 % Chapter 2 ÑSpace and Time Domain (DOM)2 % Chapter 2 � Space and Time Domain (DOM) 3 3 % ================================================================ 4 4 \chapter{Space Domain (DOM) } … … 24 24 directory routines as well as the DOM (DOMain) directory. 25 25 26 $\ $\newline % force a new lign e26 $\ $\newline % force a new lign 27 27 28 28 % ================================================================ … … 274 274 \label{DOM_size} 275 275 276 The total size of the computational domain is set by the parameters \jp{jpiglo}, 277 \jp{jpjglo} and \jp{jpk} in the $i$, $j$ and $k$ directions respectively. They are 278 given as parameters in the \mdl{par\_oce} module\footnote{When a specific 279 configuration is used (ORCA2 global ocean, etc...) the parameter are actually 280 defined in additional files introduced by \mdl{par\_oce} module via CPP 281 \textit{include} command. For example, ORCA2 parameters are set in 282 \textit{par\_ORCA\_R2.h90} file}. The use of parameters rather than variables 283 (together with dynamic allocation of arrays) was chosen because it ensured that 284 the compiler would optimize the executable code efficiently, especially on vector 285 machines (optimization may be less efficient when the problem size is unknown 286 at the time of compilation). Nevertheless, it is possible to set up the code with full 287 dynamical allocation by using the AGRIF packaged \citep{Debreu_al_CG2008}. 288 % 289 \gmcomment{ add the following ref 290 \colorbox{yellow}{(ref part of the doc)} } 291 % 292 Note that are other parameters in \mdl{par\_oce} that refer to the domain size. 293 The two parameters $jpidta$ and $jpjdta$ may be larger than $jpiglo$, $jpjglo$ 276 The total size of the computational domain is set by the parameters \np{jpiglo}, 277 \np{jpjglo} and \np{jpkdta} in the $i$, $j$ and $k$ directions respectively. They are 278 given as namelist variables in the \ngn{namcfg} namelist. 279 280 Note that are other namelist variables in the \ngn{namcfg} namelist that refer to 281 the domain size. 282 The two variables \np{jpidta} and \np{jpjdta} may be larger than \np{jpiglo}, \np{jpjglo} 294 283 when the user wants to use only a sub-region of a given configuration. This is 295 284 the "zoom" capability described in \S\ref{MISC_zoom}. In most applications of … … 300 289 301 290 302 $\ $\newline % force a new lign e291 $\ $\newline % force a new lign 303 292 304 293 % ================================================================ … … 388 377 389 378 The user has three options available in defining a horizontal grid, which involve 390 the parameter $jphgr\_mesh$ of the \mdl{par\_oce} module.379 the namelist variable \np{jphgr\_mesh} of the \ngn{namcfg} namelist. 391 380 \begin{description} 392 \item[\ jp{jphgr\_mesh}=0] The most general curvilinear orthogonal grids.381 \item[\np{jphgr\_mesh}=0] The most general curvilinear orthogonal grids. 393 382 The coordinates and their first derivatives with respect to $i$ and $j$ are provided 394 383 in a input file (\ifile{coordinates}), read in \rou{hgr\_read} subroutine of the domhgr module. 395 \item[\ jp{jphgr\_mesh}=1 to 5] A few simple analytical grids are provided (see below).384 \item[\np{jphgr\_mesh}=1 to 5] A few simple analytical grids are provided (see below). 396 385 For other analytical grids, the \mdl{domhgr} module must be modified by the user. 397 386 \end{description} 398 387 399 388 There are two simple cases of geographical grids on the sphere. With 400 \ jp{jphgr\_mesh}=1, the grid (expressed in degrees) is regular in space,401 with grid sizes specified by parameters \ pp{ppe1\_deg} and \pp{ppe2\_deg},389 \np{jphgr\_mesh}=1, the grid (expressed in degrees) is regular in space, 390 with grid sizes specified by parameters \np{ppe1\_deg} and \np{ppe2\_deg}, 402 391 respectively. Such a geographical grid can be very anisotropic at high latitudes 403 392 because of the convergence of meridians (the zonal scale factors $e_1$ 404 393 become much smaller than the meridional scale factors $e_2$). The Mercator 405 grid (\ jp{jphgr\_mesh}=4) avoids this anisotropy by refining the meridional scale394 grid (\np{jphgr\_mesh}=4) avoids this anisotropy by refining the meridional scale 406 395 factors in the same way as the zonal ones. In this case, meridional scale factors 407 396 and latitudes are calculated analytically using the formulae appropriate for 408 a Mercator projection, based on \ pp{ppe1\_deg} which is a reference grid spacing397 a Mercator projection, based on \np{ppe1\_deg} which is a reference grid spacing 409 398 at the equator (this applies even when the geographical equator is situated outside 410 399 the model domain). … … 412 401 \gmcomment{ give here the analytical expression of the Mercator mesh} 413 402 %%% 414 In these two cases (\ jp{jphgr\_mesh}=1 or 4), the grid position is defined by the415 longitude and latitude of the south-westernmost point (\ pp{ppglamt0}416 and \ pp{ppgphi0}). Note that for the Mercator grid the user need only provide403 In these two cases (\np{jphgr\_mesh}=1 or 4), the grid position is defined by the 404 longitude and latitude of the south-westernmost point (\np{ppglamt0} 405 and \np{ppgphi0}). Note that for the Mercator grid the user need only provide 417 406 an approximate starting latitude: the real latitude will be recalculated analytically, 418 407 in order to ensure that the equator corresponds to line passing through $t$- … … 420 409 421 410 Rectangular grids ignoring the spherical geometry are defined with 422 \ jp{jphgr\_mesh} = 2, 3, 5. The domain is either an $f$-plane (\jp{jphgr\_mesh} = 2,423 Coriolis factor is constant) or a beta-plane (\ jp{jphgr\_mesh} = 3, the Coriolis factor411 \np{jphgr\_mesh} = 2, 3, 5. The domain is either an $f$-plane (\np{jphgr\_mesh} = 2, 412 Coriolis factor is constant) or a beta-plane (\np{jphgr\_mesh} = 3, the Coriolis factor 424 413 is linear in the $j$-direction). The grid size is uniform in meter in each direction, 425 and given by the parameters \ pp{ppe1\_m} and \pp{ppe2\_m} respectively.414 and given by the parameters \np{ppe1\_m} and \np{ppe2\_m} respectively. 426 415 The zonal grid coordinate (\textit{glam} arrays) is in kilometers, starting at zero 427 416 with the first $t$-point. The meridional coordinate (gphi. arrays) is in kilometers, 428 417 and the second $t$-point corresponds to coordinate $gphit=0$. The input 429 parameter \pp{ppglam0} is ignored. \pp{ppgphi0} is used to set the reference418 variable \np{ppglam0} is ignored. \np{ppgphi0} is used to set the reference 430 419 latitude for computation of the Coriolis parameter. In the case of the beta plane, 431 \ pp{ppgphi0} corresponds to the center of the domain. Finally, the special case432 \ jp{jphgr\_mesh}=5 corresponds to a beta plane in a rotated domain for the420 \np{ppgphi0} corresponds to the center of the domain. Finally, the special case 421 \np{jphgr\_mesh}=5 corresponds to a beta plane in a rotated domain for the 433 422 GYRE configuration, representing a classical mid-latitude double gyre system. 434 423 The rotation allows us to maximize the jet length relative to the gyre areas … … 436 425 437 426 The choice of the grid must be consistent with the boundary conditions specified 438 by the parameter \ jp{jperio} (see {\S\ref{LBC}).427 by the parameter \np{jperio} (see {\S\ref{LBC}). 439 428 440 429 % ------------------------------------------------------------------------------------------------------------- … … 446 435 All the arrays relating to a particular ocean model configuration (grid-point 447 436 position, scale factors, masks) can be saved in files if $\np{nn\_msh} \not= 0$ 448 (namelist parameter). This can be particularly useful for plots and off-line437 (namelist variable in \ngn{namdom}). This can be particularly useful for plots and off-line 449 438 diagnostics. In some cases, the user may choose to make a local modification 450 439 of a scale factor in the code. This is the case in global configurations when … … 454 443 the output grid written when $\np{nn\_msh} \not=0$ is no more equal to the input grid. 455 444 456 $\ $\newline % force a new lign e445 $\ $\newline % force a new lign 457 446 458 447 % ================================================================ … … 467 456 %------------------------------------------------------------------------------------------------------------- 468 457 458 Variables are defined through the \ngn{namzgr} and \ngn{namdom} namelists. 469 459 In the vertical, the model mesh is determined by four things: 470 460 (1) the bathymetry given in meters ; … … 553 543 \item[\np{nn\_bathy} = 0] a flat-bottom domain is defined. The total depth $z_w (jpk)$ 554 544 is given by the coordinate transformation. The domain can either be a closed 555 basin or a periodic channel depending on the parameter \ jp{jperio}.545 basin or a periodic channel depending on the parameter \np{jperio}. 556 546 \item[\np{nn\_bathy} = -1] a domain with a bump of topography one third of the 557 547 domain width at the central latitude. This is meant for the "EEL-R5" configuration, … … 599 589 vertical scale factors. The user must provide the analytical expression of both 600 590 $z_0$ and its first derivative with respect to $k$. This is done in routine \mdl{domzgr} 601 through statement functions, using parameters provided in the \ textit{par\_oce.h90} file.602 603 It is possible to define a simple regular vertical grid by giving zero stretching (\ pp{ppacr=0}).604 In that case, the parameters \jp{jpk} (number of $w$-levels) and \ pp{pphmax}591 through statement functions, using parameters provided in the \ngn{namcfg} namelist. 592 593 It is possible to define a simple regular vertical grid by giving zero stretching (\np{ppacr=0}). 594 In that case, the parameters \jp{jpk} (number of $w$-levels) and \np{pphmax} 605 595 (total ocean depth in meters) fully define the grid. 606 596 … … 639 629 scale factors as a function of the model levels are shown in Fig.~\ref{Fig_zgr} and 640 630 given in Table \ref{Tab_orca_zgr}. Those values correspond to the parameters 641 \ pp{ppsur}, \pp{ppa0}, \pp{ppa1}, \pp{ppkth} in the parameter file \mdl{par\_oce}.631 \np{ppsur}, \np{ppa0}, \np{ppa1}, \np{ppkth} in \ngn{namcfg} namelist. 642 632 643 633 Rather than entering parameters $h_{sur}$, $h_{0}$, and $h_{1}$ directly, it is 644 634 possible to recalculate them. In that case the user sets 645 \ pp{ppsur}=\pp{ppa0}=\pp{ppa1}=\pp{pp\_to\_be\_computed}, in \mdl{par\_oce},635 \np{ppsur}=\np{ppa0}=\np{ppa1}=999999., in \ngn{namcfg} namelist, 646 636 and specifies instead the four following parameters: 647 637 \begin{itemize} 648 \item \ pp{ppacr}=$h_{cr} $: stretching factor (nondimensional). The larger649 \ pp{ppacr}, the smaller the stretching. Values from $3$ to $10$ are usual.650 \item \ pp{ppkth}=$h_{th} $: is approximately the model level at which maximum638 \item \np{ppacr}=$h_{cr} $: stretching factor (nondimensional). The larger 639 \np{ppacr}, the smaller the stretching. Values from $3$ to $10$ are usual. 640 \item \np{ppkth}=$h_{th} $: is approximately the model level at which maximum 651 641 stretching occurs (nondimensional, usually of order 1/2 or 2/3 of \jp{jpk}) 652 \item \ pp{ppdzmin}: minimum thickness for the top layer (in meters)653 \item \ pp{pphmax}: total depth of the ocean (meters).642 \item \np{ppdzmin}: minimum thickness for the top layer (in meters) 643 \item \np{pphmax}: total depth of the ocean (meters). 654 644 \end{itemize} 655 645 As an example, for the $45$ layers used in the DRAKKAR configuration those 656 parameters are: \jp{jpk}=46, \ pp{ppacr}=9, \pp{ppkth}=23.563, \pp{ppdzmin}=6m,657 \ pp{pphmax}=5750m.646 parameters are: \jp{jpk}=46, \np{ppacr}=9, \np{ppkth}=23.563, \np{ppdzmin}=6m, 647 \np{pphmax}=5750m. 658 648 659 649 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 720 710 is allowed to have either a smaller or larger thickness than $e_{3t}(jpk)$: the 721 711 maximum thickness allowed is $2*e_{3t}(jpk-1)$. This has to be kept in mind when 722 specifying the maximum depth \pp{pphmax} in partial steps: for example, with 723 \pp{pphmax}$=5750~m$ for the DRAKKAR 45 layer grid, the maximum ocean depth 712 specifying values in \ngn{namdom} namelist, as the maximum depth \np{pphmax} 713 in partial steps: for example, with 714 \np{pphmax}$=5750~m$ for the DRAKKAR 45 layer grid, the maximum ocean depth 724 715 allowed is actually $6000~m$ (the default thickness $e_{3t}(jpk-1)$ being $250~m$). 725 716 Two variables in the namdom namelist are used to define the partial step … … 740 731 \namdisplay{namzgr_sco} 741 732 %-------------------------------------------------------------------------------------------------------------- 733 Options are defined in \ngn{namzgr\_sco}. 742 734 In $s$-coordinate (\np{ln\_sco}~=~true), the depth and thickness of the model 743 735 levels are defined from the product of a depth field and either a stretching … … 905 897 %------------------------------------------------------------------------------------------ 906 898 899 Options are defined in \ngn{namtsd}. 907 900 By default, the ocean start from rest (the velocity field is set to zero) and the initialization of 908 901 temperature and salinity fields is controlled through the \np{ln\_tsd\_ini} namelist parameter. -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Chap_DYN.tex
r3764 r3989 1 1 % ================================================================ 2 % Chapter ÑOcean Dynamics (DYN)2 % Chapter � Ocean Dynamics (DYN) 3 3 % ================================================================ 4 4 \chapter{Ocean Dynamics (DYN)} … … 167 167 The vector invariant form of the momentum equations is the one most 168 168 often used in applications of the \NEMO ocean model. The flux form option 169 (see next section) has been present since version $2$. 169 (see next section) has been present since version $2$. Options are defined 170 through the \ngn{namdyn\_adv} namelist variables 170 171 Coriolis and momentum advection terms are evaluated using a leapfrog 171 172 scheme, $i.e.$ the velocity appearing in these expressions is centred in … … 184 185 %------------------------------------------------------------------------------------------------------------- 185 186 187 Options are defined through the \ngn{namdyn\_vor} namelist variables. 186 188 Four discretisations of the vorticity term (\textit{ln\_dynvor\_xxx}=true) are available: 187 189 conserving potential enstrophy of horizontally non-divergent flow (ENS scheme) ; … … 382 384 %------------------------------------------------------------------------------------------------------------- 383 385 386 Options are defined through the \ngn{namdyn\_adv} namelist variables. 384 387 In the flux form (as in the vector invariant form), the Coriolis and momentum 385 388 advection terms are evaluated using a leapfrog scheme, $i.e.$ the velocity … … 526 529 %------------------------------------------------------------------------------------------------------------- 527 530 531 Options are defined through the \ngn{namdyn\_hpg} namelist variables. 528 532 The key distinction between the different algorithms used for the hydrostatic 529 533 pressure gradient is the vertical coordinate used, since HPG is a \emph{horizontal} … … 712 716 713 717 %%% 718 Options are defined through the \ngn{namdyn\_spg} namelist variables. 714 719 The surface pressure gradient term is related to the representation of the free surface (\S\ref{PE_hor_pg}). The main distinction is between the fixed volume case (linear free surface) and the variable volume case (nonlinear free surface, \key{vvl} is defined). In the linear free surface case (\S\ref{PE_free_surface}) the vertical scale factors $e_{3}$ are fixed in time, while they are time-dependent in the nonlinear case (\S\ref{PE_free_surface}). With both linear and nonlinear free surface, external gravity waves are allowed in the equations, which imposes a very small time step when an explicit time stepping is used. Two methods are proposed to allow a longer time step for the three-dimensional equations: the filtered free surface, which is a modification of the continuous equations (see \eqref{Eq_PE_flt}), and the split-explicit free surface described below. The extra term introduced in the filtered method is calculated implicitly, so that the update of the next velocities is done in module \mdl{dynspg\_flt} and not in \mdl{dynnxt}. 715 720 … … 931 936 %------------------------------------------------------------------------------------------------------------- 932 937 938 Options are defined through the \ngn{namdyn\_ldf} namelist variables. 933 939 The options available for lateral diffusion are to use either laplacian 934 940 (rotated or not) or biharmonic operators. The coefficients may be constant … … 1060 1066 %------------------------------------------------------------------------------------------------------------- 1061 1067 1068 Options are defined through the \ngn{namzdf} namelist variables. 1062 1069 The large vertical diffusion coefficient found in the surface mixed layer together 1063 1070 with high vertical resolution implies that in the case of explicit time stepping there … … 1130 1137 %------------------------------------------------------------------------------------------------------------- 1131 1138 1139 Options are defined through the \ngn{namdom} namelist variables. 1132 1140 The general framework for dynamics time stepping is a leap-frog scheme, 1133 1141 $i.e.$ a three level centred time scheme associated with an Asselin time filter -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Chap_LBC.tex
r3294 r3989 1 1 % ================================================================ 2 % Chapter ÑLateral Boundary Condition (LBC)2 % Chapter � Lateral Boundary Condition (LBC) 3 3 % ================================================================ 4 4 \chapter{Lateral Boundary Condition (LBC) } … … 25 25 %OPA allows land and topography grid points in the computational domain due to the presence of continents or islands, and includes the use of a full or partial step representation of bottom topography. The computation is performed over the whole domain, i.e. we do not try to restrict the computation to ocean-only points. This choice has two motivations. Firstly, working on ocean only grid points overloads the code and harms the code readability. Secondly, and more importantly, it drastically reduces the vector portion of the computation, leading to a dramatic increase of CPU time requirement on vector computers. The current section describes how the masking affects the computation of the various terms of the equations with respect to the boundary condition at solid walls. The process of defining which areas are to be masked is described in \S\ref{DOM_msk}. 26 26 27 Options are defined through the \ngn{namlbc} namelist variables. 27 28 The discrete representation of a domain with complex boundaries (coastlines and 28 29 bottom topography) leads to arrays that include large portions where a computation … … 148 149 % Boundary Condition around the Model Domain 149 150 % ================================================================ 150 \section{Model Domain Boundary Condition (\ jp{jperio})}151 \section{Model Domain Boundary Condition (\np{jperio})} 151 152 \label{LBC_jperio} 152 153 … … 156 157 157 158 % ------------------------------------------------------------------------------------------------------------- 158 % Closed, cyclic, south symmetric (\ jp{jperio} = 0, 1 or 2)159 % Closed, cyclic, south symmetric (\np{jperio} = 0, 1 or 2) 159 160 % ------------------------------------------------------------------------------------------------------------- 160 \subsection{Closed, cyclic, south symmetric (\ jp{jperio} = 0, 1 or 2)}161 \subsection{Closed, cyclic, south symmetric (\np{jperio} = 0, 1 or 2)} 161 162 \label{LBC_jperio012} 162 163 163 164 The choice of closed, cyclic or symmetric model domain boundary condition is made 164 by setting \ jp{jperio} to 0, 1 or 2 in file \mdl{par\_oce}. Each time such a boundary165 by setting \np{jperio} to 0, 1 or 2 in namelist \ngn{namcfg}. Each time such a boundary 165 166 condition is needed, it is set by a call to routine \mdl{lbclnk}. The computation of 166 167 momentum and tracer trends proceeds from $i=2$ to $i=jpi-1$ and from $j=2$ to … … 295 296 domain and the overlapping rows. The number of rows to exchange (known as 296 297 the halo) is usually set to one (\jp{jpreci}=1, in \mdl{par\_oce}). The whole domain 297 dimensions are named \ jp{jpiglo}, \jp{jpjglo} and \jp{jpk}. The relationship between298 dimensions are named \np{jpiglo}, \np{jpjglo} and \jp{jpk}. The relationship between 298 299 the whole domain and a sub-domain is: 299 300 \begin{eqnarray} … … 419 420 \end{itemize} 420 421 422 Options are defined through the \ngn{namobc} namelist variables. 421 423 The package resides in the OBC directory. It is described here in four parts: the 422 424 boundary geometry (parameters to be set in \mdl{obc\_par}), the forcing data at … … 455 457 Logical flag & & & \\ 456 458 \hline 457 West & \jp{jpiwob} $>= 2$ & \jp{jpjwd}$>= 2$ & \jp{jpjwf}<= \ jp{jpjglo}-1 \\459 West & \jp{jpiwob} $>= 2$ & \jp{jpjwd}$>= 2$ & \jp{jpjwf}<= \np{jpjglo}-1 \\ 458 460 lp\_obc\_west & $i$-index of a $u$ point & $j$ of a $T$ point &$j$ of a $T$ point \\ 459 461 \hline 460 East & \jp{jpieob}$<=$\ jp{jpiglo}-2&\jp{jpjed} $>= 2$ & \jp{jpjef}$<=$ \jp{jpjglo}-1 \\462 East & \jp{jpieob}$<=$\np{jpiglo}-2&\jp{jpjed} $>= 2$ & \jp{jpjef}$<=$ \np{jpjglo}-1 \\ 461 463 lp\_obc\_east & $i$-index of a $u$ point & $j$ of a $T$ point & $j$ of a $T$ point \\ 462 464 \hline 463 South & \jp{jpjsob} $>= 2$ & \jp{jpisd} $>= 2$ & \jp{jpisf}$<=$\ jp{jpiglo}-1 \\465 South & \jp{jpjsob} $>= 2$ & \jp{jpisd} $>= 2$ & \jp{jpisf}$<=$\np{jpiglo}-1 \\ 464 466 lp\_obc\_south & $j$-index of a $v$ point & $i$ of a $T$ point & $i$ of a $T$ point \\ 465 467 \hline 466 North & \jp{jpjnob} $<=$ \ jp{jpjglo}-2& \jp{jpind} $>= 2$ & \jp{jpinf}$<=$\jp{jpiglo}-1 \\468 North & \jp{jpjnob} $<=$ \np{jpjglo}-2& \jp{jpind} $>= 2$ & \jp{jpinf}$<=$\np{jpiglo}-1 \\ 467 469 lp\_obc\_north & $j$-index of a $v$ point & $i$ of a $T$ point & $i$ of a $T$ point \\ 468 470 \hline … … 754 756 %----------------------------------------------------------------------------------------------- 755 757 758 Options are defined through the \ngn{nambdy} \ngn{nambdy\_index} 759 \ngn{nambdy\_dta} \ngn{nambdy\_dta2} namelist variables. 756 760 The BDY module is an alternative implementation of open boundary 757 761 conditions for regional configurations. It implements the Flow … … 1024 1028 %----------------------------------------------------------------------------------------------- 1025 1029 1026 To be written.... 1027 1028 1029 1030 1030 Options are defined through the \ngn{nambdy\_tide} namelist variables. 1031 To be written.... 1032 1033 1034 1035 -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Chap_LDF.tex
r3294 r3989 1 1 2 2 % ================================================================ 3 % Chapter ÑLateral Ocean Physics (LDF)3 % Chapter � Lateral Ocean Physics (LDF) 4 4 % ================================================================ 5 5 \chapter{Lateral Ocean Physics (LDF)} … … 21 21 and for tracers only, eddy induced advection on tracers). These three aspects 22 22 of the lateral diffusion are set through namelist parameters and CPP keys 23 (see the \textit{ nam\_traldf} and \textit{nam\_dynldf} below). Note23 (see the \textit{\ngn{nam\_traldf}} and \textit{\ngn{nam\_dynldf}} below). Note 24 24 that this chapter describes the default implementation of iso-neutral 25 25 tracer mixing, and Griffies's implementation, which is used if … … 104 104 105 105 Other formulations can be introduced by the user for a given configuration. 106 For example, in the ORCA2 global ocean model ( \key{orca\_r2}), the laplacian106 For example, in the ORCA2 global ocean model (see Configurations), the laplacian 107 107 viscosity operator uses \np{rn\_ahm0}~= 4.10$^4$ m$^2$/s poleward of 20$^{\circ}$ 108 108 north and south and decreases linearly to \np{rn\_aht0}~= 2.10$^3$ m$^2$/s … … 110 110 can be found in routine \rou{ldf\_dyn\_c2d\_orca} defined in \mdl{ldfdyn\_c2d}. 111 111 Similar modified horizontal variations can be found with the Antarctic or Arctic 112 sub-domain options of ORCA2 and ORCA05 (\key{antarctic} or \key{arctic} 113 defined, see \hf{ldfdyn\_antarctic} and \hf{ldfdyn\_arctic}). 112 sub-domain options of ORCA2 and ORCA05 (see \&namcfg namelist). 114 113 115 114 \subsubsection{Space Varying Mixing Coefficients (\key{traldf\_c3d} and \key{dynldf\_c3d})} … … 123 122 There is no default specification of space and time varying mixing coefficient. 124 123 The only case available is specific to the ORCA2 and ORCA05 global ocean 125 configurations (\key{orca\_r2} or \key{orca\_r05}). It provides only a tracer124 configurations. It provides only a tracer 126 125 mixing coefficient for eddy induced velocity (ORCA2) or both iso-neutral and 127 126 eddy induced velocity (ORCA05) that depends on the local growth rate of -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Chap_MISC.tex
r3294 r3989 1 1 % ================================================================ 2 % Chapter ÑMiscellaneous Topics2 % Chapter � Miscellaneous Topics 3 3 % ================================================================ 4 4 \chapter{Miscellaneous Topics} … … 33 33 Note that such modifications are so specific to a given configuration that no attempt 34 34 has been made to set them in a generic way. However, examples of how 35 they can be set up is given in the ORCA 2\deg and 0.5\deg configurations (search for 36 \key{orca\_r2} or \key{orca\_r05} in the code). 35 they can be set up is given in the ORCA 2\deg and 0.5\deg configurations. For example, 36 for details of implementation in ORCA2, search: 37 \vspace{-10pt} 38 \begin{alltt} 39 \tiny 40 \begin{verbatim} 41 IF( cp_cfg == "orca" .AND. jp_cfg == 2 ) 42 \end{verbatim} 43 \end{alltt} 37 44 38 45 % ------------------------------------------------------------------------------------------------------------- … … 83 90 84 91 \colorbox{yellow}{Add a short description of CLA staff here or in lateral boundary condition chapter?} 92 Options are defined through the \ngn{namcla} namelist variables. 85 93 86 94 %The problem is resolved here by allowing the mixing of tracers and mass/volume between non-adjacent water columns at nominated regions within the model. Momentum is not mixed. The scheme conserves total tracer content, and total volume (the latter in $z*$- or $s*$-coordinate), and maintains compatibility between the tracer and mass/volume budgets. … … 98 106 % Sub-Domain Functionality (\textit{nizoom, njzoom}, namelist parameters) 99 107 % ================================================================ 100 \section{Sub-Domain Functionality (\ jp{jpizoom}, \jp{jpjzoom})}108 \section{Sub-Domain Functionality (\np{jpizoom}, \np{jpjzoom})} 101 109 \label{MISC_zoom} 102 110 … … 119 127 In order to easily define a sub-domain over which the computation can be 120 128 performed, the dimension of all input arrays (ocean mesh, bathymetry, 121 forcing, initial state, ...) are defined as \ jp{jpidta}, \jp{jpjdta} and \jp{jpkdta}122 ( \mdl{par\_oce} module), while the computational domain is defined through123 \ jp{jpiglo}, \jp{jpjglo} and \jp{jpk} (\mdl{par\_oce} module). When running the124 model over the whole domain, the user sets \ jp{jpiglo}=\jp{jpidta} \jp{jpjglo}=\jp{jpjdta}129 forcing, initial state, ...) are defined as \np{jpidta}, \np{jpjdta} and \np{jpkdta} 130 ( in \ngn{namcfg} namelist), while the computational domain is defined through 131 \np{jpiglo}, \np{jpjglo} and \jp{jpk} (\ngn{namcfg} namelist). When running the 132 model over the whole domain, the user sets \np{jpiglo}=\np{jpidta} \np{jpjglo}=\np{jpjdta} 125 133 and \jp{jpk}=\jp{jpkdta}. When running the model over a sub-domain, the user 126 has to provide the size of the sub-domain, (\ jp{jpiglo}, \jp{jpjglo}, \jp{jpkglo}),127 and the indices of the south western corner as \ jp{jpizoom} and \jp{jpjzoom} in128 the \mdl{par\_oce} module(Fig.~\ref{Fig_LBC_zoom}).134 has to provide the size of the sub-domain, (\np{jpiglo}, \np{jpjglo}, \np{jpkglo}), 135 and the indices of the south western corner as \np{jpizoom} and \np{jpjzoom} in 136 the \ngn{namcfg} namelist (Fig.~\ref{Fig_LBC_zoom}). 129 137 130 138 Note that a third set of dimensions exist, \jp{jpi}, \jp{jpj} and \jp{jpk} which is 131 actually used to perform the computation. It is set by default to \jp{jpi}=\ jp{jpjglo}132 and \jp{jpj}=\ jp{jpjglo}, except for massively parallel computing where the139 actually used to perform the computation. It is set by default to \jp{jpi}=\np{jpjglo} 140 and \jp{jpj}=\np{jpjglo}, except for massively parallel computing where the 133 141 computational domain is laid out on local processor memories following a 2D 134 142 horizontal splitting. % (see {\S}IV.2-c) ref to the section to be updated … … 162 170 trajectory to reach it. 163 171 172 Options are defined through the \ngn{namdom} namelist variables. 164 173 The acceleration of convergence option is used when \np{nn\_acc}=1. In that case, 165 174 $\rdt=rn\_rdt$ is the time step of dynamics while $\widetilde{\rdt}=rdttra$ is the … … 284 293 285 294 \gmcomment{why not make these bullets into subsections?} 286 295 Options are defined through the \ngn{namctl} namelist variables. 287 296 288 297 $\bullet$ Vector optimisation: … … 343 352 a Successive-Over-Relaxation scheme (SOR) and a preconditioned conjugate gradient 344 353 scheme(PCG) \citep{Madec_al_OM88, Madec_PhD90}. The solver is selected trough the 345 the value of \np{nn\_solv} (namelist parameter).354 the value of \np{nn\_solv} \ngn{namsol} namelist variable. 346 355 347 356 The PCG is a very efficient method for solving elliptic equations on vector computers. -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Chap_Model_Basics_zstar.tex
r3294 r3989 1 1 % ================================================================ 2 % Chapter 1 ÑModel Basics2 % Chapter 1 � Model Basics 3 3 % ================================================================ 4 4 % ================================================================ … … 53 53 Because $z^\star$ has a time independent range, all grid cells have static increments 54 54 ds, and the sum of the ver tical increments yields the time independent ocean 55 depth % ·k ds = H.55 depth %�k ds = H. 56 56 The $z^\star$ coordinate is therefore invisible to undulations of the 57 57 free surface, since it moves along with the free surface. This proper ty means that … … 78 78 \namdisplay{nam_dynspg} 79 79 %------------------------------------------------------------------------------------------------------------ 80 Options are defined through the \ngn{nam\_dynspg} namelist variables. 80 81 The surface pressure gradient term is related to the representation of the free surface (\S\ref{PE_hor_pg}). The main distinction is between the fixed volume case (linear free surface or rigid lid) and the variable volume case (nonlinear free surface, \key{vvl} is active). In the linear free surface case (\S\ref{PE_free_surface}) and rigid lid (\S\ref{PE_rigid_lid}), the vertical scale factors $e_{3}$ are fixed in time, while in the nonlinear case (\S\ref{PE_free_surface}) they are time-dependent. With both linear and nonlinear free surface, external gravity waves are allowed in the equations, which imposes a very small time step when an explicit time stepping is used. Two methods are proposed to allow a longer time step for the three-dimensional equations: the filtered free surface, which is a modification of the continuous equations (see \eqref{Eq_PE_flt}), and the split-explicit free surface described below. The extra term introduced in the filtered method is calculated implicitly, so that the update of the next velocities is done in module \mdl{dynspg\_flt} and not in \mdl{dynnxt}. 81 82 … … 114 115 \namdisplay{namdom} 115 116 %-------------------------------------------------------------------------------------------------------------- 116 The split-explicit free surface formulation used in OPA follows the one proposed by \citet{Griffies2004}. The general idea is to solve the free surface equation with a small time step, while the three dimensional prognostic variables are solved with a longer time step that is a multiple of \np{rdtbt} (Figure III.3). 117 The split-explicit free surface formulation used in OPA follows the one proposed by \citet{Griffies2004}. The general idea is to solve the free surface equation with a small time step, while the three dimensional prognostic variables are solved with a longer time step that is a multiple of \np{rdtbt} 118 in the \ngn{namdom} namelist. 119 (Figure III.3). 117 120 118 121 %> > > > > > > > > > > > > > > > > > > > > > > > > > > > -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Chap_OBS.tex
r3294 r3989 20 20 can be used for validation or verification of model or any other data assimilation system. 21 21 22 The OBS code is called from \ np{opa.F90} for model initialisation and to calculate the model22 The OBS code is called from \mdl{nemogcm.F90} for model initialisation and to calculate the model 23 23 equivalent values for observations on the 0th timestep. The code is then called again after 24 each timestep from \ np{step.F90}. To build with the OBS code active \key{diaobs} must be24 each timestep from \mdl{step.F90}. To build with the OBS code active \key{diaobs} must be 25 25 set. 26 26 … … 66 66 67 67 \item Add the following to the NEMO namelist to run the observation 68 operator on this data. Set the \np{enactfiles} namelist parameterto the68 operator on this data. Set the \np{enactfiles} namelist variable to the 69 69 observation file name: 70 70 \end{enumerate} … … 74 74 %------------------------------------------------------------------------------------------------------------- 75 75 76 Options are defined through the \ngn{namobs} namelist variables. 76 77 The options \np{ln\_t3d} and \np{ln\_s3d} switch on the temperature and salinity 77 78 profile observation operator code. The \np{ln\_ena} switch turns on the reading … … 94 95 \label{OBS_details} 95 96 96 Here we show a more complete example namelist and also show the NetCDF headers97 Here we show a more complete example namelist \ngn{namobs} and also show the NetCDF headers 97 98 of the observation 98 99 files that may be used with the observation operator -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Chap_SBC.tex
r3795 r3989 1 1 % ================================================================ 2 % Chapter ÑSurface Boundary Condition (SBC, ICB)2 % Chapter � Surface Boundary Condition (SBC, ICB) 3 3 % ================================================================ 4 4 \chapter{Surface Boundary Condition (SBC, ICB) } … … 25 25 26 26 Five different ways to provide the first six fields to the ocean are available which 27 are controlled by namelist variables: an analytical formulation (\np{ln\_ana}~=~true),27 are controlled by namelist \ngn{namsbc} variables: an analytical formulation (\np{ln\_ana}~=~true), 28 28 a flux formulation (\np{ln\_flx}~=~true), a bulk formulae formulation (CORE 29 29 (\np{ln\_core}~=~true), CLIO (\np{ln\_clio}~=~true) or MFS … … 442 442 %-------------------------------------------------------------------------------------------------------------- 443 443 444 In some circumstances it may be useful to avoid calculating the 3D temperature, salinity and velocity fields and 445 simply read them in from a previous run. For example: 444 In some circumstances it may be useful to avoid calculating the 3D temperature, salinity and velocity fields and simply read them in from a previous run. 445 Options are defined through the \ngn{namsbc\_sas} namelist variables. 446 For example: 446 447 447 448 \begin{enumerate} … … 507 508 In this case, all the six fluxes needed by the ocean are assumed to 508 509 be uniform in space. They take constant values given in the namelist 509 namsbc{\_}anaby the variables \np{rn\_utau0}, \np{rn\_vtau0}, \np{rn\_qns0},510 \ngn{namsbc{\_}ana} by the variables \np{rn\_utau0}, \np{rn\_vtau0}, \np{rn\_qns0}, 510 511 \np{rn\_qsr0}, and \np{rn\_emp0} ($\textit{emp}=\textit{emp}_S$). The runoff is set to zero. 511 512 In addition, the wind is allowed to reach its nominal value within a given number … … 530 531 In the flux formulation (\np{ln\_flx}=true), the surface boundary 531 532 condition fields are directly read from input files. The user has to define 532 in the namelist namsbc{\_}flxthe name of the file, the name of the variable533 in the namelist \ngn{namsbc{\_}flx} the name of the file, the name of the variable 533 534 read in the file, the time frequency at which it is given (in hours), and a logical 534 535 setting whether a time interpolation to the model time step is required … … 580 581 This is the so-called DRAKKAR Forcing Set (DFS) \citep{Brodeau_al_OM09}. 581 582 583 Options are defined through the \ngn{namsbc\_core} namelist variables. 582 584 The required 8 input fields are: 583 585 … … 621 623 compute the radiative fluxes from a climatological cloud cover. 622 624 625 Options are defined through the \ngn{namsbc\_clio} namelist variables. 623 626 The required 7 input fields are: 624 627 … … 673 676 Details on the bulk formulae used can be found in \citet{Maggiore_al_PCE98} and \citet{Castellari_al_JMS1998}. 674 677 678 Options are defined through the \ngn{namsbc\_mfs} namelist variables. 675 679 The required 7 input fields must be provided on the model Grid-T and are: 676 680 \begin{itemize} … … 711 715 When PISCES biogeochemical model (\key{top} and \key{pisces}) is also used in the coupled system, 712 716 the whole carbon cycle is computed by defining \key{cpl\_carbon\_cycle}. In this case, 713 CO$_2$ fluxes will be exchanged between the atmosphere and the ice-ocean system (and need to be activated 714 in namsbc{\_}cpl). 715 716 The new namelist above allows control of various aspects of the coupling fields (particularly for 717 CO$_2$ fluxes will be exchanged between the atmosphere and the ice-ocean system (and need to be activated in \ngn{namsbc{\_}cpl} ). 718 719 The namelist above allows control of various aspects of the coupling fields (particularly for 717 720 vectors) and now allows for any coupling fields to have multiple sea ice categories (as required by LIM3 718 721 and CICE). When indicating a multi-category coupling field in namsbc{\_}cpl the number of categories will be … … 736 739 737 740 The optional atmospheric pressure can be used to force ocean and ice dynamics 738 (\np{ln\_apr\_dyn}~=~true, \textit{ namsbc} namelist ).741 (\np{ln\_apr\_dyn}~=~true, \textit{\ngn{namsbc}} namelist ). 739 742 The input atmospheric forcing defined via \np{sn\_apr} structure (\textit{namsbc\_apr} namelist) 740 743 can be interpolated in time to the model time step, and even in space when the … … 774 777 %------------------------------------------------------------------------------------------------------------- 775 778 776 Concerning the tidal potential, some parameters are available in namelist :777 778 - \ texttt{ln\_tide\_pot} activate the tidal potential forcing779 780 - \ texttt{nb\_harmo} is the number of constituent used781 782 - \ texttt{clname} is the name of constituent779 Concerning the tidal potential, some parameters are available in namelist \ngn{nam\_tide}: 780 781 - \np{ln\_tide\_pot} activate the tidal potential forcing 782 783 - \np{nb\_harmo} is the number of constituent used 784 785 - \np{clname} is the name of constituent 783 786 784 787 … … 858 861 depth (in metres) which the river should be added to. 859 862 860 Namelist options, \np{ln\_rnf\_depth}, \np{ln\_rnf\_sal} and \np{ln\_rnf\_temp} control whether863 Namelist variables in \ngn{namsbc\_rnf}, \np{ln\_rnf\_depth}, \np{ln\_rnf\_sal} and \np{ln\_rnf\_temp} control whether 861 864 the river attributes (depth, salinity and temperature) are read in and used. If these are set 862 865 as false the river is added to the surface box only, assumed to be fresh (0~psu), and/or … … 943 946 Their physical behaviour is controlled by equations as described in \citet{Martin_Adcroft_OM10} ). 944 947 (Note that the authors kindly provided a copy of their code to act as a basis for implementation in NEMO.) 945 Icebergs are initially spawned into one of ten classes which have specific mass and thickness as described by948 Icebergs are initially spawned into one of ten classes which have specific mass and thickness as described in the \ngn{namberg} namelist: 946 949 \np{rn\_initial\_mass} and \np{rn\_initial\_thickness}. 947 950 Each class has an associated scaling (\np{rn\_mass\_scaling}), which is an integer representing how many icebergs … … 1031 1034 the diurnal cycle of SWF is a scaling of the top of the atmosphere diurnal cycle 1032 1035 of incident SWF. The \cite{Bernie_al_CD07} reconstruction algorithm is available 1033 in \NEMO by setting \np{ln\_dm2dc}~=~true (a \textit{ namsbc} namelist parameter) when using1036 in \NEMO by setting \np{ln\_dm2dc}~=~true (a \textit{\ngn{namsbc}} namelist variable) when using 1034 1037 CORE bulk formulea (\np{ln\_blk\_core}~=~true) or the flux formulation (\np{ln\_flx}~=~true). 1035 1038 The reconstruction is performed in the \mdl{sbcdcy} module. The detail of the algoritm used … … 1088 1091 %------------------------------------------------------------------------------------------------------------- 1089 1092 1090 In forced mode using a flux formulation (\np{ln\_flx}~=~true), a 1093 IOptions are defined through the \ngn{namsbc\_ssr} namelist variables. 1094 n forced mode using a flux formulation (\np{ln\_flx}~=~true), a 1091 1095 feedback term \emph{must} be added to the surface heat flux $Q_{ns}^o$: 1092 1096 \begin{equation} \label{Eq_sbc_dmp_q} … … 1212 1216 in $namsbc$ namelist must be defined ${.true.}$. 1213 1217 The \mdl{sbcwave} module containing the routine \np{sbc\_wave} reads the 1214 namelist ${namsbc\_wave}$(for external data names, locations, frequency, interpolation and all1218 namelist \ngn{namsbc\_wave} (for external data names, locations, frequency, interpolation and all 1215 1219 the miscellanous options allowed by Input Data generic Interface see \S\ref{SBC_input}) 1216 1220 and a 2D field of neutral drag coefficient. Then using the routine … … 1222 1226 % Griffies doc: 1223 1227 % When running ocean-ice simulations, we are not explicitly representing land processes, such as rivers, catchment areas, snow accumulation, etc. However, to reduce model drift, it is important to balance the hydrological cycle in ocean-ice models. We thus need to prescribe some form of global normalization to the precipitation minus evaporation plus river runoff. The result of the normalization should be a global integrated zero net water input to the ocean-ice system over a chosen time scale. 1224 %How often the normalization is done is a matter of choice. In mom4p1, we choose to do so at each model time step, so that there is always a zero net input of water to the ocean-ice system. Others choose to normalize over an annual cycle, in which case the net imbalance over an annual cycle is used to alter the subsequent year Õs water budget in an attempt to damp the annual water imbalance. Note that the annual budget approach may be inappropriate with interannually varying precipitation forcing.1228 %How often the normalization is done is a matter of choice. In mom4p1, we choose to do so at each model time step, so that there is always a zero net input of water to the ocean-ice system. Others choose to normalize over an annual cycle, in which case the net imbalance over an annual cycle is used to alter the subsequent year�s water budget in an attempt to damp the annual water imbalance. Note that the annual budget approach may be inappropriate with interannually varying precipitation forcing. 1225 1229 %When running ocean-ice coupled models, it is incorrect to include the water transport between the ocean and ice models when aiming to balance the hydrological cycle. The reason is that it is the sum of the water in the ocean plus ice that should be balanced when running ocean-ice models, not the water in any one sub-component. As an extreme example to illustrate the issue, consider an ocean-ice model with zero initial sea ice. As the ocean-ice model spins up, there should be a net accumulation of water in the growing sea ice, and thus a net loss of water from the ocean. The total water contained in the ocean plus ice system is constant, but there is an exchange of water between the subcomponents. This exchange should not be part of the normalization used to balance the hydrological cycle in ocean-ice models. 1226 1230 -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Chap_STP.tex
r3294 r3989 1 1 2 2 % ================================================================ 3 % Chapter 2 ÑTime Domain (step.F90)3 % Chapter 2 � Time Domain (step.F90) 4 4 % ================================================================ 5 5 \chapter{Time Domain (STP) } … … 335 335 When restarting, if the the time step has been changed, a restart using an Euler time 336 336 stepping scheme is imposed. 337 Options are defined through the \ngn{namrun} namelist variables. 337 338 %%% 338 339 \gmcomment{ … … 358 359 %-------------------------------------------------------------------------------------------------------------- 359 360 360 361 Options are defined through the \ngn{namdom} namelist variables. 361 362 \colorbox{yellow}{add here a few word on nit000 and nitend} 362 363 -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Chap_TRA.tex
r3764 r3989 1 1 % ================================================================ 2 % Chapter 1 ÑOcean Tracers (TRA)2 % Chapter 1 � Ocean Tracers (TRA) 3 3 % ================================================================ 4 4 \chapter{Ocean Tracers (TRA)} … … 137 137 \textit{effective} velocity ($i.e.$ the sum of the eulerian and eiv velocities) which is used. 138 138 139 The choice of an advection scheme is made in the \textit{ nam\_traadv} namelist, by139 The choice of an advection scheme is made in the \textit{\ngn{nam\_traadv}} namelist, by 140 140 setting to \textit{true} one and only one of the logicals \textit{ln\_traadv\_xxx}. The 141 141 corresponding code can be found in the \textit{traadv\_xxx.F90} module, where … … 441 441 %------------------------------------------------------------------------------------------------------------- 442 442 443 Options are defined through the \ngn{namtra\_ldf} namelist variables. 443 444 The options available for lateral diffusion are a laplacian (rotated or not) 444 445 or a biharmonic operator, the latter being more scale-selective (more … … 602 603 %-------------------------------------------------------------------------------------------------------------- 603 604 605 Options are defined through the \ngn{namzdf} namelist variables. 604 606 The formulation of the vertical subgrid scale tracer physics is the same 605 607 for all the vertical coordinates, and is based on a laplacian operator. … … 757 759 %-------------------------------------------------------------------------------------------------------------- 758 760 761 Options are defined through the \ngn{namtra\_qsr} namelist variables. 759 762 When the penetrative solar radiation option is used (\np{ln\_flxqsr}=true), 760 763 the solar radiation penetrates the top few tens of meters of the ocean. If it is not used … … 879 882 Bottom Water) by a few Sverdrups \citep{Emile-Geay_Madec_OS09}. 880 883 884 Options are defined through the \ngn{namtra\_bbc} namelist variables. 881 885 The presence of geothermal heating is controlled by setting the namelist 882 886 parameter \np{ln\_trabbc} to true. Then, when \np{nn\_geoflx} is set to 1, … … 897 901 %-------------------------------------------------------------------------------------------------------------- 898 902 903 Options are defined through the \ngn{nambbl} namelist variables. 899 904 In a $z$-coordinate configuration, the bottom topography is represented by a 900 905 series of discrete steps. This is not adequate to represent gravity driven … … 1066 1071 where $\gamma$ is the inverse of a time scale, and $T_o$ and $S_o$ 1067 1072 are given temperature and salinity fields (usually a climatology). 1073 Options are defined through the \ngn{namtra\_dmp} namelist variables. 1068 1074 The restoring term is added when the namelist parameter \np{ln\_tradmp} is set to true. 1069 1075 It also requires that both \np{ln\_tsd\_init} and \np{ln\_tsd\_tradmp} are set to true … … 1128 1134 %-------------------------------------------------------------------------------------------------------------- 1129 1135 1136 Options are defined through the \ngn{namdom} namelist variables. 1130 1137 The general framework for tracer time stepping is a modified leap-frog scheme 1131 1138 \citep{Leclair_Madec_OM09}, $i.e.$ a three level centred time scheme associated … … 1205 1212 \citep{Gill1982}. 1206 1213 1214 Options are defined through the \ngn{nameos} namelist variables. 1207 1215 The default option (namelist parameter \np{nn\_eos}=0) is the \citet{JackMcD1995} 1208 1216 equation of state. Its use is highly recommended. However, for process studies, … … 1222 1230 coefficients, and $\rho_o$, the reference volumic mass, $rau0$. 1223 1231 ($\alpha$ and $\beta$ can be modified through the \np{rn\_alpha} and 1224 \np{rn\_beta} namelist parameters). Note that when $d_a$ is a function1232 \np{rn\_beta} namelist variables). Note that when $d_a$ is a function 1225 1233 of $T$ only (\np{nn\_eos}=1), the salinity is a passive tracer and can be 1226 1234 used as such. -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Chap_ZDF.tex
r3764 r3989 53 53 %-------------------------------------------------------------------------------------------------------------- 54 54 55 Options are defined through the \ngn{namzdf} namelist variables. 55 56 When \key{zdfcst} is defined, the momentum and tracer vertical eddy coefficients 56 57 are set to constant values over the whole ocean. This is the crudest way to define … … 79 80 80 81 When \key{zdfric} is defined, a local Richardson number dependent formulation 81 for the vertical momentum and tracer eddy coefficients is set. The vertical mixing 82 for the vertical momentum and tracer eddy coefficients is set through the \ngn{namzdf\_ric} 83 namelist variables.The vertical mixing 82 84 coefficients are diagnosed from the large scale variables computed by the model. 83 85 \textit{In situ} measurements have been used to link vertical turbulent activity to … … 176 178 \end{cases} 177 179 \end{align*} 178 The choice of $P_{rt}$ is controlled by the \np{nn\_pdl} namelist parameter. 180 Options are defined through the \ngn{namzdfy\_tke} namelist variables. 181 The choice of $P_{rt}$ is controlled by the \np{nn\_pdl} namelist variable. 179 182 180 183 At the sea surface, the value of $\bar{e}$ is prescribed from the wind … … 539 542 \caption{ \label{Tab_GLS} 540 543 Set of predefined GLS parameters, or equivalently predefined turbulence models available 541 with \key{zdfgls} and controlled by the \np{nn\_clos} namelist parameter.}544 with \key{zdfgls} and controlled by the \np{nn\_clos} namelist variable in \ngn{namzdf\_gls} .} 542 545 \end{center} \end{table} 543 546 %-------------------------------------------------------------------------------------------------------------- … … 581 584 582 585 The KKP scheme has been implemented by J. Chanut ... 586 Options are defined through the \ngn{namzdf\_kpp} namelist variables. 583 587 584 588 \colorbox{yellow}{Add a description of KPP here.} … … 631 635 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 632 636 637 Options are defined through the \ngn{namzdf} namelist variables. 633 638 The non-penetrative convective adjustment is used when \np{ln\_zdfnpc}=true. 634 639 It is applied at each \np{nn\_npc} time step and mixes downwards instantaneously … … 691 696 %-------------------------------------------------------------------------------------------------------------- 692 697 698 Options are defined through the \ngn{namzdf} namelist variables. 693 699 The enhanced vertical diffusion parameterisation is used when \np{ln\_zdfevd}=true. 694 700 In this case, the vertical eddy mixing coefficients are assigned very large values … … 749 755 %-------------------------------------------------------------------------------------------------------------- 750 756 757 Options are defined through the \ngn{namzdf\_ddm} namelist variables. 751 758 Double diffusion occurs when relatively warm, salty water overlies cooler, fresher 752 759 water, or vice versa. The former condition leads to salt fingering and the latter … … 830 837 %-------------------------------------------------------------------------------------------------------------- 831 838 839 Options are defined through the \ngn{nambfr} namelist variables. 832 840 Both the surface momentum flux (wind stress) and the bottom momentum 833 841 flux (bottom friction) enter the equations as a condition on the vertical … … 1136 1144 \label{ZDF_tmx_bottom} 1137 1145 1146 Options are defined through the \ngn{namzdf\_tmx} namelist variables. 1138 1147 The parameterization of tidal mixing follows the general formulation for 1139 1148 the vertical eddy diffusivity proposed by \citet{St_Laurent_al_GRL02} and -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Introduction.tex
r3625 r3989 81 81 the coefficients with \citet{Blanke1993}, \citet{Large_al_RG94}, \citet{Pacanowski_Philander_JPO81}, 82 82 or \citet{Umlauf_Burchard_JMS03} mixing schemes. 83 \vspace{1cm} 84 85 86 \noindent CPP keys and namelists are used for inputs to the code. \newline 87 88 \noindent \index{CPP keys} CPP keys \newline 89 Some CPP keys are implemented in the FORTRAN code to allow code selection at compiling step. This selection of code at compilation time reduces the reliability of the whole platform since it changes the code from one set of CPP keys to the other. It is used only when the addition/suppression of the part of code highly changes the amount of memory at run time. 90 Usual coding looks like : 91 \vspace{-10pt} 92 \begin{alltt} 93 \tiny 94 \begin{verbatim} 95 #if defined key_option1 96 This part of the FORTRAN code will be active 97 only if key_option1 is activated at compiling step 98 #endif 99 \end{verbatim} 100 \end{alltt} 101 102 103 \noindent \index{Namelist} Namelists 104 105 The namelist allows to input variables (character, logical, real and integer) into the code. There is one namelist file for each component of NEMO (dynamics, sea-ice, biogeochemistry...) containing all the FOTRAN namelists needed. The implementation in NEMO uses a two step process. For each FORTRAN namelist, two files are read: 106 \begin{enumerate} 107 \item A reference namelist ( in \textit{CONFIG/SHARED/namelist\_ref} ) is read first. This file contains all the namelist variables which are initialised to default values 108 \item A configuration namelist ( in \textit{CONFIG/CFG\_NAME/EXP00/namelist\_cfg} ) is read aferwards. This file contains only the namelist variables which are changed from default values, and overwrites those. 109 \end{enumerate} 110 A template can be found in \textit{NEMO/OPA\_SRC/module.example} 111 The effective namelist, taken in account during the run, is stored at execution time in an output\_namelist\_dyn (or \_ice or \_top) file. 112 \vspace{1cm} 113 83 114 84 115 Model outputs management and specific online diagnostics are described in chapters~\ref{DIA}. -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/nam_tide
r3294 r3989 1 1 !----------------------------------------------------------------------- 2 ! nam_tidetide parameters (#ifdef key_tide)2 &nam_tide ! tide parameters (#ifdef key_tide) 3 3 !----------------------------------------------------------------------- 4 ! ln_tide_pot = use tidal potential forcing 5 ! nb_harmo = number of constituents used 6 ! name(1) = 'M2', 'K1', etc name of constituent 7 8 &nam_tide 9 ln_tide_pot = .true. 10 nb_harmo = 11 11 clname(1) = 'M2' 4 ln_tide_pot = .true. ! use tidal potential forcing 5 clname(1) = 'M2' ! name of constituent 12 6 clname(2) = 'S2' 13 7 clname(3) = 'N2' -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namasm
r3764 r3989 2 2 &nam_asminc ! assimilation increments ('key_asminc') 3 3 !----------------------------------------------------------------------- 4 ln_bkgwri = .false. ! Logical switch for writing out background state 4 ln_bkgwri = .false. ! Logical switch for writing out background state 5 5 ln_trainc = .false. ! Logical switch for applying tracer increments 6 6 ln_dyninc = .false. ! Logical switch for applying velocity increments 7 ln_sshinc = .false. ! Logical switch for applying SSH increments 7 ln_sshinc = .false. ! Logical switch for applying SSH increments 8 8 ln_asmdin = .false. ! Logical switch for Direct Initialization (DI) 9 9 ln_asmiau = .false. ! Logical switch for Incremental Analysis Updating (IAU) … … 15 15 ln_salfix = .false. ! Logical switch for ensuring that the sa > salfixmin 16 16 salfixmin = -9999 ! Minimum salinity after applying the increments 17 n divdmp= 0 ! Number of iterations of divergence damping operator17 nn_divdmp = 0 ! Number of iterations of divergence damping operator 18 18 / -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/nambdy
r3294 r3989 2 2 &nambdy ! unstructured open boundaries ("key_bdy") 3 3 !----------------------------------------------------------------------- 4 nb_bdy = 2 ! number of open boundary sets5 ln_coords_file = .true. ,.false.! =T : read bdy coordinates from file6 cn_coords_file = 'coordinates.bdy.nc' ,''! bdy coordinates files7 ln_mask_file = .false. 8 cn_mask_file = '' 9 nn_dyn2d = 2 , 0! boundary conditions for barotropic fields10 nn_dyn2d_dta = 3 , 0! = 0, bdy data are equal to the initial state11 12 13 14 nn_dyn3d = 0 , 0! boundary conditions for baroclinic velocities15 nn_dyn3d_dta = 0 , 0! = 0, bdy data are equal to the initial state16 17 nn_tra = 1 , 1! boundary conditions for T and S18 nn_tra_dta = 1 , 1! = 0, bdy data are equal to the initial state19 20 nn_rimwidth = 10 , 5! width of the relaxation zone21 ln_vol = .false. 22 nn_volctl = 1 4 nb_bdy = 1 ! number of open boundary sets 5 ln_coords_file = .true. ! =T : read bdy coordinates from file 6 cn_coords_file = 'coordinates.bdy.nc' ! bdy coordinates files 7 ln_mask_file = .false. ! =T : read mask from file 8 cn_mask_file = '' ! name of mask file (if ln_mask_file=.TRUE.) 9 nn_dyn2d = 2 ! boundary conditions for barotropic fields 10 nn_dyn2d_dta = 3 ! = 0, bdy data are equal to the initial state 11 ! = 1, bdy data are read in 'bdydata .nc' files 12 ! = 2, use tidal harmonic forcing data from files 13 ! = 3, use external data AND tidal harmonic forcing 14 nn_dyn3d = 0 ! boundary conditions for baroclinic velocities 15 nn_dyn3d_dta = 0 ! = 0, bdy data are equal to the initial state 16 ! = 1, bdy data are read in 'bdydata .nc' files 17 nn_tra = 1 ! boundary conditions for T and S 18 nn_tra_dta = 1 ! = 0, bdy data are equal to the initial state 19 ! = 1, bdy data are read in 'bdydata .nc' files 20 nn_rimwidth = 10 ! width of the relaxation zone 21 ln_vol = .false. ! total volume correction (see nn_volctl parameter) 22 nn_volctl = 1 ! = 0, the total water flux across open boundaries is zero 23 23 / -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/nambdy_tide
r3294 r3989 1 1 !----------------------------------------------------------------------- 2 &nambdy_tide ! tidal forcing at open boundaries 2 &nambdy_tide ! tidal forcing at open boundaries 3 3 !----------------------------------------------------------------------- 4 filtide = 'bdydta/amm12_bdytide_' ! file name root of tidal forcing files 5 tide_cpt(1) ='Q1' ! names of tidal components used 6 tide_cpt(2) ='O1' ! names of tidal components used 7 tide_cpt(3) ='P1' ! names of tidal components used 8 tide_cpt(4) ='S1' ! names of tidal components used 9 tide_cpt(5) ='K1' ! names of tidal components used 10 tide_cpt(6) ='2N2' ! names of tidal components used 11 tide_cpt(7) ='MU2' ! names of tidal components used 12 tide_cpt(8) ='N2' ! names of tidal components used 13 tide_cpt(9) ='NU2' ! names of tidal components used 14 tide_cpt(10) ='M2' ! names of tidal components used 15 tide_cpt(11) ='L2' ! names of tidal components used 16 tide_cpt(12) ='T2' ! names of tidal components used 17 tide_cpt(13) ='S2' ! names of tidal components used 18 tide_cpt(14) ='K2' ! names of tidal components used 19 tide_cpt(15) ='M4' ! names of tidal components used 20 tide_speed(1) = 13.398661 ! phase speeds of tidal components (deg/hour) 21 tide_speed(2) = 13.943036 ! phase speeds of tidal components (deg/hour) 22 tide_speed(3) = 14.958932 ! phase speeds of tidal components (deg/hour) 23 tide_speed(4) = 15.000001 ! phase speeds of tidal components (deg/hour) 24 tide_speed(5) = 15.041069 ! phase speeds of tidal components (deg/hour) 25 tide_speed(6) = 27.895355 ! phase speeds of tidal components (deg/hour) 26 tide_speed(7) = 27.968210 ! phase speeds of tidal components (deg/hour) 27 tide_speed(8) = 28.439730 ! phase speeds of tidal components (deg/hour) 28 tide_speed(9) = 28.512585 ! phase speeds of tidal components (deg/hour) 29 tide_speed(10) = 28.984106 ! phase speeds of tidal components (deg/hour) 30 tide_speed(11) = 29.528479 ! phase speeds of tidal components (deg/hour) 31 tide_speed(12) = 29.958935 ! phase speeds of tidal components (deg/hour) 32 tide_speed(13) = 30.000002 ! phase speeds of tidal components (deg/hour) 33 tide_speed(14) = 30.082138 ! phase speeds of tidal components (deg/hour) 34 tide_speed(15) = 57.968212 ! phase speeds of tidal components (deg/hour) 35 ln_tide_date = .true. ! adjust tidal harmonics for start date of run 4 filtide = 'bdydta/amm12_bdytide_' ! file name root of tidal forcing files 5 ln_bdytide_2ddta = .false. 6 ln_bdytide_conj = .false. 36 7 / -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namberg
r3609 r3989 2 2 &namberg ! iceberg parameters 3 3 !----------------------------------------------------------------------- 4 ln_icebergs = . true.4 ln_icebergs = .false. 5 5 ln_bergdia = .true. ! Calculate budgets 6 6 nn_verbose_level = 1 ! Turn on more verbose output if level > 0 7 7 nn_verbose_write = 15 ! Timesteps between verbose messages 8 8 nn_sample_rate = 1 ! Timesteps between sampling for trajectory storage 9 ! Initial mass required for an iceberg of each class 9 10 rn_initial_mass = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11 11 ! Proportion of calving mass to apportion to each class 10 12 rn_distribution = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02 11 13 ! Ratio between effective and real iceberg mass (non-dim) 14 ! i.e. number of icebergs represented at a point 12 15 rn_mass_scaling = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1 13 ! Totalthickness of newly calved bergs (m)16 ! thickness of newly calved bergs (m) 14 17 rn_initial_thickness = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250. 15 18 rn_rho_bergs = 850. ! Density of icebergs … … 18 21 rn_bits_erosion_fraction = 0. ! Fraction of erosion melt flux to divert to bergy bits 19 22 rn_sicn_shift = 0. ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1) 20 ln_passive_mode = .false. ! iceberg - ocean decoupling 21 nn_test_icebergs = 10 ! Create icebergs in absence of a calving file (-1 = no) 22 rn_test_box = -61.0, -55.0, 59.0, 65.0 23 rn_speed_limit = 0. ! CFL speed limit for a berg 23 ln_passive_mode = .false. ! iceberg - ocean decoupling 24 nn_test_icebergs = 10 ! Create test icebergs of this class (-1 = no) 25 ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2) 26 rn_test_box = 108.0, 116.0, -66.0, -58.0 27 rn_speed_limit = 0. ! CFL speed limit for a berg 24 28 25 29 ! filename ! freq (hours) ! variable ! time interp. ! clim !'yearly' or ! weights ! rotation ! 26 30 ! ! (<0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 27 31 sn_icb = 'calving' , -1 , 'calvingmask', .true. , .true., 'yearly' , ' ' , ' ' 28 29 cn_dir = './' 32 33 cn_dir = './' 30 34 / -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/nambfr
r3294 r3989 7 7 rn_bfri2 = 1.e-3 ! bottom drag coefficient (non linear case) 8 8 rn_bfeb2 = 2.5e-3 ! bottom turbulent kinetic energy background (m2/s2) 9 rn_bfrz0 = 3.e-3 ! bottom roughness for loglayer bfr coeff 9 10 ln_bfr2d = .false. ! horizontal variation of the bottom friction coef (read a 2D mask file ) 10 11 rn_bfrien = 50. ! local multiplying factor of bfr (ln_bfr2d=T) 11 ln_bfrimp = . false.! implicit bottom friction (requires ln_zdfexp = .false. if true)12 ln_bfrimp = .true. ! implicit bottom friction (requires ln_zdfexp = .false. if true) 12 13 / -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namctl
r3294 r3989 12 12 nn_bench = 0 ! Bench mode (1/0): CAUTION use zero except for bench 13 13 ! (no physical validity of the results) 14 nn_timing = 0 ! timing by routine activated (=1) creates timing.output file, or not (=0) 14 15 / -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namdct
r3294 r3989 4 4 nn_dct = 15 ! time step frequency for transports computing 5 5 nn_dctwri = 15 ! time step frequency for transports writing 6 nn_secdebug = 0! 0 : no section to debug6 nn_secdebug = 112 ! 0 : no section to debug 7 7 ! -1 : debug all section 8 8 ! 0 < n : debug section number n -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namdia_harm
r3294 r3989 2 2 &nam_diaharm ! Harmonic analysis of tidal constituents ('key_diaharm') 3 3 !----------------------------------------------------------------------- 4 nit000_han =1! First time step used for harmonic analysis5 nitend_han =105! Last time step used for harmonic analysis6 nstep_han =1! Time step frequency for harmonic analysis7 nb_ana=1 ! Number of harmonics to analyse8 tname( 1)='M2' ! Name of tidal constituents4 nit000_han = 1 ! First time step used for harmonic analysis 5 nitend_han = 75 ! Last time step used for harmonic analysis 6 nstep_han = 15 ! Time step frequency for harmonic analysis 7 tname(1) = 'M2' ! Name of tidal constituents 8 tname(2) = 'K1' 9 9 / -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namdyn_hpg
r3294 r3989 2 2 &namdyn_hpg ! Hydrostatic pressure gradient option 3 3 !----------------------------------------------------------------------- 4 ln_hpg_zco = .false. ! z-coordinate - full steps 4 ln_hpg_zco = .false. ! z-coordinate - full steps 5 5 ln_hpg_zps = .true. ! z-coordinate - partial steps (interpolation) 6 6 ln_hpg_sco = .false. ! s-coordinate (standard jacobian formulation) -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namdyn_ldf
r3294 r3989 2 2 &namdyn_ldf ! lateral diffusion on momentum 3 3 !----------------------------------------------------------------------- 4 ! ! Type of the operator : 5 ln_dynldf_lap = .true. ! laplacian operator 6 ln_dynldf_bilap = .false. ! bilaplacian operator 7 ! ! Direction of action : 8 ln_dynldf_level = .false. ! iso-level 4 ! ! Type of the operator : 5 ln_dynldf_lap = .true. ! laplacian operator 6 ln_dynldf_bilap = .false. ! bilaplacian operator 7 ! ! Direction of action : 8 ln_dynldf_level = .false. ! iso-level 9 9 ln_dynldf_hor = .true. ! horizontal (geopotential) (require "key_ldfslp" in s-coord.) 10 10 ln_dynldf_iso = .false. ! iso-neutral (require "key_ldfslp") … … 12 12 rn_ahm_0_lap = 40000. ! horizontal laplacian eddy viscosity [m2/s] 13 13 rn_ahmb_0 = 0. ! background eddy viscosity for ldf_iso [m2/s] 14 rn_ahm_0_blp = 0. ! horizontal bilaplacian eddy viscosity [m4/s] 14 rn_ahm_0_blp = 0. ! horizontal bilaplacian eddy viscosity [m4/s] 15 rn_cmsmag_1 = 3. ! constant in laplacian Smagorinsky viscosity 16 rn_cmsmag_2 = 3 ! constant in bilaplacian Smagorinsky viscosity 17 rn_cmsh = 1. ! 1 or 0 , if 0 -use only shear for Smagorinsky viscosity 18 rn_ahm_m_blp = -1.e12 ! upper limit for bilap abs(ahm) < min( dx^4/128rdt, rn_ahm_m_blp) 19 rn_ahm_m_lap = 40000. ! upper limit for lap ahm < min(dx^2/16rdt, rn_ahm_m_lap) 15 20 / -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namdyn_nept
r3294 r3989 1 1 !----------------------------------------------------------------------- 2 &nam _dynnept ! Neptune effect (simplified: lateral and vertical diffusions removed)2 &namdyn_nept ! Neptune effect (simplified: lateral and vertical diffusions removed) 3 3 !----------------------------------------------------------------------- 4 4 ! Suggested lengthscale values are those of Eby & Holloway (1994) for a coarse model … … 10 10 ! Specify whether to ramp down the Neptune velocity in shallow 11 11 ! water, and if so the depth range controlling such ramping down 12 ln_neptramp = . false.! ramp down Neptune velocity in shallow water12 ln_neptramp = .true. ! ramp down Neptune velocity in shallow water 13 13 rn_htrmin = 100.0 ! min. depth of transition range 14 14 rn_htrmax = 200.0 ! max. depth of transition range -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namdyn_vor
r3764 r3989 2 2 &namdyn_vor ! option of physics/algorithm (not control by CPP keys) 3 3 !----------------------------------------------------------------------- 4 ln_dynvor_ene = .false. ! energy conserving scheme 5 ln_dynvor_ens = .false. ! enstrophy conserving scheme 6 ln_dynvor_mix = .false. ! mixed scheme 7 ln_dynvor_een = .true. ! energy & enstrophy scheme 8 ln_dynvor_con = .false. ! consistency of BC with analytical eqs. 4 ln_dynvor_ene = .false. ! enstrophy conserving scheme 5 ln_dynvor_ens = .false. ! energy conserving scheme 6 ln_dynvor_mix = .false. ! mixed scheme 7 ln_dynvor_een = .true. ! energy & enstrophy scheme 9 8 / -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namflo
r3294 r3989 2 2 &namflo ! float parameters ("key_float") 3 3 !----------------------------------------------------------------------- 4 jpnfl = 1 ! total number of floats during the run5 jpnnewflo = 0 ! number of floats for the restart6 ln_rstflo = .false. ! float restart (T) or not (F)7 nn_writefl = 75 ! frequency of writing in float output file8 nn_stockfl = 5475 ! frequency of creation of the float restart file9 ln_argo = .false. ! Argo type floats (stay at the surface each 10 days)10 ln_flork4 = .false. ! trajectories computed with a 4th order Runge-Kutta (T)11 ! or computed with Blanke' scheme (F)12 ln_ariane = .true. ! Input with Ariane tool convention(T)13 ln_ ascii = .true. ! Output with Ariane tool netcdf convention(T) or ascii file (F)4 jpnfl = 1 ! total number of floats during the run 5 jpnnewflo = 0 ! number of floats for the restart 6 ln_rstflo = .false. ! float restart (T) or not (F) 7 nn_writefl = 75 ! frequency of writing in float output file 8 nn_stockfl = 5475 ! frequency of creation of the float restart file 9 ln_argo = .false. ! Argo type floats (stay at the surface each 10 days) 10 ln_flork4 = .false. ! trajectories computed with a 4th order Runge-Kutta (T) 11 ! or computed with Blanke' scheme (F) 12 ln_ariane = .true. ! Input with Ariane tool convention(T) 13 ln_flo_ascii = .true. ! Output with Ariane tool netcdf convention(F) or ascii file (T) 14 14 / -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namlbc
r3294 r3989 4 4 rn_shlat = 2. ! shlat = 0 ! 0 < shlat < 2 ! shlat = 2 ! 2 < shlat 5 5 ! free slip ! partial slip ! no slip ! strong slip 6 ln_vorlat = .false. ! consistency of vorticity boundary condition with analytical eqs. 6 7 / -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/nammpp
r3294 r3989 5 5 ! buffer blocking send or immediate non-blocking sends, resp. 6 6 nn_buffer = 0 ! size in bytes of exported buffer ('B' case), 0 no exportation 7 ln_nnogather= .false. ! activate code to avoid mpi_allgather use at the northfold 8 jpni = 0 ! jpni number of processors following i (set automatically if < 1) 9 jpnj = 0 ! jpnj number of processors following j (set automatically if < 1) 10 jpnij = 0 ! jpnij number of local domains (set automatically if < 1) 7 11 / -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namobc
r3294 r3989 4 4 ln_obc_clim = .false. ! climatological obc data files (T) or not (F) 5 5 ln_vol_cst = .true. ! impose the total volume conservation (T) or not (F) 6 ln_obc_fla = .false. ! Flather open boundary condition 6 ln_obc_fla = .false. ! Flather open boundary condition 7 7 nn_obcdta = 1 ! = 0 the obc data are equal to the initial state 8 8 ! = 1 the obc data are read in 'obc.dta' files -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namobs
r2540 r3989 2 2 &namobs ! observation usage switch ('key_diaobs') 3 3 !----------------------------------------------------------------------- 4 ln_t3d = .false. ! Logical switch for T profile observations 5 ln_s3d = .false. ! Logical switch for S profile observations 6 ln_ena = .false. ! Logical switch for ENACT insitu data set 7 ! ! ln_cor Logical switch for Coriolis insitu data set 8 ln_profb = .false. ! Logical switch for feedback insitu data set 9 ln_sla = .false. ! Logical switch for SLA observations 4 ln_t3d = .false. ! Logical switch for T profile observations 5 ln_s3d = .false. ! Logical switch for S profile observations 6 ln_ena = .false. ! Logical switch for ENACT insitu data set 7 ! ! ln_cor Logical switch for Coriolis insitu data set 8 ln_profb = .false. ! Logical switch for feedback insitu data set 9 ln_sla = .false. ! Logical switch for SLA observations 10 10 11 ln_sladt = .false. ! Logical switch for AVISO SLA data 11 ln_sladt = .false. ! Logical switch for AVISO SLA data 12 12 13 ln_slafb = .false. ! Logical switch for feedback SLA data 14 ! ln_ssh Logical switch for SSH observations 13 ln_slafb = .false. ! Logical switch for feedback SLA data 14 ! ln_ssh Logical switch for SSH observations 15 15 16 ln_sst = . false. ! Logical switch for SST observations17 ! ln_reysst Logical switch for Reynolds observations18 ! ln_ghrsst Logical switch for GHRSST observations16 ln_sst = .true. ! Logical switch for SST observations 17 ln_reysst = .true. ! ln_reysst Logical switch for Reynolds observations 18 ln_ghrsst = .false. ! ln_ghrsst Logical switch for GHRSST observations 19 19 20 ln_sstfb = .false. ! Logical switch for feedback SST data 21 ! ln_sss Logical switch for SSS observations 22 ! ln_seaice Logical switch for Sea Ice observations 23 ! ln_vel3d Logical switch for velocity observations 24 ! ln_velavcur Logical switch for velocity daily av. cur. 25 ! ln_velhrcur Logical switch for velocity high freq. cur. 26 ! ln_velavadcp Logical switch for velocity daily av. ADCP 20 ln_sstfb = .false. ! Logical switch for feedback SST data 21 ! ln_sss Logical switch for SSS observations 22 ! ln_seaice Logical switch for Sea Ice observations 23 ! ln_vel3d Logical switch for velocity observations 24 ! ln_velavcur Logical switch for velocity daily av. cur. 25 ! ln_velhrcur Logical switch for velocity high freq. cur. 26 ! ln_velavadcp Logical switch for velocity daily av. ADCP 27 27 ! ln_velhradcp Logical switch for velocity high freq. ADCP 28 ! ln_velfb Logical switch for feedback velocity data 29 ! ln_grid_global Global distribtion of observations 30 ! ln_grid_search_lookup Logical switch for obs grid search w/lookup table 31 ! grid_search_file Grid search lookup file header 32 ! enactfiles ENACT input observation file names 33 ! coriofiles Coriolis input observation file name 34 ! ! profbfiles: Profile feedback input observation file name 28 ! ln_velfb Logical switch for feedback velocity data 29 ! ln_grid_global Global distribtion of observations 30 ! ln_grid_search_lookup Logical switch for obs grid search w/lookup table 31 ! grid_search_file Grid search lookup file header 32 ! enactfiles ENACT input observation file names 33 ! coriofiles Coriolis input observation file name 34 ! ! profbfiles: Profile feedback input observation file name 35 35 profbfiles = 'profiles_01.nc' 36 ! ln_profb_enatim Enact feedback input time setting switch 36 ! ln_profb_enatim Enact feedback input time setting switch 37 37 ! slafilesact Active SLA input observation file name 38 ! slafilespas Passive SLA input observation file name 39 ! ! slafbfiles: Feedback SLA input observation file name 38 ! slafilespas Passive SLA input observation file name 39 ! ! slafbfiles: Feedback SLA input observation file name 40 40 slafbfiles = 'sla_01.nc' 41 ! sstfiles GHRSST input observation file name 42 ! ! sstfbfiles: Feedback SST input observation file name 41 ! sstfiles GHRSST input observation file name 42 ! ! sstfbfiles: Feedback SST input observation file name 43 43 sstfbfiles = 'sst_01.nc' 'sst_02.nc' 'sst_03.nc' 'sst_04.nc' 'sst_05.nc' 44 ! seaicefiles Sea Ice input observation file name 45 ! velavcurfiles Vel. cur. daily av. input file name 46 ! velhvcurfiles Vel. cur. high freq. input file name 47 ! velavadcpfiles Vel. ADCP daily av. input file name 48 ! velhvadcpfiles Vel. ADCP high freq. input file name 49 ! velfbfiles Vel. feedback input observation file name 50 ! dobsini Initial date in window YYYYMMDD.HHMMSS 51 ! dobsend Final date in window YYYYMMDD.HHMMSS 52 ! n1dint Type of vertical interpolation method 53 ! n2dint Type of horizontal interpolation method 54 ! ln_nea Rejection of observations near land switch 55 nmsshc = 0 ! MSSH correction scheme 56 ! mdtcorr MDT correction 57 ! mdtcutoff MDT cutoff for computed correction 58 ln_altbias = .false. ! Logical switch for alt bias 59 ln_ignmis = .true. ! Logical switch for ignoring missing files 60 ! endailyavtypes ENACT daily average types 44 ! seaicefiles Sea Ice input observation file name 45 ! velavcurfiles Vel. cur. daily av. input file name 46 ! velhvcurfiles Vel. cur. high freq. input file name 47 ! velavadcpfiles Vel. ADCP daily av. input file name 48 ! velhvadcpfiles Vel. ADCP high freq. input file name 49 ! velfbfiles Vel. feedback input observation file name 50 ! dobsini Initial date in window YYYYMMDD.HHMMSS 51 ! dobsend Final date in window YYYYMMDD.HHMMSS 52 ! n1dint Type of vertical interpolation method 53 ! n2dint Type of horizontal interpolation method 54 ! ln_nea Rejection of observations near land switch 55 nmsshc = 0 ! MSSH correction scheme 56 ! mdtcorr MDT correction 57 ! mdtcutoff MDT cutoff for computed correction 58 ln_altbias = .false. ! Logical switch for alt bias 59 ln_ignmis = .true. ! Logical switch for ignoring missing files 60 ! endailyavtypes ENACT daily average types 61 61 ln_grid_global = .true. 62 62 ln_grid_search_lookup = .false. 63 / 63 / -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namptr
r3294 r3989 4 4 ln_diaptr = .false. ! Poleward heat and salt transport (T) or not (F) 5 5 ln_diaznl = .true. ! Add zonal means and meridional stream functions 6 ln_subbas = .true. ! Atlantic/Pacific/Indian basins computation (T) or not 6 ln_subbas = .true. ! Atlantic/Pacific/Indian basins computation (T) or not 7 7 ! (orca configuration only, need input basins mask file named "subbasins.nc" 8 8 ln_ptrcomp = .true. ! Add decomposition : overturning -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namrun
r3306 r3989 3 3 !----------------------------------------------------------------------- 4 4 nn_no = 0 ! job number (no more used...) 5 cn_exp = "ORCA2" ! experience name 5 cn_exp = "ORCA2" ! experience name 6 6 nn_it000 = 1 ! first time step 7 7 nn_itend = 5475 ! last time step (std 5475) 8 nn_date0 = 010101 ! date at nit_0000 (format yyyymmdd) 9 ! used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) 8 nn_date0 = 010101 ! date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) 10 9 nn_leapy = 0 ! Leap year calendar (1) or not (0) 11 10 ln_rstart = .false. ! start from rest (F) or from a restart file (T) 12 11 nn_rstctl = 0 ! restart control => activated only if ln_rstart = T 13 != 0 nn_date0 read in namelist ; nn_it000 : read in namelist14 != 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart15 != 2 nn_date0 read in restart ; nn_it000 : check consistancy between namelist and restart12 ! = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 13 ! = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart 14 ! = 2 nn_date0 read in restart ; nn_it000 : check consistancy between namelist and restart 16 15 cn_ocerst_in = "restart" ! suffix of ocean restart name (input) 17 16 cn_ocerst_out = "restart" ! suffix of ocean restart name (output) -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namsbc
r3294 r3989 2 2 &namsbc ! Surface Boundary Condition (surface module) 3 3 !----------------------------------------------------------------------- 4 nn_fsbc = 5 ! frequency of surface boundary condition computation 4 nn_fsbc = 5 ! frequency of surface boundary condition computation 5 5 ! (also = the frequency of sea-ice model call) 6 ln_ana = .false. ! analytical formulation (T => fill namsbc_ana ) 6 ln_ana = .false. ! analytical formulation (T => fill namsbc_ana ) 7 7 ln_flx = .false. ! flux formulation (T => fill namsbc_flx ) 8 ln_blk_clio = .false. ! CLIO bulk formulation (T => fill namsbc_clio) 9 ln_blk_core = .true. ! CORE bulk formulation (T => fill namsbc_core) 8 ln_blk_clio = .false. ! CLIO bulk formulation (T => fill namsbc_clio) 9 ln_blk_core = .true. ! CORE bulk formulation (T => fill namsbc_core) 10 10 ln_blk_mfs = .false. ! MFS bulk formulation (T => fill namsbc_mfs ) 11 11 ln_cpl = .false. ! Coupled formulation (T => fill namsbc_cpl ) … … 14 14 ! =1 use observed ice-cover , 15 15 ! =2 ice-model used ("key_lim3" or "key_lim2) 16 nn_ice_embd = 0 ! =0 levitating ice (no mass exchange, concentration/dilution effect) 17 ! =1 levitating ice with mass and salt exchange but no presure effect 18 ! =2 embedded sea-ice (full salt and mass exchanges and pressure) 16 19 ln_dm2dc = .false. ! daily mean to diurnal cycle on short wave 17 20 ln_rnf = .true. ! runoffs (T => fill namsbc_rnf) 18 21 ln_ssr = .true. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) 19 nn_fwb = 3 ! FreshWater Budget: =0 unchecked 20 ! =1 global mean of e-p-r set to zero at each time step 22 nn_fwb = 3 ! FreshWater Budget: =0 unchecked 23 ! =1 global mean of e-p-r set to zero at each time step 21 24 ! =2 annual global mean of e-p-r set to zero 22 25 ! =3 global emp set to zero and spread out over erp area 26 ln_wave = .false. ! Activate coupling with wave (either Stokes Drift or Drag coefficient, or both) (T => fill namsbc_wave) 23 27 ln_cdgw = .false. ! Neutral drag coefficient read from wave model (T => fill namsbc_wave) 28 ln_sdw = .false. ! Computation of 3D stokes drift (T => fill namsbc_wave) 24 29 / -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namsbc_alb
r3294 r3989 2 2 &namsbc_alb ! albedo parameters 3 3 !----------------------------------------------------------------------- 4 rn_cloud = 0.06 ! cloud correction to snow and ice albedo 4 rn_cloud = 0.06 ! cloud correction to snow and ice albedo 5 5 rn_albice = 0.53 ! albedo of melting ice in the arctic and antarctic 6 6 rn_alphd = 0.80 ! coefficients for linear interpolation used to 7 rn_alphc = 0.65 ! compute albedo between two extremes values 7 rn_alphc = 0.65 ! compute albedo between two extremes values 8 8 rn_alphdi = 0.72 ! (Pyane, 1972) 9 9 / -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namsbc_clio
r3294 r3989 1 1 !----------------------------------------------------------------------- 2 &namsbc_clio ! namsbc_clio CLIO bulk formul ea2 &namsbc_clio ! namsbc_clio CLIO bulk formulae 3 3 !----------------------------------------------------------------------- 4 4 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namsbc_core
r3294 r3989 1 1 !----------------------------------------------------------------------- 2 &namsbc_core ! namsbc_core CORE bulk formul ea2 &namsbc_core ! namsbc_core CORE bulk formulae 3 3 !----------------------------------------------------------------------- 4 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation !5 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing !6 sn_wndi = 'u_10.15JUNE2009_orca2' , 6, 'U_10_MOD', .false. , .true. , 'yearly' , '' , 'Uwnd'7 sn_wndj = 'v_10.15JUNE2009_orca2' , 6, 'V_10_MOD', .false. , .true. , 'yearly' , '' , 'Vwnd'8 sn_qsr = 'ncar_rad.15JUNE2009_orca2' , 24, 'SWDN_MOD', .false. , .true. , 'yearly' , '' , ''9 sn_qlw = 'ncar_rad.15JUNE2009_orca2' , 24, 'LWDN_MOD', .false. , .true. , 'yearly' , '' , ''10 sn_tair = 't_10.15JUNE2009_orca2' , 6, 'T_10_MOD', .false. , .true. , 'yearly' , '' , ''11 sn_humi = 'q_10.15JUNE2009_orca2' , 6, 'Q_10_MOD', .false. , .true. , 'yearly' , '' , ''12 sn_prec = 'ncar_precip.15JUNE2009_orca2', -1, 'PRC_MOD1', .false. , .true. , 'yearly' , '' , ''13 sn_snow = 'ncar_precip.15JUNE2009_orca2', -1, 'SNOW' , .false. , .true. , 'yearly' , '' , ''14 sn_tdif = 'taudif_core' , 24, 'taudif' , .false. , .true. , 'yearly' , '' , ''4 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! 5 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 6 sn_wndi = 'u_10.15JUNE2009_orca2' , 6 , 'U_10_MOD', .false. , .true. , 'yearly' , '' , 'Uwnd' 7 sn_wndj = 'v_10.15JUNE2009_orca2' , 6 , 'V_10_MOD', .false. , .true. , 'yearly' , '' , 'Vwnd' 8 sn_qsr = 'ncar_rad.15JUNE2009_orca2' , 24 , 'SWDN_MOD', .false. , .true. , 'yearly' , '' , '' 9 sn_qlw = 'ncar_rad.15JUNE2009_orca2' , 24 , 'LWDN_MOD', .false. , .true. , 'yearly' , '' , '' 10 sn_tair = 't_10.15JUNE2009_orca2' , 6 , 'T_10_MOD', .false. , .true. , 'yearly' , '' , '' 11 sn_humi = 'q_10.15JUNE2009_orca2' , 6 , 'Q_10_MOD', .false. , .true. , 'yearly' , '' , '' 12 sn_prec = 'ncar_precip.15JUNE2009_orca2', -1 , 'PRC_MOD1', .false. , .true. , 'yearly' , '' , '' 13 sn_snow = 'ncar_precip.15JUNE2009_orca2', -1 , 'SNOW' , .false. , .true. , 'yearly' , '' , '' 14 sn_tdif = 'taudif_core' , 24 , 'taudif' , .false. , .true. , 'yearly' , '' , '' 15 15 16 16 cn_dir = './' ! root directory for the location of the bulk files -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namsbc_cpl
r3294 r3989 5 5 ! ! ! categories ! reference ! orientation ! grids ! 6 6 ! send 7 sn_snd_temp = 'weighted oce and ice' , 'no' , '' , '' , '' 8 sn_snd_alb = 'weighted ice' , 'no' , '' , '' , '' 9 sn_snd_thick = 'none' , 'no' , '' , '' , ''10 sn_snd_crt = 'none' , 'no' , 'spherical' , 'eastward-northward' , 'T' 11 sn_snd_co2 = 'coupled' , 'no' , '' , '' , '' 7 sn_snd_temp = 'weighted oce and ice' , 'no' , '' , '' , '' 8 sn_snd_alb = 'weighted ice' , 'no' , '' , '' , '' 9 sn_snd_thick = 'none' , 'no' , '' , '' , '' 10 sn_snd_crt = 'none' , 'no' , 'spherical' , 'eastward-northward' , 'T' 11 sn_snd_co2 = 'coupled' , 'no' , '' , '' , '' 12 12 ! receive 13 sn_rcv_w10m = 'none' , 'no' , '' , '' , '' 14 sn_rcv_taumod = 'coupled' , 'no' , '' , '' , '' 15 sn_rcv_tau = 'oce only' , 'no' , 'cartesian' , 'eastward-northward', 'U,V' 16 sn_rcv_dqnsdt = 'coupled' , 'no' , '' , '' , '' 17 sn_rcv_qsr = 'oce and ice' , 'no' , '' , '' , '' 18 sn_rcv_qns = 'oce and ice' , 'no' , '' , '' , '' 19 sn_rcv_emp = 'conservative' , 'no' , '' , '' , '' 20 sn_rcv_rnf = 'coupled' , 'no' , '' , '' , '' 21 sn_rcv_cal = 'coupled' , 'no' , '' , '' , '' 22 sn_rcv_co2 = 'coupled' , 'no' , '' , '' , '' 13 sn_rcv_w10m = 'none' , 'no' , '' , '' , '' 14 sn_rcv_taumod = 'coupled' , 'no' , '' , '' , '' 15 sn_rcv_tau = 'oce only' , 'no' , 'cartesian' , 'eastward-northward', 'U,V' 16 sn_rcv_dqnsdt = 'coupled' , 'no' , '' , '' , '' 17 sn_rcv_qsr = 'oce and ice' , 'no' , '' , '' , '' 18 sn_rcv_qns = 'oce and ice' , 'no' , '' , '' , '' 19 sn_rcv_emp = 'conservative' , 'no' , '' , '' , '' 20 sn_rcv_rnf = 'coupled' , 'no' , '' , '' , '' 21 sn_rcv_cal = 'coupled' , 'no' , '' , '' , '' 22 sn_rcv_co2 = 'coupled' , 'no' , '' , '' , '' 23 23 / -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namsbc_flx
r3294 r3989 11 11 12 12 cn_dir = './' ! root directory for the location of the flux files 13 / 13 / -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namsbc_mfs
r3294 r3989 10 10 sn_tair = 'ecmwf' , 6 , 't2' , .true. , .false. , 'daily' ,'bicubic.nc' , '' 11 11 sn_rhm = 'ecmwf' , 6 , 'rh' , .true. , .false. , 'daily' ,'bilinear.nc', '' 12 sn_prec = ' precip' , 6 , 'precip' , .true. , .false. , 'daily' ,'bicubic', ''12 sn_prec = 'ecmwf' , 6 , 'precip' , .true. , .true. , 'daily' ,'bicubic.nc' , '' 13 13 14 14 cn_dir = './ECMWF/' ! root directory for the location of the bulk files -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namsbc_rnf
r3294 r3989 2 2 &namsbc_rnf ! runoffs namelist surface boundary condition 3 3 !----------------------------------------------------------------------- 4 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation !5 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing !6 sn_rnf = 'runoff_core_monthly', -1, 'sorunoff', .true. , .true. , 'yearly' , '' , ''7 sn_cnf = 'runoff_core_monthly', 0, 'socoefr0', .false. , .true. , 'yearly' , '' , ''8 sn_s_rnf = 'runoffs' , 24, 'rosaline', .true. , .true. , 'yearly' , '' , ''9 sn_t_rnf = 'runoffs' , 24, 'rotemper', .true. , .true. , 'yearly' , '' , ''10 sn_dep_rnf = 'runoffs' , 0, 'rodepth' , .false. , .true. , 'yearly' , '' , ''4 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! 5 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 6 sn_rnf = 'runoff_core_monthly', -1 , 'sorunoff', .true. , .true. , 'yearly' , '' , '' 7 sn_cnf = 'runoff_core_monthly', 0 , 'socoefr0', .false. , .true. , 'yearly' , '' , '' 8 sn_s_rnf = 'runoffs' , 24 , 'rosaline', .true. , .true. , 'yearly' , '' , '' 9 sn_t_rnf = 'runoffs' , 24 , 'rotemper', .true. , .true. , 'yearly' , '' , '' 10 sn_dep_rnf = 'runoffs' , 0 , 'rodepth' , .false. , .true. , 'yearly' , '' , '' 11 11 12 12 cn_dir = './' ! root directory for the location of the runoff files -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namsbc_ssr
r3294 r3989 9 9 cn_dir = './' ! root directory for the location of the runoff files 10 10 nn_sstr = 0 ! add a retroaction term in the surface heat flux (=1) or not (=0) 11 nn_sssr = 2 ! add a damping term in the surface freshwater flux (=2) 11 nn_sssr = 2 ! add a damping term in the surface freshwater flux (=2) 12 12 ! or to SSS only (=1) or no damping term (=0) 13 13 rn_dqdt = -40. ! magnitude of the retroaction on temperature [W/m2/K] 14 rn_deds = -27.7! magnitude of the damping on salinity [mm/day]14 rn_deds = -166.67 ! magnitude of the damping on salinity [mm/day] 15 15 ln_sssr_bnd = .true. ! flag to bound erp term (associated with nn_sssr=2) 16 16 rn_sssr_bnd = 4.e0 ! ABS(Max/Min) value of the damping erp term [mm/day] 17 / 17 / -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namsbc_wave
r3294 r3989 4 4 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! 5 5 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 6 sn_cdg = 'cdg_wave' , 1 , 'drag_coeff' , .true. , .false. , 'daily' , '' , '' 6 sn_cdg = 'cdg_wave' , 1 , 'drag_coeff' , .true. , .false. , 'daily' ,'' , '' 7 sn_usd = 'sdw_wave' , 1 , 'u_sd2d' , .true. , .false. , 'daily' ,'' , '' 8 sn_vsd = 'sdw_wave' , 1 , 'v_sd2d' , .true. , .false. , 'daily' ,'' , '' 9 sn_wn = 'sdw_wave' , 1 , 'wave_num' , .true. , .false. , 'daily' ,'' , '' 7 10 ! 8 11 cn_dir_cdg = './' ! root directory for the location of drag coefficient files -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namsol
r3294 r3989 1 1 !----------------------------------------------------------------------- 2 &namsol ! elliptic solver / island / free surface 2 &namsol ! elliptic solver / island / free surface 3 3 !----------------------------------------------------------------------- 4 4 nn_solv = 1 ! elliptic solver: =1 preconditioned conjugate gradient (pcg) -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namtra_adv
r3294 r3989 1 1 !----------------------------------------------------------------------- 2 &namtra_adv ! advection scheme for tracer 2 &namtra_adv ! advection scheme for tracer 3 3 !----------------------------------------------------------------------- 4 ln_traadv_cen2 = .false. ! 2nd order centered scheme 5 ln_traadv_tvd = .true. ! TVD scheme 6 ln_traadv_muscl = .false. ! MUSCL scheme 7 ln_traadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 8 ln_traadv_ubs = .false. ! UBS scheme 9 ln_traadv_qck = .false. ! QUCIKEST scheme 4 ln_traadv_cen2 = .false. ! 2nd order centered scheme 5 ln_traadv_tvd = .true. ! TVD scheme 6 ln_traadv_muscl = .false. ! MUSCL scheme 7 ln_traadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 8 ln_traadv_ubs = .false. ! UBS scheme 9 ln_traadv_qck = .false. ! QUICKEST scheme 10 ln_traadv_msc_ups= .false. ! use upstream scheme within muscl 10 11 / -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namtra_bbc
r2540 r3989 3 3 !----------------------------------------------------------------------- 4 4 ln_trabbc = .true. ! Apply a geothermal heating at the ocean bottom 5 nn_geoflx = 2 ! geothermal heat flux: = 0 no flux 5 nn_geoflx = 2 ! geothermal heat flux: = 0 no flux 6 6 ! = 1 constant flux 7 ! = 2 variable flux (read in geothermal_heating.nc in mW/m2) 7 ! = 2 variable flux (read in geothermal_heating.nc in mW/m2) 8 8 rn_geoflx_cst = 86.4e-3 ! Constant value of geothermal heat flux [W/m2] 9 9 / -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namtra_ldf
r3296 r3989 22 22 rn_ahtb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 23 23 ! (normally=0; not used with Griffies) 24 rn_slpmax = 0.01 ! slope limit 25 rn_chsmag = 1. ! multiplicative factor in Smagorinsky diffusivity 26 rn_smsh = 1. ! Smagorinsky diffusivity: = 0 - use only sheer 27 rn_aht_m = 2000. ! upper limit or stability criteria for lateral eddy diffusivity (m2/s) 24 28 / -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namtsd
r3294 r3989 2 2 &namtsd ! data : Temperature & Salinity 3 3 !----------------------------------------------------------------------- 4 ! ! file name ! frequency (hours) ! variable! time interp. ! clim !'yearly' or ! weights ! rotation !5 ! ! ! (if <0 months) ! name! (logical) ! (T/F) ! 'monthly' ! filename ! pairing !6 sn_tem = 'data_1m_potential_temperature_nomask', -1,'votemper', .true. 7 sn_sal = 'data_1m_salinity_nomask' , -1,'vosaline', .true. 4 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim !'yearly' or ! weights ! rotation ! 5 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 6 sn_tem = 'data_1m_potential_temperature_nomask', -1,'votemper', .true. , .true., 'yearly' , ' ' , ' ' 7 sn_sal = 'data_1m_salinity_nomask' , -1,'vosaline', .true. , .true., 'yearly' , '' , ' ' 8 8 ! 9 9 cn_dir = './' ! root directory for the location of the runoff files -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namzdf_kpp
r3294 r3989 2 2 &namzdf_kpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionally: 3 3 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 4 ln_kpprimix = .true. ! shear instability mixing 4 ln_kpprimix = .true. ! shear instability mixing 5 5 rn_difmiw = 1.0e-04 ! constant internal wave viscosity [m2/s] 6 6 rn_difsiw = 0.1e-04 ! constant internal wave diffusivity [m2/s] 7 7 rn_riinfty = 0.8 ! local Richardson Number limit for shear instability 8 8 rn_difri = 0.0050 ! maximum shear mixing at Rig = 0 [m2/s] 9 rn_bvsqcon = -0.01e-07 ! Brunt-Vaisala squared for maximum convection [1/s2] 10 rn_difcon = 1. ! maximum mixing in interior convection [m2/s] 9 rn_bvsqcon = -0.01e-07 ! Brunt-Vaisala squared for maximum convection [1/s2] 10 rn_difcon = 1. ! maximum mixing in interior convection [m2/s] 11 11 nn_avb = 0 ! horizontal averaged (=1) or not (=0) on avt and amv 12 12 nn_ave = 1 ! constant (=0) or profile (=1) background on avt -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namzdf_tke
r3294 r3989 7 7 rn_emin = 1.e-6 ! minimum value of tke [m2/s2] 8 8 rn_emin0 = 1.e-4 ! surface minimum value of tke [m2/s2] 9 rn_bshear = 1.e-20 ! background shear (>0) currently a numerical threshold (do not change it) 9 10 nn_mxl = 2 ! mixing length: = 0 bounded by the distance to surface and bottom 10 11 ! = 1 bounded by the local vertical scale factor -
branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namzdf_tmx
r2540 r3989 5 5 rn_n2min = 1.e-8 ! threshold of the Brunt-Vaisala frequency (s-1) 6 6 rn_tfe = 0.333 ! tidal dissipation efficiency 7 rn_me = 0.2 ! mixing efficiency 7 rn_me = 0.2 ! mixing efficiency 8 8 ln_tmx_itf = .true. ! ITF specific parameterisation 9 9 rn_tfe_itf = 1. ! ITF tidal dissipation efficiency
Note: See TracChangeset
for help on using the changeset viewer.