Skip to content
Snippets Groups Projects
Commit 248d1a6e authored by Victor Stepanenko's avatar Victor Stepanenko
Browse files

TKE budget terms output at a selected depth added

parent 985f6a82
Branches
Tags
No related merge requests found
Show whitespace changes
Inline Side-by-side
Showing
with 323 additions and 272 deletions
docs/users_guide.aux
+ 7
− 7
View file @ 248d1a6e
\relax
\@writefile{toc}{\contentsline {section}{\numberline {0.1}Compiling the model and creating a project}{2}}
\@writefile{toc}{\contentsline {section}{\numberline {0.2}Configuring the driver file}{2}}
\@writefile{lot}{\contentsline {table}{\numberline {1}{\ignorespaces The parameters of driver file\relax }}{2}}
\@writefile{toc}{\contentsline {section}{\numberline {0.1}Compiling the model and creating a project}{2}\protected@file@percent }
\@writefile{toc}{\contentsline {section}{\numberline {0.2}Configuring the driver file}{2}\protected@file@percent }
\@writefile{lot}{\contentsline {table}{\numberline {1}{\ignorespaces The parameters of driver file\relax }}{2}\protected@file@percent }
\gdef \LT@i {\LT@entry
{1}{201.44666pt}\LT@entry
{1}{296.92755pt}}
\@writefile{toc}{\contentsline {section}{\numberline {0.3}Configuring the setup file}{5}}
\@writefile{lot}{\contentsline {table}{\numberline {3}{\ignorespaces The parameters of setup file\relax }}{5}}
\@writefile{toc}{\contentsline {section}{\numberline {0.3}Configuring the setup file}{5}\protected@file@percent }
\@writefile{lot}{\contentsline {table}{\numberline {3}{\ignorespaces The parameters of setup file\relax }}{5}\protected@file@percent }
\gdef \LT@ii {\LT@entry
{1}{201.44666pt}\LT@entry
{1}{296.92755pt}}
\@writefile{toc}{\contentsline {section}{\numberline {0.4}Running the model}{8}}
\@writefile{toc}{\contentsline {section}{\numberline {0.5}Model output}{9}}
\@writefile{toc}{\contentsline {section}{\numberline {0.4}Running the model}{8}\protected@file@percent }
\@writefile{toc}{\contentsline {section}{\numberline {0.5}Model output}{9}\protected@file@percent }
docs/users_guide.log
+ 131
− 126
View file @ 248d1a6e
This is pdfTeX, Version 3.14159265-2.6-1.40.16 (TeX Live 2015/Debian) (preloaded format=pdflatex 2019.1.26) 5 APR 2020 14:22
This is pdfTeX, Version 3.14159265-2.6-1.40.20 (TeX Live 2019) (preloaded format=pdflatex 2019.11.28) 18 AUG 2020 17:00
entering extended mode
restricted \write18 enabled.
file:line:error style messages enabled.
%&-line parsing enabled.
**users_guide.tex
(./users_guide.tex
LaTeX2e <2016/02/01>
Babel <3.9q> and hyphenation patterns for 81 language(s) loaded.
(/usr/share/texlive/texmf-dist/tex/latex/base/book.cls
Document Class: book 2014/09/29 v1.4h Standard LaTeX document class
(/usr/share/texlive/texmf-dist/tex/latex/base/bk10.clo
File: bk10.clo 2014/09/29 v1.4h Standard LaTeX file (size option)
LaTeX2e <2018-12-01>
(/usr/local/texlive/2019/texmf-dist/tex/latex/base/book.cls
Document Class: book 2018/09/03 v1.4i Standard LaTeX document class
(/usr/local/texlive/2019/texmf-dist/tex/latex/base/bk10.clo
File: bk10.clo 2018/09/03 v1.4i Standard LaTeX file (size option)
)
\c@part=\count79
\c@chapter=\count80
\c@section=\count81
\c@subsection=\count82
\c@subsubsection=\count83
\c@paragraph=\count84
\c@subparagraph=\count85
\c@figure=\count86
\c@table=\count87
\c@part=\count80
\c@chapter=\count81
\c@section=\count82
\c@subsection=\count83
\c@subsubsection=\count84
\c@paragraph=\count85
\c@subparagraph=\count86
\c@figure=\count87
\c@table=\count88
\abovecaptionskip=\skip41
\belowcaptionskip=\skip42
\bibindent=\dimen102
)
(/usr/share/texlive/texmf-dist/tex/latex/amsmath/amsmath.sty
Package: amsmath 2016/03/03 v2.15a AMS math features
(/usr/local/texlive/2019/texmf-dist/tex/latex/amsmath/amsmath.sty
Package: amsmath 2018/12/01 v2.17b AMS math features
\@mathmargin=\skip43
For additional information on amsmath, use the `?' option.
(/usr/share/texlive/texmf-dist/tex/latex/amsmath/amstext.sty
(/usr/local/texlive/2019/texmf-dist/tex/latex/amsmath/amstext.sty
Package: amstext 2000/06/29 v2.01 AMS text
(/usr/share/texlive/texmf-dist/tex/latex/amsmath/amsgen.sty
(/usr/local/texlive/2019/texmf-dist/tex/latex/amsmath/amsgen.sty
File: amsgen.sty 1999/11/30 v2.0 generic functions
\@emptytoks=\toks14
\ex@=\dimen103
))
(/usr/share/texlive/texmf-dist/tex/latex/amsmath/amsbsy.sty
(/usr/local/texlive/2019/texmf-dist/tex/latex/amsmath/amsbsy.sty
Package: amsbsy 1999/11/29 v1.2d Bold Symbols
\pmbraise@=\dimen104
)
(/usr/share/texlive/texmf-dist/tex/latex/amsmath/amsopn.sty
Package: amsopn 1999/12/14 v2.01 operator names
(/usr/local/texlive/2019/texmf-dist/tex/latex/amsmath/amsopn.sty
Package: amsopn 2016/03/08 v2.02 operator names
)
\inf@bad=\count88
LaTeX Info: Redefining \frac on input line 199.
\uproot@=\count89
\leftroot@=\count90
LaTeX Info: Redefining \overline on input line 297.
\classnum@=\count91
\DOTSCASE@=\count92
LaTeX Info: Redefining \ldots on input line 394.
LaTeX Info: Redefining \dots on input line 397.
LaTeX Info: Redefining \cdots on input line 518.
\Mathstrutbox@=\box26
\strutbox@=\box27
\inf@bad=\count89
LaTeX Info: Redefining \frac on input line 223.
\uproot@=\count90
\leftroot@=\count91
LaTeX Info: Redefining \overline on input line 385.
\classnum@=\count92
\DOTSCASE@=\count93
LaTeX Info: Redefining \ldots on input line 482.
LaTeX Info: Redefining \dots on input line 485.
LaTeX Info: Redefining \cdots on input line 606.
\Mathstrutbox@=\box27
\strutbox@=\box28
\big@size=\dimen105
LaTeX Font Info: Redeclaring font encoding OML on input line 630.
LaTeX Font Info: Redeclaring font encoding OMS on input line 631.
\macc@depth=\count93
\c@MaxMatrixCols=\count94
LaTeX Font Info: Redeclaring font encoding OML on input line 729.
LaTeX Font Info: Redeclaring font encoding OMS on input line 730.
\macc@depth=\count94
\c@MaxMatrixCols=\count95
\dotsspace@=\muskip10
\c@parentequation=\count95
\dspbrk@lvl=\count96
\c@parentequation=\count96
\dspbrk@lvl=\count97
\tag@help=\toks15
\row@=\count97
\column@=\count98
\maxfields@=\count99
\row@=\count98
\column@=\count99
\maxfields@=\count100
\andhelp@=\toks16
\eqnshift@=\dimen106
\alignsep@=\dimen107
......@@ -79,63 +79,62 @@ LaTeX Font Info: Redeclaring font encoding OMS on input line 631.
\multlinegap=\skip44
\multlinetaggap=\skip45
\mathdisplay@stack=\toks18
LaTeX Info: Redefining \[ on input line 2735.
LaTeX Info: Redefining \] on input line 2736.
LaTeX Info: Redefining \[ on input line 2844.
LaTeX Info: Redefining \] on input line 2845.
)
(/usr/share/texlive/texmf-dist/tex/latex/preprint/authblk.sty
(/usr/local/texlive/2019/texmf-dist/tex/latex/preprint/authblk.sty
Package: authblk 2001/02/27 1.3 (PWD)
\affilsep=\skip46
\@affilsep=\skip47
\c@Maxaffil=\count100
\c@authors=\count101
\c@affil=\count102
\c@Maxaffil=\count101
\c@authors=\count102
\c@affil=\count103
)
(/usr/share/texlive/texmf-dist/tex/latex/listings/listings.sty
(/usr/share/texlive/texmf-dist/tex/latex/graphics/keyval.sty
(/usr/local/texlive/2019/texmf-dist/tex/latex/listings/listings.sty
(/usr/local/texlive/2019/texmf-dist/tex/latex/graphics/keyval.sty
Package: keyval 2014/10/28 v1.15 key=value parser (DPC)
\KV@toks@=\toks19
)
\lst@mode=\count103
\lst@gtempboxa=\box28
\lst@mode=\count104
\lst@gtempboxa=\box29
\lst@token=\toks20
\lst@length=\count104
\lst@length=\count105
\lst@currlwidth=\dimen112
\lst@column=\count105
\lst@pos=\count106
\lst@column=\count106
\lst@pos=\count107
\lst@lostspace=\dimen113
\lst@width=\dimen114
\lst@newlines=\count107
\lst@lineno=\count108
\lst@newlines=\count108
\lst@lineno=\count109
\lst@maxwidth=\dimen115
(/usr/share/texlive/texmf-dist/tex/latex/listings/lstmisc.sty
File: lstmisc.sty 2015/06/04 1.6 (Carsten Heinz)
\c@lstnumber=\count109
\lst@skipnumbers=\count110
\lst@framebox=\box29
(/usr/local/texlive/2019/texmf-dist/tex/latex/listings/lstmisc.sty
File: lstmisc.sty 2019/02/27 1.8b (Carsten Heinz)
\c@lstnumber=\count110
\lst@skipnumbers=\count111
\lst@framebox=\box30
)
(/usr/share/texlive/texmf-dist/tex/latex/listings/listings.cfg
File: listings.cfg 2015/06/04 1.6 listings configuration
(/usr/local/texlive/2019/texmf-dist/tex/latex/listings/listings.cfg
File: listings.cfg 2019/02/27 1.8b listings configuration
))
Package: listings 2015/06/04 1.6 (Carsten Heinz)
Package: listings 2019/02/27 1.8b (Carsten Heinz)
(/usr/share/texlive/texmf-dist/tex/latex/geometry/geometry.sty
Package: geometry 2010/09/12 v5.6 Page Geometry
(/usr/local/texlive/2019/texmf-dist/tex/latex/geometry/geometry.sty
Package: geometry 2018/04/16 v5.8 Page Geometry
(/usr/share/texlive/texmf-dist/tex/generic/oberdiek/ifpdf.sty
Package: ifpdf 2011/01/30 v2.3 Provides the ifpdf switch (HO)
Package ifpdf Info: pdfTeX in PDF mode is detected.
(/usr/local/texlive/2019/texmf-dist/tex/generic/oberdiek/ifpdf.sty
Package: ifpdf 2018/09/07 v3.3 Provides the ifpdf switch
)
(/usr/share/texlive/texmf-dist/tex/generic/oberdiek/ifvtex.sty
Package: ifvtex 2010/03/01 v1.5 Detect VTeX and its facilities (HO)
(/usr/local/texlive/2019/texmf-dist/tex/generic/oberdiek/ifvtex.sty
Package: ifvtex 2016/05/16 v1.6 Detect VTeX and its facilities (HO)
Package ifvtex Info: VTeX not detected.
)
(/usr/share/texlive/texmf-dist/tex/generic/ifxetex/ifxetex.sty
(/usr/local/texlive/2019/texmf-dist/tex/generic/ifxetex/ifxetex.sty
Package: ifxetex 2010/09/12 v0.6 Provides ifxetex conditional
)
\Gm@cnth=\count111
\Gm@cntv=\count112
\c@Gm@tempcnt=\count113
\Gm@cnth=\count112
\Gm@cntv=\count113
\c@Gm@tempcnt=\count114
\Gm@bindingoffset=\dimen116
\Gm@wd@mp=\dimen117
\Gm@odd@mp=\dimen118
......@@ -146,41 +145,43 @@ Package: ifxetex 2010/09/12 v0.6 Provides ifxetex conditional
\Gm@layoutvoffset=\dimen123
\Gm@dimlist=\toks21
)
(/usr/share/texlive/texmf-dist/tex/latex/tools/array.sty
Package: array 2014/10/28 v2.4c Tabular extension package (FMi)
(/usr/local/texlive/2019/texmf-dist/tex/latex/tools/array.sty
Package: array 2018/12/30 v2.4k Tabular extension package (FMi)
\col@sep=\dimen124
\ar@mcellbox=\box31
\extrarowheight=\dimen125
\NC@list=\toks22
\extratabsurround=\skip48
\backup@length=\skip49
\ar@cellbox=\box32
)
(/usr/share/texlive/texmf-dist/tex/latex/tools/longtable.sty
Package: longtable 2014/10/28 v4.11 Multi-page Table package (DPC)
(/usr/local/texlive/2019/texmf-dist/tex/latex/tools/longtable.sty
Package: longtable 2014/10/28 v4.11 Multi-page Table package (DPC)+ FMi change
\LTleft=\skip50
\LTright=\skip51
\LTpre=\skip52
\LTpost=\skip53
\LTchunksize=\count114
\LTchunksize=\count115
\LTcapwidth=\dimen126
\LT@head=\box30
\LT@firsthead=\box31
\LT@foot=\box32
\LT@lastfoot=\box33
\LT@cols=\count115
\LT@rows=\count116
\c@LT@tables=\count117
\c@LT@chunks=\count118
\LT@head=\box33
\LT@firsthead=\box34
\LT@foot=\box35
\LT@lastfoot=\box36
\LT@cols=\count116
\LT@rows=\count117
\c@LT@tables=\count118
\c@LT@chunks=\count119
\LT@p@ftn=\toks23
)
(/usr/share/texlive/texmf-dist/tex/latex/was/gensymb.sty
(/usr/local/texlive/2019/texmf-dist/tex/latex/was/gensymb.sty
Package: gensymb 2003/07/02 v1.0 (WaS)
)
(/usr/share/texlive/texmf-dist/tex/latex/caption/caption.sty
Package: caption 2016/02/21 v3.3-144 Customizing captions (AR)
(/usr/local/texlive/2019/texmf-dist/tex/latex/caption/caption.sty
Package: caption 2018/10/06 v3.3-154 Customizing captions (AR)
(/usr/share/texlive/texmf-dist/tex/latex/caption/caption3.sty
Package: caption3 2016/02/04 v1.7-139 caption3 kernel (AR)
Package caption3 Info: TeX engine: e-TeX on input line 67.
(/usr/local/texlive/2019/texmf-dist/tex/latex/caption/caption3.sty
Package: caption3 2018/09/12 v1.8c caption3 kernel (AR)
Package caption3 Info: TeX engine: e-TeX on input line 64.
\captionmargin=\dimen127
\captionmargin@=\dimen128
\captionwidth=\dimen129
......@@ -189,12 +190,13 @@ Package caption3 Info: TeX engine: e-TeX on input line 67.
\caption@parindent=\dimen132
\caption@hangindent=\dimen133
)
\c@ContinuedFloat=\count119
\c@caption@flags=\count120
\c@ContinuedFloat=\count121
Package caption Info: listings package is loaded.
Package caption Info: longtable package is loaded.
(/usr/share/texlive/texmf-dist/tex/latex/caption/ltcaption.sty
Package: ltcaption 2013/06/09 v1.4-94 longtable captions (AR)
(/usr/local/texlive/2019/texmf-dist/tex/latex/caption/ltcaption.sty
Package: ltcaption 2018/08/26 v1.4-95 longtable captions (AR)
))
Package geometry Warning: Over-specification in `h'-direction.
......@@ -219,7 +221,7 @@ LaTeX Font Info: Checking defaults for OMX/cmex/m/n on input line 46.
LaTeX Font Info: ... okay on input line 46.
LaTeX Font Info: Checking defaults for U/cmr/m/n on input line 46.
LaTeX Font Info: ... okay on input line 46.
\c@lstlisting=\count120
\c@lstlisting=\count122
*geometry* driver: auto-detecting
*geometry* detected driver: pdftex
......@@ -271,10 +273,10 @@ Package caption Info: End \AtBeginDocument code.
[1
{/var/lib/texmf/fonts/map/pdftex/updmap/pdftex.map}]
{/usr/local/texlive/2019/texmf-var/fonts/map/pdftex/updmap/pdftex.map}]
LaTeX Font Info: Try loading font information for OML+cmr on input line 66.
(/usr/share/texlive/texmf-dist/tex/latex/base/omlcmr.fd
(/usr/local/texlive/2019/texmf-dist/tex/latex/base/omlcmr.fd
File: omlcmr.fd 2014/09/29 v2.5h Standard LaTeX font definitions
)
LaTeX Font Info: Font shape `OML/cmr/m/n' in size <10> not available
......@@ -334,42 +336,45 @@ Overfull \hbox (14.67734pt too wide) in alignment at lines 372--416
[]
[6]
Missing character: There is no in font cmr10!
Missing character: There is no in font cmr10!
Overfull \hbox (14.67734pt too wide) in alignment at lines 416--461
[] []
[]
[7] [8]
LaTeX Font Info: Try loading font information for OMS+cmr on input line 488.
[7]
Overfull \hbox (14.67734pt too wide) in alignment at lines 461--463
[] []
[]
[8]
LaTeX Font Info: Try loading font information for OMS+cmr on input line 490.
(/usr/share/texlive/texmf-dist/tex/latex/base/omscmr.fd
(/usr/local/texlive/2019/texmf-dist/tex/latex/base/omscmr.fd
File: omscmr.fd 2014/09/29 v2.5h Standard LaTeX font definitions
)
LaTeX Font Info: Font shape `OMS/cmr/m/n' in size <10> not available
(Font) Font shape `OMS/cmsy/m/n' tried instead on input line 488.
(Font) Font shape `OMS/cmsy/m/n' tried instead on input line 490.
[9]
(./users_guide.aux) )
Here is how much of TeX's memory you used:
3924 strings out of 493029
57961 string characters out of 6136234
155919 words of memory out of 5000000
7463 multiletter control sequences out of 15000+600000
4393 strings out of 492609
62611 string characters out of 6129047
164304 words of memory out of 5000000
8317 multiletter control sequences out of 15000+600000
5995 words of font info for 22 fonts, out of 8000000 for 9000
1141 hyphenation exceptions out of 8191
36i,8n,39p,563b,1088s stack positions out of 5000i,500n,10000p,200000b,80000s
</usr/share/texlive/texmf-dist/fonts/type1/public/amsfonts
/cm/cmbx12.pfb></usr/share/texlive/texmf-dist/fonts/type1/public/amsfonts/cm/cm
mi10.pfb></usr/share/texlive/texmf-dist/fonts/type1/public/amsfonts/cm/cmr10.pf
b></usr/share/texlive/texmf-dist/fonts/type1/public/amsfonts/cm/cmr12.pfb></usr
/share/texlive/texmf-dist/fonts/type1/public/amsfonts/cm/cmr17.pfb></usr/share/
texlive/texmf-dist/fonts/type1/public/amsfonts/cm/cmr7.pfb></usr/share/texlive/
texmf-dist/fonts/type1/public/amsfonts/cm/cmsl10.pfb></usr/share/texlive/texmf-
dist/fonts/type1/public/amsfonts/cm/cmsy10.pfb></usr/share/texlive/texmf-dist/f
onts/type1/public/amsfonts/cm/cmsy7.pfb></usr/share/texlive/texmf-dist/fonts/ty
pe1/public/amsfonts/cm/cmti10.pfb>
Output written on users_guide.pdf (9 pages, 162683 bytes).
36i,8n,39p,563b,1214s stack positions out of 5000i,500n,10000p,200000b,80000s
</usr/local/texlive/2019/texmf-dist/fonts/type1/public/ams
fonts/cm/cmbx12.pfb></usr/local/texlive/2019/texmf-dist/fonts/type1/public/amsf
onts/cm/cmmi10.pfb></usr/local/texlive/2019/texmf-dist/fonts/type1/public/amsfo
nts/cm/cmr10.pfb></usr/local/texlive/2019/texmf-dist/fonts/type1/public/amsfont
s/cm/cmr12.pfb></usr/local/texlive/2019/texmf-dist/fonts/type1/public/amsfonts/
cm/cmr17.pfb></usr/local/texlive/2019/texmf-dist/fonts/type1/public/amsfonts/cm
/cmr7.pfb></usr/local/texlive/2019/texmf-dist/fonts/type1/public/amsfonts/cm/cm
sl10.pfb></usr/local/texlive/2019/texmf-dist/fonts/type1/public/amsfonts/cm/cms
y10.pfb></usr/local/texlive/2019/texmf-dist/fonts/type1/public/amsfonts/cm/cmsy
7.pfb></usr/local/texlive/2019/texmf-dist/fonts/type1/public/amsfonts/cm/cmti10
.pfb>
Output written on users_guide.pdf (9 pages, 162677 bytes).
PDF statistics:
74 PDF objects out of 1000 (max. 8388607)
52 compressed objects within 1 object stream
......
docs/users_guide.pdf
+ 0
− 0
View file @ 248d1a6e
No preview for this file type
docs/users_guide.synctex.gz
+ 0
− 0
View file @ 248d1a6e
No preview for this file type
docs/users_guide.tex
+ 3
− 1
View file @ 248d1a6e
......@@ -441,9 +441,11 @@ give significant speedup \\
\hline
\emph{time\_series} & Switch for time series output : 0 - off, 1 - on \\
\hline
\emph{zserout} & The depth of time series that are written in output (excluding time series at surface, bottom, bulk parameters, etc.), m \\
\hline
\emph{dt\_out} & Time interval for time series output, hours \\
\hline
\emph{nscreen} & Еhe period of screen output, time steps \\
\emph{nscreen} & The period of screen output, time steps \\
\hline
\emph{scale\_output} & the switch for scaling of output of turbulent characteristics: 0 - off, 1 - on \\
\hline
......
setup/river_setup.dat
+ 1
− 0
View file @ 248d1a6e
......@@ -304,6 +304,7 @@ daily 1
hourly 0
everystep 0
time_series 1
zserout 0.4
dt_out 0.16666667
nscreen 1000
scale_output 0
......
source/model/driving_params_mod.f90
+ 4
− 1
View file @ 248d1a6e
......@@ -25,7 +25,7 @@ type, public :: intpar
end type
integer(kind=iintegers), parameter :: numder = 42 ! Number derived-type model integer controls handled in unified manner
integer(kind=iintegers), parameter :: numder_real = 15 ! The same, but real controls
integer(kind=iintegers), parameter :: numder_real = 16 ! The same, but real controls
real(kind=ireals), parameter :: missing_value = -999.
......@@ -61,6 +61,7 @@ type(realpar), target :: soil_meth_prod
type(realpar), target :: c_d
type(realpar), target :: nManning
type(realpar), target :: horvisc
type(realpar), target :: zserout
real(kind=ireals) :: dttribupdate = missing_value
......@@ -227,6 +228,7 @@ if (firstcall) then
dertypepar_real(13)%p => c_d
dertypepar_real(14)%p => nManning
dertypepar_real(15)%p => horvisc
dertypepar_real(16)%p => zserout
dertypepar(1)%p => stabfunc
dertypepar(2)%p => tribheat
......@@ -287,6 +289,7 @@ if (firstcall) then
c_d%name = 'c_d'
nManning%name = 'nManning'
horvisc%name = 'horvisc'
zserout%name = 'zserout'
stabfunc%name = 'stabfunc' ;
......
source/model/lake_modules.f90
+ 1
− 0
View file @ 248d1a6e
......@@ -583,6 +583,7 @@ real(kind=ireals) :: Seps_integr_positive, Seps_integr_negative, &
& Geneps_integr, epseps_integr
real(kind=ireals) :: TKE_balance, eps_balance
real(kind=ireals) :: Eseiches
real(kind=ireals) :: TKE_budget_terms(1:5)
! Data structure for turbulent characteristics
type, public :: turb_type
......
source/model/out_mod.f90
+ 160
− 134
View file @ 248d1a6e
......@@ -781,7 +781,7 @@ use NUMERIC_PARAMS , only : &
use ARRAYS_TURB, only : row, H_mixed_layer, H_entrainment, &
& signwaveheight, w_conv_scale, u_star, Ri_bulk, maxN, i_maxN, &
& ThermThick, ReTherm, RiTherm
& ThermThick, ReTherm, RiTherm, TKE_budget_terms
use ARRAYS_WATERSTATE, only : Tw1, Ti1, lamw, salice
......@@ -878,7 +878,7 @@ real(kind=ireals), allocatable :: var_scalar (:,:,:)
real(kind=ireals), allocatable :: accum_var_scalar(:,:,:)
! Integers
integer(kind=iintegers), parameter :: nfiles = 5
integer(kind=iintegers), parameter :: nfiles = 6
integer(kind=iintegers), parameter :: numaccum = 5 ! The number of scalar variables being accumulated (averaged)
integer(kind=iintegers), allocatable :: n_unit(:,:,:)
integer(kind=iintegers), allocatable :: count_out(:,:)
......@@ -919,161 +919,182 @@ endif
if (firstcall(ix,iy)) then
write (coords_point, '(2i3)') ix, iy
i = 1
call CHECK_UNIT(lake_series_out_unit_min,lake_series_out_unit_max, &
& n_unit(1,ix,iy))
open (n_unit(1,ix,iy),file=outpath(1:len_trim(outpath))//'time_series/'// &
& n_unit(i,ix,iy))
open (n_unit(i,ix,iy),file=outpath(1:len_trim(outpath))//'time_series/'// &
& 'layers'//coords_point//'.dat', status='unknown')
write (n_unit(1,ix,iy),*)'Col. 1 - year'
write (n_unit(1,ix,iy),*)'Col. 2 - month'
write (n_unit(1,ix,iy),*)'Col. 3 - day'
write (n_unit(1,ix,iy),*)'Col. 4 - hour'
write (n_unit(1,ix,iy),*)'Col. 5 - the time from the start of integration, hours'
write (n_unit(1,ix,iy),*)'Col. 6 - water layer thickness, m'
write (n_unit(1,ix,iy),*)'Col. 7 - W mixed layer thickness, m'
write (n_unit(1,ix,iy),*)'Col. 8 - E mixed layer thickness, m'
write (n_unit(1,ix,iy),*)'Col. 9 - S mixed layer thickness, m'
write (n_unit(1,ix,iy),*)'Col. 10 - N mixed layer thickness, m'
write (n_unit(1,ix,iy),*)'Col. 11 - W lower layer thickness, m'
write (n_unit(1,ix,iy),*)'Col. 12 - E lower layer thickness, m'
write (n_unit(1,ix,iy),*)'Col. 13 - S lower layer thickness, m'
write (n_unit(1,ix,iy),*)'Col. 14 - N lower layer thickness, m'
write (n_unit(1,ix,iy),*)'Col. 15 - ice layer thickness, m'
write (n_unit(1,ix,iy),*)'Col. 16 - snow layer thickness, m'
write (n_unit(1,ix,iy),*)'Col. 17 - bottom ice thickness, m'
write (n_unit(1,ix,iy),*)'Col. 18 - reservoir volume, m**3'
write (n_unit(1,ix,iy),*)'Col. 19 - volume deficit (accumulated), m**3'
write (n_unit(i,ix,iy),*)'Col. 1 - year'
write (n_unit(i,ix,iy),*)'Col. 2 - month'
write (n_unit(i,ix,iy),*)'Col. 3 - day'
write (n_unit(i,ix,iy),*)'Col. 4 - hour'
write (n_unit(i,ix,iy),*)'Col. 5 - the time from the start of integration, hours'
write (n_unit(i,ix,iy),*)'Col. 6 - water layer thickness, m'
write (n_unit(i,ix,iy),*)'Col. 7 - W mixed layer thickness, m'
write (n_unit(i,ix,iy),*)'Col. 8 - E mixed layer thickness, m'
write (n_unit(i,ix,iy),*)'Col. 9 - S mixed layer thickness, m'
write (n_unit(i,ix,iy),*)'Col. 10 - N mixed layer thickness, m'
write (n_unit(i,ix,iy),*)'Col. 11 - W lower layer thickness, m'
write (n_unit(i,ix,iy),*)'Col. 12 - E lower layer thickness, m'
write (n_unit(i,ix,iy),*)'Col. 13 - S lower layer thickness, m'
write (n_unit(i,ix,iy),*)'Col. 14 - N lower layer thickness, m'
write (n_unit(i,ix,iy),*)'Col. 15 - ice layer thickness, m'
write (n_unit(i,ix,iy),*)'Col. 16 - snow layer thickness, m'
write (n_unit(i,ix,iy),*)'Col. 17 - bottom ice thickness, m'
write (n_unit(i,ix,iy),*)'Col. 18 - reservoir volume, m**3'
write (n_unit(i,ix,iy),*)'Col. 19 - volume deficit (accumulated), m**3'
i = i + 1
call CHECK_UNIT(lake_series_out_unit_min,lake_series_out_unit_max, &
& n_unit(2,ix,iy))
open (n_unit(2,ix,iy),file=outpath(1:len_trim(outpath))//'time_series/'// &
& n_unit(i,ix,iy))
open (n_unit(i,ix,iy),file=outpath(1:len_trim(outpath))//'time_series/'// &
& 'T_fluxes'//coords_point//'.dat',status='unknown')
write (n_unit(2,ix,iy),*)'Col. 1 - year'
write (n_unit(2,ix,iy),*)'Col. 2 - month'
write (n_unit(2,ix,iy),*)'Col. 3 - day'
write (n_unit(2,ix,iy),*)'Col. 4 - hour'
write (n_unit(2,ix,iy),*)'Col. 5 - the time from the start of integration, hours'
write (n_unit(2,ix,iy),*)'Col. 6 - surface temperature, C'
write (n_unit(2,ix,iy),*)'Col. 7 - water skin temperature, C'
write (n_unit(2,ix,iy),*)'Col. 8 - water surface temperature, C'
write (n_unit(2,ix,iy),*)'Col. 9 - mean temperature of water coloumn, C'
write (n_unit(2,ix,iy),*)'Col. 10 - maximal temperature in the water coloumn, C'
write (n_unit(2,ix,iy),*)'Col. 11 - zero-dimensional model temperature, C'
write (n_unit(2,ix,iy),*)'Col. 12 - upper ice surface temperature, C'
write (n_unit(2,ix,iy),*)'Col. 13 - upper snow surface temperature, C'
write (n_unit(2,ix,iy),*)'Col. 14 - sensible heat flux, W/m**2'
write (n_unit(2,ix,iy),*)'Col. 15 - latent heat flux, W/m**2'
write (n_unit(2,ix,iy),*)'Col. 16 - downward heat flux at the upper lake surface, W/m**2'
write (n_unit(2,ix,iy),*)'Col. 17 - downward heat flux at the lake bottom, W/m**2'
write (n_unit(2,ix,iy),*)'Col. 18 - friction velocity in the surface air layer, m/s'
write (n_unit(2,ix,iy),*)'Col. 19 - friction velocity at bottom, m/s'
write (n_unit(2,ix,iy),*)'Col. 20 - wind work at the water surface, W/m**2'
write (n_unit(2,ix,iy),*)'Col. 21 - albedo of the lake-atmosphere interface, n/d'
write (n_unit(2,ix,iy),*)'Col. 22 - shortwave radiation penetrated below surface, W/m**2'
write (n_unit(2,ix,iy),*)'Col. 23 - significant wave height, m'
write (n_unit(2,ix,iy),*)'Col. 24 - bottom ice salinity, kg/kg'
write (n_unit(2,ix,iy),*)'Col. 25 - discharge in x direction, m**3/s'
write (n_unit(2,ix,iy),*)'Col. 26 - discharge in y direction, m**3/s'
write (n_unit(2,ix,iy),*)'Cols. 27... - top soil columns temperature, m'
write (n_unit(i,ix,iy),*)'Col. 1 - year'
write (n_unit(i,ix,iy),*)'Col. 2 - month'
write (n_unit(i,ix,iy),*)'Col. 3 - day'
write (n_unit(i,ix,iy),*)'Col. 4 - hour'
write (n_unit(i,ix,iy),*)'Col. 5 - the time from the start of integration, hours'
write (n_unit(i,ix,iy),*)'Col. 6 - surface temperature, C'
write (n_unit(i,ix,iy),*)'Col. 7 - water skin temperature, C'
write (n_unit(i,ix,iy),*)'Col. 8 - water surface temperature, C'
write (n_unit(i,ix,iy),*)'Col. 9 - mean temperature of water coloumn, C'
write (n_unit(i,ix,iy),*)'Col. 10 - maximal temperature in the water coloumn, C'
write (n_unit(i,ix,iy),*)'Col. 11 - zero-dimensional model temperature, C'
write (n_unit(i,ix,iy),*)'Col. 12 - upper ice surface temperature, C'
write (n_unit(i,ix,iy),*)'Col. 13 - upper snow surface temperature, C'
write (n_unit(i,ix,iy),*)'Col. 14 - sensible heat flux, W/m**2'
write (n_unit(i,ix,iy),*)'Col. 15 - latent heat flux, W/m**2'
write (n_unit(i,ix,iy),*)'Col. 16 - downward heat flux at the upper lake surface, W/m**2'
write (n_unit(i,ix,iy),*)'Col. 17 - downward heat flux at the lake bottom, W/m**2'
write (n_unit(i,ix,iy),*)'Col. 18 - friction velocity in the surface air layer, m/s'
write (n_unit(i,ix,iy),*)'Col. 19 - friction velocity at bottom, m/s'
write (n_unit(i,ix,iy),*)'Col. 20 - wind work at the water surface, W/m**2'
write (n_unit(i,ix,iy),*)'Col. 21 - albedo of the lake-atmosphere interface, n/d'
write (n_unit(i,ix,iy),*)'Col. 22 - shortwave radiation penetrated below surface, W/m**2'
write (n_unit(i,ix,iy),*)'Col. 23 - significant wave height, m'
write (n_unit(i,ix,iy),*)'Col. 24 - bottom ice salinity, kg/kg'
write (n_unit(i,ix,iy),*)'Col. 25 - discharge in x direction, m**3/s'
write (n_unit(i,ix,iy),*)'Col. 26 - discharge in y direction, m**3/s'
write (n_unit(i,ix,iy),*)'Cols. 27... - top soil columns temperature, m'
i = i + 1
call CHECK_UNIT(lake_series_out_unit_min,lake_series_out_unit_max, &
& n_unit(3,ix,iy))
open (n_unit(3,ix,iy),file=outpath(1:len_trim(outpath))//'time_series/'// &
& n_unit(i,ix,iy))
open (n_unit(i,ix,iy),file=outpath(1:len_trim(outpath))//'time_series/'// &
& 'conv_series'//coords_point//'.dat',status='unknown')
write (n_unit(3,ix,iy),*)'Col. 1 - year'
write (n_unit(3,ix,iy),*)'Col. 2 - month'
write (n_unit(3,ix,iy),*)'Col. 3 - day'
write (n_unit(3,ix,iy),*)'Col. 4 - hour'
write (n_unit(3,ix,iy),*)'Col. 5 - the time from the start of integration, hours'
write (n_unit(3,ix,iy),*)'Col. 6 - surface temperature, C'
write (n_unit(3,ix,iy),*)'Col. 7 - surface density, kg/m**3'
write (n_unit(3,ix,iy),*)'Col. 8 - heat flux downwards at the surface, K*m/s'
write (n_unit(3,ix,iy),*)'Col. 9 - turbulent density flux at the surface, kg/m**2/s'
write (n_unit(3,ix,iy),*)'Col. 10 - inversion depth (mass flux diagnostics), m'
write (n_unit(3,ix,iy),*)'Col. 11 - mixed layer depth, m'
write (n_unit(3,ix,iy),*)'Col. 12 - entrainment depth, m'
write (n_unit(3,ix,iy),*)'Col. 13 - convective velocity scale, m/s'
write (n_unit(3,ix,iy),*)'Col. 14 - friction velocity, m/s'
write (n_unit(3,ix,iy),*)'Col. 15 - -w_star/u_star, m/s'
write (n_unit(3,ix,iy),*)'Col. 16 - bulk Richardson number, n/d'
write (n_unit(3,ix,iy),*)'Col. 17 - maximal Brunt-Vaisala frequency, s**(-1)'
write (n_unit(3,ix,iy),*)'Col. 18 - minimal thermal conductance, m**2/s'
write (n_unit(3,ix,iy),*)'Col. 19 - thermocline thickness, m'
write (n_unit(3,ix,iy),*)'Col. 20 - Reynolds number in thermocline , n/d'
write (n_unit(3,ix,iy),*)'Col. 21 - bulk Richardson number in thermocline, n/d'
write (n_unit(i,ix,iy),*)'Col. 1 - year'
write (n_unit(i,ix,iy),*)'Col. 2 - month'
write (n_unit(i,ix,iy),*)'Col. 3 - day'
write (n_unit(i,ix,iy),*)'Col. 4 - hour'
write (n_unit(i,ix,iy),*)'Col. 5 - the time from the start of integration, hours'
write (n_unit(i,ix,iy),*)'Col. 6 - surface temperature, C'
write (n_unit(i,ix,iy),*)'Col. 7 - surface density, kg/m**3'
write (n_unit(i,ix,iy),*)'Col. 8 - heat flux downwards at the surface, K*m/s'
write (n_unit(i,ix,iy),*)'Col. 9 - turbulent density flux at the surface, kg/m**2/s'
write (n_unit(i,ix,iy),*)'Col. 10 - inversion depth (mass flux diagnostics), m'
write (n_unit(i,ix,iy),*)'Col. 11 - mixed layer depth, m'
write (n_unit(i,ix,iy),*)'Col. 12 - entrainment depth, m'
write (n_unit(i,ix,iy),*)'Col. 13 - convective velocity scale, m/s'
write (n_unit(i,ix,iy),*)'Col. 14 - friction velocity, m/s'
write (n_unit(i,ix,iy),*)'Col. 15 - -w_star/u_star, m/s'
write (n_unit(i,ix,iy),*)'Col. 16 - bulk Richardson number, n/d'
write (n_unit(i,ix,iy),*)'Col. 17 - maximal Brunt-Vaisala frequency, s**(-1)'
write (n_unit(i,ix,iy),*)'Col. 18 - minimal thermal conductance, m**2/s'
write (n_unit(i,ix,iy),*)'Col. 19 - thermocline thickness, m'
write (n_unit(i,ix,iy),*)'Col. 20 - Reynolds number in thermocline , n/d'
write (n_unit(i,ix,iy),*)'Col. 21 - bulk Richardson number in thermocline, n/d'
i = i + 1
! The output of LakeMIP file 1 format
call CHECK_UNIT(lake_series_out_unit_min,lake_series_out_unit_max, &
& n_unit(4,ix,iy))
open (n_unit(4,ix,iy),file=outpath(1:len_trim(outpath))//'time_series/'// &
& n_unit(i,ix,iy))
open (n_unit(i,ix,iy),file=outpath(1:len_trim(outpath))//'time_series/'// &
& 'LakeMIP_file1'//coords_point//'.dat',status='unknown')
! write (n_unit(4,ix,iy),*)'Col. 1 - time, days'
! write (n_unit(4,ix,iy),*)'Col. 1 - year'
! write (n_unit(4,ix,iy),*)'Col. 2 - month'
! write (n_unit(4,ix,iy),*)'Col. 3 - day in a month'
! write (n_unit(4,ix,iy),*)'Col. 4 - hour'
! write (n_unit(4,ix,iy),*)'Col. 5 - minute'
! write (n_unit(4,ix,iy),*)'Col. 6 - mixed-layer temperature &
! write (n_unit(i,ix,iy),*)'Col. 1 - time, days'
! write (n_unit(i,ix,iy),*)'Col. 1 - year'
! write (n_unit(i,ix,iy),*)'Col. 2 - month'
! write (n_unit(i,ix,iy),*)'Col. 3 - day in a month'
! write (n_unit(i,ix,iy),*)'Col. 4 - hour'
! write (n_unit(i,ix,iy),*)'Col. 5 - minute'
! write (n_unit(i,ix,iy),*)'Col. 6 - mixed-layer temperature &
! &(equals to the temperature of the upper water surface), degrees Celsius'
! write (n_unit(4,ix,iy),*)'Col. 7 - mean temperature of the water column, degrees Celsius'
! write (n_unit(4,ix,iy),*)'Col. 8 - bottom temperature, degrees Celsius'
! write (n_unit(4,ix,iy),*)'Col. 9 - mixed-layer depth, meters'
! write (n_unit(4,ix,iy),*)'Col. 10 - ice thickness, meters'
! write (n_unit(4,ix,iy),*)'Col. 11 - snow thickness, meters'
! write (n_unit(4,ix,iy),*)'Col. 12 - temperature at the ice upper surface, degrees Celsius'
! write (n_unit(4,ix,iy),*)'Col. 13 - temperature at the snow upper surface, degrees Celsius'
! write (n_unit(4,ix,iy),*)'Col. 14 - sensible heat flux at the lake-atmosphere &
! write (n_unit(i,ix,iy),*)'Col. 7 - mean temperature of the water column, degrees Celsius'
! write (n_unit(i,ix,iy),*)'Col. 8 - bottom temperature, degrees Celsius'
! write (n_unit(i,ix,iy),*)'Col. 9 - mixed-layer depth, meters'
! write (n_unit(i,ix,iy),*)'Col. 10 - ice thickness, meters'
! write (n_unit(i,ix,iy),*)'Col. 11 - snow thickness, meters'
! write (n_unit(i,ix,iy),*)'Col. 12 - temperature at the ice upper surface, degrees Celsius'
! write (n_unit(i,ix,iy),*)'Col. 13 - temperature at the snow upper surface, degrees Celsius'
! write (n_unit(i,ix,iy),*)'Col. 14 - sensible heat flux at the lake-atmosphere &
! &interface, averaged over output interval, upwards, W/m**2'
! write (n_unit(4,ix,iy),*)'Col. 15 - latent heat flux at the lake-atmosphere &
! write (n_unit(i,ix,iy),*)'Col. 15 - latent heat flux at the lake-atmosphere &
! &interface, averaged over output interval, upwards, W/m**2'
! write (n_unit(4,ix,iy),*)'Col. 16 - momentum flux at the lake-atmosphere &
! write (n_unit(i,ix,iy),*)'Col. 16 - momentum flux at the lake-atmosphere &
! &interface, averaged over output interval, positive, N/m**2'
! write (n_unit(4,ix,iy),*)'Col. 17 - upward long-wave radiation flux &
! write (n_unit(i,ix,iy),*)'Col. 17 - upward long-wave radiation flux &
! &at the lake-atmosphere interface, averaged over output interval, W/m**2'
! write (n_unit(4,ix,iy),*)'Col. 18 - downward heat flux at the lake-atmosphere &
! write (n_unit(i,ix,iy),*)'Col. 18 - downward heat flux at the lake-atmosphere &
! &interface, averaged over output interval, W/m**2'
! write (n_unit(4,ix,iy),*)'Col. 19 - surface albedo, n/d'
! write (n_unit(i,ix,iy),*)'Col. 19 - surface albedo, n/d'
i = i + 1
call CHECK_UNIT(lake_series_out_unit_min,lake_series_out_unit_max, &
& n_unit(5,ix,iy))
open (n_unit(5,ix,iy),file=outpath(1:len_trim(outpath))//'time_series/'// &
& n_unit(i,ix,iy))
open (n_unit(i,ix,iy),file=outpath(1:len_trim(outpath))//'time_series/'// &
& 'methane_series'//coords_point//'.dat',status='unknown')
write (n_unit(5,ix,iy),*)'Col. 1 - year'
write (n_unit(5,ix,iy),*)'Col. 2 - month'
write (n_unit(5,ix,iy),*)'Col. 3 - day'
write (n_unit(5,ix,iy),*)'Col. 4 - hour'
write (n_unit(5,ix,iy),*)'Col. 5 - the time from the start of integration, hours'
write (n_unit(5,ix,iy),*)'Col. 6 - the talik depth, m'
write (n_unit(5,ix,iy),*)'Col. 7 - lake surface methane concentration, mol/m**3'
write (n_unit(5,ix,iy),*)'Col. 8 - lake bottom methane concentration, mol/m**3'
write (n_unit(5,ix,iy),*)'Col. 9 - soil bottom methane concentration, mol/m**3'
write (n_unit(5,ix,iy),*)'Col. 10 - lake surface oxygen concentration, mol/m**3'
write (n_unit(5,ix,iy),*)'Col. 11 - lake bottom oxygen concentration, mol/m**3'
write (n_unit(5,ix,iy),*)'Col. 12 - total methane production due to young C decomposition, mol/(m**2*s)'
write (n_unit(5,ix,iy),*)'Col. 13 - total methane production due to old C decomposition, mol/(m**2*s)'
write (n_unit(5,ix,iy),*)'Col. 14 - methane ebullition flux at the surface, mol/(m**2*s)'
write (n_unit(5,ix,iy),*)'Col. 15 - methane plant-mediated flux at the lake bottom, mol/(m**2*s)'
write (n_unit(5,ix,iy),*)'Col. 16 - methane diffusion flux at the lake bottom, upwards, mol/(m**2*s)'
write (n_unit(5,ix,iy),*)'Col. 17 - methane turbulent flux at the lake surface, upwards, mol/(m**2*s)'
write (n_unit(5,ix,iy),*)'Col. 18 - methane ebullition flux at the surface, mg/(m**2*day)'
write (n_unit(5,ix,iy),*)'Col. 19 - methane plant-mediated flux at the lake bottom, mg/(m**2*day)'
write (n_unit(5,ix,iy),*)'Col. 20 - methane diffusion flux at the lake bottom, upwards, mg/(m**2*day)'
write (n_unit(5,ix,iy),*)'Col. 21 - methane turbulent flux at the lake surface, upwards, mg/(m**2*day)'
write (n_unit(5,ix,iy),*)'Col. 22 - methane turbulent flux at the bottom of mixed layer normalized by &
write (n_unit(i,ix,iy),*)'Col. 1 - year'
write (n_unit(i,ix,iy),*)'Col. 2 - month'
write (n_unit(i,ix,iy),*)'Col. 3 - day'
write (n_unit(i,ix,iy),*)'Col. 4 - hour'
write (n_unit(i,ix,iy),*)'Col. 5 - the time from the start of integration, hours'
write (n_unit(i,ix,iy),*)'Col. 6 - the talik depth, m'
write (n_unit(i,ix,iy),*)'Col. 7 - lake surface methane concentration, mol/m**3'
write (n_unit(i,ix,iy),*)'Col. 8 - lake bottom methane concentration, mol/m**3'
write (n_unit(i,ix,iy),*)'Col. 9 - soil bottom methane concentration, mol/m**3'
write (n_unit(i,ix,iy),*)'Col. 10 - lake surface oxygen concentration, mol/m**3'
write (n_unit(i,ix,iy),*)'Col. 11 - lake bottom oxygen concentration, mol/m**3'
write (n_unit(i,ix,iy),*)'Col. 12 - total methane production due to young C decomposition, mol/(m**2*s)'
write (n_unit(i,ix,iy),*)'Col. 13 - total methane production due to old C decomposition, mol/(m**2*s)'
write (n_unit(i,ix,iy),*)'Col. 14 - methane ebullition flux at the surface, mol/(m**2*s)'
write (n_unit(i,ix,iy),*)'Col. 15 - methane plant-mediated flux at the lake bottom, mol/(m**2*s)'
write (n_unit(i,ix,iy),*)'Col. 16 - methane diffusion flux at the lake bottom, upwards, mol/(m**2*s)'
write (n_unit(i,ix,iy),*)'Col. 17 - methane turbulent flux at the lake surface, upwards, mol/(m**2*s)'
write (n_unit(i,ix,iy),*)'Col. 18 - methane ebullition flux at the surface, mg/(m**2*day)'
write (n_unit(i,ix,iy),*)'Col. 19 - methane plant-mediated flux at the lake bottom, mg/(m**2*day)'
write (n_unit(i,ix,iy),*)'Col. 20 - methane diffusion flux at the lake bottom, upwards, mg/(m**2*day)'
write (n_unit(i,ix,iy),*)'Col. 21 - methane turbulent flux at the lake surface, upwards, mg/(m**2*day)'
write (n_unit(i,ix,iy),*)'Col. 22 - methane turbulent flux at the bottom of mixed layer normalized by &
& surface area, upwards, mg/(m**2*day)'
write (n_unit(5,ix,iy),*)'Col. 23 - methane flux from sediments in the mixed layer normalized by &
write (n_unit(i,ix,iy),*)'Col. 23 - methane flux from sediments in the mixed layer normalized by &
& surface area, upwards, mg/(m**2*day)'
write (n_unit(5,ix,iy),*)'Col. 24 - methane bubble flux at the bottom of the mixed layer normalized by &
write (n_unit(i,ix,iy),*)'Col. 24 - methane bubble flux at the bottom of the mixed layer normalized by &
& surface area, upwards, mg/(m**2*day)'
write (n_unit(5,ix,iy),*)'Col. 25 - total methane oxidation in water normalized by &
write (n_unit(i,ix,iy),*)'Col. 25 - total methane oxidation in water normalized by &
& surface area, mg/(m**2*day)'
write (n_unit(5,ix,iy),*)'Col. 26 - methane oxidation in mixed layer normalized by &
write (n_unit(i,ix,iy),*)'Col. 26 - methane oxidation in mixed layer normalized by &
& surface area, mg/(m**2*day)'
write (n_unit(5,ix,iy),*)'Col. 27 - co2 turbulent flux at the lake surface, upwards, mol/(m**2*s)'
write (n_unit(5,ix,iy),*)'Col. 28 - co2 ebullition flux at the surface, mol/(m**2*s)'
write (n_unit(5,ix,iy),*)'Col. 29 - oxygen turbulent flux at the lake surface, upwards, mol/(m**2*s)'
write (n_unit(5,ix,iy),*)'Col. 30 - oxygen ebullition flux at the surface, mol/(m**2*s)'
write (n_unit(5,ix,iy),*)'Col. 31 - methane flux through outlet, normalized by surface area, mol/(m**2*s)'
write (n_unit(5,ix,iy),*)'Col. 32..32+ns-1 - methane ebullition flux at the surface from different soil columns, mg/(m**2*day)'
write (n_unit(i,ix,iy),*)'Col. 27 - co2 turbulent flux at the lake surface, upwards, mol/(m**2*s)'
write (n_unit(i,ix,iy),*)'Col. 28 - co2 ebullition flux at the surface, mol/(m**2*s)'
write (n_unit(i,ix,iy),*)'Col. 29 - oxygen turbulent flux at the lake surface, upwards, mol/(m**2*s)'
write (n_unit(i,ix,iy),*)'Col. 30 - oxygen ebullition flux at the surface, mol/(m**2*s)'
write (n_unit(i,ix,iy),*)'Col. 31 - methane flux through outlet, normalized by surface area, mol/(m**2*s)'
write (n_unit(i,ix,iy),*)'Col. 32..32+ns-1 - methane ebullition flux at the surface from different soil columns, mg/(m**2*day)'
i = i + 1
call CHECK_UNIT(lake_series_out_unit_min,lake_series_out_unit_max, &
& n_unit(i,ix,iy))
open (n_unit(i,ix,iy),file=outpath(1:len_trim(outpath))//'time_series/'// &
& 'TKE_budget_terms_series'//coords_point//'.dat',status='unknown')
write (n_unit(i,ix,iy),*)'Col. 1 - year'
write (n_unit(i,ix,iy),*)'Col. 2 - month'
write (n_unit(i,ix,iy),*)'Col. 3 - day'
write (n_unit(i,ix,iy),*)'Col. 4 - hour'
write (n_unit(i,ix,iy),*)'Col. 5 - the time from the start of integration, hours'
write (n_unit(i,ix,iy),*)'Col. 6 - TKE tendency, m**2/s**3'
write (n_unit(i,ix,iy),*)'Col. 7 - TKE turbulent transport, m**2/s**3'
write (n_unit(i,ix,iy),*)'Col. 8 - TKE shear production, m**2/s**3'
write (n_unit(i,ix,iy),*)'Col. 9 - TKE buoyancy production/sink, m**2/s**3'
write (n_unit(i,ix,iy),*)'Col. 10 - TKE dissipation rate, m**2/s**3'
count_out(ix,iy) = int(time/(dt_out%par*hour_sec))
endif
......@@ -1198,6 +1219,11 @@ if (int(time/(dt_out%par*hour_sec))>count_out(ix,iy) .or. &
! & - fdiff2(1:2)*mfs, & ! two-meth
! & - fdiff_lake_surf2(1:2)*mfs ! two-meth
format_char = '(3i7,f8.2,f13.2,5E15.5)'
write (n_unit(6,ix,iy),format_char) year,month,day,hour, &
& time/hour_sec, &
& TKE_budget_terms(1:5)
count_out(ix,iy) = count_out(ix,iy) + 1
endif
......
source/model/turb_mod.f90
+ 16
− 3
View file @ 248d1a6e
......@@ -4,6 +4,7 @@ use NUMERICS, only : PROGONKA, IND_STAB_FACT_DB
use NUMERIC_PARAMS, only : vector_length, small_value
use INOUT, only : CHECK_UNIT
use LAKE_DATATYPES, only : ireals, iintegers
use DRIVING_PARAMS, only : missing_value
contains
SUBROUTINE K_EPSILON(ix, iy, nx, ny, year, month, day, &
......@@ -29,7 +30,8 @@ use DRIVING_PARAMS , only : &
& kwe, &
& M, &
& omp, &
& eos, lindens
& eos, lindens, &
& zserout
use ATMOS, only : &
& uwind, vwind
......@@ -77,11 +79,12 @@ use ARRAYS_TURB , only : &
& H_entrainment, &
& Buoyancy0, signwaveheight, &
& knum, &
& Eseiches
& Eseiches, &
& TKE_budget_terms
use ARRAYS_GRID, only : &
& ddz, ddz2, ddz05, ddz052, ddz054, &
& dzeta_int, dzeta_05int
& dzeta_int, dzeta_05int, z_half
use ARRAYS_BATHYM, only : &
& dhw, dhw0, &
......@@ -1012,6 +1015,16 @@ allocate (rhotemp(1:M+1), rhosal(1:M+1))
Gen(i) = Gen(i)*k2(i)
enddo
! Components of TKE budget for output
if (zserout%par /= missing_value) then
i = minloc(abs(z_half(:) - zserout%par),1)
TKE_budget_terms(1) = ( E2(i) - E1(i) )/dt
TKE_budget_terms(2) = TKE_turb_trans(i)
TKE_budget_terms(3) = Gen(i)
TKE_budget_terms(4) = S(i)
TKE_budget_terms(5) = eps1(i)
endif
S_integr_positive = 0._ireals
S_integr_negative = 0._ireals
Gen_integr = 0._ireals
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment