Changeset 11537 for NEMO/trunk/doc/latex/NEMO/subfiles/chap_DYN.tex
- Timestamp:
- 2019-09-12T10:24:48+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/doc/latex/NEMO/subfiles/chap_DYN.tex
r11435 r11537 192 192 193 193 Options are defined through the \nam{dyn\_vor} namelist variables. 194 Four discretisations of the vorticity term (\texttt{ln\_dynvor\_xxx}\forcode{ =.true.}) are available:194 Four discretisations of the vorticity term (\texttt{ln\_dynvor\_xxx}\forcode{=.true.}) are available: 195 195 conserving potential enstrophy of horizontally non-divergent flow (ENS scheme); 196 196 conserving horizontal kinetic energy (ENE scheme); … … 200 200 (EEN scheme) (see \autoref{subsec:C_vorEEN}). 201 201 In the case of ENS, ENE or MIX schemes the land sea mask may be slightly modified to ensure the consistency of 202 vorticity term with analytical equations (\np{ln\_dynvor\_con}\forcode{ =.true.}).202 vorticity term with analytical equations (\np{ln\_dynvor\_con}\forcode{=.true.}). 203 203 The vorticity terms are all computed in dedicated routines that can be found in the \mdl{dynvor} module. 204 204 … … 206 206 % enstrophy conserving scheme 207 207 %------------------------------------------------------------- 208 \subsubsection[Enstrophy conserving scheme (\forcode{ln_dynvor_ens =.true.})]209 {Enstrophy conserving scheme (\protect\np{ln\_dynvor\_ens}\forcode{ =.true.})}208 \subsubsection[Enstrophy conserving scheme (\forcode{ln_dynvor_ens=.true.})] 209 {Enstrophy conserving scheme (\protect\np{ln\_dynvor\_ens}\forcode{=.true.})} 210 210 \label{subsec:DYN_vor_ens} 211 211 … … 230 230 % energy conserving scheme 231 231 %------------------------------------------------------------- 232 \subsubsection[Energy conserving scheme (\forcode{ln_dynvor_ene =.true.})]233 {Energy conserving scheme (\protect\np{ln\_dynvor\_ene}\forcode{ =.true.})}232 \subsubsection[Energy conserving scheme (\forcode{ln_dynvor_ene=.true.})] 233 {Energy conserving scheme (\protect\np{ln\_dynvor\_ene}\forcode{=.true.})} 234 234 \label{subsec:DYN_vor_ene} 235 235 … … 251 251 % mix energy/enstrophy conserving scheme 252 252 %------------------------------------------------------------- 253 \subsubsection[Mixed energy/enstrophy conserving scheme (\forcode{ln_dynvor_mix =.true.})]254 {Mixed energy/enstrophy conserving scheme (\protect\np{ln\_dynvor\_mix}\forcode{ =.true.})}253 \subsubsection[Mixed energy/enstrophy conserving scheme (\forcode{ln_dynvor_mix=.true.})] 254 {Mixed energy/enstrophy conserving scheme (\protect\np{ln\_dynvor\_mix}\forcode{=.true.})} 255 255 \label{subsec:DYN_vor_mix} 256 256 … … 277 277 % energy and enstrophy conserving scheme 278 278 %------------------------------------------------------------- 279 \subsubsection[Energy and enstrophy conserving scheme (\forcode{ln_dynvor_een =.true.})]280 {Energy and enstrophy conserving scheme (\protect\np{ln\_dynvor\_een}\forcode{ =.true.})}279 \subsubsection[Energy and enstrophy conserving scheme (\forcode{ln_dynvor_een=.true.})] 280 {Energy and enstrophy conserving scheme (\protect\np{ln\_dynvor\_een}\forcode{=.true.})} 281 281 \label{subsec:DYN_vor_een} 282 282 … … 328 328 A key point in \autoref{eq:een_e3f} is how the averaging in the \textbf{i}- and \textbf{j}- directions is made. 329 329 It uses the sum of masked t-point vertical scale factor divided either by the sum of the four t-point masks 330 (\np{nn\_een\_e3f}\forcode{ = 1}), or just by $4$ (\np{nn\_een\_e3f}\forcode{ =.true.}).330 (\np{nn\_een\_e3f}\forcode{=1}), or just by $4$ (\np{nn\_een\_e3f}\forcode{=.true.}). 331 331 The latter case preserves the continuity of $e_{3f}$ when one or more of the neighbouring $e_{3t}$ tends to zero and 332 332 extends by continuity the value of $e_{3f}$ into the land areas. … … 410 410 \right. 411 411 \] 412 When \np{ln\_dynzad\_zts}\forcode{ =.true.},412 When \np{ln\_dynzad\_zts}\forcode{=.true.}, 413 413 a split-explicit time stepping with 5 sub-timesteps is used on the vertical advection term. 414 414 This option can be useful when the value of the timestep is limited by vertical advection \citep{lemarie.debreu.ea_OM15}. … … 495 495 % 2nd order centred scheme 496 496 %------------------------------------------------------------- 497 \subsubsection[CEN2: $2^{nd}$ order centred scheme (\forcode{ln_dynadv_cen2 =.true.})]498 {CEN2: $2^{nd}$ order centred scheme (\protect\np{ln\_dynadv\_cen2}\forcode{ =.true.})}497 \subsubsection[CEN2: $2^{nd}$ order centred scheme (\forcode{ln_dynadv_cen2=.true.})] 498 {CEN2: $2^{nd}$ order centred scheme (\protect\np{ln\_dynadv\_cen2}\forcode{=.true.})} 499 499 \label{subsec:DYN_adv_cen2} 500 500 … … 519 519 % UBS scheme 520 520 %------------------------------------------------------------- 521 \subsubsection[UBS: Upstream Biased Scheme (\forcode{ln_dynadv_ubs =.true.})]522 {UBS: Upstream Biased Scheme (\protect\np{ln\_dynadv\_ubs}\forcode{ =.true.})}521 \subsubsection[UBS: Upstream Biased Scheme (\forcode{ln_dynadv_ubs=.true.})] 522 {UBS: Upstream Biased Scheme (\protect\np{ln\_dynadv\_ubs}\forcode{=.true.})} 523 523 \label{subsec:DYN_adv_ubs} 524 524 … … 542 542 But the amplitudes of the false extrema are significantly reduced over those in the centred second order method. 543 543 As the scheme already includes a diffusion component, it can be used without explicit lateral diffusion on momentum 544 (\ie\ \np{ln\_dynldf\_lap}\forcode{ = }\np{ln\_dynldf\_bilap}\forcode{ =.false.}),544 (\ie\ \np{ln\_dynldf\_lap}\forcode{=}\np{ln\_dynldf\_bilap}\forcode{=.false.}), 545 545 and it is recommended to do so. 546 546 … … 596 596 % z-coordinate with full step 597 597 %-------------------------------------------------------------------------------------------------------------- 598 \subsection[Full step $Z$-coordinate (\forcode{ln_dynhpg_zco =.true.})]599 {Full step $Z$-coordinate (\protect\np{ln\_dynhpg\_zco}\forcode{ =.true.})}598 \subsection[Full step $Z$-coordinate (\forcode{ln_dynhpg_zco=.true.})] 599 {Full step $Z$-coordinate (\protect\np{ln\_dynhpg\_zco}\forcode{=.true.})} 600 600 \label{subsec:DYN_hpg_zco} 601 601 … … 642 642 % z-coordinate with partial step 643 643 %-------------------------------------------------------------------------------------------------------------- 644 \subsection[Partial step $Z$-coordinate (\forcode{ln_dynhpg_zps =.true.})]645 {Partial step $Z$-coordinate (\protect\np{ln\_dynhpg\_zps}\forcode{ =.true.})}644 \subsection[Partial step $Z$-coordinate (\forcode{ln_dynhpg_zps=.true.})] 645 {Partial step $Z$-coordinate (\protect\np{ln\_dynhpg\_zps}\forcode{=.true.})} 646 646 \label{subsec:DYN_hpg_zps} 647 647 … … 672 672 density Jacobian with cubic polynomial method is currently disabled whilst known bugs are under investigation. 673 673 674 $\bullet$ Traditional coding (see for example \citet{madec.delecluse.ea_JPO96}: (\np{ln\_dynhpg\_sco}\forcode{ =.true.})674 $\bullet$ Traditional coding (see for example \citet{madec.delecluse.ea_JPO96}: (\np{ln\_dynhpg\_sco}\forcode{=.true.}) 675 675 \begin{equation} 676 676 \label{eq:dynhpg_sco} … … 690 690 ($e_{3w}$). 691 691 692 $\bullet$ Traditional coding with adaptation for ice shelf cavities (\np{ln\_dynhpg\_isf}\forcode{ =.true.}).693 This scheme need the activation of ice shelf cavities (\np{ln\_isfcav}\forcode{ =.true.}).694 695 $\bullet$ Pressure Jacobian scheme (prj) (a research paper in preparation) (\np{ln\_dynhpg\_prj}\forcode{ =.true.})692 $\bullet$ Traditional coding with adaptation for ice shelf cavities (\np{ln\_dynhpg\_isf}\forcode{=.true.}). 693 This scheme need the activation of ice shelf cavities (\np{ln\_isfcav}\forcode{=.true.}). 694 695 $\bullet$ Pressure Jacobian scheme (prj) (a research paper in preparation) (\np{ln\_dynhpg\_prj}\forcode{=.true.}) 696 696 697 697 $\bullet$ Density Jacobian with cubic polynomial scheme (DJC) \citep{shchepetkin.mcwilliams_OM05} 698 (\np{ln\_dynhpg\_djc}\forcode{ =.true.}) (currently disabled; under development)698 (\np{ln\_dynhpg\_djc}\forcode{=.true.}) (currently disabled; under development) 699 699 700 700 Note that expression \autoref{eq:dynhpg_sco} is commonly used when the variable volume formulation is activated 701 701 (\texttt{vvl?}) because in that case, even with a flat bottom, 702 702 the coordinate surfaces are not horizontal but follow the free surface \citep{levier.treguier.ea_rpt07}. 703 The pressure jacobian scheme (\np{ln\_dynhpg\_prj}\forcode{ =.true.}) is available as704 an improved option to \np{ln\_dynhpg\_sco}\forcode{ =.true.} when \texttt{vvl?} is active.703 The pressure jacobian scheme (\np{ln\_dynhpg\_prj}\forcode{=.true.}) is available as 704 an improved option to \np{ln\_dynhpg\_sco}\forcode{=.true.} when \texttt{vvl?} is active. 705 705 The pressure Jacobian scheme uses a constrained cubic spline to 706 706 reconstruct the density profile across the water column. … … 713 713 \label{subsec:DYN_hpg_isf} 714 714 Beneath an ice shelf, the total pressure gradient is the sum of the pressure gradient due to the ice shelf load and 715 the pressure gradient due to the ocean load (\np{ln\_dynhpg\_isf}\forcode{ =.true.}).\\715 the pressure gradient due to the ocean load (\np{ln\_dynhpg\_isf}\forcode{=.true.}).\\ 716 716 717 717 The main hypothesis to compute the ice shelf load is that the ice shelf is in an isostatic equilibrium. … … 728 728 % Time-scheme 729 729 %-------------------------------------------------------------------------------------------------------------- 730 \subsection[Time-scheme (\forcode{ln_dynhpg_imp = .{true,false}.})]731 {Time-scheme (\protect\np{ln\_dynhpg\_imp}\forcode{ = .\{true,false\}}.)}730 \subsection[Time-scheme (\forcode{ln_dynhpg_imp={.true.,.false.}})] 731 {Time-scheme (\protect\np{ln\_dynhpg\_imp}\forcode{=.true.,.false.})} 732 732 \label{subsec:DYN_hpg_imp} 733 733 … … 745 745 rather than at the central time level $t$ only, as in the standard leapfrog scheme. 746 746 747 $\bullet$ leapfrog scheme (\np{ln\_dynhpg\_imp}\forcode{ =.true.}):747 $\bullet$ leapfrog scheme (\np{ln\_dynhpg\_imp}\forcode{=.true.}): 748 748 749 749 \begin{equation} … … 753 753 \end{equation} 754 754 755 $\bullet$ semi-implicit scheme (\np{ln\_dynhpg\_imp}\forcode{ =.true.}):755 $\bullet$ semi-implicit scheme (\np{ln\_dynhpg\_imp}\forcode{=.true.}): 756 756 \begin{equation} 757 757 \label{eq:dynhpg_imp} … … 771 771 such as the stability limits associated with advection or diffusion. 772 772 773 In practice, the semi-implicit scheme is used when \np{ln\_dynhpg\_imp}\forcode{ =.true.}.773 In practice, the semi-implicit scheme is used when \np{ln\_dynhpg\_imp}\forcode{=.true.}. 774 774 In this case, we choose to apply the time filter to temperature and salinity used in the equation of state, 775 775 instead of applying it to the hydrostatic pressure or to the density, … … 829 829 % Explicit free surface formulation 830 830 %-------------------------------------------------------------------------------------------------------------- 831 \subsection[Explicit free surface (\texttt{ln\_dynspg\_exp}\forcode{ =.true.})]832 {Explicit free surface (\protect\np{ln\_dynspg\_exp}\forcode{ =.true.})}831 \subsection[Explicit free surface (\texttt{ln\_dynspg\_exp}\forcode{=.true.})] 832 {Explicit free surface (\protect\np{ln\_dynspg\_exp}\forcode{=.true.})} 833 833 \label{subsec:DYN_spg_exp} 834 834 … … 856 856 % Split-explict free surface formulation 857 857 %-------------------------------------------------------------------------------------------------------------- 858 \subsection[Split-explicit free surface (\texttt{ln\_dynspg\_ts}\forcode{ =.true.})]859 {Split-explicit free surface (\protect\np{ln\_dynspg\_ts}\forcode{ =.true.})}858 \subsection[Split-explicit free surface (\texttt{ln\_dynspg\_ts}\forcode{=.true.})] 859 {Split-explicit free surface (\protect\np{ln\_dynspg\_ts}\forcode{=.true.})} 860 860 \label{subsec:DYN_spg_ts} 861 861 %------------------------------------------namsplit----------------------------------------------------------- … … 871 871 The size of the small time step, $\rdt_e$ (the external mode or barotropic time step) is provided through 872 872 the \np{nn\_baro} namelist parameter as: $\rdt_e = \rdt / nn\_baro$. 873 This parameter can be optionally defined automatically (\np{ln\_bt\_nn\_auto}\forcode{ =.true.}) considering that873 This parameter can be optionally defined automatically (\np{ln\_bt\_nn\_auto}\forcode{=.true.}) considering that 874 874 the stability of the barotropic system is essentially controled by external waves propagation. 875 875 Maximum Courant number is in that case time independent, and easily computed online from the input bathymetry. … … 916 916 The former are used to obtain time filtered quantities at $t+\rdt$ while 917 917 the latter are used to obtain time averaged transports to advect tracers. 918 a) Forward time integration: \protect\np{ln\_bt\_fw}\forcode{ =.true.},919 \protect\np{ln\_bt\_av}\forcode{ =.true.}.920 b) Centred time integration: \protect\np{ln\_bt\_fw}\forcode{ =.false.},921 \protect\np{ln\_bt\_av}\forcode{ =.true.}.918 a) Forward time integration: \protect\np{ln\_bt\_fw}\forcode{=.true.}, 919 \protect\np{ln\_bt\_av}\forcode{=.true.}. 920 b) Centred time integration: \protect\np{ln\_bt\_fw}\forcode{=.false.}, 921 \protect\np{ln\_bt\_av}\forcode{=.true.}. 922 922 c) Forward time integration with no time filtering (POM-like scheme): 923 \protect\np{ln\_bt\_fw}\forcode{ = .true.}, \protect\np{ln\_bt\_av}\forcode{ =.false.}.923 \protect\np{ln\_bt\_fw}\forcode{=.true.}, \protect\np{ln\_bt\_av}\forcode{=.false.}. 924 924 } 925 925 \end{center} … … 927 927 %> > > > > > > > > > > > > > > > > > > > > > > > > > > > 928 928 929 In the default case (\np{ln\_bt\_fw}\forcode{ =.true.}),929 In the default case (\np{ln\_bt\_fw}\forcode{=.true.}), 930 930 the external mode is integrated between \textit{now} and \textit{after} baroclinic time-steps 931 931 (\autoref{fig:DYN_dynspg_ts}a). 932 932 To avoid aliasing of fast barotropic motions into three dimensional equations, 933 time filtering is eventually applied on barotropic quantities (\np{ln\_bt\_av}\forcode{ =.true.}).933 time filtering is eventually applied on barotropic quantities (\np{ln\_bt\_av}\forcode{=.true.}). 934 934 In that case, the integration is extended slightly beyond \textit{after} time step to 935 935 provide time filtered quantities. … … 938 938 asselin filtering is not applied to barotropic quantities.\\ 939 939 Alternatively, one can choose to integrate barotropic equations starting from \textit{before} time step 940 (\np{ln\_bt\_fw}\forcode{ =.false.}).940 (\np{ln\_bt\_fw}\forcode{=.false.}). 941 941 Although more computationaly expensive ( \np{nn\_baro} additional iterations are indeed necessary), 942 942 the baroclinic to barotropic forcing term given at \textit{now} time step become centred in … … 963 963 964 964 One can eventually choose to feedback instantaneous values by not using any time filter 965 (\np{ln\_bt\_av}\forcode{ =.false.}).965 (\np{ln\_bt\_av}\forcode{=.false.}). 966 966 In that case, external mode equations are continuous in time, 967 967 \ie\ they are not re-initialized when starting a new sub-stepping sequence. … … 1165 1165 1166 1166 % ================================================================ 1167 \subsection[Iso-level laplacian (\forcode{ln_dynldf_lap =.true.})]1168 {Iso-level laplacian operator (\protect\np{ln\_dynldf\_lap}\forcode{ =.true.})}1167 \subsection[Iso-level laplacian (\forcode{ln_dynldf_lap=.true.})] 1168 {Iso-level laplacian operator (\protect\np{ln\_dynldf\_lap}\forcode{=.true.})} 1169 1169 \label{subsec:DYN_ldf_lap} 1170 1170 … … 1191 1191 % Rotated laplacian operator 1192 1192 %-------------------------------------------------------------------------------------------------------------- 1193 \subsection[Rotated laplacian (\forcode{ln_dynldf_iso =.true.})]1194 {Rotated laplacian operator (\protect\np{ln\_dynldf\_iso}\forcode{ =.true.})}1193 \subsection[Rotated laplacian (\forcode{ln_dynldf_iso=.true.})] 1194 {Rotated laplacian operator (\protect\np{ln\_dynldf\_iso}\forcode{=.true.})} 1195 1195 \label{subsec:DYN_ldf_iso} 1196 1196 1197 1197 A rotation of the lateral momentum diffusion operator is needed in several cases: 1198 for iso-neutral diffusion in the $z$-coordinate (\np{ln\_dynldf\_iso}\forcode{ =.true.}) and1199 for either iso-neutral (\np{ln\_dynldf\_iso}\forcode{ =.true.}) or1200 geopotential (\np{ln\_dynldf\_hor}\forcode{ =.true.}) diffusion in the $s$-coordinate.1198 for iso-neutral diffusion in the $z$-coordinate (\np{ln\_dynldf\_iso}\forcode{=.true.}) and 1199 for either iso-neutral (\np{ln\_dynldf\_iso}\forcode{=.true.}) or 1200 geopotential (\np{ln\_dynldf\_hor}\forcode{=.true.}) diffusion in the $s$-coordinate. 1201 1201 In the partial step case, coordinates are horizontal except at the deepest level and 1202 no rotation is performed when \np{ln\_dynldf\_hor}\forcode{ =.true.}.1202 no rotation is performed when \np{ln\_dynldf\_hor}\forcode{=.true.}. 1203 1203 The diffusion operator is defined simply as the divergence of down gradient momentum fluxes on 1204 1204 each momentum component. … … 1250 1250 % Iso-level bilaplacian operator 1251 1251 %-------------------------------------------------------------------------------------------------------------- 1252 \subsection[Iso-level bilaplacian (\forcode{ln_dynldf_bilap =.true.})]1253 {Iso-level bilaplacian operator (\protect\np{ln\_dynldf\_bilap}\forcode{ =.true.})}1252 \subsection[Iso-level bilaplacian (\forcode{ln_dynldf_bilap=.true.})] 1253 {Iso-level bilaplacian operator (\protect\np{ln\_dynldf\_bilap}\forcode{=.true.})} 1254 1254 \label{subsec:DYN_ldf_bilap} 1255 1255 … … 1277 1277 Two time stepping schemes can be used for the vertical diffusion term: 1278 1278 $(a)$ a forward time differencing scheme 1279 (\np{ln\_zdfexp}\forcode{ =.true.}) using a time splitting technique (\np{nn\_zdfexp} $>$ 1) or1280 $(b)$ a backward (or implicit) time differencing scheme (\np{ln\_zdfexp}\forcode{ =.false.})1279 (\np{ln\_zdfexp}\forcode{=.true.}) using a time splitting technique (\np{nn\_zdfexp} $>$ 1) or 1280 $(b)$ a backward (or implicit) time differencing scheme (\np{ln\_zdfexp}\forcode{=.false.}) 1281 1281 (see \autoref{chap:STP}). 1282 1282 Note that namelist variables \np{ln\_zdfexp} and \np{nn\_zdfexp} apply to both tracers and dynamics. … … 1328 1328 three other forcings may enter the dynamical equations by affecting the surface pressure gradient. 1329 1329 1330 (1) When \np{ln\_apr\_dyn}\forcode{ =.true.} (see \autoref{sec:SBC_apr}),1330 (1) When \np{ln\_apr\_dyn}\forcode{=.true.} (see \autoref{sec:SBC_apr}), 1331 1331 the atmospheric pressure is taken into account when computing the surface pressure gradient. 1332 1332 1333 (2) When \np{ln\_tide\_pot}\forcode{ = .true.} and \np{ln\_tide}\forcode{ =.true.} (see \autoref{sec:SBC_tide}),1333 (2) When \np{ln\_tide\_pot}\forcode{=.true.} and \np{ln\_tide}\forcode{=.true.} (see \autoref{sec:SBC_tide}), 1334 1334 the tidal potential is taken into account when computing the surface pressure gradient. 1335 1335 1336 (3) When \np{nn\_ice\_embd}\forcode{ =2} and LIM or CICE is used1336 (3) When \np{nn\_ice\_embd}\forcode{=2} and LIM or CICE is used 1337 1337 (\ie\ when the sea-ice is embedded in the ocean), 1338 1338 the snow-ice mass is taken into account when computing the surface pressure gradient. … … 1409 1409 1410 1410 The flux across each $u$-face of a tracer cell is multiplied by a factor zuwdmask (an array which depends on ji and jj). 1411 If the user sets \np{ln\_wd\_dl\_ramp}\forcode{ =.false.} then zuwdmask is 1 when the1411 If the user sets \np{ln\_wd\_dl\_ramp}\forcode{=.false.} then zuwdmask is 1 when the 1412 1412 flux is from a cell with water depth greater than \np{rn\_wdmin1} and 0 otherwise. If the user sets 1413 \np{ln\_wd\_dl\_ramp}\forcode{ =.true.} the flux across the face is ramped down as the water depth decreases1413 \np{ln\_wd\_dl\_ramp}\forcode{=.true.} the flux across the face is ramped down as the water depth decreases 1414 1414 from 2 * \np{rn\_wdmin1} to \np{rn\_wdmin1}. The use of this ramp reduced grid-scale noise in idealised test cases. 1415 1415 … … 1428 1428 fields (tracers independent of $x$, $y$ and $z$). Our scheme conserves constant tracers because 1429 1429 the velocities used at the tracer cell faces on the baroclinic timesteps are carefully calculated by dynspg\_ts 1430 to equal their mean value during the barotropic steps. If the user sets \np{ln\_wd\_dl\_bc}\forcode{ =.true.}, the1430 to equal their mean value during the barotropic steps. If the user sets \np{ln\_wd\_dl\_bc}\forcode{=.true.}, the 1431 1431 baroclinic velocities are also multiplied by a suitably weighted average of zuwdmask. 1432 1432 … … 1655 1655 1656 1656 $\bullet$ vector invariant form or linear free surface 1657 (\np{ln\_dynhpg\_vec}\forcode{ =.true.} ; \texttt{vvl?} not defined):1657 (\np{ln\_dynhpg\_vec}\forcode{=.true.} ; \texttt{vvl?} not defined): 1658 1658 \[ 1659 1659 % \label{eq:dynnxt_vec} … … 1667 1667 1668 1668 $\bullet$ flux form and nonlinear free surface 1669 (\np{ln\_dynhpg\_vec}\forcode{ =.false.} ; \texttt{vvl?} defined):1669 (\np{ln\_dynhpg\_vec}\forcode{=.false.} ; \texttt{vvl?} defined): 1670 1670 \[ 1671 1671 % \label{eq:dynnxt_flux} … … 1681 1681 the subscript $f$ denotes filtered values and $\gamma$ is the Asselin coefficient. 1682 1682 $\gamma$ is initialized as \np{nn\_atfp} (namelist parameter). 1683 Its default value is \np{nn\_atfp}\forcode{ =10.e-3}.1683 Its default value is \np{nn\_atfp}\forcode{=10.e-3}. 1684 1684 In both cases, the modified Asselin filter is not applied since perfect conservation is not an issue for 1685 1685 the momentum equations.
Note: See TracChangeset
for help on using the changeset viewer.