Changeset 9019 for branches/2017/dev_merge_2017/DOC/TexFiles/Chapters/Chap_DOM.tex

Timestamp:

2017-12-13T15:58:53+01:00 (6 years ago)

Author:

timgraham

Message:

Merge of dev_CNRS_2017 into branch

File:

• branches/2017/dev_merge_2017/DOC/TexFiles/Chapters/Chap_DOM.tex (modified) (6 diffs)

branches/2017/dev_merge_2017/DOC/TexFiles/Chapters/Chap_DOM.tex

@@ -268 +268 @@
 
 The total size of the computational domain is set by the parameters \np{jpiglo}, 
-\np{jpjglo} and \np{jpkglo} in the $i$, $j$ and $k$ directions respectively. They are 
-given as namelist variables in the \ngn{namcfg} namelist. 
+\np{jpjglo} and \np{jpkglo} in the $i$, $j$ and $k$ directions respectively. 
 %%%
 %%%
@@ -278 +277 @@
 
 $\ $\newline    % force a new line
-
-%%%
-\sfcomment {Hereafter I want to create new subsection 4.2: "fields needed by opa engine or something like this"
-and add list of fields :
-case 1: read in domain.nc
-case 2: defined in userdef\_hrg\/zgr.F90
-longitude, latitude, domaine size 
-number of points
-factor scales (e1, e2, e3)
-coriolis
-k\_top, k\_bottom (first and last ocean level)
-periodicity
-}
-%%%
 
 % ================================================================
@@ -303 +288 @@
 The grid-points are located at integer or integer and a half values of as indicated 
 in Table~\ref{Tab_cell}. The associated scale factors are defined using the  
-analytical first derivative of the transformation \eqref{Eq_scale_factors}. These 
-definitions are done in two modules given by example, \mdl{userdef\_hgr} and \mdl{userdef\_zgr}, which 
-provide the horizontal and vertical meshes, respectively. Otherwise all needed fields can be read in file \np{cn\_domcfg} specified in \ngn{namcfg}.
+analytical first derivative of the transformation \eqref{Eq_scale_factors}. 
+Necessary fields for configuration definition are: \\
+Geographic position :
+
+longitude : glamt , glamu , glamv and glamf  (at T, U, V and F point)
+
+latitude : gphit , gphiu , gphiv and gphif (at T, U, V and F point)\\
+Coriolis parameter (if domain not on the sphere): 
+
+ ff\_f  and  ff\_t (at T and F point)\\
+Scale factors : 
  
-The needed fields for domain are: 
-
-geographic position :
-
-longitude : glamt , glamu , glamv and glamf  (at T, U, V and F point)
-
-latitude : gphit , gphiu , gphiv and gphif (at T, U, V and F point)
-
-Coriolis parameter (if domain not on the sphere):  ff\_f  and  ff\_t (at T and F point)
-
-Scale factors : e1t, e1u, e1v and e1f (on i direction),
-
-   e2t, e2u, e2v and e2f (on j direction)
-
-   and ie1e2u\_v, e1e2u , e1e2v   
-
-%%%
-\sfcomment {
-say something about ie1e2u\_v, e1e2u , e1e2v
-
-and add list of fields :
-case 1: read in domain.nc
-case 2: defined in userdef\_hrg\/zgr.F90
-longitude, latitude, domaine size
-number of points
-factor scales (e1, e2, e3)
-coriolis
-k\_top, k\_bottom (first and last ocean level)
-periodicity
-----
-        int ORCA ;
-   int ORCA\_index ;
-   int jpiglo ; j, k
-   int jperio ;
-   int ln_zco ; zps, sco
-   int ln_isfcav ;
-   double glamt(t, y, x) ; u,v,f
-   double gphit(t, y, x) ; u,v,f
-   double e1t(t, y, x) ; u,v,w,
-   double e2t(t, y, x) ; u,v,w
-   double ff\_f(t, y, x) ;  double ff\_t(t, y, x) ;
-   double e3t\_1d(t, z) ;
-   double e3w\_1d(t, z) ;
-   double e3t\_0(t, z, y, x) ; u0, v0 , w0
-----
-}
-
+ e1t, e1u, e1v and e1f (on i direction),
+
+ e2t, e2u, e2v and e2f (on j direction)
+
+ and ie1e2u\_v, e1e2u , e1e2v   
+ 
+e1e2u , e1e2v are u and v surfaces (if gridsize reduction in some straits)\\
+ie1e2u\_v is a flag to flag set u and  v surfaces are neither read nor computed.\\
+ 
+These fields can be read in an domain input file which name is setted in \np{cn\_domcfg} parameter specified in \ngn{namcfg}.
+\namdisplay{namcfg}
+or they can be defined in an analytical way in MY\_SRC directory of the configuration.
+For Reference Configurations of NEMO input domain files are supplied by NEMO System Team. For analytical definition of input fields two routines are supplied: \mdl{userdef\_hgr} and \mdl{userdef\_zgr}. They are an example of GYRE configuration parameters, and they are available in NEMO/OPA\_SRC/USR directory, they provide the horizontal and vertical mesh. 
 % -------------------------------------------------------------------------------------------------------------
 %        Needed fields 
@@ -446 +405 @@
 \label{DOM_hgr_msh_choice}
 
-The user has three options available in defining a horizontal grid, which involve 
-the namelist variable \np{jphgr\_mesh} of the \ngn{namcfg} namelist. 
-\begin{description}
-\item[\np{jphgr\_mesh}=0]  The most general curvilinear orthogonal grids.
-The coordinates and their first derivatives with respect to $i$ and $j$ are provided
-in a input file (\ifile{coordinates}), read in \rou{hgr\_read} subroutine of the domhgr module.
-\item[\np{jphgr\_mesh}=1 to 5] A few simple analytical grids are provided (see below). 
-For other analytical grids, the \mdl{domhgr} module must be modified by the user. 
-\end{description}
-
-There are two simple cases of geographical grids on the sphere. With 
-\np{jphgr\_mesh}=1, the grid (expressed in degrees) is regular in space, 
-with grid sizes specified by parameters \np{ppe1\_deg} and \np{ppe2\_deg}, 
-respectively. Such a geographical grid can be very anisotropic at high latitudes 
-because of the convergence of meridians (the zonal scale factors $e_1$ 
-become much smaller than the meridional scale factors $e_2$). The Mercator 
-grid (\np{jphgr\_mesh}=4) avoids this anisotropy by refining the meridional scale 
-factors in the same way as the zonal ones. In this case, meridional scale factors 
-and latitudes are calculated analytically using the formulae appropriate for 
-a Mercator projection, based on \np{ppe1\_deg} which is a reference grid spacing 
-at the equator (this applies even when the geographical equator is situated outside 
-the model domain). 
-%%%
-\gmcomment{ give here the analytical expression of the Mercator mesh}
-%%%
-In these two cases (\np{jphgr\_mesh}=1 or 4), the grid position is defined by the 
-longitude and latitude of the south-westernmost point (\np{ppglamt0} 
-and \np{ppgphi0}). Note that for the Mercator grid the user need only provide 
-an approximate starting latitude: the real latitude will be recalculated analytically, 
-in order to ensure that the equator corresponds to line passing through $t$- 
-and $u$-points.  
-
-Rectangular grids ignoring the spherical geometry are defined with 
-\np{jphgr\_mesh} = 2, 3, 5. The domain is either an $f$-plane (\np{jphgr\_mesh} = 2, 
-Coriolis factor is constant) or a beta-plane (\np{jphgr\_mesh} = 3, the Coriolis factor 
-is linear in the $j$-direction). The grid size is uniform in meter in each direction, 
-and given by the parameters \np{ppe1\_m} and \np{ppe2\_m} respectively. 
-The zonal grid coordinate (\textit{glam} arrays) is in kilometers, starting at zero 
-with the first $t$-point. The meridional coordinate (gphi. arrays) is in kilometers, 
-and the second $t$-point corresponds to coordinate $gphit=0$. The input 
-variable \np{ppglam0} is ignored. \np{ppgphi0} is used to set the reference 
-latitude for computation of the Coriolis parameter. In the case of the beta plane, 
-\np{ppgphi0} corresponds to the center of the domain. Finally, the special case 
-\np{jphgr\_mesh}=5 corresponds to a beta plane in a rotated domain for the 
-GYRE configuration, representing a classical mid-latitude double gyre system. 
-The rotation allows us to maximize the jet length relative to the gyre areas 
-(and the number of grid points). 
-
-The choice of the grid must be consistent with the boundary conditions specified 
-by \np{jperio}, a parameter found in \ngn{namcfg} namelist (see {\S\ref{LBC}).
 
 % -------------------------------------------------------------------------------------------------------------
@@ -684 +593 @@
 (Fig.~\ref{Fig_zgr}).
 
-If the ice shelf cavities are opened (\np{ln\_isfcav}=~true~}), the definition of $z_0$ is the same. 
+If the ice shelf cavities are opened (\np{ln\_isfcav}=~true~), the definition of $z_0$ is the same. 
 However, definition of $e_3^0$ at $t$- and $w$-points is respectively changed to:
 \begin{equation} \label{DOM_zgr_ana}
@@ -737 +646 @@
 \begin{table}     \begin{center} \begin{tabular}{c||r|r|r|r}
 \hline
-\textbf{LEVEL}& \textbf{gdept}& \textbf{gdepw}& \textbf{e3t }& \textbf{e3w  } \\ \hline
+\textbf{LEVEL}& \textbf{gdept\_1d}& \textbf{gdepw\_1d}& \textbf{e3t\_1d }& \textbf{e3w\_1d  } \\ \hline
 1  &  \textbf{  5.00}   &       0.00 & \textbf{ 10.00} &  10.00 \\   \hline
 2  &  \textbf{15.00} &    10.00 &   \textbf{ 10.00} &  10.00 \\   \hline