Changeset 3989
- Timestamp:
- 2013-07-24T11:48:35+02:00 (11 years ago)
- Location:
- branches/2013/dev_r3853_CNRS9_ConfSetting
- Files:
-
- 2 deleted
- 64 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 -
branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/CONFIG/SHARED/README_other_configurations_namelist_namcfg
r3973 r3989 152 152 !----------------------------------------------------------------------- 153 153 cp_cfg = "orca" ! name of the configuration 154 cp_cfz = "antarctic" ! name of the zoom of configuration 154 155 jp_cfg = 05 ! resolution of the configuration 155 156 jpidta = 722 ! 1st lateral dimension ( >= jpi ) … … 188 189 !----------------------------------------------------------------------- 189 190 cp_cfg = "orca" ! name of the configuration 191 cp_cfz = "arctic" ! name of the zoom of configuration 190 192 jp_cfg = 05 ! resolution of the configuration 191 193 jpidta = 722 ! 1st lateral dimension ( >= jpi ) … … 216 218 ppacr2 = 999999.0 ! 217 219 / 220 221 ORCA2 - Antarctic zoom 222 ====================== 223 !----------------------------------------------------------------------- 224 &namcfg ! parameters of the configuration 225 !----------------------------------------------------------------------- 226 cp_cfg = "orca" ! name of the configuration 227 cp_cfz = "antarctic" ! name of the zoom of configuration 228 jp_cfg = 2 ! resolution of the configuration 229 jpidta = 182 ! 1st lateral dimension ( >= jpi ) 230 jpjdta = 149 ! 2nd " " ( >= jpj ) 231 jpkdta = 31 ! number of levels ( >= jpk ) 232 jpiglo = 182 ! 1st dimension of global domain --> i =jpidta 233 jpjglo = 50 ! 2nd - - --> j =jpjdta 234 jpizoom = 1 ! left bottom (i,j) indices of the zoom 235 jpjzoom = 1 ! in data domain indices 236 jperio = 1 ! lateral cond. type (between 0 and 6) 237 jphgr_msh = 0 ! type of horizontal mesh 238 ppglam0 = 999999.0 ! longitude of first raw and column T-point (jphgr_msh = 1) 239 ppgphi0 = 999999.0 ! latitude of first raw and column T-point (jphgr_msh = 1) 240 ppe1_deg = 999999.0 ! zonal grid-spacing (degrees) 241 ppe2_deg = 999999.0 ! meridional grid-spacing (degrees) 242 ppe1_m = 999999.0 ! zonal grid-spacing (degrees) 243 ppe2_m = 999999.0 ! meridional grid-spacing (degrees) 244 ppsur = -4762.96143546300 ! ORCA r4, r2 and r05 coefficients 245 ppa0 = 255.58049070440 ! (default coefficients) 246 ppa1 = 245.58132232490 ! 247 ppkth = 21.43336197938 ! 248 ppacr = 3.0 ! 249 ppdzmin = 999999. ! Minimum vertical spacing 250 pphmax = 999999. ! Maximum depth 251 ldbletanh = .FALSE. ! Use/do not use double tanf function for vertical coordinates 252 ppa2 = 999999. ! Double tanh function parameters 253 ppkth2 = 999999. ! 254 ppacr2 = 999999. ! 255 / 256 257 258 ORCA2 - Arctic zoom 259 =================== 260 !----------------------------------------------------------------------- 261 &namcfg ! parameters of the configuration 262 !----------------------------------------------------------------------- 263 cp_cfg = "orca" ! name of the configuration 264 cp_cfz = "arctic" ! name of the zoom of configuration 265 jp_cfg = 2 ! resolution of the configuration 266 jpidta = 182 ! 1st lateral dimension ( >= jpi ) 267 jpjdta = 149 ! 2nd " " ( >= jpj ) 268 jpkdta = 31 ! number of levels ( >= jpk ) 269 jpiglo = 142 ! 1st dimension of global domain --> i =jpidta 270 jpjglo = 53 ! 2nd - - --> j =jpjdta 271 jpizoom = 21 ! left bottom (i,j) indices of the zoom 272 jpjzoom = 97 ! in data domain indices 273 jperio = 3 ! lateral cond. type (between 0 and 6) 274 jphgr_msh = 0 ! type of horizontal mesh 275 ppglam0 = 999999.0 ! longitude of first raw and column T-point (jphgr_msh = 1) 276 ppgphi0 = 999999.0 ! latitude of first raw and column T-point (jphgr_msh = 1) 277 ppe1_deg = 999999.0 ! zonal grid-spacing (degrees) 278 ppe2_deg = 999999.0 ! meridional grid-spacing (degrees) 279 ppe1_m = 999999.0 ! zonal grid-spacing (degrees) 280 ppe2_m = 999999.0 ! meridional grid-spacing (degrees) 281 ppsur = -4762.96143546300 ! ORCA r4, r2 and r05 coefficients 282 ppa0 = 255.58049070440 ! (default coefficients) 283 ppa1 = 245.58132232490 ! 284 ppkth = 21.43336197938 ! 285 ppacr = 3.0 ! 286 ppdzmin = 999999. ! Minimum vertical spacing 287 pphmax = 999999. ! Maximum depth 288 ldbletanh = .FALSE. ! Use/do not use double tanf function for vertical coordinates 289 ppa2 = 999999. ! Double tanh function parameters 290 ppkth2 = 999999. ! 291 ppacr2 = 999999. ! 292 / 293 294 218 295 ORCA025 - 75 vertical levels 219 296 ======= … … 428 505 ppacr2 = 999999.0 ! 429 506 / 507 508 C1D - 1D configuration. Add key_c1d in active cpp keys 509 !----------------------------------------------------------------------- 510 &namcfg ! parameters of the configuration 511 !----------------------------------------------------------------------- 512 cp_cfg = "orca" ! name of the configuration 513 jp_cfg = 2 ! resolution of the configuration 514 jpidta = 182 ! 1st lateral dimension ( >= jpi ) 515 jpjdta = 149 ! 2nd " " ( >= jpj ) 516 jpkdta = 31 ! number of levels ( >= jpk ) 517 jpiglo = 3 ! 1st dimension of global domain --> i =jpidta 518 jpjglo = 3 ! 2nd - - --> j =jpjdta 519 ! Choose postion of the 1D column: 520 ! jpizoom = 61, jpjzoom = 133 (160W,75N) 521 ! jpizoom = 61, jpjzoom = 110 (160W,50N) 522 ! jpizoom = 61, jpjzoom = 97 (160W,30N) 523 ! jpizoom = 61, jpjzoom = 86 (160W,10N) 524 ! jpizoom = 61, jpjzoom = 49 (160W,30S) 525 ! jpizoom = 61, jpjzoom = 27 (160W,60S) 526 ! jpizoom = 61, jpjzoom = 7 (160W,75S) 527 ! jpizoom = 110,jpjzoom = 97 (64W,31.5N) BATS site 528 jpizoom = 1 ! left bottom (i,j) indices of the zoom 529 jpjzoom = 1 ! in data domain indices 530 jperio = 0 ! lateral cond. type (between 0 and 6) 531 jphgr_msh = 0 ! type of horizontal mesh 532 ppglam0 = 999999.0 ! longitude of first raw and column T-point (jphgr_msh = 1) 533 ppgphi0 = 999999.0 ! latitude of first raw and column T-point (jphgr_msh = 1) 534 ppe1_deg = 999999.0 ! zonal grid-spacing (degrees) 535 ppe2_deg = 999999.0 ! meridional grid-spacing (degrees) 536 ppe1_m = 999999.0 ! zonal grid-spacing (degrees) 537 ppe2_m = 999999.0 ! meridional grid-spacing (degrees) 538 ppsur = -4762.96143546300 ! ORCA r4, r2 and r05 coefficients 539 ppa0 = 255.58049070440 ! (default coefficients) 540 ppa1 = 245.58132232490 ! 541 ppkth = 21.43336197938 ! 542 ppacr = 3.0 ! 543 ppdzmin = 999999. ! Minimum vertical spacing 544 pphmax = 999999. ! Maximum depth 545 ldbletanh = .FALSE. ! Use/do not use double tanf function for vertical coordinates 546 ppa2 = 999999. ! Double tanh function parameters 547 ppkth2 = 999999. ! 548 ppacr2 = 999999. ! 549 / -
branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/CONFIG/SHARED/namelist_ref
r3973 r3989 48 48 ! 49 49 !----------------------------------------------------------------------- 50 &namcfg ! parameters of the configuration (to be changed in your namelist_cfg in CONFIG/"YOUR_CONFIG"/EXP00 51 ! For ORCA2 Antartic zoom, use &namcfg from ORCA2_LIM/EXP00/namelist_cfg changing 52 ! jpjglo = 50, jperio = 1, 53 ! For ORCA2 Arctic zoom, use &namcfg from ORCA2_LIM/EXP00/namelist_cfg changing 54 ! jpiglo = 142, jpjglo = jpjdta-97+1, jpizoom = 21, jpjzoom = 97, jperio = 3 55 ! For 1D configuration, use &namcfg from ORCA2_LIM/EXP00/namelist_cfg changing 56 ! jpiglo = 3, jpjglo = 3, jperio = 0 and choose postion of the 1D column: 57 ! jpizoom = 61, jpjzoom = 133 (160W,75N) 58 ! jpizoom = 61, jpjzoom = 110 (160W,50N) 59 ! jpizoom = 61, jpjzoom = 97 (160W,30N) 60 ! jpizoom = 61, jpjzoom = 86 (160W,10N) 61 ! jpizoom = 61, jpjzoom = 49 (160W,30S) 62 ! jpizoom = 61, jpjzoom = 27 (160W,60S) 63 ! jpizoom = 61, jpjzoom = 7 (160W,75S) 64 ! jpizoom = 110,jpjzoom = 97 (64W,31.5N) BATS site 65 ! 50 &namcfg ! default parameters of the configuration 66 51 !----------------------------------------------------------------------- 67 52 cp_cfg = "default" ! name of the configuration 53 cp_cfz = '' ! name of the zoom of configuration 68 54 jp_cfg = 0 ! resolution of the configuration 69 55 jpidta = 10 ! 1st lateral dimension ( >= jpi ) -
branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/DOM/dom_oce.F90
r3973 r3989 72 72 LOGICAL, PUBLIC :: lzoom_s = .FALSE. !: South zoom type flag 73 73 LOGICAL, PUBLIC :: lzoom_n = .FALSE. !: North zoom type flag 74 LOGICAL, PUBLIC :: lzoom_arct = .FALSE. !: ORCA arctic zoom flag75 LOGICAL, PUBLIC :: lzoom_anta = .FALSE. !: ORCA antarctic zoom flag76 74 77 75 ! !!! domain parameters linked to mpp -
branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/DOM/domcfg.F90
r3973 r3989 168 168 SELECT CASE ( jp_cfg ) 169 169 CASE ( 2 ) ! ORCA_R2 configuration 170 IF( jpiglo == 142 .AND. jpjglo == 53 .AND. & 171 & jpizoom == 21 .AND. jpjzoom == 97 ) lzoom_arct = .TRUE. 172 IF( jpiglo == jpidta .AND. jpjglo == 50 .AND. & 173 & jpizoom == 1 .AND. jpjzoom == 1 ) lzoom_anta = .TRUE. 170 IF( cp_cfz == "arctic" .AND. jpiglo == 142 .AND. jpjglo == 53 .AND. & 171 & jpizoom == 21 .AND. jpjzoom == 97 ) THEN 172 IF(lwp) WRITE(numout,*) ' ORCA configuration: arctic zoom ' 173 ENDIF 174 IF( cp_cfz == "antarctic" .AND. jpiglo == jpidta .AND. jpjglo == 50 .AND. & 175 & jpizoom == 1 .AND. jpjzoom == 1 ) THEN 176 IF(lwp) WRITE(numout,*) ' ORCA configuration: antarctic zoom ' 177 ENDIF 174 178 ! 175 179 CASE ( 05 ) ! ORCA_R05 configuration 176 IF( jpiglo == 562 .AND. jpjglo == 202 .AND. & 177 & jpizoom == 81 .AND. jpjzoom == 301 ) lzoom_arct = .TRUE. 178 IF( jpiglo == jpidta .AND. jpjglo == 187 .AND. & 179 & jpizoom == 1 .AND. jpjzoom == 1 ) lzoom_anta = .TRUE. 180 IF( cp_cfz == "arctic" .AND. jpiglo == 562 .AND. jpjglo == 202 .AND. & 181 & jpizoom == 81 .AND. jpjzoom == 301 ) THEN 182 IF(lwp) WRITE(numout,*) ' ORCA configuration: arctic zoom ' 183 ENDIF 184 IF( cp_cfz == "antarctic" .AND. jpiglo == jpidta .AND. jpjglo == 187 .AND. & 185 & jpizoom == 1 .AND. jpjzoom == 1 ) THEN 186 IF(lwp) WRITE(numout,*) ' ORCA configuration: antarctic zoom ' 187 ENDIF 180 188 END SELECT 181 !182 IF(lwp) WRITE(numout,*) ' ORCA configuration: antarctic/arctic zoom flags : '183 IF(lwp) WRITE(numout,*) ' lzoom_arct = ', lzoom_arct, ' (T= arctic zoom, F=global)'184 IF(lwp) WRITE(numout,*) ' lzoom_anta = ', lzoom_anta, ' (T=antarctic zoom, F=global)'185 189 ! 186 190 ENDIF -
branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90
r3901 r3989 564 564 ! Configuration specific domain modifications 565 565 ! (here, ORCA arctic configuration: suppress Med Sea) 566 IF( cp_cfg == "orca" .AND. lzoom_arct) THEN566 IF( cp_cfg == "orca" .AND. cp_cfz == "arctic" ) THEN 567 567 SELECT CASE ( jp_cfg ) 568 568 ! ! ======================= -
branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_c2d.h90
r3294 r3989 160 160 IF(lwp) WRITE(numout,*) ' orca ocean configuration' 161 161 162 #if defined key_antarctic 163 # include "ldfdyn_antarctic.h90" 164 #elif defined key_arctic 165 # include "ldfdyn_arctic.h90" 166 #else 167 ! Read 2d integer array to specify western boundary increase in the 168 ! ===================== equatorial strip (20N-20S) defined at t-points 169 170 CALL ctl_opn( inum, 'ahmcoef', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp ) 171 READ(inum,9101) clexp, iim, ijm 172 READ(inum,'(/)') 173 ifreq = 40 174 il1 = 1 175 DO jn = 1, jpidta/ifreq+1 162 IF( cp_cfg == "orca" .AND. cp_cfz == "antarctic" ) THEN 163 ! 164 ! 1.2 Modify ahm 165 ! -------------- 166 IF(lwp)WRITE(numout,*) ' inildf: Antarctic ocean' 167 IF(lwp)WRITE(numout,*) ' no tropics, no reduction of ahm' 168 IF(lwp)WRITE(numout,*) ' north boundary increase' 169 170 ahm1(:,:) = ahm0 171 ahm2(:,:) = ahm0 172 173 ijpt0=max(1,min(49 -njmpp+1,jpj)) 174 ijpt1=max(0,min(49-njmpp+1,jpj-1)) 175 DO jj=ijpt0,ijpt1 176 ahm2(:,jj)=ahm0*2. 177 ahm1(:,jj)=ahm0*2. 178 END DO 179 ijpt0=max(1,min(48 -njmpp+1,jpj)) 180 ijpt1=max(0,min(48-njmpp+1,jpj-1)) 181 DO jj=ijpt0,ijpt1 182 ahm2(:,jj)=ahm0*1.9 183 ahm1(:,jj)=ahm0*1.75 184 END DO 185 ijpt0=max(1,min(47 -njmpp+1,jpj)) 186 ijpt1=max(0,min(47-njmpp+1,jpj-1)) 187 DO jj=ijpt0,ijpt1 188 ahm2(:,jj)=ahm0*1.5 189 ahm1(:,jj)=ahm0*1.25 190 END DO 191 ijpt0=max(1,min(46 -njmpp+1,jpj)) 192 ijpt1=max(0,min(46-njmpp+1,jpj-1)) 193 DO jj=ijpt0,ijpt1 194 ahm2(:,jj)=ahm0*1.1 195 END DO 196 197 ELSE IF( cp_cfg == "orca" .AND. cp_cfz == "arctic" ) THEN 198 ! 1.2 Modify ahm 199 ! -------------- 200 IF(lwp)WRITE(numout,*) ' inildf: Arctic ocean' 201 IF(lwp)WRITE(numout,*) ' no tropics, no reduction of ahm' 202 IF(lwp)WRITE(numout,*) ' south and west boundary increase' 203 204 205 ahm1(:,:) = ahm0 206 ahm2(:,:) = ahm0 207 208 ijpt0=max(1,min(98-jpjzoom+1-njmpp+1,jpj)) 209 ijpt1=max(0,min(98-jpjzoom+1-njmpp+1,jpj-1)) 210 DO jj=ijpt0,ijpt1 211 ahm2(:,jj)=ahm0*2. 212 ahm1(:,jj)=ahm0*2. 213 END DO 214 ijpt0=max(1,min(99-jpjzoom+1-njmpp+1,jpj)) 215 ijpt1=max(0,min(99-jpjzoom+1-njmpp+1,jpj-1)) 216 DO jj=ijpt0,ijpt1 217 ahm2(:,jj)=ahm0*1.9 218 ahm1(:,jj)=ahm0*1.75 219 END DO 220 ijpt0=max(1,min(100-jpjzoom+1-njmpp+1,jpj)) 221 ijpt1=max(0,min(100-jpjzoom+1-njmpp+1,jpj-1)) 222 DO jj=ijpt0,ijpt1 223 ahm2(:,jj)=ahm0*1.5 224 ahm1(:,jj)=ahm0*1.25 225 END DO 226 ijpt0=max(1,min(101-jpjzoom+1-njmpp+1,jpj)) 227 ijpt1=max(0,min(101-jpjzoom+1-njmpp+1,jpj-1)) 228 DO jj=ijpt0,ijpt1 229 ahm2(:,jj)=ahm0*1.1 230 END DO 231 ELSE 232 ! Read 2d integer array to specify western boundary increase in the 233 ! ===================== equatorial strip (20N-20S) defined at t-points 234 235 CALL ctl_opn( inum, 'ahmcoef', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp ) 236 READ(inum,9101) clexp, iim, ijm 176 237 READ(inum,'(/)') 177 il2 = MIN( jpidta, il1+ifreq-1 ) 178 READ(inum,9201) ( ii, ji = il1, il2, 5 ) 179 READ(inum,'(/)') 180 DO jj = jpjdta, 1, -1 181 READ(inum,9202) ij, ( idata(ji,jj), ji = il1, il2 ) 182 END DO 183 il1 = il1 + ifreq 184 END DO 185 186 DO jj = 1, nlcj 187 DO ji = 1, nlci 188 icof(ji,jj) = idata( mig(ji), mjg(jj) ) 189 END DO 190 END DO 191 DO jj = nlcj+1, jpj 192 DO ji = 1, nlci 193 icof(ji,jj) = icof(ji,nlcj) 194 END DO 195 END DO 196 DO jj = 1, jpj 197 DO ji = nlci+1, jpi 198 icof(ji,jj) = icof(nlci,jj) 199 END DO 200 END DO 201 202 9101 FORMAT(1x,a15,2i8) 203 9201 FORMAT(3x,13(i3,12x)) 204 9202 FORMAT(i3,41i3) 205 206 207 ! Set ahm1 and ahm2 ( T- and F- points) (used for laplacian operator) 208 ! ================= 209 ! define ahm1 and ahm2 at the right grid point position 210 ! (USER: modify ahm1 and ahm2 following your desiderata) 211 212 213 ! Decrease ahm to zahmeq m2/s in the tropics 214 ! (from 90 to 20 degre: ahm = constant 215 ! from 20 to 2.5 degre: ahm = decrease in (1-cos)/2 216 ! from 2.5 to 0 degre: ahm = constant 217 ! symmetric in the south hemisphere) 218 219 zahmeq = aht0 220 221 DO jj = 1, jpj 222 DO ji = 1, jpi 223 IF( ABS( gphif(ji,jj) ) >= 20. ) THEN 224 ahm2(ji,jj) = ahm0 225 ELSEIF( ABS( gphif(ji,jj) ) <= 2.5 ) THEN 226 ahm2(ji,jj) = zahmeq 227 ELSE 228 ahm2(ji,jj) = zahmeq + (ahm0-zahmeq)/2. & 229 * ( 1. - COS( rad * ( ABS(gphif(ji,jj))-2.5 ) * 180. / 17.5 ) ) 230 ENDIF 231 IF( ABS( gphit(ji,jj) ) >= 20. ) THEN 232 ahm1(ji,jj) = ahm0 233 ELSEIF( ABS( gphit(ji,jj) ) <= 2.5 ) THEN 234 ahm1(ji,jj) = zahmeq 235 ELSE 236 ahm1(ji,jj) = zahmeq + (ahm0-zahmeq)/2. & 237 * ( 1. - COS( rad * ( ABS(gphit(ji,jj))-2.5 ) * 180. / 17.5 ) ) 238 ENDIF 239 END DO 240 END DO 241 242 ! increase along western boundaries of equatorial strip 243 ! t-point 244 DO jj = 1, jpjm1 245 DO ji = 1, jpim1 246 zcoft = FLOAT( icof(ji,jj) ) / 100. 247 ahm1(ji,jj) = zcoft * ahm0 + (1.-zcoft) * ahm1(ji,jj) 248 END DO 249 END DO 250 ! f-point 251 icof(:,:) = icof(:,:) * tmask(:,:,1) 252 DO jj = 1, jpjm1 253 DO ji = 1, jpim1 ! NO vector opt. 254 zmsk = tmask(ji,jj+1,1) + tmask(ji+1,jj+1,1) + tmask(ji,jj,1) + tmask(ji,jj+1,1) 255 IF( zmsk == 0. ) THEN 256 zcoff = 1. 257 ELSE 258 zcoff = FLOAT( icof(ji,jj+1) + icof(ji+1,jj+1) + icof(ji,jj) + icof(ji,jj+1) ) & 238 ifreq = 40 239 il1 = 1 240 DO jn = 1, jpidta/ifreq+1 241 READ(inum,'(/)') 242 il2 = MIN( jpidta, il1+ifreq-1 ) 243 READ(inum,9201) ( ii, ji = il1, il2, 5 ) 244 READ(inum,'(/)') 245 DO jj = jpjdta, 1, -1 246 READ(inum,9202) ij, ( idata(ji,jj), ji = il1, il2 ) 247 END DO 248 il1 = il1 + ifreq 249 END DO 250 251 DO jj = 1, nlcj 252 DO ji = 1, nlci 253 icof(ji,jj) = idata( mig(ji), mjg(jj) ) 254 END DO 255 END DO 256 DO jj = nlcj+1, jpj 257 DO ji = 1, nlci 258 icof(ji,jj) = icof(ji,nlcj) 259 END DO 260 END DO 261 DO jj = 1, jpj 262 DO ji = nlci+1, jpi 263 icof(ji,jj) = icof(nlci,jj) 264 END DO 265 END DO 266 267 9101 FORMAT(1x,a15,2i8) 268 9201 FORMAT(3x,13(i3,12x)) 269 9202 FORMAT(i3,41i3) 270 271 272 ! Set ahm1 and ahm2 ( T- and F- points) (used for laplacian operator) 273 ! ================= 274 ! define ahm1 and ahm2 at the right grid point position 275 ! (USER: modify ahm1 and ahm2 following your desiderata) 276 277 278 ! Decrease ahm to zahmeq m2/s in the tropics 279 ! (from 90 to 20 degre: ahm = constant 280 ! from 20 to 2.5 degre: ahm = decrease in (1-cos)/2 281 ! from 2.5 to 0 degre: ahm = constant 282 ! symmetric in the south hemisphere) 283 284 zahmeq = aht0 285 286 DO jj = 1, jpj 287 DO ji = 1, jpi 288 IF( ABS( gphif(ji,jj) ) >= 20. ) THEN 289 ahm2(ji,jj) = ahm0 290 ELSEIF( ABS( gphif(ji,jj) ) <= 2.5 ) THEN 291 ahm2(ji,jj) = zahmeq 292 ELSE 293 ahm2(ji,jj) = zahmeq + (ahm0-zahmeq)/2. & 294 * ( 1. - COS( rad * ( ABS(gphif(ji,jj))-2.5 ) * 180. / 17.5 ) ) 295 ENDIF 296 IF( ABS( gphit(ji,jj) ) >= 20. ) THEN 297 ahm1(ji,jj) = ahm0 298 ELSEIF( ABS( gphit(ji,jj) ) <= 2.5 ) THEN 299 ahm1(ji,jj) = zahmeq 300 ELSE 301 ahm1(ji,jj) = zahmeq + (ahm0-zahmeq)/2. & 302 * ( 1. - COS( rad * ( ABS(gphit(ji,jj))-2.5 ) * 180. / 17.5 ) ) 303 ENDIF 304 END DO 305 END DO 306 307 ! increase along western boundaries of equatorial strip 308 ! t-point 309 DO jj = 1, jpjm1 310 DO ji = 1, jpim1 311 zcoft = FLOAT( icof(ji,jj) ) / 100. 312 ahm1(ji,jj) = zcoft * ahm0 + (1.-zcoft) * ahm1(ji,jj) 313 END DO 314 END DO 315 ! f-point 316 icof(:,:) = icof(:,:) * tmask(:,:,1) 317 DO jj = 1, jpjm1 318 DO ji = 1, jpim1 ! NO vector opt. 319 zmsk = tmask(ji,jj+1,1) + tmask(ji+1,jj+1,1) + tmask(ji,jj,1) + tmask(ji,jj+1,1) 320 IF( zmsk == 0. ) THEN 321 zcoff = 1. 322 ELSE 323 zcoff = FLOAT( icof(ji,jj+1) + icof(ji+1,jj+1) + icof(ji,jj) + icof(ji,jj+1) ) & 259 324 / (zmsk * 100.) 260 ENDIF261 ahm2(ji,jj) = zcoff * ahm0 + (1.-zcoff) * ahm2(ji,jj)262 END DO263 END DO264 #endif 325 ENDIF 326 ahm2(ji,jj) = zcoff * ahm0 + (1.-zcoff) * ahm2(ji,jj) 327 END DO 328 END DO 329 ENDIF 265 330 266 331 ! Lateral boundary conditions on ( ahm1, ahm2 ) … … 323 388 IF(lwp) WRITE(numout,*) ' orca_r1 configuration' 324 389 325 #if defined key_antarctic 326 # include "ldfdyn_antarctic.h90" 327 #elif defined key_arctic 328 # include "ldfdyn_arctic.h90" 329 #else 330 ! Read 2d integer array to specify western boundary increase in the 331 ! ===================== equatorial strip (20N-20S) defined at t-points 332 333 CALL ctl_opn( inum, 'ahmcoef', 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', & 334 & 1, numout, lwp ) 335 REWIND inum 336 READ(inum,9101) clexp, iim, ijm 337 READ(inum,'(/)') 338 ifreq = 40 339 il1 = 1 340 DO jn = 1, jpidta/ifreq+1 390 IF( cp_cfg == "orca" .AND. cp_cfz == "antarctic" ) THEN 391 ! 392 ! 1.2 Modify ahm 393 ! -------------- 394 IF(lwp)WRITE(numout,*) ' inildf: Antarctic ocean' 395 IF(lwp)WRITE(numout,*) ' no tropics, no reduction of ahm' 396 IF(lwp)WRITE(numout,*) ' north boundary increase' 397 398 ahm1(:,:) = ahm0 399 ahm2(:,:) = ahm0 400 401 ijpt0=max(1,min(49 -njmpp+1,jpj)) 402 ijpt1=max(0,min(49-njmpp+1,jpj-1)) 403 DO jj=ijpt0,ijpt1 404 ahm2(:,jj)=ahm0*2. 405 ahm1(:,jj)=ahm0*2. 406 END DO 407 ijpt0=max(1,min(48 -njmpp+1,jpj)) 408 ijpt1=max(0,min(48-njmpp+1,jpj-1)) 409 DO jj=ijpt0,ijpt1 410 ahm2(:,jj)=ahm0*1.9 411 ahm1(:,jj)=ahm0*1.75 412 END DO 413 ijpt0=max(1,min(47 -njmpp+1,jpj)) 414 ijpt1=max(0,min(47-njmpp+1,jpj-1)) 415 DO jj=ijpt0,ijpt1 416 ahm2(:,jj)=ahm0*1.5 417 ahm1(:,jj)=ahm0*1.25 418 END DO 419 ijpt0=max(1,min(46 -njmpp+1,jpj)) 420 ijpt1=max(0,min(46-njmpp+1,jpj-1)) 421 DO jj=ijpt0,ijpt1 422 ahm2(:,jj)=ahm0*1.1 423 END DO 424 425 ELSE IF( cp_cfg == "orca" .AND. cp_cfz == "arctic" ) THEN 426 ! 1.2 Modify ahm 427 ! -------------- 428 IF(lwp)WRITE(numout,*) ' inildf: Arctic ocean' 429 IF(lwp)WRITE(numout,*) ' no tropics, no reduction of ahm' 430 IF(lwp)WRITE(numout,*) ' south and west boundary increase' 431 432 433 ahm1(:,:) = ahm0 434 ahm2(:,:) = ahm0 435 436 ijpt0=max(1,min(98-jpjzoom+1-njmpp+1,jpj)) 437 ijpt1=max(0,min(98-jpjzoom+1-njmpp+1,jpj-1)) 438 DO jj=ijpt0,ijpt1 439 ahm2(:,jj)=ahm0*2. 440 ahm1(:,jj)=ahm0*2. 441 END DO 442 ijpt0=max(1,min(99-jpjzoom+1-njmpp+1,jpj)) 443 ijpt1=max(0,min(99-jpjzoom+1-njmpp+1,jpj-1)) 444 DO jj=ijpt0,ijpt1 445 ahm2(:,jj)=ahm0*1.9 446 ahm1(:,jj)=ahm0*1.75 447 END DO 448 ijpt0=max(1,min(100-jpjzoom+1-njmpp+1,jpj)) 449 ijpt1=max(0,min(100-jpjzoom+1-njmpp+1,jpj-1)) 450 DO jj=ijpt0,ijpt1 451 ahm2(:,jj)=ahm0*1.5 452 ahm1(:,jj)=ahm0*1.25 453 END DO 454 ijpt0=max(1,min(101-jpjzoom+1-njmpp+1,jpj)) 455 ijpt1=max(0,min(101-jpjzoom+1-njmpp+1,jpj-1)) 456 DO jj=ijpt0,ijpt1 457 ahm2(:,jj)=ahm0*1.1 458 END DO 459 ELSE 460 461 ! Read 2d integer array to specify western boundary increase in the 462 ! ===================== equatorial strip (20N-20S) defined at t-points 463 464 CALL ctl_opn( inum, 'ahmcoef', 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', & 465 & 1, numout, lwp ) 466 REWIND inum 467 READ(inum,9101) clexp, iim, ijm 341 468 READ(inum,'(/)') 342 il2 = MIN( jpidta, il1+ifreq-1 ) 343 READ(inum,9201) ( ii, ji = il1, il2, 5 ) 344 READ(inum,'(/)') 345 DO jj = jpjdta, 1, -1 346 READ(inum,9202) ij, ( idata(ji,jj), ji = il1, il2 ) 347 END DO 348 il1 = il1 + ifreq 349 END DO 350 351 DO jj = 1, nlcj 352 DO ji = 1, nlci 353 icof(ji,jj) = idata( mig(ji), mjg(jj) ) 354 END DO 355 END DO 356 DO jj = nlcj+1, jpj 357 DO ji = 1, nlci 358 icof(ji,jj) = icof(ji,nlcj) 359 END DO 360 END DO 361 DO jj = 1, jpj 362 DO ji = nlci+1, jpi 363 icof(ji,jj) = icof(nlci,jj) 364 END DO 365 END DO 366 367 9101 FORMAT(1x,a15,2i8) 368 9201 FORMAT(3x,13(i3,12x)) 369 9202 FORMAT(i3,41i3) 370 371 372 ! Set ahm1 and ahm2 ( T- and F- points) (used for laplacian operator) 373 ! ================= 374 ! define ahm1 and ahm2 at the right grid point position 375 ! (USER: modify ahm1 and ahm2 following your desiderata) 376 377 378 ! Decrease ahm to zahmeq m2/s in the tropics 379 ! (from 90 to 20 degrees: ahm = scaled by local metrics 380 ! from 20 to 2.5 degrees: ahm = decrease in (1-cos)/2 381 ! from 2.5 to 0 degrees: ahm = constant 382 ! symmetric in the south hemisphere) 383 384 zahmeq = aht0 385 zam20s = ahm0*COS( rad * 20. ) 386 387 DO jj = 1, jpj 388 DO ji = 1, jpi 389 IF( ABS( gphif(ji,jj) ) >= 20. ) THEN 390 ! leave as set in ldf_dyn_c2d 391 ELSEIF( ABS( gphif(ji,jj) ) <= 2.5 ) THEN 392 ahm2(ji,jj) = zahmeq 393 ELSE 394 ahm2(ji,jj) = zahmeq + (zam20s-zahmeq)/2. & 395 * ( 1. - COS( rad * ( ABS(gphif(ji,jj))-2.5 ) * 180. / 17.5 ) ) 396 ENDIF 397 IF( ABS( gphit(ji,jj) ) >= 20. ) THEN 398 ! leave as set in ldf_dyn_c2d 399 ELSEIF( ABS( gphit(ji,jj) ) <= 2.5 ) THEN 400 ahm1(ji,jj) = zahmeq 401 ELSE 402 ahm1(ji,jj) = zahmeq + (zam20s-zahmeq)/2. & 403 * ( 1. - COS( rad * ( ABS(gphit(ji,jj))-2.5 ) * 180. / 17.5 ) ) 404 ENDIF 405 END DO 406 END DO 407 408 ! increase along western boundaries of equatorial strip 409 ! t-point 410 DO jj = 1, jpjm1 411 DO ji = 1, jpim1 412 IF( ABS( gphit(ji,jj) ) < 20. ) THEN 413 zcoft = FLOAT( icof(ji,jj) ) / 100. 414 ahm1(ji,jj) = zcoft * ahm0 + (1.-zcoft) * ahm1(ji,jj) 415 ENDIF 416 END DO 417 END DO 418 ! f-point 419 icof(:,:) = icof(:,:) * tmask(:,:,1) 420 DO jj = 1, jpjm1 421 DO ji = 1, jpim1 422 IF( ABS( gphif(ji,jj) ) < 20. ) THEN 423 zmsk = tmask(ji,jj+1,1) + tmask(ji+1,jj+1,1) + tmask(ji,jj,1) + tmask(ji,jj+1,1) 424 IF( zmsk == 0. ) THEN 425 zcoff = 1. 426 ELSE 427 zcoff = FLOAT( icof(ji,jj+1) + icof(ji+1,jj+1) + icof(ji,jj) + icof(ji,jj+1) ) & 428 / (zmsk * 100.) 429 ENDIF 430 ahm2(ji,jj) = zcoff * ahm0 + (1.-zcoff) * ahm2(ji,jj) 431 ENDIF 432 END DO 433 END DO 434 #endif 469 ifreq = 40 470 il1 = 1 471 DO jn = 1, jpidta/ifreq+1 472 READ(inum,'(/)') 473 il2 = MIN( jpidta, il1+ifreq-1 ) 474 READ(inum,9201) ( ii, ji = il1, il2, 5 ) 475 READ(inum,'(/)') 476 DO jj = jpjdta, 1, -1 477 READ(inum,9202) ij, ( idata(ji,jj), ji = il1, il2 ) 478 END DO 479 il1 = il1 + ifreq 480 END DO 481 482 DO jj = 1, nlcj 483 DO ji = 1, nlci 484 icof(ji,jj) = idata( mig(ji), mjg(jj) ) 485 END DO 486 END DO 487 DO jj = nlcj+1, jpj 488 DO ji = 1, nlci 489 icof(ji,jj) = icof(ji,nlcj) 490 END DO 491 END DO 492 DO jj = 1, jpj 493 DO ji = nlci+1, jpi 494 icof(ji,jj) = icof(nlci,jj) 495 END DO 496 END DO 497 498 9101 FORMAT(1x,a15,2i8) 499 9201 FORMAT(3x,13(i3,12x)) 500 9202 FORMAT(i3,41i3) 501 502 503 ! Set ahm1 and ahm2 ( T- and F- points) (used for laplacian operator) 504 ! ================= 505 ! define ahm1 and ahm2 at the right grid point position 506 ! (USER: modify ahm1 and ahm2 following your desiderata) 507 508 509 ! Decrease ahm to zahmeq m2/s in the tropics 510 ! (from 90 to 20 degrees: ahm = scaled by local metrics 511 ! from 20 to 2.5 degrees: ahm = decrease in (1-cos)/2 512 ! from 2.5 to 0 degrees: ahm = constant 513 ! symmetric in the south hemisphere) 514 515 zahmeq = aht0 516 zam20s = ahm0*COS( rad * 20. ) 517 518 DO jj = 1, jpj 519 DO ji = 1, jpi 520 IF( ABS( gphif(ji,jj) ) >= 20. ) THEN 521 ! leave as set in ldf_dyn_c2d 522 ELSEIF( ABS( gphif(ji,jj) ) <= 2.5 ) THEN 523 ahm2(ji,jj) = zahmeq 524 ELSE 525 ahm2(ji,jj) = zahmeq + (zam20s-zahmeq)/2. & 526 * ( 1. - COS( rad * ( ABS(gphif(ji,jj))-2.5 ) * 180. / 17.5 ) ) 527 ENDIF 528 IF( ABS( gphit(ji,jj) ) >= 20. ) THEN 529 ! leave as set in ldf_dyn_c2d 530 ELSEIF( ABS( gphit(ji,jj) ) <= 2.5 ) THEN 531 ahm1(ji,jj) = zahmeq 532 ELSE 533 ahm1(ji,jj) = zahmeq + (zam20s-zahmeq)/2. & 534 * ( 1. - COS( rad * ( ABS(gphit(ji,jj))-2.5 ) * 180. / 17.5 ) ) 535 ENDIF 536 END DO 537 END DO 538 539 ! increase along western boundaries of equatorial strip 540 ! t-point 541 DO jj = 1, jpjm1 542 DO ji = 1, jpim1 543 IF( ABS( gphit(ji,jj) ) < 20. ) THEN 544 zcoft = FLOAT( icof(ji,jj) ) / 100. 545 ahm1(ji,jj) = zcoft * ahm0 + (1.-zcoft) * ahm1(ji,jj) 546 ENDIF 547 END DO 548 END DO 549 ! f-point 550 icof(:,:) = icof(:,:) * tmask(:,:,1) 551 DO jj = 1, jpjm1 552 DO ji = 1, jpim1 553 IF( ABS( gphif(ji,jj) ) < 20. ) THEN 554 zmsk = tmask(ji,jj+1,1) + tmask(ji+1,jj+1,1) + tmask(ji,jj,1) + tmask(ji,jj+1,1) 555 IF( zmsk == 0. ) THEN 556 zcoff = 1. 557 ELSE 558 zcoff = FLOAT( icof(ji,jj+1) + icof(ji+1,jj+1) + icof(ji,jj) + icof(ji,jj+1) ) & 559 / (zmsk * 100.) 560 ENDIF 561 ahm2(ji,jj) = zcoff * ahm0 + (1.-zcoff) * ahm2(ji,jj) 562 ENDIF 563 END DO 564 END DO 565 ENDIF 435 566 436 567 ! Lateral boundary conditions on ( ahm1, ahm2 ) -
branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/TRA/tradmp.F90
r3901 r3989 303 303 304 304 ! ! ==================================================== 305 IF( lzoom_arct .AND. lzoom_anta ) THEN ! ORCA configuration : arctic zoomor antarctic zoom305 IF( cp_cfz == "arctic" . OR. cp_cfz == "antarctic" ) THEN ! ORCA configuration : arctic or antarctic zoom 306 306 ! ! ==================================================== 307 307 IF(lwp) WRITE(numout,*) 308 IF(lwp .AND. lzoom_arct) WRITE(numout,*) ' dtacof_zoom : ORCA Arctic zoom'309 IF(lwp .AND. lzoom_arct ) WRITE(numout,*) 'dtacof_zoom : ORCA Antarctic zoom'308 IF(lwp .AND. cp_cfz == "arctic" ) WRITE(numout,*) ' dtacof_zoom : ORCA Arctic zoom' 309 IF(lwp .AND. cp_cfz == "antarctic" ) WRITE(numout,*) ' dtacof_zoom : ORCA Antarctic zoom' 310 310 IF(lwp) WRITE(numout,*) 311 311 ! -
branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/nemogcm.F90
r3973 r3989 228 228 & nn_isplt, nn_jsplt, nn_jctls, nn_jctle, & 229 229 & nn_bench, nn_timing 230 NAMELIST/namcfg/ cp_cfg, jp_cfg, jpidta, jpjdta, jpkdta, jpiglo, jpjglo, &230 NAMELIST/namcfg/ cp_cfg, cp_cfz, jp_cfg, jpidta, jpjdta, jpkdta, jpiglo, jpjglo, & 231 231 & jpizoom, jpjzoom, jperio, jphgr_msh, & 232 232 & ppglam0, ppgphi0, ppe1_deg, ppe2_deg, ppe1_m, ppe2_m, & … … 482 482 WRITE(numout,*) ' Namelist namcfg' 483 483 WRITE(numout,*) ' configuration name cp_cfg = ', TRIM(cp_cfg) 484 WRITE(numout,*) ' configuration zoom name cp_cfz = ', TRIM(cp_cfz) 484 485 WRITE(numout,*) ' configuration resolution jp_cfg = ', jp_cfg 485 486 WRITE(numout,*) ' 1st lateral dimension ( >= jpi ) jpidta = ', jpidta -
branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/par_oce.F90
r3973 r3989 29 29 !!---------------------------------------------------------------------- 30 30 CHARACTER(lc) :: cp_cfg !: name of the configuration 31 CHARACTER(lc) :: cp_cfz !: name of the zoom of configuration 31 32 INTEGER :: jp_cfg !: resolution of the configuration 32 33
Note: See TracChangeset
for help on using the changeset viewer.