Based on a flux formulation of the advection term
TracMath macro processor has detected an error. Please fix the problem before continuing.
The command:
'/usr/bin/pdflatex -interaction=nonstopmode 0b578335815282d81f3d3a0b1ac74e1e506bcb61.tex'
failed with the following output:
"This is pdfTeX, Version 3.1415926-2.5-1.40.14 (TeX Live 2013)\n restricted \\write18 enabled.\nentering extended mode\n(./0b578335815282d81f3d3a0b1ac74e1e506bcb61.tex\nLaTeX2e <2011/06/27>\nBabel <3.9h> and hyphenation patterns for 78 languages loaded.\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/article.cls\nDocument Class: article 2007/10/19 v1.4h Standard LaTeX document class\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/size10.clo))\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/inputenc.sty\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/utf8.def\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/t1enc.dfu)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/ot1enc.dfu)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/omsenc.dfu)))\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/cmap/cmap.sty)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/cm-super/type1ec.sty\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/t1cmr.fd))\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/fontenc.sty\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/t1enc.def)<<t1.cmap\n>>) (/localdata/software/texlive/2013/texmf-dist/tex/latex/amsmath/amsmath.sty\nFor additional information on amsmath, use the `?' option.\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsmath/amstext.sty\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsmath/amsgen.sty))\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsmath/amsbsy.sty)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsmath/amsopn.sty))\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amscls/amsthm.sty)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsfonts/amssymb.sty\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsfonts/amsfonts.sty))\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/tools/bm.sty)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/preview/preview.sty\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/preview/prtightpage.def)\n)\nNo file 0b578335815282d81f3d3a0b1ac74e1e506bcb61.aux.\nPreview: Fontsize 10pt\nPreview: PDFoutput 1\n<<ot1.cmap>><<oml.cmap>><<oms.cmap>><<omx.cmap>>\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsfonts/umsa.fd)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsfonts/umsb.fd)\n! Misplaced alignment tab character &.\nl.14 \\[ \\frac{\\partial \\psi}{\\partial t} &\n = & -\\frac{\\partial ( u \\psi )}{\\...\n\n! Misplaced alignment tab character &.\nl.14 \\[ \\frac{\\partial \\psi}{\\partial t} & = &\n -\\frac{\\partial ( u \\psi )}{\\...\n\nPreview: Tightpage -32891 -32891 32891 32891\n[1{/localdata/software/texlive/2013/texmf-var/fonts/map/pdftex/updmap/pdftex.ma\np}] (./0b578335815282d81f3d3a0b1ac74e1e506bcb61.aux) )\n(see the transcript file for additional information)</localdata/software/texliv\ne/2013/texmf-dist/fonts/type1/public/amsfonts/cm/cmmi10.pfb></localdata/softwar\ne/texlive/2013/texmf-dist/fonts/type1/public/amsfonts/cm/cmr10.pfb></localdata/\nsoftware/texlive/2013/texmf-dist/fonts/type1/public/amsfonts/cm/cmsy10.pfb>\nOutput written on 0b578335815282d81f3d3a0b1ac74e1e506bcb61.pdf (1 page, 28774 b\nytes).\nTranscript written on 0b578335815282d81f3d3a0b1ac74e1e506bcb61.log.\n"
the one dimensional advection equation can written in following semidiscrete form:
![\[ \frac{\partial \psi_{i} }{\partial t} = - \frac{F_{i+\frac{1}{2}}(u\psi) - F_{i-\frac{1}{2}}(u\psi)}{\Delta x}, \]](/trac/tracmath/53e7c3bc53ddcd1dc1c1c9a6e1e640e3eb1c5892.png)
where
![\[F_{i\pm\frac{1}{2}} \]](/trac/tracmath/59709fd6caa67dfd553c086346471bd94accb061.png)
denotes the fluxes staggered a half grid length related to the advected quantity.
Wicker and Skamarock (2002) dicretized the 6th and 5th order fluxes as follows:
TracMath macro processor has detected an error. Please fix the problem before continuing.
The command:
'/usr/bin/pdflatex -interaction=nonstopmode 1ebc1fd2eb680917f6c99414e3b9eadd8cd36d48.tex'
failed with the following output:
"This is pdfTeX, Version 3.1415926-2.5-1.40.14 (TeX Live 2013)\n restricted \\write18 enabled.\nentering extended mode\n(./1ebc1fd2eb680917f6c99414e3b9eadd8cd36d48.tex\nLaTeX2e <2011/06/27>\nBabel <3.9h> and hyphenation patterns for 78 languages loaded.\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/article.cls\nDocument Class: article 2007/10/19 v1.4h Standard LaTeX document class\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/size10.clo))\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/inputenc.sty\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/utf8.def\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/t1enc.dfu)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/ot1enc.dfu)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/omsenc.dfu)))\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/cmap/cmap.sty)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/cm-super/type1ec.sty\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/t1cmr.fd))\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/fontenc.sty\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/t1enc.def)<<t1.cmap\n>>) (/localdata/software/texlive/2013/texmf-dist/tex/latex/amsmath/amsmath.sty\nFor additional information on amsmath, use the `?' option.\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsmath/amstext.sty\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsmath/amsgen.sty))\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsmath/amsbsy.sty)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsmath/amsopn.sty))\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amscls/amsthm.sty)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsfonts/amssymb.sty\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsfonts/amsfonts.sty))\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/tools/bm.sty)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/preview/preview.sty\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/preview/prtightpage.def)\n)\nNo file 1ebc1fd2eb680917f6c99414e3b9eadd8cd36d48.aux.\nPreview: Fontsize 10pt\nPreview: PDFoutput 1\n<<ot1.cmap>><<oml.cmap>><<oms.cmap>><<omx.cmap>>\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsfonts/umsa.fd)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsfonts/umsb.fd)\n! Misplaced alignment tab character &.\nl.14 \\[ F_{i-\\frac{1}{2}}^{6} &\n =& \\frac{u_{i-\\frac{1}{2}}}{60} \\left[ 37(\\ps...\n\n! Misplaced alignment tab character &.\nl.14 \\[ F_{i-\\frac{1}{2}}^{6} &=&\n \\frac{u_{i-\\frac{1}{2}}}{60} \\left[ 37(\\ps...\n\nPreview: Tightpage -32891 -32891 32891 32891\n[1{/localdata/software/texlive/2013/texmf-var/fonts/map/pdftex/updmap/pdftex.ma\np}] (./1ebc1fd2eb680917f6c99414e3b9eadd8cd36d48.aux) )\n(see the transcript file for additional information)</localdata/software/texliv\ne/2013/texmf-dist/fonts/type1/public/amsfonts/cm/cmmi10.pfb></localdata/softwar\ne/texlive/2013/texmf-dist/fonts/type1/public/amsfonts/cm/cmmi7.pfb></localdata/\nsoftware/texlive/2013/texmf-dist/fonts/type1/public/amsfonts/cm/cmr10.pfb></loc\naldata/software/texlive/2013/texmf-dist/fonts/type1/public/amsfonts/cm/cmr5.pfb\n></localdata/software/texlive/2013/texmf-dist/fonts/type1/public/amsfonts/cm/cm\nr7.pfb></localdata/software/texlive/2013/texmf-dist/fonts/type1/public/amsfonts\n/cm/cmsy10.pfb></localdata/software/texlive/2013/texmf-dist/fonts/type1/public/\namsfonts/cm/cmsy7.pfb>\nOutput written on 1ebc1fd2eb680917f6c99414e3b9eadd8cd36d48.pdf (1 page, 59949 b\nytes).\nTranscript written on 1ebc1fd2eb680917f6c99414e3b9eadd8cd36d48.log.\n" TracMath macro processor has detected an error. Please fix the problem before continuing.
The command:
'/usr/bin/pdflatex -interaction=nonstopmode daad285d6ef38475045284963456f3c0bece2f45.tex'
failed with the following output:
"This is pdfTeX, Version 3.1415926-2.5-1.40.14 (TeX Live 2013)\n restricted \\write18 enabled.\nentering extended mode\n(./daad285d6ef38475045284963456f3c0bece2f45.tex\nLaTeX2e <2011/06/27>\nBabel <3.9h> and hyphenation patterns for 78 languages loaded.\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/article.cls\nDocument Class: article 2007/10/19 v1.4h Standard LaTeX document class\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/size10.clo))\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/inputenc.sty\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/utf8.def\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/t1enc.dfu)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/ot1enc.dfu)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/omsenc.dfu)))\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/cmap/cmap.sty)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/cm-super/type1ec.sty\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/t1cmr.fd))\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/fontenc.sty\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/t1enc.def)<<t1.cmap\n>>) (/localdata/software/texlive/2013/texmf-dist/tex/latex/amsmath/amsmath.sty\nFor additional information on amsmath, use the `?' option.\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsmath/amstext.sty\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsmath/amsgen.sty))\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsmath/amsbsy.sty)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsmath/amsopn.sty))\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amscls/amsthm.sty)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsfonts/amssymb.sty\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsfonts/amsfonts.sty))\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/tools/bm.sty)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/preview/preview.sty\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/preview/prtightpage.def)\n)\nNo file daad285d6ef38475045284963456f3c0bece2f45.aux.\nPreview: Fontsize 10pt\nPreview: PDFoutput 1\n<<ot1.cmap>><<oml.cmap>><<oms.cmap>><<omx.cmap>>\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsfonts/umsa.fd)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsfonts/umsb.fd)\n! Misplaced alignment tab character &.\nl.14 \\[ F_{i-\\frac{1}{2}}^{5} &\n =& F_{i-\\frac{1}{2}}^{6} - \\frac{|u_{i-\\frac{...\n\n! Misplaced alignment tab character &.\nl.14 \\[ F_{i-\\frac{1}{2}}^{5} &=&\n F_{i-\\frac{1}{2}}^{6} - \\frac{|u_{i-\\frac{...\n\nPreview: Tightpage -32891 -32891 32891 32891\n[1{/localdata/software/texlive/2013/texmf-var/fonts/map/pdftex/updmap/pdftex.ma\np}] (./daad285d6ef38475045284963456f3c0bece2f45.aux) )\n(see the transcript file for additional information)</localdata/software/texliv\ne/2013/texmf-dist/fonts/type1/public/amsfonts/cm/cmmi10.pfb></localdata/softwar\ne/texlive/2013/texmf-dist/fonts/type1/public/amsfonts/cm/cmmi7.pfb></localdata/\nsoftware/texlive/2013/texmf-dist/fonts/type1/public/amsfonts/cm/cmr10.pfb></loc\naldata/software/texlive/2013/texmf-dist/fonts/type1/public/amsfonts/cm/cmr5.pfb\n></localdata/software/texlive/2013/texmf-dist/fonts/type1/public/amsfonts/cm/cm\nr7.pfb></localdata/software/texlive/2013/texmf-dist/fonts/type1/public/amsfonts\n/cm/cmsy10.pfb></localdata/software/texlive/2013/texmf-dist/fonts/type1/public/\namsfonts/cm/cmsy7.pfb>\nOutput written on daad285d6ef38475045284963456f3c0bece2f45.pdf (1 page, 60101 b\nytes).\nTranscript written on daad285d6ef38475045284963456f3c0bece2f45.log.\n"
The 5th order upwind discretization (WS5) consists of the centered non dissipative 6th (WS6) order flux and an artificial added numerical dissipation term. This term is necessary to stabilize the numerical solution, because higher order centered fluxes exhibits worse stability properties.The absolute value assures a dissipative effect also for u < 0.
Numerical properties
A semidiscrete fourier transformation for the spatial derivatives maps the one dimensional advection equation in fourier space as follows (Baldauf, 2008):
![\[ \frac{\partial \hat{\psi}_{\kappa}}{\partial t} = - \frac{i}{\Delta t} C_{r} \kappa_{eff}\,\hat{\psi}_{\kappa}, \]](/trac/tracmath/ec07d3fd04956bdcbd3d16f239aee2e41b3911f9.png)
where
TracMath macro processor has detected an error. Please fix the problem before continuing.
The command:
'/usr/bin/pdflatex -interaction=nonstopmode 894de7f9c6fe4938fb0f6054fe301352c88e1442.tex'
failed with the following output:
"This is pdfTeX, Version 3.1415926-2.5-1.40.14 (TeX Live 2013)\n restricted \\write18 enabled.\nentering extended mode\n(./894de7f9c6fe4938fb0f6054fe301352c88e1442.tex\nLaTeX2e <2011/06/27>\nBabel <3.9h> and hyphenation patterns for 78 languages loaded.\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/article.cls\nDocument Class: article 2007/10/19 v1.4h Standard LaTeX document class\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/size10.clo))\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/inputenc.sty\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/utf8.def\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/t1enc.dfu)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/ot1enc.dfu)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/omsenc.dfu)))\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/cmap/cmap.sty)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/cm-super/type1ec.sty\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/t1cmr.fd))\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/fontenc.sty\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/base/t1enc.def)<<t1.cmap\n>>) (/localdata/software/texlive/2013/texmf-dist/tex/latex/amsmath/amsmath.sty\nFor additional information on amsmath, use the `?' option.\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsmath/amstext.sty\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsmath/amsgen.sty))\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsmath/amsbsy.sty)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsmath/amsopn.sty))\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amscls/amsthm.sty)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsfonts/amssymb.sty\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsfonts/amsfonts.sty))\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/tools/bm.sty)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/preview/preview.sty\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/preview/prtightpage.def)\n)\nNo file 894de7f9c6fe4938fb0f6054fe301352c88e1442.aux.\nPreview: Fontsize 10pt\nPreview: PDFoutput 1\n<<ot1.cmap>><<oml.cmap>><<oms.cmap>><<omx.cmap>>\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsfonts/umsa.fd)\n(/localdata/software/texlive/2013/texmf-dist/tex/latex/amsfonts/umsb.fd)\n! Missing $ inserted.\n<inserted text> \n $\nl.16 \\end{preview}\n \n! Display math should end with $$.\n<to be read again> \n \\endgroup \nl.16 \\end{preview}\n \nPreview: Tightpage -32891 -32891 32891 32891\n[1{/localdata/software/texlive/2013/texmf-var/fonts/map/pdftex/updmap/pdftex.ma\np}]\n\n! LaTeX Error: \\begin{preview} on input line 13 ended by \\end{document}.\n\nSee the LaTeX manual or LaTeX Companion for explanation.\nType H <return> for immediate help.\n ... \n \nl.17 \\end{document}\n \n(./894de7f9c6fe4938fb0f6054fe301352c88e1442.aux) )\n(\\end occurred inside a group at level 1)\n\n### semi simple group (level 1) entered at line 13 (\\begingroup)\n### bottom level\n(see the transcript file for additional information)</localdata/software/texliv\ne/2013/texmf-dist/fonts/type1/public/amsfonts/cm/cmmi10.pfb></localdata/softwar\ne/texlive/2013/texmf-dist/fonts/type1/public/amsfonts/cm/cmr10.pfb>\nOutput written on 894de7f9c6fe4938fb0f6054fe301352c88e1442.pdf (1 page, 20063 b\nytes).\nTranscript written on 894de7f9c6fe4938fb0f6054fe301352c88e1442.log.\n"
denotes the fourier transformed of
![\[
\psi. \quad C_{r} = \frac{u \Delta t}{\Delta x}
\]](/trac/tracmath/02c57de02e3b0f92e728c76137f216f5c5539f4a.png)
the Courant number which characterizes stability properties and i is a complex number.
![\[ \kappa_{eff} \]](/trac/tracmath/dfaab09dc4e9502b473e25f29befb06ffe8a4337.png)
is the effective wavenumber of a mode in fourier space which characterizes the modified wavenumber through the discretization. The real part of the effective wavenumber describes the dispersion error, the imaginary part the dissipation error.
Fig. 1 shows the dispersion and dissipation errors of WS3, WS4, WS5, WS6 and the 2nd order scheme of Piascek and Williams (1970) (PW) presented as a function of the dimensionless wavenumber
![\[ \kappa\Delta x. \]](/trac/tracmath/3176d2c5f49bf69ac9813da2bffeec12f1861f03.png)
The dispersion error of the upwind schemes and the dispersion error of the next higher even ordered scheme are identical. Generally the dispersion error decreases with increasing order of the dicretization. However, no scheme is able to resolute wavelength with the range near 2-Delta x.
The centered even ordered schemes holds no dissipation errors. The numerical dissipation is more local with increasing order, so the maximal affected wavelength by the dissipation term is round about 8 gridlength with WS5, whereat with WS3 16 gridlength waves are still affected. Accordingly to a maximum of the amplification factor at
![\[ \kappa\Delta x = 1.69 \]](/trac/tracmath/ee4efc9a95945c1eae05036801f1a81786d996b3.png)
in conjunction with the used Runge-Kutta method (Baldauf, 2008), the 5th order dissipation is more than sufficient to avoid instabilities.
The maximal stable Courant-number is Cr = 1.4 (Wicker and Skamarock, 2002), (Baldauf, 2008).
Note: A stable numerical solution can only guaranteed with the 3 rd oder Runge-Kutta method.
Boundaries
Due to the large stencil of WS5 additional ghost layers are necessary on each lateral boundary of each processor subdomain to avoid local data dependencies. Therefor the exchange of ghost layers is adapted to a dynamic number of ghost layers.
For the bottom and top boundaries a successive degradation from WS5 to WS3 to a 2nd order scheme is required to avoid unphysical fluxes which would arise from the bottom and top ghost layers. The used 2nd order scheme is based on a flux discretization to ensure consistency with the WS-schemes. So the skew symmetric PW-scheme cannot be used. Furthermore an additional numerical dissipation term of 2nd order, based on (Shchepetkin and McWilliams?, 1998),is required for the 2nd order scheme to perform a numerical stable switching of advection schemes of different order.
These successive degradation is also done for the lateral radiation boundary condition at the outflow and near topography (Note: Topography is not implemented at the moment.)
Statical evaluation of turbulent fluxes
The evaluation of turbulent fluxes should be consistent to the discretization in the prognostic equations, else some unphysical effects will occur. For example, the computation of the fluxes as variances and covariances will induce some conspicuous kinks in vertical heat and momentum fluxes near the surface, while the temperature and velocity profiles show no conspicuity. For computing turbulent fluxes as appearing in the prognostic equations the computed fluxes in the advection routines are buffered and reused also for the statistics. For getting only the turbulent, not the mean signal and to remove the influence of Galilei transformation, the centered flux Fi+1/2 has to be multiplied with a factor
![\[ \frac{u_{i+\frac{1}{2}} - 2 \overline u}{u_{i+\frac{1}{2}} - u_{i, Galilei}} \]](/trac/tracmath/1df3fb54042fd16698dd4d9ce32eccd53c9e5f0f.png)
and the dissipative flux with a factor
![\[ \frac{|u_{i+\frac{1}{2}} - 2 \overline u|}{|u_{i+\frac{1}{2}} - u_{i, Galilei}|}, \]](/trac/tracmath/7a1569ed9bba32d310e74dc0bca866e8e3ebbd37.png)
where u denotes the respective velocity component. Furthermore the turbulent fluxes are evaluated on each Runge-Kutta substep and weighted with the respective Runge-Kutta coefficients, to remove dependencies of the Runge-Kutta substep. The interpretation of the turbulent fluxes as variances and covariances is no longer valid when using WS5. For other advection schemes like the PW-scheme the interpretation of turbulent fluxes as co/variances is still valid, because the discretization is alike the computation of the co/variances.
