Changeset 3989


Ignore:
Timestamp:
2013-07-24T11:48:35+02:00 (7 years ago)
Author:
clevy
Message:

Configuration setting/Step3 and doc, see ticket:#1074

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  
    177177\newcommand{\rou} [1] {\textit{#1}\index{Routines!#1}}            %module (routine) 
    178178\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  
    180181\newcommand{\jp} [1] {\textit{#1}\index{Model parameters!#1}}        %model parameter (jp) 
    181182\newcommand{\pp} [1] {\textit{#1}\index{Model parameters!#1}}        %namelist parameter (pp) 
     
    245246% ================================================================ 
    246247% ================================================================ 
     248 
    247249\begin{document} 
    248250 
     
    291293\include{./TexFiles/Chapters/Chap_ZDF}       % Vertical diffusion 
    292294 
    293 \include{./TexFiles/Chapters/Chap_DIA}       % Miscellaneous topics 
     295\include{./TexFiles/Chapters/Chap_DIA}       % Outputs and Diagnostics 
    294296 
    295297\include{./TexFiles/Chapters/Chap_OBS}                  % Observation operator 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Annex_ISO.tex

    r3297 r3989  
    77\minitoc 
    88\pagebreak 
    9 \section{Choice of namelist parameters} 
     9\section{Choice of \ngn{namtra\_ldf} namelist parameters} 
    1010%-----------------------------------------nam_traldf------------------------------------------------------ 
    1111\namdisplay{namtra_ldf} 
     
    2626\np{rn\_aeiv\_0}. If 2D-varying coefficients are set with 
    2727\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}, where 
     28scale factor according to \eqref{Eq_title} \footnote{Except in global ORCA 
     29  $0.5^{\circ}$ runs with \key{traldf\_eiv}, where 
    3030  $A_l$ is set like $A_e$ but with a minimum vale of 
    3131  $100\;\mathrm{m}^2\;\mathrm{s}^{-1}$}. In idealised setups with 
    3232\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 
     33is set in the global configurations with \key{traldf\_c2d}, a horizontally varying $A_e$ is 
    3534instead set from the Held-Larichev parameterisation\footnote{In this 
    3635  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  
    1818assimilation code.  The code can also output model background fields which are used 
    1919as 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 options 
     20\textit{\ngn{nam\_asminc} }.  There is a brief description of all the namelist options 
    2121provided.  To build the ASM code \key{asminc} must be set. 
    2222 
     
    125125\label{ASM_details} 
    126126 
    127 Here we show an example namelist and the header of an example assimilation  
     127Here we show an example \ngn{namasm} namelist and the header of an example assimilation  
    128128increments file on the ORCA2 grid. 
    129129 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Chap_CFG.tex

    r3764 r3989  
    11% ================================================================ 
    2 % Chapter Ñ Configurations 
     2% Chapter Configurations 
    33% ================================================================ 
    44\chapter{Configurations} 
     
    1616 
    1717 
    18 The purpose of this part of the manual is to introduce the \NEMO predefined configuration.  
     18The purpose of this part of the manual is to introduce the \NEMO reference configurations.  
    1919These configurations are offered as means to explore various numerical and physical options,  
    2020thus allowing the user to verify that the code is performing in a manner consistent with that  
    2121we are running. This form of verification is critical as one adopts the code for his or her particular  
    2222research 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  
     23in the code, though by no means are all options exercised in the reference configurations. 
     24 
     25Configuration is defined mainly through the \ngn{namcfg} namelist variables: 
     26%------------------------------------------namcfg---------------------------------------------------- 
     27\namdisplay{namcfg} 
     28%------------------------------------------------------------------------------------------------------------- 
    3129 
    3230% ================================================================ 
    3331% 1D model configuration 
    3432% ================================================================ 
    35 \section{Water column model: 1D model (C1D) (\key{c1d})} 
     33\section{Water column model: 1D model (C1D) (\key{c1d}) } 
    3634\label{CFG_c1d} 
    3735 
    3836The 1D model option simulates a stand alone water column within the 3D \NEMO system.  
    3937It 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.  
     38or 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} ).  
    4140The 1D model is a very useful tool   
    4241\textit{(a)} to learn about the physics and numerical treatment of vertical mixing processes ;  
     
    4847 
    4948The methodology is based on the use of the zoom functionality over the smallest possible  
    50 domain : a 3x3 domain centred on the grid point of interest (see \S\ref{MISC_zoom}),  
     49domain : a 3x3 domain centered on the grid point of interest,  
    5150with some extra routines. There is no need to define a new mesh, bathymetry,  
    5251initial 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}  
     52The chosen grid point is set in \textit{\ngn{namcfg}} namelist by setting the \np{jpizoom} and \np{jpjzoom}  
    5453parameters to the indices of the location of the chosen grid point. 
    5554 
     
    7675% ORCA family configurations 
    7776% ================================================================ 
    78 \section{ORCA family: global ocean with tripolar grid (\key{orca\_rX})} 
     77\section{ORCA family: global ocean with tripolar grid } 
    7978\label{CFG_orca} 
    8079 
     
    8281the LIM sea-ice model (ORCA-LIM) and possibly with PISCES biogeochemical model  
    8382(ORCA-LIM-PISCES), using various resolutions. 
     83The appropriate \textit{\&namcfg} namelist is available in \textit{CONFIG/ORCA2\_LIM/EXP00/namelist\_cfg}  
     84for ORCA2 and in \textit{CONFIG/SHARED/README\_other\_configurations\_namelist\_namcfg}  
     85for other resolutions 
    8486 
    8587 
     
    147149The NEMO system is provided with five built-in ORCA configurations which differ in the  
    148150horizontal 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 set the grid size and configuration  
    151 name parameters  (Tab.~\ref{Tab_ORCA}). 
     151expressed in degrees. Each of configuration is set through the \textit{\ngn{namcfg}} namelist,  
     152which sets the grid size and configuration  
     153name parameters  (Tab. \ref{Tab_ORCA}). 
    152154. 
    153155 
     
    155157\begin{table}[!t]     \begin{center} 
    156158\begin{tabular}{p{4cm} c c c c} 
    157 CPP key                        & \jp{jp\_cfg} &  \jp{jpiglo} & \jp{jpiglo} &       \\   
     159Horizontal Grid                         & \np{jp\_cfg} &  \np{jpiglo} & \np{jpjglo} &       \\   
    158160\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     &        \\ 
    164166%\key{orca\_r8}       &        8         &      2882    &   2042     &        \\ 
    165167%\key{orca\_r12}     &      12         &      4322    &   3062      &       \\ 
     
    168170\caption{ \label{Tab_ORCA}    
    169171Set 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). } 
     172In all cases, the name of the configuration is set to "orca" ($i.e.$ \np{cp\_cfg}~=~orca). } 
    171173\end{center} 
    172174\end{table} 
     
    197199in the upper 150m (see Tab.~\ref{Tab_orca_zgr} and Fig.~\ref{Fig_zgr}).  
    198200The bottom topography and the coastlines are derived from the global atlas of Smith and Sandwell (1997).  
    199 The default forcing employ the boundary forcing from \citet{Large_Yeager_Rep04} (see \S\ref{SBC_blk_core}),  
     201The default forcing uses the boundary forcing from \citet{Large_Yeager_Rep04} (see \S\ref{SBC_blk_core}),  
    200202which was developed for the purpose of running global coupled ocean-ice simulations  
    201203without an interactive atmosphere. This \citet{Large_Yeager_Rep04} dataset is available  
     
    205207 
    206208ORCA\_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},  
     209current area ( \key{agrif}  defined) and,  by setting the appropriate variables in  
     210\textit{\&namcfg}, see \textit{CONFIG/SHARED/namelist\_ref} 
    208211a regional Arctic or peri-Antarctic configuration is extracted from an ORCA\_R2 or R05 configurations 
    209212using sponge layers at open boundaries.  
     
    212215%       GYRE family: double gyre basin 
    213216% ------------------------------------------------------------------------------------------------------------- 
    214 \section{GYRE family: double gyre basin (\key{gyre})} 
     217\section{GYRE family: double gyre basin } 
    215218\label{CFG_gyre} 
    216219 
    217 The GYRE configuration \citep{Levy_al_OM10} have been built to simulated  
    218 the seasonal cycle of a double-gyre box model. It consist in an idealized domain  
     220The GYRE configuration \citep{Levy_al_OM10} has been built to simulate 
     221the seasonal cycle of a double-gyre box model. It consists in an idealized domain  
    219222similar to that used in the studies of \citet{Drijfhout_JPO94} and \citet{Hazeleger_Drijfhout_JPO98,  
    220223Hazeleger_Drijfhout_JPO99, Hazeleger_Drijfhout_JGR00, Hazeleger_Drijfhout_JPO00},  
     
    242245uniformly applied to the whole domain. 
    243246 
    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}). 
     247The GYRE configuration is set through the \textit{\&namcfg} namelist defined in the reference  
     248configuration \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   \\ 
     252Obviously, the namelist parameters have to be adjusted to the chosen resolution, see the Configurations  
     253pages on the NEMO web site (Using NEMO\/Configurations) . 
     254In the vertical, GYRE uses the default 30 ocean levels (\pp{jpk}=31) (Fig.~\ref{Fig_zgr}). 
    250255 
    251256The GYRE configuration is also used in benchmark test as it is very simple to increase  
     
    270275 
    271276\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 
    275280\end{description} 
    276  
     281The appropriate \textit{\&namcfg} namelists are available in   
     282\textit{CONFIG/SHARED/README\_other\_configurations\_namelist\_namcfg} 
    277283% ------------------------------------------------------------------------------------------------------------- 
    278284%       AMM configuration 
    279285% ------------------------------------------------------------------------------------------------------------- 
    280 \section{AMM: atlantic margin configuration (\key{amm\_12km})} 
     286\section{AMM: atlantic margin configuration } 
    281287\label{MISC_config_AMM} 
    282288 
    283289The AMM, Atlantic Margins Model, is a regional model covering the 
    284290Northwest 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. 
     291approximately 12km horizontal resolution. The appropriate  
     292\textit{\&namcfg} namelist  is available in \textit{CONFIG/AMM12/EXP00/namelist\_cfg}. 
     293It is used to build the correct dimensions of the AMM domain. 
    287294 
    288295This configuration tests several features of NEMO functionality specific 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Chap_DIA.tex

    r3764 r3989  
    11% ================================================================ 
    2 % Chapter Ñ I/O & Diagnostics 
     2% Chapter I/O & Diagnostics 
    33% ================================================================ 
    44\chapter{Ouput and Diagnostics (IOM, DIA, TRD, FLO)} 
     
    609609set via an equivalent and identically named namelist to \textit{namnc4}  
    610610in \np{xmlio\_server.def}. Typically this namelist serves the mean files 
    611 whilst the \np{ namnc4} in the main namelist file continues to serve the 
     611whilst the \ngn{ namnc4} in the main namelist file continues to serve the 
    612612restart files. This duplication is unfortunate but appropriate since, if 
    613613using io\_servers, the domain sizes of the individual files produced by the 
     
    631631trend of the dynamics and/or temperature and salinity time evolution equations  
    632632is 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  
     633of each $dyn\cdots.F90$ and/or $tra\cdots.F90$ routines). Options are defined by 
     634\ngn{namtrd} namelist variables. These trends are then  
    634635used in \mdl{trdmod} (see TRD directory) every \textit{nn\_trd } time-steps. 
    635636 
     
    675676The on-line computation of floats advected either by the three dimensional velocity  
    676677field 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  
     678introduced in the system during the CLIPPER project. Options are defined by \ngn{namflo} 
     679namelis variables. The algorithm used is based  
    678680either on the work of \cite{Blanke_Raynaud_JPO97} (default option), or on a $4^th$ 
    679681Runge-Hutta algorithm (\np{ln\_flork4}=true). Note that the \cite{Blanke_Raynaud_JPO97}  
     
    687689 
    688690 
    689 In case of Ariane convention, input filename is \np{"init\_float\_ariane"}. Its format is: 
     691In case of Ariane convention, input filename is \np{init\_float\_ariane}. Its format is: 
    690692 
    691693\texttt{ I J K nisobfl itrash itrash } 
     
    709711 
    710712 
    711 In the other case ( longitude and latitude ), input filename is \np{"init\_float"}. Its format is: 
     713In the other case ( longitude and latitude ), input filename is init\_float. Its format is: 
    712714 
    713715\texttt{ Long Lat depth nisobfl ngrpfl itrash} 
     
    732734 
    733735\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} 
     736When initial positions are read in a restart file ( \np{ln\_rstflo}= .TRUE. ),  \np{jpnflnewflo} 
    735737can be added in the initialization file.  
    736738 
     
    740742is the frequency of creation of the float restart file. 
    741743 
    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 specification  
     744Output data can be written in ascii files (\np{ln\_flo\_ascii} = .TRUE. ). In that case,  
     745output filename is trajec\_float. 
     746 
     747Another 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  
    748750   to put in files description section: 
    749751 
     
    768770 
    769771 
    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. 
    771773 
    772774 
     
    789791%---------------------------------------------------------------------------------------------------------- 
    790792 
    791 Concerning the on-line Harmonic analysis, some parameters are available in namelist: 
     793Concerning the on-line Harmonic analysis, some parameters are available in namelist 
     794\ngn{namdia\_harm} : 
    792795 
    793796- \texttt{nit000\_han} is the first time step used for harmonic analysis 
     
    841844 
    842845 
    843 Namelist parameters control how frequently the flows are summed and the time scales over which 
    844  they are averaged, as well as the level of output for debugging: 
     846Namelist variables in \ngn{namdct} control how frequently the flows are summed  
     847and the time scales over which they are averaged, as well as the level of output for debugging: 
    845848 
    846849%------------------------------------------namdct---------------------------------------------------- 
     
    992995 
    993996The 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 setting  
    995 \np{ln\_diaptr} to true (see the \textit{namptr} namelist below).   
     997the meriodional stream function can be computed on-line in \mdl{diaptr}  
     998\np{ln\_diaptr} to true (see the \textit{\ngn{namptr} } namelist below).   
    996999When \np{ln\_subbas}~=~true, transports and stream function are computed  
    9971000for 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  
    11% ================================================================ 
    2 % Chapter 2 Ñ Space and Time Domain (DOM) 
     2% Chapter 2 Space and Time Domain (DOM) 
    33% ================================================================ 
    44\chapter{Space Domain (DOM) } 
     
    2424directory routines as well as the DOM (DOMain) directory.  
    2525 
    26 $\ $\newline    % force a new ligne 
     26$\ $\newline    % force a new lign 
    2727 
    2828% ================================================================ 
     
    274274\label{DOM_size} 
    275275 
    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$  
     276The 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  
     278given as namelist variables in the \ngn{namcfg} namelist.  
     279 
     280Note that are other namelist variables in the \ngn{namcfg} namelist that refer to 
     281 the domain size.  
     282The two variables \np{jpidta} and \np{jpjdta} may be larger than \np{jpiglo}, \np{jpjglo} 
    294283when the user wants to use only a sub-region of a given configuration. This is  
    295284the "zoom" capability described in \S\ref{MISC_zoom}. In most applications of  
     
    300289 
    301290 
    302 $\ $\newline    % force a new ligne 
     291$\ $\newline    % force a new lign 
    303292 
    304293% ================================================================ 
     
    388377 
    389378The user has three options available in defining a horizontal grid, which involve  
    390 the parameter $jphgr\_mesh$ of the \mdl{par\_oce} module.  
     379the namelist variable \np{jphgr\_mesh} of the \ngn{namcfg} namelist.  
    391380\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. 
    393382The coordinates and their first derivatives with respect to $i$ and $j$ are provided 
    394383in 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).  
    396385For other analytical grids, the \mdl{domhgr} module must be modified by the user.  
    397386\end{description} 
    398387 
    399388There 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,  
     390with grid sizes specified by parameters \np{ppe1\_deg} and \np{ppe2\_deg},  
    402391respectively. Such a geographical grid can be very anisotropic at high latitudes  
    403392because of the convergence of meridians (the zonal scale factors $e_1$  
    404393become much smaller than the meridional scale factors $e_2$). The Mercator  
    405 grid (\jp{jphgr\_mesh}=4) avoids this anisotropy by refining the meridional scale  
     394grid (\np{jphgr\_mesh}=4) avoids this anisotropy by refining the meridional scale  
    406395factors in the same way as the zonal ones. In this case, meridional scale factors  
    407396and latitudes are calculated analytically using the formulae appropriate for  
    408 a Mercator projection, based on \pp{ppe1\_deg} which is a reference grid spacing  
     397a Mercator projection, based on \np{ppe1\_deg} which is a reference grid spacing  
    409398at the equator (this applies even when the geographical equator is situated outside  
    410399the model domain).  
     
    412401\gmcomment{ give here the analytical expression of the Mercator mesh} 
    413402%%% 
    414 In these two cases (\jp{jphgr\_mesh}=1 or 4), the grid position is defined by the  
    415 longitude and latitude of the south-westernmost point (\pp{ppglamt0}  
    416 and \pp{ppgphi0}). Note that for the Mercator grid the user need only provide  
     403In these two cases (\np{jphgr\_mesh}=1 or 4), the grid position is defined by the  
     404longitude and latitude of the south-westernmost point (\np{ppglamt0}  
     405and \np{ppgphi0}). Note that for the Mercator grid the user need only provide  
    417406an approximate starting latitude: the real latitude will be recalculated analytically,  
    418407in order to ensure that the equator corresponds to line passing through $t$-  
     
    420409 
    421410Rectangular 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 factor  
     411\np{jphgr\_mesh} = 2, 3, 5. The domain is either an $f$-plane (\np{jphgr\_mesh} = 2,  
     412Coriolis factor is constant) or a beta-plane (\np{jphgr\_mesh} = 3, the Coriolis factor  
    424413is 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.  
     414and given by the parameters \np{ppe1\_m} and \np{ppe2\_m} respectively.  
    426415The zonal grid coordinate (\textit{glam} arrays) is in kilometers, starting at zero  
    427416with the first $t$-point. The meridional coordinate (gphi. arrays) is in kilometers,  
    428417and the second $t$-point corresponds to coordinate $gphit=0$. The input  
    429 parameter \pp{ppglam0} is ignored. \pp{ppgphi0} is used to set the reference  
     418variable \np{ppglam0} is ignored. \np{ppgphi0} is used to set the reference  
    430419latitude 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 case  
    432 \jp{jphgr\_mesh}=5 corresponds to a beta plane in a rotated domain for the  
     420\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  
    433422GYRE configuration, representing a classical mid-latitude double gyre system.  
    434423The rotation allows us to maximize the jet length relative to the gyre areas  
     
    436425 
    437426The choice of the grid must be consistent with the boundary conditions specified  
    438 by the parameter \jp{jperio} (see {\S\ref{LBC}). 
     427by the parameter \np{jperio} (see {\S\ref{LBC}). 
    439428 
    440429% ------------------------------------------------------------------------------------------------------------- 
     
    446435All the arrays relating to a particular ocean model configuration (grid-point  
    447436position, 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-line  
     437(namelist variable in \ngn{namdom}). This can be particularly useful for plots and off-line  
    449438diagnostics. In some cases, the user may choose to make a local modification  
    450439of a scale factor in the code. This is the case in global configurations when  
     
    454443the output grid written when $\np{nn\_msh} \not=0$ is no more equal to the input grid. 
    455444 
    456 $\ $\newline    % force a new ligne 
     445$\ $\newline    % force a new lign 
    457446 
    458447% ================================================================ 
     
    467456%------------------------------------------------------------------------------------------------------------- 
    468457 
     458Variables are defined through the \ngn{namzgr} and \ngn{namdom} namelists. 
    469459In the vertical, the model mesh is determined by four things:  
    470460(1) the bathymetry given in meters ;  
     
    553543\item[\np{nn\_bathy} = 0] a flat-bottom domain is defined. The total depth $z_w (jpk)$  
    554544is given by the coordinate transformation. The domain can either be a closed  
    555 basin or a periodic channel depending on the parameter \jp{jperio}.  
     545basin or a periodic channel depending on the parameter \np{jperio}.  
    556546\item[\np{nn\_bathy} = -1] a domain with a bump of topography one third of the  
    557547domain width at the central latitude. This is meant for the "EEL-R5" configuration,  
     
    599589vertical scale factors. The user must provide the analytical expression of both  
    600590$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}  
     591through statement functions, using parameters provided in the \ngn{namcfg} namelist.  
     592 
     593It is possible to define a simple regular vertical grid by giving zero stretching (\np{ppacr=0}).  
     594In that case, the parameters \jp{jpk} (number of $w$-levels) and \np{pphmax}  
    605595(total ocean depth in meters) fully define the grid.  
    606596 
     
    639629scale factors as a function of the model levels are shown in Fig.~\ref{Fig_zgr} and  
    640630given 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.  
    642632 
    643633Rather than entering parameters $h_{sur}$, $h_{0}$, and $h_{1}$ directly, it is  
    644634possible 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,  
    646636and specifies instead the four following parameters: 
    647637\begin{itemize} 
    648 \item    \pp{ppacr}=$h_{cr} $: stretching factor (nondimensional). The larger  
    649 \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 maximum  
     638\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  
    651641stretching 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). 
    654644\end{itemize} 
    655645As 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. 
     646parameters are: \jp{jpk}=46, \np{ppacr}=9, \np{ppkth}=23.563, \np{ppdzmin}=6m,  
     647\np{pphmax}=5750m. 
    658648 
    659649%>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     
    720710is allowed to have either a smaller or larger thickness than $e_{3t}(jpk)$: the  
    721711maximum 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  
     712specifying values in \ngn{namdom} namelist, as the maximum depth \np{pphmax}  
     713in partial steps: for example, with  
     714\np{pphmax}$=5750~m$ for the DRAKKAR 45 layer grid, the maximum ocean depth  
    724715allowed is actually $6000~m$ (the default thickness $e_{3t}(jpk-1)$ being $250~m$).  
    725716Two variables in the namdom namelist are used to define the partial step  
     
    740731\namdisplay{namzgr_sco}  
    741732%-------------------------------------------------------------------------------------------------------------- 
     733Options are defined in \ngn{namzgr\_sco}. 
    742734In $s$-coordinate (\np{ln\_sco}~=~true), the depth and thickness of the model  
    743735levels are defined from the product of a depth field and either a stretching  
     
    905897%------------------------------------------------------------------------------------------ 
    906898 
     899Options are defined in \ngn{namtsd}. 
    907900By default, the ocean start from rest (the velocity field is set to zero) and the initialization of  
    908901temperature 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  
    11% ================================================================ 
    2 % Chapter Ñ Ocean Dynamics (DYN) 
     2% Chapter Ocean Dynamics (DYN) 
    33% ================================================================ 
    44\chapter{Ocean Dynamics (DYN)} 
     
    167167The vector invariant form of the momentum equations is the one most  
    168168often 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 
     170through the \ngn{namdyn\_adv} namelist variables 
    170171Coriolis and momentum advection terms are evaluated using a leapfrog  
    171172scheme, $i.e.$ the velocity appearing in these expressions is centred in  
     
    184185%------------------------------------------------------------------------------------------------------------- 
    185186 
     187Options are defined through the \ngn{namdyn\_vor} namelist variables. 
    186188Four discretisations of the vorticity term (\textit{ln\_dynvor\_xxx}=true) are available:  
    187189conserving potential enstrophy of horizontally non-divergent flow (ENS scheme) ;  
     
    382384%------------------------------------------------------------------------------------------------------------- 
    383385 
     386Options are defined through the \ngn{namdyn\_adv} namelist variables. 
    384387In the flux form (as in the vector invariant form), the Coriolis and momentum  
    385388advection terms are evaluated using a leapfrog scheme, $i.e.$ the velocity  
     
    526529%------------------------------------------------------------------------------------------------------------- 
    527530 
     531Options are defined through the \ngn{namdyn\_hpg} namelist variables. 
    528532The key distinction between the different algorithms used for the hydrostatic  
    529533pressure gradient is the vertical coordinate used, since HPG is a \emph{horizontal}  
     
    712716 
    713717%%% 
     718Options are defined through the \ngn{namdyn\_spg} namelist variables. 
    714719The 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}. 
    715720 
     
    931936%------------------------------------------------------------------------------------------------------------- 
    932937 
     938Options are defined through the \ngn{namdyn\_ldf} namelist variables. 
    933939The options available for lateral diffusion are to use either laplacian  
    934940(rotated or not) or biharmonic operators. The coefficients may be constant  
     
    10601066%------------------------------------------------------------------------------------------------------------- 
    10611067 
     1068Options are defined through the \ngn{namzdf} namelist variables. 
    10621069The large vertical diffusion coefficient found in the surface mixed layer together  
    10631070with high vertical resolution implies that in the case of explicit time stepping there  
     
    11301137%------------------------------------------------------------------------------------------------------------- 
    11311138 
     1139Options are defined through the \ngn{namdom} namelist variables. 
    11321140The general framework for dynamics time stepping is a leap-frog scheme,  
    11331141$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  
    11% ================================================================ 
    2 % Chapter Ñ Lateral Boundary Condition (LBC)  
     2% Chapter Lateral Boundary Condition (LBC)  
    33% ================================================================ 
    44\chapter{Lateral Boundary Condition (LBC) } 
     
    2525%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}. 
    2626 
     27Options are defined through the \ngn{namlbc} namelist variables. 
    2728The discrete representation of a domain with complex boundaries (coastlines and  
    2829bottom topography) leads to arrays that include large portions where a computation  
     
    148149% Boundary Condition around the Model Domain 
    149150% ================================================================ 
    150 \section{Model Domain Boundary Condition (\jp{jperio})} 
     151\section{Model Domain Boundary Condition (\np{jperio})} 
    151152\label{LBC_jperio} 
    152153 
     
    156157 
    157158% ------------------------------------------------------------------------------------------------------------- 
    158 %        Closed, cyclic, south symmetric (\jp{jperio} = 0, 1 or 2)  
     159%        Closed, cyclic, south symmetric (\np{jperio} = 0, 1 or 2)  
    159160% ------------------------------------------------------------------------------------------------------------- 
    160 \subsection{Closed, cyclic, south symmetric (\jp{jperio} = 0, 1 or 2)} 
     161\subsection{Closed, cyclic, south symmetric (\np{jperio} = 0, 1 or 2)} 
    161162\label{LBC_jperio012} 
    162163 
    163164The 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 boundary  
     165by setting \np{jperio} to 0, 1 or 2 in namelist \ngn{namcfg}. Each time such a boundary  
    165166condition is needed, it is set by a call to routine \mdl{lbclnk}. The computation of  
    166167momentum and tracer trends proceeds from $i=2$ to $i=jpi-1$ and from $j=2$ to  
     
    295296 domain and the overlapping rows. The number of rows to exchange (known as  
    296297 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 between  
     298 dimensions are named \np{jpiglo}, \np{jpjglo} and \jp{jpk}. The relationship between  
    298299 the whole domain and a sub-domain is: 
    299300\begin{eqnarray}  
     
    419420\end{itemize} 
    420421 
     422Options are defined through the \ngn{namobc} namelist variables. 
    421423The package resides in the OBC directory. It is described here in four parts: the  
    422424boundary geometry (parameters to be set in \mdl{obc\_par}), the forcing data at  
     
    455457Logical flag  &                 &                            &                     \\ 
    456458\hline 
    457 West          & \jp{jpiwob} $>= 2$         &  \jp{jpjwd}$>= 2$          &  \jp{jpjwf}<= \jp{jpjglo}-1 \\ 
     459West          & \jp{jpiwob} $>= 2$         &  \jp{jpjwd}$>= 2$          &  \jp{jpjwf}<= \np{jpjglo}-1 \\ 
    458460lp\_obc\_west & $i$-index of a $u$ point   & $j$ of a $T$ point   &$j$ of a $T$ point \\ 
    459461\hline 
    460 East            & \jp{jpieob}$<=$\jp{jpiglo}-2&\jp{jpjed} $>= 2$         & \jp{jpjef}$<=$ \jp{jpjglo}-1 \\ 
     462East            & \jp{jpieob}$<=$\np{jpiglo}-2&\jp{jpjed} $>= 2$         & \jp{jpjef}$<=$ \np{jpjglo}-1 \\ 
    461463 lp\_obc\_east  & $i$-index of a $u$ point    & $j$ of a $T$ point & $j$ of a $T$ point \\ 
    462464\hline 
    463 South           & \jp{jpjsob} $>= 2$         & \jp{jpisd} $>= 2$          & \jp{jpisf}$<=$\jp{jpiglo}-1 \\ 
     465South           & \jp{jpjsob} $>= 2$         & \jp{jpisd} $>= 2$          & \jp{jpisf}$<=$\np{jpiglo}-1 \\ 
    464466lp\_obc\_south  & $j$-index of a $v$ point   & $i$ of a $T$ point   & $i$ of a $T$ point \\ 
    465467\hline 
    466 North           & \jp{jpjnob} $<=$ \jp{jpjglo}-2& \jp{jpind} $>= 2$        & \jp{jpinf}$<=$\jp{jpiglo}-1 \\ 
     468North           & \jp{jpjnob} $<=$ \np{jpjglo}-2& \jp{jpind} $>= 2$        & \jp{jpinf}$<=$\np{jpiglo}-1 \\ 
    467469lp\_obc\_north  & $j$-index of a $v$ point      & $i$  of a $T$ point & $i$ of a $T$ point \\ 
    468470\hline 
     
    754756%----------------------------------------------------------------------------------------------- 
    755757 
     758Options are defined through the \ngn{nambdy} \ngn{nambdy\_index}  
     759\ngn{nambdy\_dta} \ngn{nambdy\_dta2} namelist variables. 
    756760The BDY module is an alternative implementation of open boundary 
    757761conditions for regional configurations. It implements the Flow 
     
    10241028%----------------------------------------------------------------------------------------------- 
    10251029 
    1026 To be written.... 
    1027  
    1028  
    1029  
    1030  
     1030Options 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  
    11 
    22% ================================================================ 
    3 % Chapter Ñ Lateral Ocean Physics (LDF) 
     3% Chapter Lateral Ocean Physics (LDF) 
    44% ================================================================ 
    55\chapter{Lateral Ocean Physics (LDF)} 
     
    2121and for tracers only, eddy induced advection on tracers). These three aspects  
    2222of the lateral diffusion are set through namelist parameters and CPP keys  
    23 (see the \textit{nam\_traldf} and \textit{nam\_dynldf} below). Note 
     23(see the \textit{\ngn{nam\_traldf}} and \textit{\ngn{nam\_dynldf}} below). Note 
    2424that this chapter describes the default implementation of iso-neutral 
    2525tracer mixing, and Griffies's implementation, which is used if 
     
    104104 
    105105Other formulations can be introduced by the user for a given configuration.  
    106 For example, in the ORCA2 global ocean model (\key{orca\_r2}), the laplacian  
     106For example, in the ORCA2 global ocean model (see Configurations), the laplacian  
    107107viscosity operator uses \np{rn\_ahm0}~= 4.10$^4$ m$^2$/s poleward of 20$^{\circ}$  
    108108north and south and decreases linearly to \np{rn\_aht0}~= 2.10$^3$ m$^2$/s  
     
    110110can be found in routine \rou{ldf\_dyn\_c2d\_orca} defined in \mdl{ldfdyn\_c2d}.  
    111111Similar 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}). 
     112sub-domain options of ORCA2 and ORCA05 (see \&namcfg namelist). 
    114113 
    115114\subsubsection{Space Varying Mixing Coefficients (\key{traldf\_c3d} and \key{dynldf\_c3d})} 
     
    123122There is no default specification of space and time varying mixing coefficient.  
    124123The 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 tracer  
     124configurations. It provides only a tracer  
    126125mixing coefficient for eddy induced velocity (ORCA2) or both iso-neutral and  
    127126eddy 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  
    11% ================================================================ 
    2 % Chapter Ñ Miscellaneous Topics 
     2% Chapter Miscellaneous Topics 
    33% ================================================================ 
    44\chapter{Miscellaneous Topics} 
     
    3333Note that such modifications are so specific to a given configuration that no attempt  
    3434has 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). 
     35they can be set up is given in the ORCA 2\deg and 0.5\deg configurations. For example,  
     36for details of implementation in ORCA2, search: 
     37\vspace{-10pt}   
     38\begin{alltt}   
     39\tiny     
     40\begin{verbatim} 
     41IF( cp_cfg == "orca" .AND. jp_cfg == 2 ) 
     42\end{verbatim}   
     43\end{alltt} 
    3744 
    3845% ------------------------------------------------------------------------------------------------------------- 
     
    8390 
    8491\colorbox{yellow}{Add a short description of CLA staff here or in lateral boundary condition chapter?} 
     92Options are defined through the  \ngn{namcla} namelist variables. 
    8593 
    8694%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.   
     
    98106% Sub-Domain Functionality (\textit{nizoom, njzoom}, namelist parameters) 
    99107% ================================================================ 
    100 \section{Sub-Domain Functionality (\jp{jpizoom}, \jp{jpjzoom})} 
     108\section{Sub-Domain Functionality (\np{jpizoom}, \np{jpjzoom})} 
    101109\label{MISC_zoom} 
    102110 
     
    119127In order to easily define a sub-domain over which the computation can be  
    120128performed, 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 through  
    123 \jp{jpiglo}, \jp{jpjglo} and \jp{jpk} (\mdl{par\_oce} module). When running the  
    124 model over the whole domain, the user sets \jp{jpiglo}=\jp{jpidta} \jp{jpjglo}=\jp{jpjdta}  
     129forcing, 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  
     132model over the whole domain, the user sets \np{jpiglo}=\np{jpidta} \np{jpjglo}=\np{jpjdta}  
    125133and \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} in  
    128 the \mdl{par\_oce} module (Fig.~\ref{Fig_LBC_zoom}).  
     134has to provide the size of the sub-domain, (\np{jpiglo}, \np{jpjglo}, \np{jpkglo}),  
     135and the indices of the south western corner as \np{jpizoom} and \np{jpjzoom} in  
     136the  \ngn{namcfg} namelist (Fig.~\ref{Fig_LBC_zoom}).  
    129137 
    130138Note 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 the  
     139actually used to perform the computation. It is set by default to \jp{jpi}=\np{jpjglo}  
     140and \jp{jpj}=\np{jpjglo}, except for massively parallel computing where the  
    133141computational domain is laid out on local processor memories following a 2D  
    134142horizontal splitting. % (see {\S}IV.2-c) ref to the section to be updated 
     
    162170trajectory to reach it. 
    163171 
     172Options are defined through the  \ngn{namdom} namelist variables. 
    164173The acceleration of convergence option is used when \np{nn\_acc}=1. In that case,  
    165174$\rdt=rn\_rdt$ is the time step of dynamics while $\widetilde{\rdt}=rdttra$ is the  
     
    284293 
    285294 \gmcomment{why not make these bullets into subsections?} 
    286  
     295Options are defined through the  \ngn{namctl} namelist variables. 
    287296 
    288297$\bullet$ Vector optimisation: 
     
    343352a Successive-Over-Relaxation scheme (SOR) and a preconditioned conjugate gradient  
    344353scheme(PCG) \citep{Madec_al_OM88, Madec_PhD90}. The solver is selected trough the 
    345 the value of \np{nn\_solv} (namelist parameter).  
     354the value of \np{nn\_solv}   \ngn{namsol} namelist variable.  
    346355 
    347356The 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  
    11% ================================================================ 
    2 % Chapter 1 Ñ Model Basics 
     2% Chapter 1 Model Basics 
    33% ================================================================ 
    44% ================================================================ 
     
    5353Because $z^\star$ has a time independent range, all grid cells have static increments  
    5454ds, and the sum of the ver tical increments yields the time independent ocean  
    55 depth %·k ds = H.  
     55depth %k ds = H.  
    5656The $z^\star$ coordinate is therefore invisible to undulations of the  
    5757free surface, since it moves along with the free surface. This proper ty means that  
     
    7878\namdisplay{nam_dynspg}  
    7979%------------------------------------------------------------------------------------------------------------ 
     80Options are defined through the  \ngn{nam\_dynspg} namelist variables. 
    8081The 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}. 
    8182 
     
    114115\namdisplay{namdom}  
    115116%-------------------------------------------------------------------------------------------------------------- 
    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).  
     117The 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} 
     118in the  \ngn{namdom} namelist.  
     119(Figure III.3).  
    117120 
    118121%>   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   >   > 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Chap_OBS.tex

    r3294 r3989  
    2020can be used for validation or verification of model or  any other data assimilation system. 
    2121 
    22 The OBS code is called from \np{opa.F90} for model initialisation and to calculate the model 
     22The OBS code is called from \mdl{nemogcm.F90} for model initialisation and to calculate the model 
    2323equivalent 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 be 
     24each timestep from \mdl{step.F90}. To build with the OBS code active \key{diaobs} must be 
    2525set. 
    2626 
     
    6666 
    6767\item Add the following to the NEMO namelist to run the observation 
    68 operator on this data. Set the \np{enactfiles} namelist parameter to the 
     68operator on this data. Set the \np{enactfiles} namelist variable to the 
    6969observation  file name: 
    7070\end{enumerate} 
     
    7474%------------------------------------------------------------------------------------------------------------- 
    7575 
     76Options are defined through the  \ngn{namobs} namelist variables. 
    7677The options \np{ln\_t3d} and \np{ln\_s3d} switch on the temperature and salinity 
    7778profile observation operator code. The \np{ln\_ena} switch turns on the reading 
     
    9495\label{OBS_details} 
    9596 
    96 Here we show a more complete example namelist and also show the NetCDF headers 
     97Here we show a more complete example namelist  \ngn{namobs} and also show the NetCDF headers 
    9798of the observation 
    9899files that may be used with the observation operator 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Chap_SBC.tex

    r3795 r3989  
    11% ================================================================ 
    2 % Chapter Ñ Surface Boundary Condition (SBC, ICB)  
     2% Chapter Surface Boundary Condition (SBC, ICB)  
    33% ================================================================ 
    44\chapter{Surface Boundary Condition (SBC, ICB) } 
     
    2525 
    2626Five 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),  
     27are controlled by namelist \ngn{namsbc} variables: an analytical formulation (\np{ln\_ana}~=~true),  
    2828a flux formulation (\np{ln\_flx}~=~true), a bulk formulae formulation (CORE  
    2929(\np{ln\_core}~=~true), CLIO (\np{ln\_clio}~=~true) or MFS 
     
    442442%-------------------------------------------------------------------------------------------------------------- 
    443443 
    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: 
     444In 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.   
     445Options are defined through the  \ngn{namsbc\_sas} namelist variables. 
     446For example: 
    446447 
    447448\begin{enumerate} 
     
    507508In this case, all the six fluxes needed by the ocean are assumed to  
    508509be uniform in space. They take constant values given in the namelist  
    509 namsbc{\_}ana by 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},  
    510511\np{rn\_qsr0}, and \np{rn\_emp0} ($\textit{emp}=\textit{emp}_S$). The runoff is set to zero.  
    511512In addition, the wind is allowed to reach its nominal value within a given number  
     
    530531In the flux formulation (\np{ln\_flx}=true), the surface boundary  
    531532condition fields are directly read from input files. The user has to define  
    532 in the namelist namsbc{\_}flx the name of the file, the name of the variable  
     533in the namelist \ngn{namsbc{\_}flx} the name of the file, the name of the variable  
    533534read in the file, the time frequency at which it is given (in hours), and a logical  
    534535setting whether a time interpolation to the model time step is required  
     
    580581This is the so-called DRAKKAR Forcing Set (DFS) \citep{Brodeau_al_OM09}.  
    581582 
     583Options are defined through the  \ngn{namsbc\_core} namelist variables. 
    582584The required 8 input fields are: 
    583585 
     
    621623compute the radiative fluxes from a climatological cloud cover.  
    622624 
     625Options are defined through the  \ngn{namsbc\_clio} namelist variables. 
    623626The required 7 input fields are: 
    624627 
     
    673676Details on the bulk formulae used can be found in \citet{Maggiore_al_PCE98} and \citet{Castellari_al_JMS1998}. 
    674677 
     678Options are defined through the  \ngn{namsbc\_mfs} namelist variables. 
    675679The required 7 input fields must be provided on the model Grid-T and  are: 
    676680\begin{itemize} 
     
    711715When PISCES biogeochemical model (\key{top} and \key{pisces}) is also used in the coupled system,  
    712716the 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 
     717CO$_2$ fluxes will be exchanged between the atmosphere and the ice-ocean system (and need to be activated in \ngn{namsbc{\_}cpl} ). 
     718 
     719The namelist above allows control of various aspects of the coupling fields (particularly for 
    717720vectors) and now allows for any coupling fields to have multiple sea ice categories (as required by LIM3 
    718721and CICE).  When indicating a multi-category coupling field in namsbc{\_}cpl the number of categories will be 
     
    736739 
    737740The 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 ). 
    739742The input atmospheric forcing defined via \np{sn\_apr} structure (\textit{namsbc\_apr} namelist)  
    740743can be interpolated in time to the model time step, and even in space when the  
     
    774777%------------------------------------------------------------------------------------------------------------- 
    775778 
    776 Concerning the tidal potential, some parameters are available in namelist: 
    777  
    778 - \texttt{ln\_tide\_pot} activate the tidal potential forcing 
    779  
    780 - \texttt{nb\_harmo} is the number of constituent used 
    781  
    782 - \texttt{clname} is the name of constituent 
     779Concerning 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 
    783786 
    784787 
     
    858861depth (in metres) which the river should be added to. 
    859862 
    860 Namelist options, \np{ln\_rnf\_depth}, \np{ln\_rnf\_sal} and \np{ln\_rnf\_temp} control whether  
     863Namelist variables in \ngn{namsbc\_rnf}, \np{ln\_rnf\_depth}, \np{ln\_rnf\_sal} and \np{ln\_rnf\_temp} control whether  
    861864the river attributes (depth, salinity and temperature) are read in and used.  If these are set  
    862865as false the river is added to the surface box only, assumed to be fresh (0~psu), and/or  
     
    943946Their physical behaviour is controlled by equations as described in  \citet{Martin_Adcroft_OM10} ). 
    944947(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 by 
     948Icebergs are initially spawned into one of ten classes which have specific mass and thickness as described in the \ngn{namberg} namelist:  
    946949\np{rn\_initial\_mass} and \np{rn\_initial\_thickness}. 
    947950Each class has an associated scaling (\np{rn\_mass\_scaling}), which is an integer representing how many icebergs  
     
    10311034the diurnal cycle of SWF is a scaling of the top of the atmosphere diurnal cycle  
    10321035of 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 using  
     1036in \NEMO by setting \np{ln\_dm2dc}~=~true (a \textit{\ngn{namsbc}} namelist variable) when using  
    10341037CORE bulk formulea (\np{ln\_blk\_core}~=~true) or the flux formulation (\np{ln\_flx}~=~true).  
    10351038The reconstruction is performed in the \mdl{sbcdcy} module. The detail of the algoritm used  
     
    10881091%------------------------------------------------------------------------------------------------------------- 
    10891092 
    1090 In forced mode using a flux formulation (\np{ln\_flx}~=~true), a  
     1093IOptions are defined through the  \ngn{namsbc\_ssr} namelist variables. 
     1094n forced mode using a flux formulation (\np{ln\_flx}~=~true), a  
    10911095feedback term \emph{must} be added to the surface heat flux $Q_{ns}^o$: 
    10921096\begin{equation} \label{Eq_sbc_dmp_q} 
     
    12121216 in $namsbc$ namelist must be defined ${.true.}$.  
    12131217The \mdl{sbcwave} module containing the routine \np{sbc\_wave} reads the 
    1214 namelist ${namsbc\_wave}$ (for external data names, locations, frequency, interpolation and all  
     1218namelist \ngn{namsbc\_wave} (for external data names, locations, frequency, interpolation and all  
    12151219the miscellanous options allowed by Input Data generic Interface see \S\ref{SBC_input})  
    12161220and a 2D field of neutral drag coefficient. Then using the routine  
     
    12221226% Griffies doc: 
    12231227% 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 years 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.  
    12251229%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.  
    12261230 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Chap_STP.tex

    r3294 r3989  
    11 
    22% ================================================================ 
    3 % Chapter 2 Ñ Time Domain (step.F90) 
     3% Chapter 2 Time Domain (step.F90) 
    44% ================================================================ 
    55\chapter{Time Domain (STP) } 
     
    335335When restarting, if the the time step has been changed, a restart using an Euler time  
    336336stepping scheme is imposed.  
     337Options are defined through the  \ngn{namrun} namelist variables. 
    337338%%% 
    338339\gmcomment{ 
     
    358359%-------------------------------------------------------------------------------------------------------------- 
    359360 
    360  
     361Options are defined through the  \ngn{namdom} namelist variables. 
    361362 \colorbox{yellow}{add here a few word on nit000 and nitend} 
    362363 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Chap_TRA.tex

    r3764 r3989  
    11% ================================================================ 
    2 % Chapter 1 Ñ Ocean Tracers (TRA) 
     2% Chapter 1 Ocean Tracers (TRA) 
    33% ================================================================ 
    44\chapter{Ocean Tracers (TRA)} 
     
    137137\textit{effective} velocity ($i.e.$ the sum of the eulerian and eiv velocities) which is used. 
    138138 
    139 The choice of an advection scheme is made in the \textit{nam\_traadv} namelist, by  
     139The choice of an advection scheme is made in the \textit{\ngn{nam\_traadv}} namelist, by  
    140140setting to \textit{true} one and only one of the logicals \textit{ln\_traadv\_xxx}. The  
    141141corresponding code can be found in the \textit{traadv\_xxx.F90} module, where  
     
    441441%------------------------------------------------------------------------------------------------------------- 
    442442  
     443Options are defined through the  \ngn{namtra\_ldf} namelist variables. 
    443444The options available for lateral diffusion are a laplacian (rotated or not)  
    444445or a biharmonic operator, the latter being more scale-selective (more  
     
    602603%-------------------------------------------------------------------------------------------------------------- 
    603604 
     605Options are defined through the  \ngn{namzdf} namelist variables. 
    604606The formulation of the vertical subgrid scale tracer physics is the same  
    605607for all the vertical coordinates, and is based on a laplacian operator.  
     
    757759%-------------------------------------------------------------------------------------------------------------- 
    758760 
     761Options are defined through the  \ngn{namtra\_qsr} namelist variables. 
    759762When the penetrative solar radiation option is used (\np{ln\_flxqsr}=true),  
    760763the solar radiation penetrates the top few tens of meters of the ocean. If it is not used  
     
    879882Bottom Water) by a few Sverdrups  \citep{Emile-Geay_Madec_OS09}.  
    880883 
     884Options are defined through the  \ngn{namtra\_bbc} namelist variables. 
    881885The presence of geothermal heating is controlled by setting the namelist  
    882886parameter  \np{ln\_trabbc} to true. Then, when \np{nn\_geoflx} is set to 1,  
     
    897901%-------------------------------------------------------------------------------------------------------------- 
    898902 
     903Options are defined through the  \ngn{nambbl} namelist variables. 
    899904In a $z$-coordinate configuration, the bottom topography is represented by a  
    900905series of discrete steps. This is not adequate to represent gravity driven  
     
    10661071where $\gamma$ is the inverse of a time scale, and $T_o$ and $S_o$  
    10671072are given temperature and salinity fields (usually a climatology).  
     1073Options are defined through the  \ngn{namtra\_dmp} namelist variables. 
    10681074The restoring term is added when the namelist parameter \np{ln\_tradmp} is set to true.  
    10691075It also requires that both \np{ln\_tsd\_init} and \np{ln\_tsd\_tradmp} are set to true 
     
    11281134%-------------------------------------------------------------------------------------------------------------- 
    11291135 
     1136Options are defined through the  \ngn{namdom} namelist variables. 
    11301137The general framework for tracer time stepping is a modified leap-frog scheme  
    11311138\citep{Leclair_Madec_OM09}, $i.e.$ a three level centred time scheme associated  
     
    12051212\citep{Gill1982}. 
    12061213 
     1214Options are defined through the  \ngn{nameos} namelist variables. 
    12071215The default option (namelist parameter \np{nn\_eos}=0) is the \citet{JackMcD1995}  
    12081216equation of state. Its use is highly recommended. However, for process studies,  
     
    12221230coefficients, and $\rho_o$, the reference volumic mass, $rau0$.  
    12231231($\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 function  
     1232\np{rn\_beta} namelist variables). Note that when $d_a$ is a function  
    12251233of $T$ only (\np{nn\_eos}=1), the salinity is a passive tracer and can be  
    12261234used as such. 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Chap_ZDF.tex

    r3764 r3989  
    5353%-------------------------------------------------------------------------------------------------------------- 
    5454 
     55Options are defined through the  \ngn{namzdf} namelist variables. 
    5556When \key{zdfcst} is defined, the momentum and tracer vertical eddy coefficients  
    5657are set to constant values over the whole ocean. This is the crudest way to define  
     
    7980 
    8081When \key{zdfric} is defined, a local Richardson number dependent formulation  
    81 for the vertical momentum and tracer eddy coefficients is set. The vertical mixing  
     82for the vertical momentum and tracer eddy coefficients is set through the  \ngn{namzdf\_ric}  
     83namelist variables.The vertical mixing  
    8284coefficients are diagnosed from the large scale variables computed by the model.  
    8385\textit{In situ} measurements have been used to link vertical turbulent activity to  
     
    176178            \end{cases} 
    177179\end{align*} 
    178 The choice of $P_{rt}$ is controlled by the \np{nn\_pdl} namelist parameter. 
     180Options are defined through the  \ngn{namzdfy\_tke} namelist variables. 
     181The choice of $P_{rt}$ is controlled by the \np{nn\_pdl} namelist variable. 
    179182 
    180183At the sea surface, the value of $\bar{e}$ is prescribed from the wind  
     
    539542\caption{   \label{Tab_GLS}  
    540543Set of predefined GLS parameters, or equivalently predefined turbulence models available  
    541 with \key{zdfgls} and controlled by the \np{nn\_clos} namelist parameter.} 
     544with \key{zdfgls} and controlled by the \np{nn\_clos} namelist variable in \ngn{namzdf\_gls} .} 
    542545\end{center}   \end{table} 
    543546%-------------------------------------------------------------------------------------------------------------- 
     
    581584 
    582585The KKP scheme has been implemented by J. Chanut ... 
     586Options are defined through the  \ngn{namzdf\_kpp} namelist variables. 
    583587 
    584588\colorbox{yellow}{Add a description of KPP here.} 
     
    631635%>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    632636 
     637Options are defined through the  \ngn{namzdf} namelist variables. 
    633638The non-penetrative convective adjustment is used when \np{ln\_zdfnpc}=true.  
    634639It is applied at each \np{nn\_npc} time step and mixes downwards instantaneously  
     
    691696%-------------------------------------------------------------------------------------------------------------- 
    692697 
     698Options are defined through the  \ngn{namzdf} namelist variables. 
    693699The enhanced vertical diffusion parameterisation is used when \np{ln\_zdfevd}=true.  
    694700In this case, the vertical eddy mixing coefficients are assigned very large values  
     
    749755%-------------------------------------------------------------------------------------------------------------- 
    750756 
     757Options are defined through the  \ngn{namzdf\_ddm} namelist variables. 
    751758Double diffusion occurs when relatively warm, salty water overlies cooler, fresher  
    752759water, or vice versa. The former condition leads to salt fingering and the latter  
     
    830837%-------------------------------------------------------------------------------------------------------------- 
    831838 
     839Options are defined through the  \ngn{nambfr} namelist variables. 
    832840Both the surface momentum flux (wind stress) and the bottom momentum  
    833841flux (bottom friction) enter the equations as a condition on the vertical  
     
    11361144\label{ZDF_tmx_bottom} 
    11371145 
     1146Options are defined through the  \ngn{namzdf\_tmx} namelist variables. 
    11381147The parameterization of tidal mixing follows the general formulation for  
    11391148the vertical eddy diffusivity proposed by \citet{St_Laurent_al_GRL02} and  
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Chapters/Introduction.tex

    r3625 r3989  
    8181the coefficients with \citet{Blanke1993}, \citet{Large_al_RG94}, \citet{Pacanowski_Philander_JPO81},  
    8282or \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 
     89Some 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. 
     90Usual 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 
     105The 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} 
     110A template can be found in \textit{NEMO/OPA\_SRC/module.example} 
     111The 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 
    83114 
    84115Model 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  
    11!----------------------------------------------------------------------- 
    2 !       nam_tide       tide parameters (#ifdef key_tide) 
     2&nam_tide      !   tide parameters (#ifdef key_tide) 
    33!----------------------------------------------------------------------- 
    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 
    126   clname(2)     =   'S2' 
    137   clname(3)     =   'N2' 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namasm

    r3764 r3989  
    22&nam_asminc   !   assimilation increments                               ('key_asminc') 
    33!----------------------------------------------------------------------- 
    4     ln_bkgwri = .false.    !  Logical switch for writing out background state  
     4    ln_bkgwri = .false.    !  Logical switch for writing out background state 
    55    ln_trainc = .false.    !  Logical switch for applying tracer increments 
    66    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 
    88    ln_asmdin = .false.    !  Logical switch for Direct Initialization (DI) 
    99    ln_asmiau = .false.    !  Logical switch for Incremental Analysis Updating (IAU) 
     
    1515    ln_salfix = .false.    !  Logical switch for ensuring that the sa > salfixmin 
    1616    salfixmin = -9999      !  Minimum salinity after applying the increments 
    17     ndivdmp  = 0          !  Number of iterations of divergence damping operator 
     17    nn_divdmp = 0          !  Number of iterations of divergence damping operator 
    1818/ 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/nambdy

    r3294 r3989  
    22&nambdy        !  unstructured open boundaries                          ("key_bdy") 
    33!----------------------------------------------------------------------- 
    4     nb_bdy = 2                               !  number of open boundary sets        
    5     ln_coords_file = .true.,.false.          !  =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, 0                    !  boundary conditions for barotropic fields 
    10     nn_dyn2d_dta  =  3, 0                    !  = 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, 0                    !  boundary conditions for baroclinic velocities 
    15     nn_dyn3d_dta  =  0, 0                    !  = 0, bdy data are equal to the initial state 
    16                                              !  = 1, bdy data are read in 'bdydata   .nc' files 
    17     nn_tra        =  1, 1                    !  boundary conditions for T and S 
    18     nn_tra_dta    =  1, 1                    !  = 0, bdy data are equal to the initial state 
    19                                              !  = 1, bdy data are read in 'bdydata   .nc' files 
    20     nn_rimwidth  = 10, 5                     !  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 
     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 
    2323/ 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/nambdy_tide

    r3294 r3989  
    11!----------------------------------------------------------------------- 
    2 &nambdy_tide     ! tidal forcing at open boundaries               
     2&nambdy_tide     ! tidal forcing at open boundaries 
    33!----------------------------------------------------------------------- 
    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. 
    367/ 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namberg

    r3609 r3989  
    22&namberg       !   iceberg parameters 
    33!----------------------------------------------------------------------- 
    4       ln_icebergs              = .true. 
     4      ln_icebergs              = .false. 
    55      ln_bergdia               = .true.               ! Calculate budgets 
    66      nn_verbose_level         = 1                    ! Turn on more verbose output if level > 0 
    77      nn_verbose_write         = 15                   ! Timesteps between verbose messages 
    88      nn_sample_rate           = 1                    ! Timesteps between sampling for trajectory storage 
     9                                                      ! Initial mass required for an iceberg of each class 
    910      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   
    1012      rn_distribution          = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02 
    1113                                                      ! Ratio between effective and real iceberg mass (non-dim) 
     14                                                      ! i.e. number of icebergs represented at a point          
    1215      rn_mass_scaling          = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1 
    13                                                       ! Total thickness of newly calved bergs (m) 
     16                                                      ! thickness of newly calved bergs (m) 
    1417      rn_initial_thickness     = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250. 
    1518      rn_rho_bergs             = 850.                 ! Density of icebergs 
     
    1821      rn_bits_erosion_fraction = 0.                   ! Fraction of erosion melt flux to divert to bergy bits 
    1922      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    
    2428 
    2529               ! filename ! freq (hours) ! variable ! time interp. ! clim  !'yearly' or ! weights  ! rotation ! 
    2630               !          ! (<0  months) !   name   !  (logical)   ! (T/F) ! 'monthly'  ! filename ! pairing  ! 
    2731      sn_icb =  'calving' ,     -1       , 'calvingmask',  .true.      , .true., 'yearly'   , ' '      , ' ' 
    28  
    29       cn_dir = './' 
     32    
     33      cn_dir = './'  
    3034/ 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/nambfr

    r3294 r3989  
    77   rn_bfri2    =    1.e-3  !  bottom drag coefficient (non linear case) 
    88   rn_bfeb2    =    2.5e-3 !  bottom turbulent kinetic energy background  (m2/s2) 
     9   rn_bfrz0    =    3.e-3  ! bottom roughness for loglayer bfr coeff  
    910   ln_bfr2d    = .false.   !  horizontal variation of the bottom friction coef (read a 2D mask file ) 
    1011   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) 
    1213/ 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namctl

    r3294 r3989  
    1212   nn_bench    =    0      !  Bench mode (1/0): CAUTION use zero except for bench 
    1313                           !     (no physical validity of the results) 
     14   nn_timing   =    0      !  timing by routine activated (=1) creates timing.output file, or not (=0) 
    1415/ 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namdct

    r3294 r3989  
    44    nn_dct      = 15       !  time step frequency for transports computing 
    55    nn_dctwri   = 15       !  time step frequency for transports writing 
    6     nn_secdebug =  0       !      0 : no section to debug 
     6    nn_secdebug = 112      !      0 : no section to debug 
    77                           !     -1 : debug all section 
    88                           !  0 < n : debug section number n 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namdia_harm

    r3294 r3989  
    22&nam_diaharm   !   Harmonic analysis of tidal constituents ('key_diaharm') 
    33!----------------------------------------------------------------------- 
    4     nit000_han=1           ! First time step used for harmonic analysis 
    5     nitend_han=105         ! Last time step used for harmonic analysis 
    6     nstep_han=1            ! Time step frequency for harmonic analysis 
    7     nb_ana=1               ! Number of harmonics to analyse 
    8     tname(1)='M2'          ! Name of tidal constituents 
     4    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' 
    99/ 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namdyn_hpg

    r3294 r3989  
    22&namdyn_hpg    !   Hydrostatic pressure gradient option 
    33!----------------------------------------------------------------------- 
    4    ln_hpg_zco  = .false.   !  z-coordinate - full steps                    
     4   ln_hpg_zco  = .false.   !  z-coordinate - full steps 
    55   ln_hpg_zps  = .true.    !  z-coordinate - partial steps (interpolation) 
    66   ln_hpg_sco  = .false.   !  s-coordinate (standard jacobian formulation) 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namdyn_ldf

    r3294 r3989  
    22&namdyn_ldf    !   lateral diffusion on momentum 
    33!----------------------------------------------------------------------- 
    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 
    99   ln_dynldf_hor    =  .true.   !  horizontal (geopotential)            (require "key_ldfslp" in s-coord.) 
    1010   ln_dynldf_iso    =  .false.  !  iso-neutral                          (require "key_ldfslp") 
     
    1212   rn_ahm_0_lap     = 40000.    !  horizontal laplacian eddy viscosity   [m2/s] 
    1313   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) 
    1520/ 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namdyn_nept

    r3294 r3989  
    11!----------------------------------------------------------------------- 
    2 &nam_dynnept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
     2&namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    33!----------------------------------------------------------------------- 
    44   ! Suggested lengthscale values are those of Eby & Holloway (1994) for a coarse model 
     
    1010   ! Specify whether to ramp down the Neptune velocity in shallow 
    1111   ! water, and if so the depth range controlling such ramping down 
    12    ln_neptramp       = .false.  ! ramp down Neptune velocity in shallow water 
     12   ln_neptramp       = .true.   ! ramp down Neptune velocity in shallow water 
    1313   rn_htrmin         =  100.0   ! min. depth of transition range 
    1414   rn_htrmax         =  200.0   ! max. depth of transition range 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namdyn_vor

    r3764 r3989  
    22&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    33!----------------------------------------------------------------------- 
    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 
    98/ 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namflo

    r3294 r3989  
    22&namflo       !   float parameters                                      ("key_float") 
    33!----------------------------------------------------------------------- 
    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_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) 
    1414/ 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namlbc

    r3294 r3989  
    44   rn_shlat    =    2.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat 
    55                           !  free slip  !   partial slip  !   no slip   ! strong slip 
     6   ln_vorlat   = .false.   !  consistency of vorticity boundary condition with analytical eqs. 
    67/ 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/nammpp

    r3294 r3989  
    55                           !  buffer blocking send or immediate non-blocking sends, resp. 
    66   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) 
    711/ 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namobc

    r3294 r3989  
    44   ln_obc_clim = .false.   !  climatological obc data files (T) or not (F) 
    55   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 
    77   nn_obcdta   =    1      !  = 0 the obc data are equal to the initial state 
    88                           !  = 1 the obc data are read in 'obc.dta' files 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namobs

    r2540 r3989  
    22&namobs       !  observation usage switch                               ('key_diaobs') 
    33!----------------------------------------------------------------------- 
    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 
    1010 
    11    ln_sladt   = .false.    ! Logical switch for AVISO SLA data               
     11   ln_sladt   = .false.    ! Logical switch for AVISO SLA data 
    1212 
    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 
    1515 
    16    ln_sst     = .false.    ! Logical switch for SST observations               
    17                            !     ln_reysst               Logical switch for Reynolds observations        
    18                            !     ln_ghrsst               Logical switch for GHRSST observations           
     16   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       
    1919 
    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 
    2727                           !     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 
    3535   profbfiles = 'profiles_01.nc' 
    36                            !     ln_profb_enatim         Enact feedback input time setting switch     
     36                           !     ln_profb_enatim         Enact feedback input time setting switch 
    3737                           !     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 
    4040   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 
    4343   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 
    6161   ln_grid_global = .true. 
    6262   ln_grid_search_lookup = .false. 
    63 /  
     63/ 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namptr

    r3294 r3989  
    44   ln_diaptr  = .false.    !  Poleward heat and salt transport (T) or not (F) 
    55   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 
    77                           !  (orca configuration only, need input basins mask file named "subbasins.nc" 
    88   ln_ptrcomp = .true.     !  Add decomposition : overturning 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namrun

    r3306 r3989  
    33!----------------------------------------------------------------------- 
    44   nn_no       =       0   !  job number (no more used...) 
    5    cn_exp      =  "ORCA2"  !  experience name  
     5   cn_exp      =  "ORCA2"  !  experience name 
    66   nn_it000    =       1   !  first time step 
    77   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) 
    109   nn_leapy    =       0   !  Leap year calendar (1) or not (0) 
    1110   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T) 
    1211   nn_rstctl   =       0   !  restart control => activated only if ln_rstart = T 
    13        !  = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 
    14        !  = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart 
    15        !  = 2 nn_date0 read in restart  ; nn_it000 : check consistancy between namelist and restart 
     12                           !    = 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 
    1615   cn_ocerst_in  = "restart"   !  suffix of ocean restart name (input) 
    1716   cn_ocerst_out = "restart"   !  suffix of ocean restart name (output) 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namsbc

    r3294 r3989  
    22&namsbc        !   Surface Boundary Condition (surface module) 
    33!----------------------------------------------------------------------- 
    4    nn_fsbc     = 5         !  frequency of surface boundary condition computation  
     4   nn_fsbc     = 5         !  frequency of surface boundary condition computation 
    55                           !     (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 ) 
    77   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) 
    1010   ln_blk_mfs  = .false.   !  MFS bulk formulation                      (T => fill namsbc_mfs ) 
    1111   ln_cpl      = .false.   !  Coupled formulation                       (T => fill namsbc_cpl ) 
     
    1414                           !  =1 use observed ice-cover      , 
    1515                           !  =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) 
    1619   ln_dm2dc    = .false.   !  daily mean to diurnal cycle on short wave 
    1720   ln_rnf      = .true.    !  runoffs                                   (T => fill namsbc_rnf) 
    1821   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 
    2124                           !     =2 annual global mean of e-p-r set to zero 
    2225                           !     =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) 
    2327   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) 
    2429/ 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namsbc_alb

    r3294 r3989  
    22&namsbc_alb    !   albedo parameters 
    33!----------------------------------------------------------------------- 
    4    rn_cloud    =    0.06   !  cloud correction to snow and ice albedo  
     4   rn_cloud    =    0.06   !  cloud correction to snow and ice albedo 
    55   rn_albice   =    0.53   !  albedo of melting ice in the arctic and antarctic 
    66   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 
    88   rn_alphdi   =    0.72   !  (Pyane, 1972) 
    99/ 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namsbc_clio

    r3294 r3989  
    11!----------------------------------------------------------------------- 
    2 &namsbc_clio   !   namsbc_clio  CLIO bulk formulea 
     2&namsbc_clio   !   namsbc_clio  CLIO bulk formulae 
    33!----------------------------------------------------------------------- 
    44!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namsbc_core

    r3294 r3989  
    11!----------------------------------------------------------------------- 
    2 &namsbc_core   !   namsbc_core  CORE bulk formulea 
     2&namsbc_core   !   namsbc_core  CORE bulk formulae 
    33!----------------------------------------------------------------------- 
    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'  , ''       , '' 
    1515 
    1616   cn_dir      = './'      !  root directory for the location of the bulk files 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namsbc_cpl

    r3294 r3989  
    55!                    !                       ! categories !  reference  !    orientation       ! grids  ! 
    66! 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'    ,     ''      ,         ''           ,   ''         
     7sn_snd_temp   =       'weighted oce and ice' ,    'no'    ,     ''      ,         ''           ,   '' 
     8sn_snd_alb    =       'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   '' 
     9sn_snd_thick  =       'none'                 ,    'no'   ,     ''      ,         ''           ,   '' 
     10sn_snd_crt    =       'none'                 ,    'no'    , 'spherical' , 'eastward-northward' ,  'T' 
     11sn_snd_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
    1212! 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'    ,     ''      ,         ''          ,   ''     
     13sn_rcv_w10m   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     14sn_rcv_taumod =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     15sn_rcv_tau    =       'oce only'             ,    'no'    , 'cartesian' , 'eastward-northward',  'U,V' 
     16sn_rcv_dqnsdt =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     17sn_rcv_qsr    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
     18sn_rcv_qns    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
     19sn_rcv_emp    =       'conservative'         ,    'no'    ,     ''      ,         ''          ,   '' 
     20sn_rcv_rnf    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     21sn_rcv_cal    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     22sn_rcv_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    2323/ 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namsbc_flx

    r3294 r3989  
    1111 
    1212   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  
    1010   sn_tair     =   'ecmwf'   ,        6          , 't2'      ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , '' 
    1111   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' , '' 
    1313 
    1414   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  
    22&namsbc_rnf    !   runoffs namelist surface boundary condition 
    33!----------------------------------------------------------------------- 
    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'  , ''       , '' 
    1111 
    1212   cn_dir       = './'      !  root directory for the location of the runoff files 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namsbc_ssr

    r3294 r3989  
    99   cn_dir      = './'      !  root directory for the location of the runoff files 
    1010   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) 
    1212                           !  or to SSS only (=1) or no damping term (=0) 
    1313   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] 
    1515   ln_sssr_bnd =   .true.  !  flag to bound erp term (associated with nn_sssr=2) 
    1616   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  
    44!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
    55!              !             !  (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'  ,''         , '' 
    710! 
    811   cn_dir_cdg  = './'  !  root directory for the location of drag coefficient files 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namsol

    r3294 r3989  
    11!----------------------------------------------------------------------- 
    2 &namsol        !   elliptic solver / island / free surface  
     2&namsol        !   elliptic solver / island / free surface 
    33!----------------------------------------------------------------------- 
    44   nn_solv     =      1    !  elliptic solver: =1 preconditioned conjugate gradient (pcg) 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namtra_adv

    r3294 r3989  
    11!----------------------------------------------------------------------- 
    2 &namtra_adv    !   advection scheme for tracer  
     2&namtra_adv    !   advection scheme for tracer 
    33!----------------------------------------------------------------------- 
    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 
    1011/ 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namtra_bbc

    r2540 r3989  
    33!----------------------------------------------------------------------- 
    44   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 
    66                           !     = 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) 
    88   rn_geoflx_cst = 86.4e-3 !  Constant value of geothermal heat flux [W/m2] 
    99/ 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namtra_ldf

    r3296 r3989  
    2222   rn_ahtb_0        =     0.    !  background eddy diffusivity for ldf_iso [m2/s] 
    2323   !                                           (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) 
    2428/ 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namtsd

    r3294 r3989  
    22&namtsd    !   data : Temperature  & Salinity 
    33!----------------------------------------------------------------------- 
    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'   , ''       , ' ' 
     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'   , ''       , ' ' 
    88   ! 
    99   cn_dir        = './'     !  root directory for the location of the runoff files 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namzdf_kpp

    r3294 r3989  
    22&namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally: 
    33!------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 
    4    ln_kpprimix = .true.    !  shear instability mixing  
     4   ln_kpprimix = .true.    !  shear instability mixing 
    55   rn_difmiw   =  1.0e-04  !  constant internal wave viscosity [m2/s] 
    66   rn_difsiw   =  0.1e-04  !  constant internal wave diffusivity [m2/s] 
    77   rn_riinfty  =  0.8      !  local Richardson Number limit for shear instability 
    88   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] 
    1111   nn_avb      =  0        !  horizontal averaged (=1) or not (=0) on avt and amv 
    1212   nn_ave      =  1        !  constant (=0) or profile (=1) background on avt 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namzdf_tke

    r3294 r3989  
    77   rn_emin     =   1.e-6   !  minimum value of tke [m2/s2] 
    88   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) 
    910   nn_mxl      =   2       !  mixing length: = 0 bounded by the distance to surface and bottom 
    1011                           !                 = 1 bounded by the local vertical scale factor 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/DOC/TexFiles/Namelist/namzdf_tmx

    r2540 r3989  
    55   rn_n2min    = 1.e-8     !  threshold of the Brunt-Vaisala frequency (s-1) 
    66   rn_tfe      = 0.333     !  tidal dissipation efficiency 
    7    rn_me       = 0.2       !  mixing efficiency  
     7   rn_me       = 0.2       !  mixing efficiency 
    88   ln_tmx_itf  = .true.    !  ITF specific parameterisation 
    99   rn_tfe_itf  = 1.        !  ITF tidal dissipation efficiency 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/CONFIG/SHARED/README_other_configurations_namelist_namcfg

    r3973 r3989  
    152152!----------------------------------------------------------------------- 
    153153   cp_cfg      =  "orca"               !  name of the configuration 
     154   cp_cfz      =  "antarctic"          !  name of the zoom of configuration 
    154155   jp_cfg      =      05               !  resolution of the configuration 
    155156   jpidta      =     722               !  1st lateral dimension ( >= jpi ) 
     
    188189!----------------------------------------------------------------------- 
    189190   cp_cfg      =  "orca"               !  name of the configuration 
     191   cp_cfz      =  "arctic"             !  name of the zoom of configuration 
    190192   jp_cfg      =      05               !  resolution of the configuration 
    191193   jpidta      =     722               !  1st lateral dimension ( >= jpi ) 
     
    216218   ppacr2      =  999999.0             ! 
    217219/ 
     220 
     221ORCA2 - 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 
     258ORCA2 - 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 
    218295ORCA025 - 75 vertical levels 
    219296======= 
     
    428505   ppacr2      =  999999.0             ! 
    429506/ 
     507 
     508C1D - 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  
    4848! 
    4949!----------------------------------------------------------------------- 
    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       
    6651!----------------------------------------------------------------------- 
    6752   cp_cfg      =  "default"            !  name of the configuration 
     53   cp_cfz      =         ''            !  name of the zoom of configuration 
    6854   jp_cfg      =       0               !  resolution of the configuration 
    6955   jpidta      =      10               !  1st lateral dimension ( >= jpi ) 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/DOM/dom_oce.F90

    r3973 r3989  
    7272   LOGICAL, PUBLIC ::   lzoom_s    =  .FALSE.   !: South zoom type flag 
    7373   LOGICAL, PUBLIC ::   lzoom_n    =  .FALSE.   !: North zoom type flag 
    74    LOGICAL, PUBLIC ::   lzoom_arct =  .FALSE.   !: ORCA    arctic zoom flag 
    75    LOGICAL, PUBLIC ::   lzoom_anta =  .FALSE.   !: ORCA antarctic zoom flag 
    7674 
    7775   !                                     !!! domain parameters linked to mpp 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/DOM/domcfg.F90

    r3973 r3989  
    168168         SELECT CASE ( jp_cfg ) 
    169169         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 
    174178            !                              
    175179         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 
    180188         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)' 
    185189         ! 
    186190      ENDIF 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

    r3901 r3989  
    564564      ! Configuration specific domain modifications 
    565565      ! (here, ORCA arctic configuration: suppress Med Sea) 
    566       IF( cp_cfg == "orca" .AND. lzoom_arct ) THEN 
     566      IF( cp_cfg == "orca" .AND. cp_cfz == "arctic" ) THEN 
    567567         SELECT CASE ( jp_cfg ) 
    568568         !                                        ! ======================= 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_c2d.h90

    r3294 r3989  
    160160      IF(lwp) WRITE(numout,*) '        orca ocean configuration' 
    161161 
    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 
    176237         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 
     2679101     FORMAT(1x,a15,2i8) 
     2689201     FORMAT(3x,13(i3,12x)) 
     2699202     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) )   & 
    259324                     / (zmsk * 100.) 
    260             ENDIF 
    261             ahm2(ji,jj) = zcoff * ahm0 + (1.-zcoff) * ahm2(ji,jj) 
    262          END DO 
    263       END DO 
    264 #endif 
     325               ENDIF 
     326               ahm2(ji,jj) = zcoff * ahm0 + (1.-zcoff) * ahm2(ji,jj) 
     327            END DO 
     328         END DO 
     329      ENDIF 
    265330       
    266331      ! Lateral boundary conditions on ( ahm1, ahm2 ) 
     
    323388      IF(lwp) WRITE(numout,*) '        orca_r1 configuration' 
    324389 
    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 
    341468         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 
     4989101     FORMAT(1x,a15,2i8) 
     4999201     FORMAT(3x,13(i3,12x)) 
     5009202     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 
    435566       
    436567      ! Lateral boundary conditions on ( ahm1, ahm2 ) 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/TRA/tradmp.F90

    r3901 r3989  
    303303 
    304304      !                                           ! ==================================================== 
    305       IF( lzoom_arct .AND. lzoom_anta ) THEN      !  ORCA configuration : arctic zoom or antarctic zoom 
     305      IF( cp_cfz == "arctic" . OR. cp_cfz == "antarctic" ) THEN   !  ORCA configuration : arctic or antarctic zoom 
    306306         !                                        ! ==================================================== 
    307307         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' 
    310310         IF(lwp) WRITE(numout,*) 
    311311         ! 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/nemogcm.F90

    r3973 r3989  
    228228         &             nn_isplt, nn_jsplt, nn_jctls, nn_jctle,   & 
    229229         &             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, & 
    231231         &             jpizoom, jpjzoom, jperio, jphgr_msh, & 
    232232         &             ppglam0, ppgphi0, ppe1_deg, ppe2_deg, ppe1_m, ppe2_m, & 
     
    482482         WRITE(numout,*) '   Namelist namcfg' 
    483483         WRITE(numout,*) '      configuration name              cp_cfg      = ', TRIM(cp_cfg) 
     484         WRITE(numout,*) '      configuration zoom name         cp_cfz      = ', TRIM(cp_cfz) 
    484485         WRITE(numout,*) '      configuration resolution        jp_cfg      = ', jp_cfg 
    485486         WRITE(numout,*) '      1st lateral dimension ( >= jpi ) jpidta     = ', jpidta 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/par_oce.F90

    r3973 r3989  
    2929   !!---------------------------------------------------------------------- 
    3030   CHARACTER(lc) ::   cp_cfg           !: name of the configuration 
     31   CHARACTER(lc) ::   cp_cfz           !: name of the zoom of configuration 
    3132   INTEGER       ::   jp_cfg           !: resolution of the configuration 
    3233 
Note: See TracChangeset for help on using the changeset viewer.