- Timestamp:
- 2017-02-10T18:16:26+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_merge_2016/DOC/TexFiles/Chapters/Chap_DOM.tex
r6997 r7669 18 18 19 19 \newpage 20 $\ $\newline % force a new li gne20 $\ $\newline % force a new line 21 21 22 22 Having defined the continuous equations in Chap.~\ref{PE} and chosen a time … … 26 26 directory routines as well as the DOM (DOMain) directory. 27 27 28 $\ $\newline % force a new li gn28 $\ $\newline % force a new line 29 29 30 30 % ================================================================ … … 268 268 269 269 The total size of the computational domain is set by the parameters \np{jpiglo}, 270 \np{jpjglo} and \np{jpk dta} in the $i$, $j$ and $k$ directions respectively. They are270 \np{jpjglo} and \np{jpkglo} in the $i$, $j$ and $k$ directions respectively. They are 271 271 given as namelist variables in the \ngn{namcfg} namelist. 272 273 Note that are other namelist variables in the \ngn{namcfg} namelist that refer to 274 the domain size. 275 The two variables \np{jpidta} and \np{jpjdta} may be larger than \np{jpiglo}, \np{jpjglo} 276 when the user wants to use only a sub-region of a given configuration. This is 277 the "zoom" capability described in \S\ref{MISC_zoom}. In most applications of 278 the model, $jpidta=jpiglo$, $jpjdta=jpjglo$, and $jpizoom=jpjzoom=1$. Parameters 279 $jpi$ and $jpj$ refer to the size of each processor subdomain when the code is 280 run in parallel using domain decomposition (\key{mpp\_mpi} defined, see 281 \S\ref{LBC_mpp}). 282 283 284 $\ $\newline % force a new lign 272 %%% 273 %%% 274 %%% 275 Parameters $jpi$ and $jpj$ refer to the size of each processor subdomain when the code is 276 run in parallel using domain decomposition (\key{mpp\_mpi} defined, see \S\ref{LBC_mpp}). 277 278 279 $\ $\newline % force a new line 280 281 %%% 282 \sfcomment {Hereafter I want to create new subsection 4.2: "fields needed by opa engine or something like this" 283 and add list of fields : 284 case 1: read in domain.nc 285 case 2: defined in userdef\_hrg\/zgr.F90 286 longitude, latitude, domaine size 287 number of points 288 factor scales (e1, e2, e3) 289 coriolis 290 k\_top, k\_bottom (first and last ocean level) 291 periodicity 292 } 293 %%% 294 295 % ================================================================ 296 % Domain: List of fields needed 297 % ================================================================ 298 \section [Domain: Needed fields] 299 {Domain: Needed fields} 300 \label{DOM_fields} 301 The ocean mesh ($i.e.$ the position of all the scalar and vector points) is defined 302 by the transformation that gives $(\lambda,\varphi,z)$ as a function of $(i,j,k)$. 303 The grid-points are located at integer or integer and a half values of as indicated 304 in Table~\ref{Tab_cell}. The associated scale factors are defined using the 305 analytical first derivative of the transformation \eqref{Eq_scale_factors}. These 306 definitions are done in two modules given by example, \mdl{userdef\_hgr} and \mdl{userdef\_zgr}, which 307 provide the horizontal and vertical meshes, respectively. Otherwise all needed fields can be read in file \np{cn\_domcfg} specified in \ngn{namcfg}. 308 309 The needed fields for domain are: 310 311 geographic position : 312 313 longitude : glamt , glamu , glamv and glamf (at T, U, V and F point) 314 315 latitude : gphit , gphiu , gphiv and gphif (at T, U, V and F point) 316 317 Coriolis parameter (if domain not on the sphere): ff\_f and ff\_t (at T and F point) 318 319 Scale factors : e1t, e1u, e1v and e1f (on i direction), 320 321 e2t, e2u, e2v and e2f (on j direction) 322 323 and ie1e2u\_v, e1e2u , e1e2v 324 325 %%% 326 \sfcomment { 327 say something about ie1e2u\_v, e1e2u , e1e2v 328 329 and add list of fields : 330 case 1: read in domain.nc 331 case 2: defined in userdef\_hrg\/zgr.F90 332 longitude, latitude, domaine size 333 number of points 334 factor scales (e1, e2, e3) 335 coriolis 336 k\_top, k\_bottom (first and last ocean level) 337 periodicity 338 ---- 339 int ORCA ; 340 int ORCA\_index ; 341 int jpiglo ; j, k 342 int jperio ; 343 int ln_zco ; zps, sco 344 int ln_isfcav ; 345 double glamt(t, y, x) ; u,v,f 346 double gphit(t, y, x) ; u,v,f 347 double e1t(t, y, x) ; u,v,w, 348 double e2t(t, y, x) ; u,v,w 349 double ff\_f(t, y, x) ; double ff\_t(t, y, x) ; 350 double e3t\_1d(t, z) ; 351 double e3w\_1d(t, z) ; 352 double e3t\_0(t, z, y, x) ; u0, v0 , w0 353 ---- 354 } 355 356 % ------------------------------------------------------------------------------------------------------------- 357 % Needed fields 358 % ------------------------------------------------------------------------------------------------------------- 359 %\subsection{List of needed fields to build DOMAIN} 360 %\label{DOM_fields_list} 361 285 362 286 363 % ================================================================ … … 436 513 the output grid written when $\np{nn\_msh} \not=0$ is no more equal to the input grid. 437 514 438 $\ $\newline % force a new li gn515 $\ $\newline % force a new line 439 516 440 517 % ================================================================
Note: See TracChangeset
for help on using the changeset viewer.