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

Initial soil temperature profile in the deepest soil column is added to setup file

parent a70c1372
Branches
Tags
No related merge requests found
Show whitespace changes
Inline Side-by-side
Showing
with 180 additions and 144 deletions
docs/users_guide.aux
+ 7
− 7
View file @ 8046e9d5
\relax
\@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 }
\@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}}
\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}\protected@file@percent }
\@writefile{lot}{\contentsline {table}{\numberline {3}{\ignorespaces The parameters of setup file\relax }}{5}\protected@file@percent }
\@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}}
\gdef \LT@ii {\LT@entry
{1}{201.44666pt}\LT@entry
{1}{296.92755pt}}
\@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 }
\@writefile{toc}{\contentsline {section}{\numberline {0.4}Running the model}{9}}
\@writefile{toc}{\contentsline {section}{\numberline {0.5}Model output}{9}}
docs/users_guide.log
+ 126
− 127
View file @ 8046e9d5
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
This is pdfTeX, Version 3.14159265-2.6-1.40.18 (TeX Live 2017/Debian) (preloaded format=pdflatex 2020.8.10) 12 OCT 2020 09:46
entering extended mode
restricted \write18 enabled.
file:line:error style messages enabled.
%&-line parsing enabled.
**users_guide.tex
(./users_guide.tex
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)
LaTeX2e <2017-04-15>
Babel <3.18> and hyphenation patterns for 84 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)
)
\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
\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
\abovecaptionskip=\skip41
\belowcaptionskip=\skip42
\bibindent=\dimen102
)
(/usr/local/texlive/2019/texmf-dist/tex/latex/amsmath/amsmath.sty
Package: amsmath 2018/12/01 v2.17b AMS math features
(/usr/share/texlive/texmf-dist/tex/latex/amsmath/amsmath.sty
Package: amsmath 2017/09/02 v2.17a AMS math features
\@mathmargin=\skip43
For additional information on amsmath, use the `?' option.
(/usr/local/texlive/2019/texmf-dist/tex/latex/amsmath/amstext.sty
(/usr/share/texlive/texmf-dist/tex/latex/amsmath/amstext.sty
Package: amstext 2000/06/29 v2.01 AMS text
(/usr/local/texlive/2019/texmf-dist/tex/latex/amsmath/amsgen.sty
(/usr/share/texlive/texmf-dist/tex/latex/amsmath/amsgen.sty
File: amsgen.sty 1999/11/30 v2.0 generic functions
\@emptytoks=\toks14
\ex@=\dimen103
))
(/usr/local/texlive/2019/texmf-dist/tex/latex/amsmath/amsbsy.sty
(/usr/share/texlive/texmf-dist/tex/latex/amsmath/amsbsy.sty
Package: amsbsy 1999/11/29 v1.2d Bold Symbols
\pmbraise@=\dimen104
)
(/usr/local/texlive/2019/texmf-dist/tex/latex/amsmath/amsopn.sty
(/usr/share/texlive/texmf-dist/tex/latex/amsmath/amsopn.sty
Package: amsopn 2016/03/08 v2.02 operator names
)
\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
\inf@bad=\count88
LaTeX Info: Redefining \frac on input line 213.
\uproot@=\count89
\leftroot@=\count90
LaTeX Info: Redefining \overline on input line 375.
\classnum@=\count91
\DOTSCASE@=\count92
LaTeX Info: Redefining \ldots on input line 472.
LaTeX Info: Redefining \dots on input line 475.
LaTeX Info: Redefining \cdots on input line 596.
\Mathstrutbox@=\box26
\strutbox@=\box27
\big@size=\dimen105
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
LaTeX Font Info: Redeclaring font encoding OML on input line 712.
LaTeX Font Info: Redeclaring font encoding OMS on input line 713.
\macc@depth=\count93
\c@MaxMatrixCols=\count94
\dotsspace@=\muskip10
\c@parentequation=\count96
\dspbrk@lvl=\count97
\c@parentequation=\count95
\dspbrk@lvl=\count96
\tag@help=\toks15
\row@=\count98
\column@=\count99
\maxfields@=\count100
\row@=\count97
\column@=\count98
\maxfields@=\count99
\andhelp@=\toks16
\eqnshift@=\dimen106
\alignsep@=\dimen107
......@@ -79,62 +79,62 @@ LaTeX Font Info: Redeclaring font encoding OMS on input line 730.
\multlinegap=\skip44
\multlinetaggap=\skip45
\mathdisplay@stack=\toks18
LaTeX Info: Redefining \[ on input line 2844.
LaTeX Info: Redefining \] on input line 2845.
LaTeX Info: Redefining \[ on input line 2817.
LaTeX Info: Redefining \] on input line 2818.
)
(/usr/local/texlive/2019/texmf-dist/tex/latex/preprint/authblk.sty
(/usr/share/texlive/texmf-dist/tex/latex/preprint/authblk.sty
Package: authblk 2001/02/27 1.3 (PWD)
\affilsep=\skip46
\@affilsep=\skip47
\c@Maxaffil=\count101
\c@authors=\count102
\c@affil=\count103
\c@Maxaffil=\count100
\c@authors=\count101
\c@affil=\count102
)
(/usr/local/texlive/2019/texmf-dist/tex/latex/listings/listings.sty
(/usr/local/texlive/2019/texmf-dist/tex/latex/graphics/keyval.sty
(/usr/share/texlive/texmf-dist/tex/latex/listings/listings.sty
(/usr/share/texlive/texmf-dist/tex/latex/graphics/keyval.sty
Package: keyval 2014/10/28 v1.15 key=value parser (DPC)
\KV@toks@=\toks19
)
\lst@mode=\count104
\lst@gtempboxa=\box29
\lst@mode=\count103
\lst@gtempboxa=\box28
\lst@token=\toks20
\lst@length=\count105
\lst@length=\count104
\lst@currlwidth=\dimen112
\lst@column=\count106
\lst@pos=\count107
\lst@column=\count105
\lst@pos=\count106
\lst@lostspace=\dimen113
\lst@width=\dimen114
\lst@newlines=\count108
\lst@lineno=\count109
\lst@newlines=\count107
\lst@lineno=\count108
\lst@maxwidth=\dimen115
(/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/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/listings.cfg
File: listings.cfg 2019/02/27 1.8b listings configuration
(/usr/share/texlive/texmf-dist/tex/latex/listings/listings.cfg
File: listings.cfg 2015/06/04 1.6 listings configuration
))
Package: listings 2019/02/27 1.8b (Carsten Heinz)
Package: listings 2015/06/04 1.6 (Carsten Heinz)
(/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/latex/geometry/geometry.sty
Package: geometry 2010/09/12 v5.6 Page Geometry
(/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/ifpdf.sty
Package: ifpdf 2017/03/15 v3.2 Provides the ifpdf switch
)
(/usr/local/texlive/2019/texmf-dist/tex/generic/oberdiek/ifvtex.sty
(/usr/share/texlive/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/local/texlive/2019/texmf-dist/tex/generic/ifxetex/ifxetex.sty
(/usr/share/texlive/texmf-dist/tex/generic/ifxetex/ifxetex.sty
Package: ifxetex 2010/09/12 v0.6 Provides ifxetex conditional
)
\Gm@cnth=\count112
\Gm@cntv=\count113
\c@Gm@tempcnt=\count114
\Gm@cnth=\count111
\Gm@cntv=\count112
\c@Gm@tempcnt=\count113
\Gm@bindingoffset=\dimen116
\Gm@wd@mp=\dimen117
\Gm@odd@mp=\dimen118
......@@ -145,43 +145,41 @@ Package: ifxetex 2010/09/12 v0.6 Provides ifxetex conditional
\Gm@layoutvoffset=\dimen123
\Gm@dimlist=\toks21
)
(/usr/local/texlive/2019/texmf-dist/tex/latex/tools/array.sty
Package: array 2018/12/30 v2.4k Tabular extension package (FMi)
(/usr/share/texlive/texmf-dist/tex/latex/tools/array.sty
Package: array 2016/10/06 v2.4d Tabular extension package (FMi)
\col@sep=\dimen124
\ar@mcellbox=\box31
\extrarowheight=\dimen125
\NC@list=\toks22
\extratabsurround=\skip48
\backup@length=\skip49
\ar@cellbox=\box32
)
(/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
(/usr/share/texlive/texmf-dist/tex/latex/tools/longtable.sty
Package: longtable 2014/10/28 v4.11 Multi-page Table package (DPC)
\LTleft=\skip50
\LTright=\skip51
\LTpre=\skip52
\LTpost=\skip53
\LTchunksize=\count115
\LTchunksize=\count114
\LTcapwidth=\dimen126
\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@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@p@ftn=\toks23
)
(/usr/local/texlive/2019/texmf-dist/tex/latex/was/gensymb.sty
(/usr/share/texlive/texmf-dist/tex/latex/was/gensymb.sty
Package: gensymb 2003/07/02 v1.0 (WaS)
)
(/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/caption.sty
Package: caption 2016/02/21 v3.3-144 Customizing captions (AR)
(/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.
(/usr/share/texlive/texmf-dist/tex/latex/caption/caption3.sty
Package: caption3 2016/05/22 v1.7-166 caption3 kernel (AR)
Package caption3 Info: TeX engine: e-TeX on input line 67.
\captionmargin=\dimen127
\captionmargin@=\dimen128
\captionwidth=\dimen129
......@@ -190,13 +188,12 @@ Package caption3 Info: TeX engine: e-TeX on input line 64.
\caption@parindent=\dimen132
\caption@hangindent=\dimen133
)
\c@caption@flags=\count120
\c@ContinuedFloat=\count121
\c@ContinuedFloat=\count119
Package caption Info: listings package is loaded.
Package caption Info: longtable package is loaded.
(/usr/local/texlive/2019/texmf-dist/tex/latex/caption/ltcaption.sty
Package: ltcaption 2018/08/26 v1.4-95 longtable captions (AR)
(/usr/share/texlive/texmf-dist/tex/latex/caption/ltcaption.sty
Package: ltcaption 2013/06/09 v1.4-94 longtable captions (AR)
))
Package geometry Warning: Over-specification in `h'-direction.
......@@ -221,7 +218,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=\count122
\c@lstlisting=\count120
*geometry* driver: auto-detecting
*geometry* detected driver: pdftex
......@@ -273,10 +270,10 @@ Package caption Info: End \AtBeginDocument code.
[1
{/usr/local/texlive/2019/texmf-var/fonts/map/pdftex/updmap/pdftex.map}]
{/var/lib/texmf/fonts/map/pdftex/updmap/pdftex.map}]
LaTeX Font Info: Try loading font information for OML+cmr on input line 66.
(/usr/local/texlive/2019/texmf-dist/tex/latex/base/omlcmr.fd
(/usr/share/texlive/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
......@@ -341,40 +338,42 @@ Overfull \hbox (14.67734pt too wide) in alignment at lines 416--461
[]
[7]
Overfull \hbox (14.67734pt too wide) in alignment at lines 461--463
Overfull \hbox (14.67734pt too wide) in alignment at lines 461--465
[] []
[]
Underfull \vbox (badness 10000) has occurred while \output is active []
[8]
LaTeX Font Info: Try loading font information for OMS+cmr on input line 490.
LaTeX Font Info: Try loading font information for OMS+cmr on input line 492.
(/usr/local/texlive/2019/texmf-dist/tex/latex/base/omscmr.fd
(/usr/share/texlive/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 490.
[9]
(./users_guide.aux) )
(Font) Font shape `OMS/cmsy/m/n' tried instead on input line 492.
[9] (./users_guide.aux) )
Here is how much of TeX's memory you used:
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
3926 strings out of 492982
58072 string characters out of 6134896
156563 words of memory out of 5000000
7495 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,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).
36i,8n,39p,563b,1088s stack positions out of 5000i,500n,10000p,200000b,80000s
</usr/share/texlive/texmf-dist/fonts/type1/public/ams
fonts/cm/cmbx12.pfb></usr/share/texlive/texmf-dist/fonts/type1/public/amsfonts/
cm/cmmi10.pfb></usr/share/texlive/texmf-dist/fonts/type1/public/amsfonts/cm/cmr
10.pfb></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/s
hare/texlive/texmf-dist/fonts/type1/public/amsfonts/cm/cmr7.pfb></usr/share/tex
live/texmf-dist/fonts/type1/public/amsfonts/cm/cmsl10.pfb></usr/share/texlive/t
exmf-dist/fonts/type1/public/amsfonts/cm/cmsy10.pfb></usr/share/texlive/texmf-d
ist/fonts/type1/public/amsfonts/cm/cmsy7.pfb></usr/share/texlive/texmf-dist/fon
ts/type1/public/amsfonts/cm/cmti10.pfb>
Output written on users_guide.pdf (9 pages, 162451 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 @ 8046e9d5
No preview for this file type
docs/users_guide.synctex.gz
+ 0
− 0
View file @ 8046e9d5
No preview for this file type
docs/users_guide.tex
+ 3
− 1
View file @ 8046e9d5
......@@ -403,7 +403,9 @@ give significant speedup \\
\multicolumn{2}{|c|}{\large{Initial conditions}} \\
\hline
\emph{T\_profile} & Initial temperature (\celsius, col.2), salinity ($kg/kg$, col.3), $CH_4$ ($mol/m^3$, col.4), $CO_2$ ($mol/m^3$, col.5), $O_2$ ($mol/m^3$, col.6) profiles, at depth levels ($m$, col.1) \\
\emph{T\_profile} & Initial temperature (\celsius, col.2), salinity ($kg/kg$, col.3), $CH_4$ ($mol/m^3$, col.4), $CO_2$ ($mol/m^3$, col.5), $O_2$ ($mol/m^3$, col.6) profiles, at depth levels ($m$, col.1); used if \emph{init\_T}=3 \\
\hline
\emph{T\_soilprofile} & Initial soil temperature (\celsius, col.2) profiles, at depth levels ($m$, col.1); used if \emph{init\_T}=3 \\
\hline
\multicolumn{2}{|c|}{\large{Tributaries and effluents}} \\
......
setup/Mojai2016-2017_setup.dat
+ 4
− 0
View file @ 8046e9d5
......@@ -199,6 +199,9 @@ T_profile 16
16. 8.4 0. 5.0E-3 0.20 0.01 1.3E-2
17. 8.1 0. 5.0E-3 0.30 0.01 1.3E-2
17.9 7.9 0. 5.0E-3 0.30 0.01 1.3E-2
T_soilprofile 2
0. 7.9
10. 4.
#.
#----------------------------------------------------------------------------------------
# BOUNDARY CONDITIONS: TRIBUTARIES AND EFFLUENTS
......@@ -267,6 +270,7 @@ everystep 0
time_series 1
dt_out 1.
nscreen 1000
zserout -999.
scale_output 0
accum_begin 2003060100
accum_end 2004060100
......
source/model/driving_params_mod.f90
+ 20
− 1
View file @ 8046e9d5
......@@ -97,7 +97,7 @@ type, public :: initprof
real(kind=ireals), allocatable :: Tinitprof(:), Sinitprof(:)
real(kind=ireals), allocatable :: ch4initprof(:), co2initprof(:), o2initprof(:), phosphinitprof(:)
end type initprof
type(initprof) :: ip
type(initprof) :: ip, isp
! Group of 1D arrays
type, public :: grarr1
......@@ -140,6 +140,7 @@ type(intpar), target :: saltice
type(intpar), target :: ifbubble
type(intpar), target :: carbon_model
type(intpar), target :: Tinitlength
type(intpar), target :: Tsoilinitlength
type(intpar), target :: dyn_pgrad
type(intpar), target :: botfric
type(intpar), target :: outflpar
......@@ -384,6 +385,10 @@ if (firstcall) then
write(*,*) 'The parameter T_profile is missing in setup file'
all_par_present = .false.
endif
if (.not.Tsoilinitlength%ok) then
write(*,*) 'The parameter T_soilprofile is missing in setup file'
all_par_present = .false.
endif
if (.not.rtemp%ok) then
write(*,*) 'The parameter rtemp is missing in setup file'
all_par_present = .false.
......@@ -536,6 +541,20 @@ lineread: if (line(n1:n1) /= '#') then
deallocate(work)
Tinitlength%ok = .true.
CASE ('T_soilprofile')
Tsoilinitlength%par &
& = igetvarval (n1,n2,line,'Init. T-soilprof length ')
isp%lenprof = Tsoilinitlength%par
allocate (isp%Tinitprof (1:isp%lenprof) )
allocate (isp%zTinitprof(1:isp%lenprof) )
ncol = 2
allocate (work(Tsoilinitlength%par,ncol))
call READPROFILE(nunit,Tsoilinitlength%par,ncol,work)
isp%zTinitprof (:) = work (:,1)
isp%Tinitprof (:) = work (:,2)
deallocate(work)
Tsoilinitlength%ok = .true.
! The group of output controls
CASE ('accum_begin')
j = igetvarval(n1,n2,line,'Accumulation period begin')
......
source/model/init_var_mod.f90
+ 19
− 7
View file @ 8046e9d5
......@@ -8,7 +8,7 @@
& init_T, skin, zero_model, &
& h10, l10, hs10, ls10, tempair, Ts0, Tb0, Tm, Tbb0, &
& Sals0, Salb0, us0, vs0, h_ML0, &
& rosoil, rosdry, por, depth_area, ip, &
& rosoil, rosdry, por, depth_area, ip, isp, &
& flag_snow, flag_snow_init, itop, nstep, itherm, &
& h1, l1, hs1, ls1, &
......@@ -72,7 +72,7 @@
real(kind=ireals), intent(in) :: h_ML0
real(kind=ireals), intent(in) :: rosoil(1:ns), rosdry(1:ns), por(1:ns)
real(kind=ireals), intent(in) :: depth_area(1:ndatamax,1:2) ! Data for lake bathymetry
type(initprof), intent(in) :: ip
type(initprof), intent(in) :: ip, isp
real(kind=ireals), intent(in) :: zsoil(1:ns)
real(kind=ireals), intent(in) :: pressure
real(kind=ireals), intent(in) :: dzeta_int(1:M+1)
......@@ -199,7 +199,8 @@
endif
Sals1(1:ns,1:nsoilcols) = Salb0*row0/( rosdry(1)*(1 - por(1)) ) ! assuming, that salinity in ground is the same
Sals1(1:ns,1:nsoilcols) = Salb0*row0/( rosdry(1)*(1 - por(1)) ) ! assuming, that salinity in ground
! is the same
! as one in near bottom layer of water
allocate(pressoil(1:ns))
......@@ -238,9 +239,20 @@
h_talik = 0.
! Linear temperature profile in soil
if (init_T == 1 .or. init_T == 2) then
do i = 1, ns
Tsoil1(i,nsoilcols) = Tb0 + float(i-1)/float(ns-1)*(Tbb0-Tb0)
enddo
! The initial soil temperature profile is given from the input file
elseif (init_T == 3) then
call LININTERPOL (isp%zTinitprof,isp%Tinitprof,isp%lenprof, &
& zsoil,Tsoil1(1,nsoilcols),ns,flag)
if (.not.flag) then
print*, 'The error while interpolating the initial &
&soil temperature profile: terminating program'
STOP
endif
endif
endif
......
source/model/lake.f90
+ 1
− 1
View file @ 8046e9d5
......@@ -650,7 +650,7 @@ if (init(ix,iy) == 0_iintegers) then
& init_T, skin%par, zero_model%par, &
& h10, l10, hs10, ls10, tempair, Ts0, Tb0, Tm, Tbb0, &
& Sals0, Salb0, us0, vs0, h_ML0, &
& rosoil, rosdry, por, depth_area, ip, &
& rosoil, rosdry, por, depth_area, ip, isp, &
& flag_snow, flag_snow_init, itop, nstep, itherm, &
& h1, l1, hs1, ls1, &
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment