Commit 7a912492 authored by Robert Schweppe's avatar Robert Schweppe

Merge branch 'develop' into 'master'

Creating version 0.1

See merge request !10
parents fdf8b286 0cecfaab
File deleted
......@@ -264,6 +264,11 @@ mpr_flex
mpr_flex_dummy
*.release*
*.debug*
wrfinput.nc
*wrfinput.nc
wrfinput_125.nc
wrfinput_25.nc
.idea/misc.xml
.idea/modules.xml
.idea/mpr_flex.iml
.idea/codeStyles/
*.f90.new
<?xml version="1.0" encoding="UTF-8"?>
<module type="PYTHON_MODULE" version="4">
<component name="NewModuleRootManager">
<content url="file://$MODULE_DIR$" />
<orderEntry type="inheritedJdk" />
<orderEntry type="sourceFolder" forTests="false" />
</component>
<component name="TestRunnerService">
<option name="PROJECT_TEST_RUNNER" value="Unittests" />
</component>
</module>
\ No newline at end of file
<component name="ProjectCodeStyleConfiguration">
<code_scheme name="Project" version="173">
<option name="ALIGN_MULTILINE_PARAMETERS_IN_CALLS" value="true" />
</code_scheme>
</component>
\ No newline at end of file
<component name="ProjectCodeStyleConfiguration">
<state>
<option name="PREFERRED_PROJECT_CODE_STYLE" value="Default" />
</state>
</component>
\ No newline at end of file
<component name="ProjectDictionaryState">
<dictionary name="ottor" />
</component>
\ No newline at end of file
<?xml version="1.0" encoding="UTF-8"?>
<project version="4">
<component name="ProjectRootManager" version="2" project-jdk-name="Python 3.6 (playground)" project-jdk-type="Python SDK" />
</project>
\ No newline at end of file
<?xml version="1.0" encoding="UTF-8"?>
<project version="4">
<component name="ProjectModuleManager">
<modules>
<module fileurl="file://$PROJECT_DIR$/.idea/20180205_mpr_playground.iml" filepath="$PROJECT_DIR$/.idea/20180205_mpr_playground.iml" />
</modules>
</component>
</project>
\ No newline at end of file
<?xml version="1.0" encoding="UTF-8"?>
<project version="4">
<component name="VcsDirectoryMappings">
<mapping directory="$PROJECT_DIR$" vcs="Git" />
</component>
</project>
\ No newline at end of file
This diff is collapsed.
project: MPR
src_dir: ./src/
output_dir: ./doc
project_download: https://git.ufz.de/chs/MPR.git
website: http://www.ufz.de/index.php?en=40114
summary: MPR as a flexible standalone tool.
author: CHS
email: robert.schweppe@ufz.de
media_dir: ./media
docmark: !
predocmark: >
media_dir: ./media
docmark_alt: #
predocmark_alt: <
coloured_edges: true
display: public
protected
private
source: true
graph: true
search: true
warn: false
MPR project description
@note
lorem ipsum
@endnote
@warning
@todo
@bug
......@@ -101,7 +101,7 @@ MAKEDPATH := $(CONFIGPATH) # where is the make.d.sh script
CHECKPATH := test # path for $(CHECKPATH)/test* and $(CHECKPATH)/check* directories if target is check
DOXCONFIG := ./doxygen.config # the doxygen config file
#
PROGNAME := mpr_flex # Name of executable
PROGNAME := mpr # Name of executable
LIBNAME := #libminpack.a # Name of library
#
# Options
......
# Set aliases so that one can, for example, say ifort to invoke standard intel11 on eve
ifneq (,$(findstring $(compiler),gnu gfortran gcc gfortran49 gcc49))
icompiler := gnu61
icompiler := gnu81
endif
ifneq (,$(findstring $(compiler),nag nag60 nagfor))
icompiler := nag61
icompiler := nag62
endif
......@@ -53,7 +53,7 @@ ifeq ($(release),debug)
# Ordinary assignment of pointer function fname reference to pointer variable vname
# Last statement of DO loop body is an unconditional jump
# Array constructor has a polymorphic element but the constructor value will not be polymorphic
F90FLAGS += -C -C=dangling -g -nan -O0 -strict95 # -w=ques -g90
F90FLAGS += -C=all -g -nan -O0 -strict95 # -w=ques -g90
ifneq ($(openmp),true)
F90FLAGS += -gline
endif
......@@ -88,7 +88,7 @@ ifeq ($(netcdf),netcdf3)
NCDIR :=
else
NCDIR := /usr/local
NCFDIR := /usr/local/netcdf-fortran-4.4.3-nagfor
NCFDIR := /usr/local/netcdf-fortran-4.4.4-nagfor
CURLLIB := /usr/lib
SZLIB := /usr/local/lib
HDF5LIB := /usr/local/lib
......
......@@ -2,11 +2,11 @@
ifneq (,$(findstring $(compiler),gnu gfortran gcc))
icompiler := gnu48
endif
ifneq (,$(findstring $(compiler),intel ifort ifort13))
icompiler := intel13
ifneq (,$(findstring $(compiler),intel ifort ifort13 ifort18))
icompiler := intel18
endif
ifneq (,$(findstring $(compiler),nag nagfor))
icompiler := nag60
icompiler := nag62
endif
ifneq (,$(findstring $(compiler),pgi pgfortran pgfortran159))
icompiler := pgi159
......
# -*- Makefile -*-
#
# Setup file for Intel compiler 13.1.0 at eve.ufz.de
#
# LICENSE
# This file is part of the UFZ makefile project.
#
# The UFZ makefile project is free software: you can redistribute it and/or modify
# it under the terms of the GNU Lesser General Public License as published by
# the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# The UFZ makefile project is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU Lesser General Public License for more details.
#
# You should have received a copy of the GNU Lesser General Public License
# along with the UFZ makefile project. If not, see <http://www.gnu.org/licenses/>.
#
# Copyright 2012 Matthias Cuntz
# The Makefile assumes the standard bin, include and lib directories
# i.e. if SOMEDIR = /path/to/lib, the make file will define the two dirs
# SOMEINC ?= SOMEDIR/include
# SOMELIB ?= SOMEDIR/lib
# Define subirectories if another structure
# Paths
# INTEL := /opt/intel/composer_xe_2013.1.117
# INTELDIR := /usr/local/intel/13.0.1
INTEL := /software/intel/parallel-studio/2018_u2-1/compilers_and_libraries_2018.2.199/linux
INTELDIR := /software/intel/parallel-studio/2018_u2-1/compilers_and_libraries_2018.2.199/linux
INTELLIB := $(INTELDIR)/compiler/lib/intel64_lin
INTELBIN := $(INTEL)/bin/intel64
GNUDIR := /usr/local/gcc/4.8.1-3_CentOS6/
OPENSSLLIB := /usr/local/openssl/1.0.2-1/usr/lib
# Compiling
F90 := $(INTELBIN)/ifort
FC := $(F90)
CC := $(INTELBIN)/icc
CPP := /usr/bin/cpp
ifeq ($(release),debug)
F90FLAGS += -warn all -g -debug -traceback -fp-stack-check -O0 -debug -gxx-name=/usr/local/gcc/6.2.0-1/bin/g++
F90FLAGS += -check all
# F90FLAGS += -check all,noarg_temp_created
FCFLAGS += -g -debug -traceback -fp-stack-check -O0 -debug -gxx-name=/usr/local/gcc/6.2.0-1/bin/g++
CFLAGS += -g -debug -traceback -fp-stack-check -O0 -debug -Wall -Wcheck
else
# -fast (=-ipo -O3 -static) Multifile interprocedure optimization
# -vec-report1 to see vectorized loops; -vec-report2 to see also non-vectorized loops
F90FLAGS += -O3 -qoverride-limits -gxx-name=/usr/local/gcc/6.2.0-1/bin/g++
FCFLAGS += -O3 -qoverride-limits -gxx-name=/usr/local/gcc/6.2.0-1/bin/g++
CFLAGS += -O3
endif
F90FLAGS += -nofixed -assume byterecl -cpp -fp-model source -m64 -assume realloc-lhs # precise -> source: suppress warning, computation identical
MODFLAG := -module # space significant
FCFLAGS += -assume byterecl -cpp -fp-model source -m64 -fixed
CFLAGS += -m64 -std=c99
DEFINES += -DINTEL -DINTEL_COMPILER
LDFLAGS += -gxx-name=/usr/local/gcc/4.8.1-3_CentOS6/bin/g++
# Linking
ifeq ($(istatic),static)
LIBS += -static-intel
endif
LIBS += -L$(INTELLIB) -limf -lm -lsvml
ifneq ($(istatic),static)
LIBS += -lintlc
endif
RPATH += -Wl,-rpath,$(INTELLIB) -Wl,-rpath,$(OPENSSLLIB)
iLDPATH = $(GNUDIR)/lib64:/usr/local/cloog/0.18.0-2/lib:/usr/local/isl/0.11.1-2/lib:/usr/local/mpc/1.0.1-3/lib:/usr/local/mpfr/3.1.2-2/lib:/usr/local/gmp/5.1.2-1/lib
ifneq ($(LDPATH),)
LDPATH += :$(iLDPATH)
else
LDPATH := $(iLDPATH)
endif
# IMSL
IMSLDIR := /usr/local/imsl/imsl/fnl700/rdhin111e64
# MKL
MKLDIR := $(INTEL)/mkl
MKLINC := $(MKLDIR)/include/intel64/lp64
MKLLIB := $(INTELDIR)/mkl/lib/intel64
MKL95DIR := $(MKLDIR)
MKL95INC := $(MKL95DIR)/include/intel64/lp64
MKL95LIB := $(MKL95DIR)/lib/intel64
# NETCDF
ifeq ($(netcdf),netcdf3)
NCDIR :=
else
ZLIB := /usr/local/zlib/1.2.11-1_intel_18.0.2-1/lib
SZLIB := /usr/local/szip/2.1.1-1_intel_18.0.2-1/lib
HDF5LIB := /usr/local/hdf5/1.10.1-1_intel_18.0.2-1/lib
NCDIR := /usr/local/netcdf/4.6.1-1_intel_18.0.2-1/
NCFDIR := /usr/local/netcdf-fortran/4.4.4-1_intel_18.0.2-1/
endif
# PROJ
PROJ4DIR :=
FPROJDIR :=
# LAPACK
LAPACKDIR := /usr/local/lapack/3.5.0-1_gcc_4.8.1_CentOS6
GFORTRANDIR := /usr/local/gcc/4.8.1-3_CentOS6/
GFORTRANLIB := $(GFORTRANDIR)/lib64
# MPI
MPIDIR := /usr/local/openmpi/intel/1.8.4-2
# Documentation
DOXYGENDIR := /usr/local/doxygen/1.8.7-1_gcc_4.8.1_CentOS6/bin
DOTDIR := /usr/bin
TEXDIR := /usr/local/texlive/2011/bin/x86_64-linux
PERLDIR := /usr/bin
iiLDPATH := /usr/local/flex/2.5.39-1_gcc_4.8.1_CentOS6/lib:/usr/local/bison/3.0.2-1_gcc_4.8.1_CentOS6/lib
ifneq ($(LDPATH),)
LDPATH += :$(iiLDPATH)
else
LDPATH := $(iiLDPATH)
endif
# Check some dependices, e.g. IMSL needs intel11 on eve
ifneq (,$(findstring $(system),eve))
ifneq (,$(findstring $(imsl),vendor imsl))
ifneq ($(icompiler),intel11)
$(error Error: IMSL needs intel11.0.075, set 'compiler=intel11')
endif
ifeq ($(imsl),vendor)
ifeq (,$(findstring $(mkl),mkl mkl95))
$(error Error: IMSL vendor needs MKL, set 'mkl=mkl' or 'mkl=mkl95')
endif
endif
endif
endif
# Lapack on Eve needs libgfortran
ifeq ($(lapack),true)
ifneq (exists, $(shell if [ -d "$(GFORTRANDIR)" ] ; then echo 'exists' ; fi))
$(error Error: GFORTRAN path '$(GFORTRANDIR)' not found.)
endif
GFORTRANLIB ?= $(GFORTRANDIR)/lib
iLIBS += -L$(GFORTRANLIB) -lgfortran
RPATH += -Wl,-rpath,$(GFORTRANLIB)
endif
# -*- Makefile -*-
#
# Setup file for NAG compiler 6.0 at eve.ufz.de
# Setup file for NAG compiler 6.2 at eve.ufz.de
#
# LICENSE
# This file is part of the UFZ makefile project.
......@@ -28,7 +28,7 @@
# Define subdirectories if another structure
# Paths
NAGDIR := /opt/nag/6.0
NAGDIR := /opt/nag/6.2
NAGLIB := $(NAGDIR)/lib
NAGBIN := $(NAGDIR)/bin
......@@ -45,9 +45,9 @@ ifeq ($(release),debug)
#F90FLAGS += -C=all -g -gline -nan -O0 -strict95 -g90
# -gline not supported in NAG v5.3 under OpenMP
ifeq ($(openmp),true)
F90FLAGS += -C -C=dangling -g -nan -O0 -strict95 -g90
F90FLAGS += -C=all -g -nan -O0 -strict95 -g90
else
F90FLAGS += -C -C=dangling -g -gline -nan -O0 -strict95 -g90
F90FLAGS += -C=all -g -gline -nan -O0 -strict95 -g90
endif
FCFLAGS += $(F90FLAGS)
CFLAGS += -pedantic -Wall -W -O -g -Wno-maybe-uninitialized
......
......@@ -43,9 +43,9 @@ CC := $(GNUBIN)/gcc
CPP := /usr/bin/cpp
ifeq ($(release),debug)
ifeq ($(openmp),true)
F90FLAGS += -C -C=dangling -g -nan -O0 -strict95 -g90
F90FLAGS += -C=all -g -nan -O0 -strict95 -g90
else
F90FLAGS += -C -C=dangling -g -gline -nan -O0 -strict95 -g90
F90FLAGS += -C=all -g -gline -nan -O0 -strict95 -g90
endif
FCFLAGS += $(F90FLAGS)
#FCFLAGS += -C -C=dangling -g -nan -O0 -fixed
......
&mainconfig
nInput_Fields = 3
nEff_Params = 3
nDims = 4
out_filename = "cosby_-1.nc"
out_filename = "cosby_-1_fine.nc"
dim_name_alias(:,1) = "XLONG", "lon", "west_east", "x"
dim_name_alias(:,2) = "XLAT", "lat", "south_north", "y"
dim_name_alias(:,3) = "soil_layers_stag", "horizons", "z"
dim_name_alias(:,3) = "soil_layers_stag", "horizons", "z", "sum_z"
/
&Input_Fields
filenames(1) = "./test/BD1.nc"
filenames(2) = "./test/SAND1.nc"
filenames(3) = "./test/CLAY1.nc"
varnames(1) = "Bd_mu"
varnames(2) = "SAND"
varnames(3) = "CLAY"
/
&Data_Arrays
names(1) = "Bd_mu"
from_file(1) = "./test/BD.nc"
to_file(1) = .false.
names(2) = "SAND"
from_file(2) = "./test/SAND.nc"
to_file(2) = .false.
names(3) = "CLAY"
from_file(3) = "./test/CLAY.nc"
to_file(3) = .false.
names(4) = "CLAY_deeper"
from_file(4) = "./test/CLAY_deeper.nc"
to_file(4) = .false.
&Eff_Params
eff_param_names(1) = "smcmax"
!transfer_funcs(1) = "linear_3arg_thresh"
transfer_funcs(1) = "linear_2arg"
!input_field_names(1:3,1) = "SAND", "CLAY", "Bd_mu"
input_field_names(1:2,1) = "SAND", "CLAY"
names(5) = "smcmax"
!transfer_funcs(5) = "linear_3arg_thresh"
transfer_funcs(5) = "linear_2arg"
!from_data_arrays(1:3,5) = "SAND", "CLAY", "Bd_mu"
from_data_arrays(1:2,5) = "SAND", "CLAY"
! this is C-ordered dimension
param_dim_names(1:3,1) = "soil_layers_stag", "south_north", "west_east"
upscale_ops(1:3,1) = -1.0, -1.0, -1.0
!limits(1:2,1) = 0.01, 1.0
eff_param_names(2) = "dksat"
transfer_funcs(2) = "linear_2arg_exp"
input_field_names(1:2,2) = "SAND", "CLAY"
param_dim_names(1:3,2) = "soil_layers_stag", "south_north", "west_east"
upscale_ops(1:3,2) = -1.0, -1.0, -1.0
!limits(1,2) = 0.0
eff_param_names(3) = "bexp"
transfer_funcs(3) = "linear_2arg"
input_field_names(1:2,3) = "SAND", "CLAY"
param_dim_names(1:3,3) = "soil_layers_stag", "south_north", "west_east"
upscale_ops(1:3,3) = -1.0, -1.0, -1.0
!limits(1,3) = 0.0
target_dim_names(1:3,5) = "soil_layers_stag", "south_north", "west_east"
upscale_ops(1:3,5) = "-1.0", "-1.0", "-1.0"
!limits(1:2,5) = 0.01, 1.0
names(6) = "dksat"
transfer_funcs(6) = "linear_2arg_exp"
from_data_arrays(1:2,6) = "SAND", "CLAY"
target_dim_names(1:3,6) = "soil_layers_stag", "south_north", "west_east"
upscale_ops(1:3,6) = "-1.0", "-1.0", "-1.0"
!limits(1,6) = 0.0
names(7) = "bexp"
transfer_funcs(7) = "linear_2arg"
from_data_arrays(1:2,7) = "SAND", "CLAY"
target_dim_names(1:3,7) = "soil_layers_stag", "south_north", "west_east"
upscale_ops(1:3,7) = "0.0", "0.0", "0.0"
!limits(1,7) = 0.0
! omitting transfer function and concatenating input fields
names(8) = "clay_concat"
from_data_arrays(1:2,8) = "CLAY", "CLAY_deeper"
target_dim_names(1:3,8) = "sum_z", "south_north", "west_east"
upscale_ops(1:3,8) = "sum", "1.0", "1.0"
! omitting upscale function and transfering to derived fields
names(9) = "derived_field"
transfer_funcs(9) = "linear_2arg"
from_data_arrays(1:2,9) = "smcmax", "bexp"
/
&Dimensions
dim_names(1:4) = "soil_layers_stag", "horizons", "south_north", "west_east"
dim_reference(1:4) = "end", "end", "center", "center"
dim_names(1:5) = "soil_layers_stag", "horizons", "south_north", "west_east", "sum_z"
dim_reference(1:5) = "end", "end", "center", "center", "end"
!dim_start(3:4) = 43.1875, -100.3125
!dim_count(3:4) = 7, 7
dim_step(3:4) = 0.125, 0.125
dim_step(3:4) = 0.03125, 0.03125
dim_bound(1:2) = 0.0, 0.0
dim_vector(:,1) = 0.1, 0.3, 0.6, 1.0
!target_dim_files(1:4) = "target_dim.nc", "target_dim.nc", "target_dim.nc", "target_dim.nc"
dim_bound(5) = 0.0,
dim_vector(:,5) = 4.0
!target_dim_files(1:4) = target_dim.nc, target_dim.nc, target_dim.nc, target_dim.nc
/
......@@ -2,28 +2,33 @@
! global_parameter([lower_bound, upper_bound, value], index_transfer_func_parameter, index_effective_parameter)
! Zacharias porosity
!global_parameter(1:3,1,1) = 0.6462, 0.9506, 0.76 ! PTF_zach_lower_constant
!global_parameter(1:3,2,1) = 0.0001, 0.0029, 0.0009 ! PTF_zach_lower_clay
!global_parameter(1:3,3,1) =-0.3727, -0.1871, -0.264 ! PTF_zach_lower_Db
!global_parameter(1:3,4,1) = 0.5358, 1.1232, 0.89 ! PTF_zach_higher_constant
!global_parameter(1:3,5,1) =-0.0055, 0.0049, -0.001 ! PTF_zach_higher_clay
!global_parameter(1:3,6,1) =-0.5513, -0.0913, -0.324 ! PTF_zach_higher_Db
!global_parameter(1:3,7,1) = 66.4, 66.6, 66.5 ! PTF_zach_thres
!global_parameter(1:3,1,5) = 0.6462, 0.9506, 0.76 ! PTF_zach_lower_constant
!global_parameter(1:3,2,5) = 0.0001, 0.0029, 0.0009 ! PTF_zach_lower_clay
!global_parameter(1:3,3,5) =-0.3727, -0.1871, -0.264 ! PTF_zach_lower_Db
!global_parameter(1:3,4,5) = 0.5358, 1.1232, 0.89 ! PTF_zach_higher_constant
!global_parameter(1:3,5,5) =-0.0055, 0.0049, -0.001 ! PTF_zach_higher_clay
!global_parameter(1:3,6,5) =-0.5513, -0.0913, -0.324 ! PTF_zach_higher_Db
!global_parameter(1:3,7,5) = 66.4, 66.6, 66.5 ! PTF_zach_thres
! Cosby porosity
global_parameter(1:3,1,1) = 42.27, 58.73, 50.5 ! PTF_porosity_constant
global_parameter(1:3,2,1) = -0.149, -0.135, -0.142 ! PTF_porosity_sand
global_parameter(1:3,3,1) =-0.1175, 0.0435, -0.037 ! PTF_porosity_clay
global_parameter(1:3,1,5) = 42.27, 58.73, 50.5 ! PTF_porosity_constant
global_parameter(1:3,2,5) = -0.149, -0.135, -0.142 ! PTF_porosity_sand
global_parameter(1:3,3,5) =-0.1175, 0.0435, -0.037 ! PTF_porosity_clay
! Cosby Ks
global_parameter(1:3,1,2) = -1.03, -0.17, -0.6 ! PTF_Ks_constant
global_parameter(1:3,2,2) = 0.0060, 0.0260, 0.0126 ! PTF_Ks_sand
global_parameter(1:3,3,2) =-0.0030, -0.0130, -0.0064 ! PTF_Ks_clay
global_parameter(1:3,4,2) = 0.000007055555, 0.000007055555, 0.000007055555 ! PTF_Ks_curveSlope conversion factor
global_parameter(1:3,1,6) = -1.03, -0.17, -0.6 ! PTF_Ks_constant
global_parameter(1:3,2,6) = 0.0060, 0.0260, 0.0126 ! PTF_Ks_sand
global_parameter(1:3,3,6) =-0.0030, -0.0130, -0.0064 ! PTF_Ks_clay
global_parameter(1:3,4,6) = 0.000007055555, 0.000007055555, 0.000007055555 ! PTF_Ks_curveSlope conversion factor
! inch/h to m/s
! Cosby b
global_parameter(1:3,1,3) = 2.18, 4.02, 3.1 ! PTF_b_constant
global_parameter(1:3,2,3) = -0.01, 0.004, -0.003 ! PTF_b_sand
global_parameter(1:3,3,3) =-0.0335, 0.0549, 0.0157 ! PTF_b_clay
global_parameter(1:3,1,7) = 2.18, 4.02, 3.1 ! PTF_b_constant
global_parameter(1:3,2,7) = -0.01, 0.004, -0.003 ! PTF_b_sand
global_parameter(1:3,3,7) =-0.0335, 0.0549, 0.0157 ! PTF_b_clay
! Cosby b
global_parameter(1:3,1,9) = 2.18, 4.02, 3.1 ! PTF_b_constant
global_parameter(1:3,2,9) = -0.01, 0.004, -0.003 ! PTF_b_sand
global_parameter(1:3,3,9) =-0.0335, 0.0549, 0.0157 ! PTF_b_clay
/
! namelist using pedotransfer function of Gupta and Larson 1979 (WRR)
! doi: 10.1029/WR015i006p01633
!
! sand, silt, and clay are defined as follows:
! sand = 2 - 0.05 mm, silt = 0.05 - 0.002 mm, and clay < 0.002 mm
!
&Input_Fields
filenames(1) = "./test/sand.nc"
filenames(2) = "./test/silt.nc"
filenames(3) = "./test/clay.nc"
filenames(4) = "./test/om.nc"
filenames(5) = "./test/bd.nc"
varnames(1) = "sand"
varnames(2) = "silt"
varnames(3) = "clay"
varnames(4) = "om"
varnames(5) = "bd"
/
&Eff_Params
eff_param_names(1) = "theta"
transfer_funcs(1) = "linear_5arg_no_intercept"
input_field_names(1:5,1) = "sand", "silt", "clay", "om", "bd"
! ! this is C-ordered dimension
! param_dim_names(1:3,1) = "soil_layers_stag", "south_north", "west_east"
upscale_ops(1:3,1) = -1.0, -1.0, -1.0
!limits(1:2,1) = 0.01, 1.0
/
&global_parameters
! global_parameter([lower_bound, upper_bound, value], index_transfer_func_parameter, index_effective_parameter)
! water content at given matric potential kPa
! matric potential of -0.04
global_parameter(1:3,1,1) = 7.053e-3, 7.053e-3, 7.053e-3 ! PTF_watercontent_sand
global_parameter(1:3,2,1) = 10.242e-3, 10.242e-3, 10.242e-3 ! PTF_watercontent_silt
global_parameter(1:3,3,1) = 10.070e-3, 10.070e-3, 10.070e-3 ! PTF_watercontent_clay
global_parameter(1:3,4,1) = 6.333e-3, 6.333e-3, 6.333e-3 ! PTF_watercontent_organic_matter
global_parameter(1:3,5,1) = -32.120e-3, -32.120e-3, -32.120e-3 ! PTF_watercontent_bulk_density
! ! matric potential of -0.07
! global_parameter(1:3,1,1) = 5.678e-3, 5.678e-3, 5.678e-3 ! PTF_watercontent_sand
! global_parameter(1:3,2,1) = 9.228e-3, 9.228e-3, 9.228e-3 ! PTF_watercontent_silt
! global_parameter(1:3,3,1) = 9.135e-3, 9.135e-3, 9.135e-3 ! PTF_watercontent_clay
! global_parameter(1:3,4,1) = 6.103e-3, 6.103e-3, 6.103e-3 ! PTF_watercontent_organic_matter
! global_parameter(1:3,5,1) = -26.960e-3, -26.960e-3, -26.960e-3 ! PTF_watercontent_bulk_density
! ! matric potential of -0.10
! global_parameter(1:3,1,1) = 5.018e-3, 5.018e-3, 5.018e-3 ! PTF_watercontent_sand
! global_parameter(1:3,2,1) = 8.548e-3, 8.548e-3, 8.548e-3 ! PTF_watercontent_silt
! global_parameter(1:3,3,1) = 8.833e-3, 8.833e-3, 8.833e-3 ! PTF_watercontent_clay
! global_parameter(1:3,4,1) = 4.966e-3, 4.966e-3, 4.966e-3 ! PTF_watercontent_organic_matter
! global_parameter(1:3,5,1) = -24.230e-3, -24.230e-3, -24.230e-3 ! PTF_watercontent_bulk_density
! ! matric potential of -0.20
! global_parameter(1:3,1,1) = 3.890e-e, 3.890e-e, 3.890e-e ! PTF_watercontent_sand
! global_parameter(1:3,2,1) = 7.066e-3, 7.066e-3, 7.066e-3 ! PTF_watercontent_silt
! global_parameter(1:3,3,1) = 8.408e-3, 8.408e-3, 8.408e-3 ! PTF_watercontent_clay
! global_parameter(1:3,4,1) = 2.817e-3, 2.817e-3, 2.817e-3 ! PTF_watercontent_organic_matter
! global_parameter(1:3,5,1) = -18.780e-3, -18.780e-3, -18.780e-3 ! PTF_watercontent_bulk_density
! ! matric potential of -0.33
! global_parameter(1:3,1,1) = 3.075e-3, 3.075e-3, 3.075e-3 ! PTF_watercontent_sand
! global_parameter(1:3,2,1) = 5.886e-3, 5.886e-3, 5.886e-3 ! PTF_watercontent_silt
! global_parameter(1:3,3,1) = 8.039e-3, 8.039e-3, 8.039e-3 ! PTF_watercontent_clay
! global_parameter(1:3,4,1) = 2.208e-3, 2.208e-3, 2.208e-3 ! PTF_watercontent_organic_matter
! global_parameter(1:3,5,1) = -14.340e-3, -14.340e-3, -14.340e-3 ! PTF_watercontent_bulk_density
! ! matric potential of -0.60
! global_parameter(1:3,1,1) = 2.181e-3, 2.181e-3, 2.181e-3 ! PTF_watercontent_sand
! global_parameter(1:3,2,1) = 4.557e-3, 4.557e-3, 4.557e-3 ! PTF_watercontent_silt
! global_parameter(1:3,3,1) = 7.557e-3, 7.557e-3, 7.557e-3 ! PTF_watercontent_clay
! global_parameter(1:3,4,1) = 2.191e-3, 2.191e-3, 2.191e-3 ! PTF_watercontent_organic_matter
! global_parameter(1:3,5,1) = -9.276e-3, -9.276e-3, -9.276e-3 ! PTF_watercontent_bulk_density
! ! matric potential of -1.00
! global_parameter(1:3,1,1) = 1.563e-3, 1.563e-3, 1.563e-3 ! PTF_watercontent_sand
! global_parameter(1:3,2,1) = 3.620e-3, 3.620e-3, 3.620e-3 ! PTF_watercontent_silt
! global_parameter(1:3,3,1) = 7.154e-3, 7.154e-3, 7.154e-3 ! PTF_watercontent_clay
! global_parameter(1:3,4,1) = 2.388e-3, 2.388e-3, 2.388e-3 ! PTF_watercontent_organic_matter
! global_parameter(1:3,5,1) = -5.754e-3, -5.754e-3, -5.754e-3 ! PTF_watercontent_bulk_density
! ! matric potential of -2.00
! global_parameter(1:3,1,1) = 0.932e-3, 0.932e-3, 0.932e-3 ! PTF_watercontent_sand
! global_parameter(1:3,2,1) = 2.643e-3, 2.643e-3, 2.643e-3 ! PTF_watercontent_silt
! global_parameter(1:3,3,1) = 6.636e-3, 6.636e-3, 6.636e-3 ! PTF_watercontent_clay
! global_parameter(1:3,4,1) = 2.717e-3, 2.717e-3, 2.717e-3 ! PTF_watercontent_organic_matter
! global_parameter(1:3,5,1) = -2.214e-3, -2.214e-3, -2.214e-3 ! PTF_watercontent_bulk_density
! ! matric potential of -4.00
! global_parameter(1:3,1,1) = 0.483e-3, 0.483e-3, 0.483e-3 ! PTF_watercontent_sand
! global_parameter(1:3,2,1) = 1.943e-3, 1.943e-3, 1.943e-3 ! PTF_watercontent_silt
! global_parameter(1:3,3,1) = 6.128e-3, 6.128e-3, 6.128e-3 ! PTF_watercontent_clay
! global_parameter(1:3,4,1) = 2.925e-3, 2.925e-3, 2.925e-3 ! PTF_watercontent_organic_matter
! global_parameter(1:3,5,1) = -0.204e-3, -0.204e-3, -0.204e-3 ! PTF_watercontent_bulk_density
! ! matric potential of -7.00
! global_parameter(1:3,1,1) = 0.214e-3, 0.214e-3, 0.214e-3 ! PTF_watercontent_sand
! global_parameter(1:3,2,1) = 1.538e-3, 1.538e-3, 1.538e-3 ! PTF_watercontent_silt
! global_parameter(1:3,3,1) = 5.908e-3, 5.908e-3, 5.908e-3 ! PTF_watercontent_clay
! global_parameter(1:3,4,1) = 2.855e-3, 2.855e-3, 2.855e-3 ! PTF_watercontent_organic_matter
! global_parameter(1:3,5,1) = 1.530e-3, 1.530e-3, 1.530e-3 ! PTF_watercontent_bulk_density
! ! matric potential of -10.00
! global_parameter(1:3,1,1) = 0.076e-3, 0.076e-3, 0.076e-3 ! PTF_watercontent_sand