def create_connections(self):
# Route snow melt to surface
cmf.SimpleTindexSnowMelt(self.cell.snow,
self.cell.surfacewater,
rate=7)
# Infiltration
cmf.SimpleInfiltration(self.soil, self.cell.surfacewater, W0=0.8)
# Route infiltration / saturation excess to outlet
cmf.WaterBalanceFlux(self.cell.surfacewater, self.outlet)
# Parameterize soil water capacity
self.soil.soil.porosity = 0.2
C = self.soil.get_capacity()
# Parameterize water stress function
self.cell.set_uptakestress(cmf.VolumeStress(0.2 * C, 0 * C))
cmf.TurcET(self.soil, self.cell.transpiration)
# Route water from soil to gw
cmf.PowerLawConnection(self.soil,
self.gw,
Q0=self.mm_to_m3(50),
V0=0.5 * C,
beta=4)
# Route water from gw to outlet
cmf.LinearStorageConnection(self.gw,
self.outlet,
residencetime=20,
residual=0 * C)
def create_project(W0=0.9):
"""
Creates the cmf project with a single cell and one layer
:return:
"""
p = cmf.project()
c = p.NewCell(0, 0, 0, 1000, with_surfacewater=True)
l = c.add_layer(0.1)
c.set_rainfall(l.get_capacity())
si_con = cmf.SimpleInfiltration(l, c.surfacewater, W0)
return p, si_con
def create_connections(self, p: Parameters):
"""
Creates the connections and parameterizes the storages of the model
"""
# Infiltration
cmf.SimpleInfiltration(self.soil, self.cell.surfacewater, W0=p.infiltration_w0)
# Route infiltration / saturation excess to outlet
cmf.waterbalance_connection(self.cell.surfacewater, self.outlet)
capacity = self.set_soil_capacity(p)
cmf.timeseriesETpot(self.soil, self.cell.transpiration, self.data.ETpot)
# Parameterize infiltration capacity
self.soil.soil.Ksat = p.infiltration_capacity / 1000
# Parameterize water stress function
self.cell.set_uptakestress(cmf.VolumeStress(
p.ETV1 * capacity, 0)
)
def setparameters(self,
tr_soil_GW = 12.36870481,
tr_soil_fulda = 12.,
tr_surf = 3.560855356,
tr_GW_l = 829.7188064,
tr_GW_u_fulda = 270.05035,
tr_GW_u_GW_l = 270.,
tr_fulda = 2.264612944,
V0_soil = 280.0850875,
beta_soil_GW=1.158865311,
beta_fulda = 1.1,
ETV1=2.575261852,
fETV0=0.014808919,
meltrate = 4.464735097,
snow_melt_temp = 4.51938545,
# Qd_max = 0.250552812,
# TW_threshold = 10.,
LAI = 2.992013336,
CanopyClosure = 5.,
Ksat = 0.02
): # this list has to be identical with the one above
"""
sets the parameters, all parameterized connections will be created anew
"""
# Get all definitions from init method
p = self.project
c = p[0]
outlet = self.outlet
fulda = self.fulda
# trinkwasser = self.trinkwasser
# Adjustment of the evapotranspiration
c.set_uptakestress(cmf.VolumeStress(ETV1,ETV1 * fETV0))
# Flux from the surfaces to the river
cmf.kinematic_wave(c.surfacewater,fulda,tr_surf)
# flux from surfaces to the soil (infiltration)
cmf.SimpleInfiltration(c.layers[0], c.surfacewater)
# change the saturated conductivity of the soil
c.layers[0].soil.Ksat = Ksat
# Flux from soil to river (interflow)
cmf.kinematic_wave(c.layers[0],fulda,tr_soil_fulda/V0_soil, V0 = V0_soil)
# flux from the soil to the upper groundwater (percolation)
cmf.kinematic_wave(c.layers[0], c.layers[1],tr_soil_GW, exponent=beta_soil_GW)
# flux from the upper groundwater to the river (baseflow)
cmf.kinematic_wave(c.layers[1], fulda, tr_GW_u_fulda)
# flux from upper to lower groundwater (percolation)
cmf.kinematic_wave(c.layers[1], c.layers[2],tr_GW_u_GW_l)#, exponent=beta_GW_u_GW_l)
# flux from the lower groundwater to river (baseflow)
cmf.kinematic_wave(c.layers[2], fulda, tr_GW_l)
# Flux from the lower groundwater to the drinking water outlet
# the fourths argument is the amount that is now allowed to be slurped
# # out of the lower groundwater
# cmf.TechnicalFlux(c.layers[2],trinkwasser,Qd_max,TW_threshold,cmf.day)
#
# # Flux from drinking water to the river
# cmf.waterbalance_connection(trinkwasser, fulda)
# flux from the river to the outlet
cmf.kinematic_wave(fulda, outlet, tr_fulda, exponent = beta_fulda)
# set snowmelt temperature
cmf.Weather.set_snow_threshold(snow_melt_temp)
# Snowmelt at the surfaces
snowmelt_surf = cmf.SimpleTindexSnowMelt(c.snow,c.surfacewater,c,rate=meltrate)
# Splits the rainfall in interzeption and throughfall
cmf.Rainfall(c.canopy,c, False, True)
cmf.Rainfall(c.surfacewater,c, True, False)
# Makes a overflow for the interception storage
cmf.RutterInterception(c.canopy,c.surfacewater,c)
# Transpiration on the plants is added
cmf.CanopyStorageEvaporation(c.canopy,c.evaporation,c)
# Sets the parameters for the interception
c.vegetation.LAI= LAI
# Defines how much throughfall there is (in %)
c.vegetation.CanopyClosure = CanopyClosure