New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 2381 for branches/nemo_v3_3_beta/DOC/TexFiles/Chapters/Chap_CFG.tex – NEMO

Ignore:
Timestamp:
2010-11-13T12:09:27+01:00 (13 years ago)
Author:
gm
Message:

v3.3beta - DOC: add ORCA1 in CFG

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/nemo_v3_3_beta/DOC/TexFiles/Chapters/Chap_CFG.tex

    r2376 r2381  
    4141The 1D model is a very useful tool   
    4242\textit{(a)} to learn about the physics and numerical treatment of vertical mixing processes ;  
    43 \textit{(b)} to investigate suitable parameterisations of unresolved turbulence (wind steering,  
    44 langmuir circulation, skin layers, ...) ;  
     43\textit{(b)} to investigate suitable parameterisations of unresolved turbulence (surface wave 
     44breaking, Langmuir circulation, ...) ;  
    4545\textit{(c)} to compare the behaviour of different vertical mixing schemes  ;  
    4646\textit{(d)} to perform sensitivity studies on the vertical diffusion at a particular point of an ocean domain ;  
     
    5151with some extra routines. There is no need to define a new mesh, bathymetry,  
    5252initial state or forcing, since the 1D model will use those of the configuration it is a zoom of.  
    53 The chosen grid point is set in par\_oce.F90 module by setting the jpizoom and jpjzoom  
     53The chosen grid point is set in par\_oce.F90 module by setting the \jp{jpizoom} and \jp{jpjzoom}  
    5454parameters to the indices of the location of the chosen grid point. 
    5555 
     56The 1D model has some specifies. First, all the horizontal derivatives are assumed to be zero.  
     57Therefore a simplified \rou{step} routine is used (\rou{step\_c1d}) in which both lateral tendancy  
     58terms and lateral physics are not called, and the vertical velocity is zero (so far, no attempt at 
     59introducing a Ekman pumping velocity has been made). 
     60Second, the two components of the velocity are moved on a $T$-point.  
     61This requires a specific treatment of the Coriolis term (see \rou{dyncor\_c1d}) and of the  
     62dynamic time stepping (\rou{dynnxt\_c1d}). 
     63All the relevant modules can be found in the NEMOGCM/NEMO/C1D\_SRC directory of  
     64the \NEMO distribution. 
     65 
     66% to be added:  a test case on the yearlong Ocean Weather Station (OWS) Papa dataset of Martin (1985) 
    5667 
    5768% ================================================================ 
     
    6475the LIM sea-ice model (ORCA-LIM) and possibly with PISCES biogeochemical model  
    6576(ORCA-LIM-PISCES), using various resolutions. 
     77 
     78 
     79%>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     80\begin{figure}[!t]   \begin{center} 
     81\includegraphics[width=0.98\textwidth]{./TexFiles/Figures/Fig_ORCA_NH_mesh.pdf} 
     82\caption{  \label{Fig_MISC_ORCA_msh}      
     83ORCA mesh conception. The departure from an isotropic Mercator grid start poleward of 20\deg N. 
     84The two "north pole" are the foci of a series of embedded ellipses (blue curves)  
     85which are determined analytically and form the i-lines of the ORCA mesh (pseudo latitudes).  
     86Then, following \citet{Madec_Imbard_CD96}, the normal to the series of ellipses (red curves) is computed  
     87which provide the j-lines of the mesh (pseudo longitudes).  } 
     88\end{center}   \end{figure} 
     89%>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    6690 
    6791% ------------------------------------------------------------------------------------------------------------- 
     
    80104either the mesh lines or the scale factors, or even the scale factor derivatives over the whole  
    81105ocean domain, as the mesh is not a composite mesh.  
     106%>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     107\begin{figure}[!tbp]  \begin{center} 
     108\includegraphics[width=1.0\textwidth]{./TexFiles/Figures/Fig_ORCA_NH_msh05_e1_e2.pdf} 
     109\includegraphics[width=0.80\textwidth]{./TexFiles/Figures/Fig_ORCA_aniso.pdf} 
     110\caption {  \label{Fig_MISC_ORCA_e1e2} 
     111\textit{Top}: Horizontal scale factors ($e_1$, $e_2$) and  
     112\textit{Bottom}: ratio of anisotropy ($e_1 / e_2$) 
     113for ORCA 0.5\deg ~mesh. South of 20\deg N a Mercator grid is used ($e_1 = e_2$)  
     114so that the anisotropy ratio is 1. Poleward of 20\deg N, the two "north pole"  
     115introduce a weak anisotropy over the ocean areas ($< 1.2$) except in vicinity of Victoria Island  
     116(Canadian Arctic Archipelago). } 
     117\end{center}   \end{figure} 
     118%>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     119 
    82120 
    83121The method is applied to Mercator grid ($i.e.$ same zonal and meridional grid spacing) poleward  
     
    93131in the Canadian Archipelago.  
    94132 
    95  
    96 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    97 \begin{figure}[!t]   \begin{center} 
    98 \includegraphics[width=0.98\textwidth]{./TexFiles/Figures/Fig_ORCA_NH_mesh.pdf} 
    99 \caption{  \label{Fig_MISC_ORCA_msh}      
    100 ORCA mesh conception. The departure from an isotropic Mercator grid start poleward of 20\deg N. 
    101 The two "north pole" are the foci of a series of embedded ellipses (blue curves)  
    102 which are determined analytically and form the i-lines of the ORCA mesh (pseudo latitudes).  
    103 Then, following \citet{Madec_Imbard_CD96}, the normal to the series of ellipses (red curves) is computed  
    104 which provide the j-lines of the mesh (pseudo longitudes).  } 
    105 \end{center}   \end{figure} 
    106 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    107  
    108  
    109 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    110 \begin{figure}[!tbp]  \begin{center} 
    111 \includegraphics[width=1.0\textwidth]{./TexFiles/Figures/Fig_ORCA_NH_msh05_e1_e2.pdf} 
    112 \includegraphics[width=0.80\textwidth]{./TexFiles/Figures/Fig_ORCA_aniso.pdf} 
    113 \caption {  \label{Fig_MISC_ORCA_e1e2} 
    114 \textit{Top}: Horizontal scale factors ($e_1$, $e_2$) and  
    115 \textit{Bottom}: ratio of anisotropy ($e_1 / e_2$) 
    116 for ORCA 0.5\deg ~mesh. South of 20\deg N a Mercator grid is used ($e_1 = e_2$)  
    117 so that the anisotropy ratio is 1. Poleward of 20\deg N, the two "north pole"  
    118 introduce a weak anisotropy over the ocean areas ($< 1.2$) except in vicinity of Victoria Island  
    119 (Canadian Arctic Archipelago). } 
    120 \end{center}   \end{figure} 
    121 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    122  
    123  
    124  
    125133% ------------------------------------------------------------------------------------------------------------- 
    126134%       ORCA-LIM(-PISCES) configurations 
    127135% ------------------------------------------------------------------------------------------------------------- 
    128 \subsection{ORCA-LIM(-PISCES) configurations} 
    129 \label{CFG_orca_grid} 
    130  
    131  
    132 The NEMO system is provided with four built-in ORCA configurations which differ in the  
    133 horizontal resolution 
    134 \footnote{the value of the resolution is given by the resolution at the Equator expressed in degrees.}  
    135 used: 
    136 \begin{description} 
    137 \item[\key{orca\_r4}]  \jp{cp\_cfg}~=~orca ; \jp{jp\_cfg}~=~4 
    138 \item[\key{orca\_r2}]  \jp{cp\_cfg}~=~orca ; \jp{jp\_cfg}~=~2 
    139 %\item[\key{orca\_r1}]  \jp{cp\_cfg}~=~orca ; \jp{jp\_cfg}~=~1 
    140 \item[\key{orca\_r05}]  \jp{cp\_cfg}~=~orca ; \jp{jp\_cfg}~=~05 
    141 \item[\key{orca\_r025}]  \jp{cp\_cfg}~=~orca ; \jp{jp\_cfg}~=~025 
    142 \end{description} 
     136\subsection{ORCA pre-defined resolution} 
     137\label{CFG_orca_resolution} 
     138 
     139 
     140The NEMO system is provided with five built-in ORCA configurations which differ in the  
     141horizontal resolution. The value of the resolution is given by the resolution at the Equator  
     142expressed in degrees. Each of configuration is set through a CPP key with set the grid size  
     143and configuration name parameters  (Tab.~\ref{Tab_ORCA}). 
     144. 
     145 
     146%--------------------------------------------------TABLE-------------------------------------------------- 
     147\begin{table}[!b]     \begin{center} 
     148\begin{tabular}{p{4cm} c c c c} 
     149CPP key                        & \jp{jp\_cfg} &  \jp{jpiglo} & \jp{jpiglo} &       \\   
     150\hline  \hline 
     151\key{orca\_r4}        &        4         &         92     &      76      &       \\ 
     152\key{orca\_r2}       &        2         &       182     &    149      &        \\ 
     153\key{orca\_r1}       &        1         &       362     &     292     &        \\ 
     154\key{orca\_r05}     &        05       &       722     &     511     &        \\ 
     155\key{orca\_r025}   &        025     &      1442    &   1021     &        \\ 
     156%\key{orca\_r8}       &        8         &      2882    &   2042     &        \\ 
     157%\key{orca\_r12}     &      12         &      4322    &   3062      &       \\ 
     158\hline   \hline 
     159\end{tabular} 
     160\caption{ \label{Tab_ORCA}    
     161Set of predefined parameters for ORCA family configurations. 
     162In all cases, the name of the configuration is set to "orca" ($i.e.$ \jp{cp\_cfg}~=~orca). } 
     163\end{center} 
     164\end{table} 
     165%-------------------------------------------------------------------------------------------------------------- 
     166 
    143167 
    144168The ORCA\_R2 configuration has the following specificity : starting from a 2\deg~ORCA mesh,  
     
    151175and R1 strongly increases the mesh anisotropy there. 
    152176 
    153 The ORCA\_R05 configuration and higher ones does not incorporate any regional refinements.   
     177The ORCA\_R05 and higher global configurations do not incorporate any regional refinements.   
     178 
     179For ORCA\_R1 and R025, setting the configuration key to 75 allows to use 75 vertical levels,  
     180otherwise 46 are used. In the other ORCA configurations, 31 levels are used  
     181(see Tab.~\ref{Tab_orca_zgr} and Fig.~\ref{Fig_zgr}). 
    154182 
    155183Only the ORCA\_R2 is provided with all its input files in the \NEMO distribution.  
     
    159187 
    160188This version of ORCA\_R2 has 31 levels in the vertical, with the highest resolution (10m)  
    161 in the upper 150m. The bottom topography and the coastlines are derived  
    162 from the global atlas of Smith and Sandwell (1997). The default forcing employ the boundary  
    163 forcing from \citet{Large_Yeager_Rep04} (see \S\ref{SBC_blk_core}),  
     189in the upper 150m (see Tab.~\ref{Tab_orca_zgr} and Fig.~\ref{Fig_zgr}).  
     190The bottom topography and the coastlines are derived from the global atlas of Smith and Sandwell (1997).  
     191The default forcing employ the boundary forcing from \citet{Large_Yeager_Rep04} (see \S\ref{SBC_blk_core}),  
    164192which was developed for the purpose of running global coupled ocean-ice simulations  
    165193without an interactive atmosphere. This \citet{Large_Yeager_Rep04} dataset is available  
     
    169197 
    170198The vertical resolution can be increased by a factor of 10 by defining the \key{orca\_lev10} CPP key.  
    171 It can also be run with an AGRIF zoom over the Agulhas current area ( \key{agrif}  defined).  
    172 Also available are to keys, \key{arctic} and \key{antarctic}, which allows to run a regional Arctic  
    173 or peri-Antarctic configuration extracted from an ORCA configuration. (This does not work with ORCA\_R4 and R1). 
    174  
    175  
    176 %--------------------------------------------------TABLE-------------------------------------------------- 
    177 \begin{table}[htbp]     \begin{center} 
    178 \begin{tabular}{ccccc} 
    179 key                         & \jp{jp\_cfg} &  \jp{jpiglo} & \jp{jpiglo} &       \\   
    180 \hline  \hline 
    181 \key{orca\_r4}        &        4         &         92     &      76      &       \\ 
    182 \key{orca\_r2}       &        2         &       182     &    149      &        \\ 
    183 %\key{orca\_r1}       &        1         &       362     &     511     &        \\ 
    184 \key{orca\_r05}     &        05       &       722     &     261     &        \\ 
    185 \key{orca\_r025}   &        025     &      1442    &   1021     &        \\ 
    186 %\key{orca\_r8}       &        8         &      2882    &   2042     &        \\ 
    187 %\key{orca\_r12}     &      12         &      4322    &   3062      &       \\ 
    188 \hline 
    189 \hline 
    190 \end{tabular} 
    191 \caption{ \label{Tab_ORCA}    
    192 Set of predefined ORCA parameters. } 
    193 \end{center} 
    194 \end{table} 
    195 %-------------------------------------------------------------------------------------------------------------- 
    196  
     199ORCA\_R2 can also be run with an AGRIF zoom over the Agulhas current area ( \key{agrif}  defined)  
     200and,  by setting the key \key{arctic} or \key{antarctic}, a regional Arctic or peri-Antarctic configuration  
     201is extracted from an ORCA\_R2 or R05 configurations.  
    197202 
    198203% ------------------------------------------------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.