HydPy 6.2.dev1__cp313-cp313-win_amd64.whl

hydpy/cythons/autogen/c_evap_aet_hbv96.pyx

@@ -0,0 +1,584 @@
+#!python
+# distutils: define_macros=NPY_NO_DEPRECATED_API=NPY_1_7_API_VERSION
+# cython: language_level=3
+# cython: cpow=True
+# cython: boundscheck=False
+# cython: wraparound=False
+# cython: initializedcheck=False
+# cython: cdivision=True
+from typing import Optional
+import numpy
+cimport numpy
+from libc.math cimport exp, fabs, log, sin, cos, tan, tanh, asin, acos, atan, isnan, isinf
+from libc.math cimport NAN as nan
+from libc.math cimport INFINITY as inf
+import cython
+from cpython.mem cimport PyMem_Malloc
+from cpython.mem cimport PyMem_Realloc
+from cpython.mem cimport PyMem_Free
+from hydpy.cythons.autogen cimport configutils
+from hydpy.cythons.autogen cimport interfaceutils
+from hydpy.cythons.autogen cimport interputils
+from hydpy.cythons.autogen import pointerutils
+from hydpy.cythons.autogen cimport pointerutils
+from hydpy.cythons.autogen cimport quadutils
+from hydpy.cythons.autogen cimport rootutils
+from hydpy.cythons.autogen cimport smoothutils
+from hydpy.cythons.autogen cimport masterinterface
+
+
+cdef void do_nothing(Model model)  noexcept nogil:
+    pass
+
+cpdef get_wrapper():
+    cdef CallbackWrapper wrapper = CallbackWrapper()
+    wrapper.callback = do_nothing
+    return wrapper
+
+@cython.final
+cdef class Parameters:
+    pass
+@cython.final
+cdef class ControlParameters:
+    pass
+@cython.final
+cdef class Sequences:
+    pass
+@cython.final
+cdef class FactorSequences:
+    cpdef inline void load_data(self, numpy.int64_t idx)  noexcept nogil:
+        cdef numpy.int64_t jdx0
+        cdef numpy.int64_t k
+        if self._airtemperature_diskflag_reading:
+            k = 0
+            for jdx0 in range(self._airtemperature_length_0):
+                self.airtemperature[jdx0] = self._airtemperature_ncarray[k]
+                k += 1
+        elif self._airtemperature_ramflag:
+            for jdx0 in range(self._airtemperature_length_0):
+                self.airtemperature[jdx0] = self._airtemperature_array[idx, jdx0]
+        if self._interceptedwater_diskflag_reading:
+            k = 0
+            for jdx0 in range(self._interceptedwater_length_0):
+                self.interceptedwater[jdx0] = self._interceptedwater_ncarray[k]
+                k += 1
+        elif self._interceptedwater_ramflag:
+            for jdx0 in range(self._interceptedwater_length_0):
+                self.interceptedwater[jdx0] = self._interceptedwater_array[idx, jdx0]
+        if self._soilwater_diskflag_reading:
+            k = 0
+            for jdx0 in range(self._soilwater_length_0):
+                self.soilwater[jdx0] = self._soilwater_ncarray[k]
+                k += 1
+        elif self._soilwater_ramflag:
+            for jdx0 in range(self._soilwater_length_0):
+                self.soilwater[jdx0] = self._soilwater_array[idx, jdx0]
+        if self._snowcover_diskflag_reading:
+            k = 0
+            for jdx0 in range(self._snowcover_length_0):
+                self.snowcover[jdx0] = self._snowcover_ncarray[k]
+                k += 1
+        elif self._snowcover_ramflag:
+            for jdx0 in range(self._snowcover_length_0):
+                self.snowcover[jdx0] = self._snowcover_array[idx, jdx0]
+    cpdef inline void save_data(self, numpy.int64_t idx)  noexcept nogil:
+        cdef numpy.int64_t jdx0
+        cdef numpy.int64_t k
+        if self._airtemperature_diskflag_writing:
+            k = 0
+            for jdx0 in range(self._airtemperature_length_0):
+                self._airtemperature_ncarray[k] = self.airtemperature[jdx0]
+                k += 1
+        if self._airtemperature_ramflag:
+            for jdx0 in range(self._airtemperature_length_0):
+                self._airtemperature_array[idx, jdx0] = self.airtemperature[jdx0]
+        if self._interceptedwater_diskflag_writing:
+            k = 0
+            for jdx0 in range(self._interceptedwater_length_0):
+                self._interceptedwater_ncarray[k] = self.interceptedwater[jdx0]
+                k += 1
+        if self._interceptedwater_ramflag:
+            for jdx0 in range(self._interceptedwater_length_0):
+                self._interceptedwater_array[idx, jdx0] = self.interceptedwater[jdx0]
+        if self._soilwater_diskflag_writing:
+            k = 0
+            for jdx0 in range(self._soilwater_length_0):
+                self._soilwater_ncarray[k] = self.soilwater[jdx0]
+                k += 1
+        if self._soilwater_ramflag:
+            for jdx0 in range(self._soilwater_length_0):
+                self._soilwater_array[idx, jdx0] = self.soilwater[jdx0]
+        if self._snowcover_diskflag_writing:
+            k = 0
+            for jdx0 in range(self._snowcover_length_0):
+                self._snowcover_ncarray[k] = self.snowcover[jdx0]
+                k += 1
+        if self._snowcover_ramflag:
+            for jdx0 in range(self._snowcover_length_0):
+                self._snowcover_array[idx, jdx0] = self.snowcover[jdx0]
+    cpdef inline set_pointeroutput(self, str name, pointerutils.PDouble value):
+        pass
+    cpdef inline void update_outputs(self) noexcept nogil:
+        pass
+@cython.final
+cdef class FluxSequences:
+    cpdef inline void load_data(self, numpy.int64_t idx)  noexcept nogil:
+        cdef numpy.int64_t jdx0
+        cdef numpy.int64_t k
+        if self._potentialinterceptionevaporation_diskflag_reading:
+            k = 0
+            for jdx0 in range(self._potentialinterceptionevaporation_length_0):
+                self.potentialinterceptionevaporation[jdx0] = self._potentialinterceptionevaporation_ncarray[k]
+                k += 1
+        elif self._potentialinterceptionevaporation_ramflag:
+            for jdx0 in range(self._potentialinterceptionevaporation_length_0):
+                self.potentialinterceptionevaporation[jdx0] = self._potentialinterceptionevaporation_array[idx, jdx0]
+        if self._potentialsoilevapotranspiration_diskflag_reading:
+            k = 0
+            for jdx0 in range(self._potentialsoilevapotranspiration_length_0):
+                self.potentialsoilevapotranspiration[jdx0] = self._potentialsoilevapotranspiration_ncarray[k]
+                k += 1
+        elif self._potentialsoilevapotranspiration_ramflag:
+            for jdx0 in range(self._potentialsoilevapotranspiration_length_0):
+                self.potentialsoilevapotranspiration[jdx0] = self._potentialsoilevapotranspiration_array[idx, jdx0]
+        if self._potentialwaterevaporation_diskflag_reading:
+            k = 0
+            for jdx0 in range(self._potentialwaterevaporation_length_0):
+                self.potentialwaterevaporation[jdx0] = self._potentialwaterevaporation_ncarray[k]
+                k += 1
+        elif self._potentialwaterevaporation_ramflag:
+            for jdx0 in range(self._potentialwaterevaporation_length_0):
+                self.potentialwaterevaporation[jdx0] = self._potentialwaterevaporation_array[idx, jdx0]
+        if self._waterevaporation_diskflag_reading:
+            k = 0
+            for jdx0 in range(self._waterevaporation_length_0):
+                self.waterevaporation[jdx0] = self._waterevaporation_ncarray[k]
+                k += 1
+        elif self._waterevaporation_ramflag:
+            for jdx0 in range(self._waterevaporation_length_0):
+                self.waterevaporation[jdx0] = self._waterevaporation_array[idx, jdx0]
+        if self._interceptionevaporation_diskflag_reading:
+            k = 0
+            for jdx0 in range(self._interceptionevaporation_length_0):
+                self.interceptionevaporation[jdx0] = self._interceptionevaporation_ncarray[k]
+                k += 1
+        elif self._interceptionevaporation_ramflag:
+            for jdx0 in range(self._interceptionevaporation_length_0):
+                self.interceptionevaporation[jdx0] = self._interceptionevaporation_array[idx, jdx0]
+        if self._soilevapotranspiration_diskflag_reading:
+            k = 0
+            for jdx0 in range(self._soilevapotranspiration_length_0):
+                self.soilevapotranspiration[jdx0] = self._soilevapotranspiration_ncarray[k]
+                k += 1
+        elif self._soilevapotranspiration_ramflag:
+            for jdx0 in range(self._soilevapotranspiration_length_0):
+                self.soilevapotranspiration[jdx0] = self._soilevapotranspiration_array[idx, jdx0]
+    cpdef inline void save_data(self, numpy.int64_t idx)  noexcept nogil:
+        cdef numpy.int64_t jdx0
+        cdef numpy.int64_t k
+        if self._potentialinterceptionevaporation_diskflag_writing:
+            k = 0
+            for jdx0 in range(self._potentialinterceptionevaporation_length_0):
+                self._potentialinterceptionevaporation_ncarray[k] = self.potentialinterceptionevaporation[jdx0]
+                k += 1
+        if self._potentialinterceptionevaporation_ramflag:
+            for jdx0 in range(self._potentialinterceptionevaporation_length_0):
+                self._potentialinterceptionevaporation_array[idx, jdx0] = self.potentialinterceptionevaporation[jdx0]
+        if self._potentialsoilevapotranspiration_diskflag_writing:
+            k = 0
+            for jdx0 in range(self._potentialsoilevapotranspiration_length_0):
+                self._potentialsoilevapotranspiration_ncarray[k] = self.potentialsoilevapotranspiration[jdx0]
+                k += 1
+        if self._potentialsoilevapotranspiration_ramflag:
+            for jdx0 in range(self._potentialsoilevapotranspiration_length_0):
+                self._potentialsoilevapotranspiration_array[idx, jdx0] = self.potentialsoilevapotranspiration[jdx0]
+        if self._potentialwaterevaporation_diskflag_writing:
+            k = 0
+            for jdx0 in range(self._potentialwaterevaporation_length_0):
+                self._potentialwaterevaporation_ncarray[k] = self.potentialwaterevaporation[jdx0]
+                k += 1
+        if self._potentialwaterevaporation_ramflag:
+            for jdx0 in range(self._potentialwaterevaporation_length_0):
+                self._potentialwaterevaporation_array[idx, jdx0] = self.potentialwaterevaporation[jdx0]
+        if self._waterevaporation_diskflag_writing:
+            k = 0
+            for jdx0 in range(self._waterevaporation_length_0):
+                self._waterevaporation_ncarray[k] = self.waterevaporation[jdx0]
+                k += 1
+        if self._waterevaporation_ramflag:
+            for jdx0 in range(self._waterevaporation_length_0):
+                self._waterevaporation_array[idx, jdx0] = self.waterevaporation[jdx0]
+        if self._interceptionevaporation_diskflag_writing:
+            k = 0
+            for jdx0 in range(self._interceptionevaporation_length_0):
+                self._interceptionevaporation_ncarray[k] = self.interceptionevaporation[jdx0]
+                k += 1
+        if self._interceptionevaporation_ramflag:
+            for jdx0 in range(self._interceptionevaporation_length_0):
+                self._interceptionevaporation_array[idx, jdx0] = self.interceptionevaporation[jdx0]
+        if self._soilevapotranspiration_diskflag_writing:
+            k = 0
+            for jdx0 in range(self._soilevapotranspiration_length_0):
+                self._soilevapotranspiration_ncarray[k] = self.soilevapotranspiration[jdx0]
+                k += 1
+        if self._soilevapotranspiration_ramflag:
+            for jdx0 in range(self._soilevapotranspiration_length_0):
+                self._soilevapotranspiration_array[idx, jdx0] = self.soilevapotranspiration[jdx0]
+    cpdef inline set_pointeroutput(self, str name, pointerutils.PDouble value):
+        pass
+    cpdef inline void update_outputs(self) noexcept nogil:
+        pass
+@cython.final
+cdef class Model(masterinterface.MasterInterface):
+    def __init__(self):
+        super().__init__()
+        self.intercmodel = None
+        self.intercmodel_is_mainmodel = False
+        self.petmodel = None
+        self.petmodel_is_mainmodel = False
+        self.snowcovermodel = None
+        self.snowcovermodel_is_mainmodel = False
+        self.soilwatermodel = None
+        self.soilwatermodel_is_mainmodel = False
+        self.tempmodel = None
+        self.tempmodel_is_mainmodel = False
+    def get_intercmodel(self) -> masterinterface.MasterInterface | None:
+        return self.intercmodel
+    def set_intercmodel(self, intercmodel: masterinterface.MasterInterface | None) -> None:
+        self.intercmodel = intercmodel
+    def get_petmodel(self) -> masterinterface.MasterInterface | None:
+        return self.petmodel
+    def set_petmodel(self, petmodel: masterinterface.MasterInterface | None) -> None:
+        self.petmodel = petmodel
+    def get_snowcovermodel(self) -> masterinterface.MasterInterface | None:
+        return self.snowcovermodel
+    def set_snowcovermodel(self, snowcovermodel: masterinterface.MasterInterface | None) -> None:
+        self.snowcovermodel = snowcovermodel
+    def get_soilwatermodel(self) -> masterinterface.MasterInterface | None:
+        return self.soilwatermodel
+    def set_soilwatermodel(self, soilwatermodel: masterinterface.MasterInterface | None) -> None:
+        self.soilwatermodel = soilwatermodel
+    def get_tempmodel(self) -> masterinterface.MasterInterface | None:
+        return self.tempmodel
+    def set_tempmodel(self, tempmodel: masterinterface.MasterInterface | None) -> None:
+        self.tempmodel = tempmodel
+    cpdef inline void simulate(self, numpy.int64_t idx)  noexcept nogil:
+        self.idx_sim = idx
+        self.reset_reuseflags()
+        self.load_data(idx)
+        self.run()
+        self.update_outputs()
+    cpdef void reset_reuseflags(self) noexcept nogil:
+        if (self.intercmodel is not None) and not self.intercmodel_is_mainmodel:
+            self.intercmodel.reset_reuseflags()
+        if (self.petmodel is not None) and not self.petmodel_is_mainmodel:
+            self.petmodel.reset_reuseflags()
+        if (self.snowcovermodel is not None) and not self.snowcovermodel_is_mainmodel:
+            self.snowcovermodel.reset_reuseflags()
+        if (self.soilwatermodel is not None) and not self.soilwatermodel_is_mainmodel:
+            self.soilwatermodel.reset_reuseflags()
+        if (self.tempmodel is not None) and not self.tempmodel_is_mainmodel:
+            self.tempmodel.reset_reuseflags()
+    cpdef void load_data(self, numpy.int64_t idx) noexcept nogil:
+        self.idx_sim = idx
+        if (self.intercmodel is not None) and not self.intercmodel_is_mainmodel:
+            self.intercmodel.load_data(idx)
+        if (self.petmodel is not None) and not self.petmodel_is_mainmodel:
+            self.petmodel.load_data(idx)
+        if (self.snowcovermodel is not None) and not self.snowcovermodel_is_mainmodel:
+            self.snowcovermodel.load_data(idx)
+        if (self.soilwatermodel is not None) and not self.soilwatermodel_is_mainmodel:
+            self.soilwatermodel.load_data(idx)
+        if (self.tempmodel is not None) and not self.tempmodel_is_mainmodel:
+            self.tempmodel.load_data(idx)
+    cpdef void save_data(self, numpy.int64_t idx) noexcept nogil:
+        self.idx_sim = idx
+        self.sequences.factors.save_data(idx)
+        self.sequences.fluxes.save_data(idx)
+        if (self.intercmodel is not None) and not self.intercmodel_is_mainmodel:
+            self.intercmodel.save_data(idx)
+        if (self.petmodel is not None) and not self.petmodel_is_mainmodel:
+            self.petmodel.save_data(idx)
+        if (self.snowcovermodel is not None) and not self.snowcovermodel_is_mainmodel:
+            self.snowcovermodel.save_data(idx)
+        if (self.soilwatermodel is not None) and not self.soilwatermodel_is_mainmodel:
+            self.soilwatermodel.save_data(idx)
+        if (self.tempmodel is not None) and not self.tempmodel_is_mainmodel:
+            self.tempmodel.save_data(idx)
+    cpdef void new2old(self) noexcept nogil:
+        if (self.intercmodel is not None) and not self.intercmodel_is_mainmodel:
+            self.intercmodel.new2old()
+        if (self.petmodel is not None) and not self.petmodel_is_mainmodel:
+            self.petmodel.new2old()
+        if (self.snowcovermodel is not None) and not self.snowcovermodel_is_mainmodel:
+            self.snowcovermodel.new2old()
+        if (self.soilwatermodel is not None) and not self.soilwatermodel_is_mainmodel:
+            self.soilwatermodel.new2old()
+        if (self.tempmodel is not None) and not self.tempmodel_is_mainmodel:
+            self.tempmodel.new2old()
+    cpdef inline void run(self) noexcept nogil:
+        self.determine_interceptionevaporation_v1()
+        self.determine_soilevapotranspiration_v1()
+        self.determine_waterevaporation_v1()
+    cpdef inline void update_inlets(self) noexcept nogil:
+        pass
+    cpdef inline void update_outlets(self) noexcept nogil:
+        pass
+    cpdef inline void update_receivers(self, numpy.int64_t idx) noexcept nogil:
+        self.idx_sim = idx
+        pass
+    cpdef inline void update_senders(self, numpy.int64_t idx) noexcept nogil:
+        self.idx_sim = idx
+        pass
+    cpdef inline void update_outputs(self) noexcept nogil:
+        pass
+    cpdef double get_waterevaporation_v1(self, numpy.int64_t k) noexcept nogil:
+        return self.sequences.fluxes.waterevaporation[k]
+    cpdef double get_interceptionevaporation_v1(self, numpy.int64_t k) noexcept nogil:
+        return self.sequences.fluxes.interceptionevaporation[k]
+    cpdef double get_soilevapotranspiration_v1(self, numpy.int64_t k) noexcept nogil:
+        return self.sequences.fluxes.soilevapotranspiration[k]
+    cpdef inline void calc_airtemperature_v1(self) noexcept nogil:
+        if self.tempmodel_typeid == 1:
+            self.calc_airtemperature_tempmodel_v1(                (<masterinterface.MasterInterface>self.tempmodel)            )
+        elif self.tempmodel_typeid == 2:
+            self.calc_airtemperature_tempmodel_v2(                (<masterinterface.MasterInterface>self.tempmodel)            )
+    cpdef inline void calc_potentialinterceptionevaporation_petmodel_v1(self, masterinterface.MasterInterface submodel) noexcept nogil:
+        cdef numpy.int64_t k
+        submodel.determine_potentialevapotranspiration()
+        for k in range(self.parameters.control.nmbhru):
+            self.sequences.fluxes.potentialinterceptionevaporation[k] = (                submodel.get_potentialevapotranspiration(k)            )
+    cpdef inline void calc_potentialinterceptionevaporation_petmodel_v2(self, masterinterface.MasterInterface submodel) noexcept nogil:
+        cdef numpy.int64_t k
+        submodel.determine_potentialinterceptionevaporation()
+        for k in range(self.parameters.control.nmbhru):
+            self.sequences.fluxes.potentialinterceptionevaporation[k] = (                submodel.get_potentialinterceptionevaporation(k)            )
+    cpdef inline void calc_potentialinterceptionevaporation_v3(self) noexcept nogil:
+        if self.petmodel_typeid == 1:
+            self.calc_potentialinterceptionevaporation_petmodel_v1(                (<masterinterface.MasterInterface>self.petmodel)            )
+        elif self.petmodel_typeid == 2:
+            self.calc_potentialinterceptionevaporation_petmodel_v2(                (<masterinterface.MasterInterface>self.petmodel)            )
+    cpdef inline void calc_potentialwaterevaporation_petmodel_v1(self, masterinterface.MasterInterface submodel) noexcept nogil:
+        cdef numpy.int64_t k
+        for k in range(self.parameters.control.nmbhru):
+            self.sequences.fluxes.potentialwaterevaporation[k] = submodel.get_potentialevapotranspiration(                k            )
+    cpdef inline void calc_potentialwaterevaporation_petmodel_v2(self, masterinterface.MasterInterface submodel) noexcept nogil:
+        cdef numpy.int64_t k
+        submodel.determine_potentialwaterevaporation()
+        for k in range(self.parameters.control.nmbhru):
+            self.sequences.fluxes.potentialwaterevaporation[k] = submodel.get_potentialwaterevaporation(k)
+    cpdef inline void calc_potentialwaterevaporation_v1(self) noexcept nogil:
+        if self.petmodel_typeid == 1:
+            self.calc_potentialwaterevaporation_petmodel_v1(                (<masterinterface.MasterInterface>self.petmodel)            )
+        elif self.petmodel_typeid == 2:
+            self.calc_potentialwaterevaporation_petmodel_v2(                (<masterinterface.MasterInterface>self.petmodel)            )
+    cpdef inline void calc_waterevaporation_v1(self) noexcept nogil:
+        cdef numpy.int64_t k
+        for k in range(self.parameters.control.nmbhru):
+            if self.parameters.control.water[k] and self.sequences.factors.airtemperature[k] > self.parameters.control.temperaturethresholdice[k]:
+                self.sequences.fluxes.waterevaporation[k] = self.sequences.fluxes.potentialwaterevaporation[k]
+            else:
+                self.sequences.fluxes.waterevaporation[k] = 0.0
+    cpdef inline void calc_interceptedwater_v1(self) noexcept nogil:
+        if self.intercmodel_typeid == 1:
+            self.calc_interceptedwater_intercmodel_v1(                (<masterinterface.MasterInterface>self.intercmodel)            )
+    cpdef inline void calc_interceptionevaporation_v1(self) noexcept nogil:
+        cdef numpy.int64_t k
+        for k in range(self.parameters.control.nmbhru):
+            if self.parameters.control.interception[k]:
+                self.sequences.fluxes.interceptionevaporation[k] = min(                    self.sequences.fluxes.potentialinterceptionevaporation[k], self.sequences.factors.interceptedwater[k]                )
+            else:
+                self.sequences.fluxes.interceptionevaporation[k] = 0.0
+    cpdef inline void calc_soilwater_v1(self) noexcept nogil:
+        if self.soilwatermodel_typeid == 1:
+            self.calc_soilwater_soilwatermodel_v1(                (<masterinterface.MasterInterface>self.soilwatermodel)            )
+    cpdef inline void calc_snowcover_v1(self) noexcept nogil:
+        if self.snowcovermodel_typeid == 1:
+            self.calc_snowcover_snowcovermodel_v1(                (<masterinterface.MasterInterface>self.snowcovermodel)            )
+    cpdef inline void calc_potentialsoilevapotranspiration_petmodel_v1(self, masterinterface.MasterInterface submodel) noexcept nogil:
+        cdef numpy.int64_t k
+        for k in range(self.parameters.control.nmbhru):
+            self.sequences.fluxes.potentialsoilevapotranspiration[k] = (                submodel.get_potentialevapotranspiration(k)            )
+    cpdef inline void calc_potentialsoilevapotranspiration_petmodel_v2(self, masterinterface.MasterInterface submodel) noexcept nogil:
+        cdef numpy.int64_t k
+        submodel.determine_potentialsoilevapotranspiration()
+        for k in range(self.parameters.control.nmbhru):
+            self.sequences.fluxes.potentialsoilevapotranspiration[k] = (                submodel.get_potentialsoilevapotranspiration(k)            )
+    cpdef inline void calc_potentialsoilevapotranspiration_v2(self) noexcept nogil:
+        if self.petmodel_typeid == 1:
+            self.calc_potentialsoilevapotranspiration_petmodel_v1(                (<masterinterface.MasterInterface>self.petmodel)            )
+        elif self.petmodel_typeid == 2:
+            self.calc_potentialsoilevapotranspiration_petmodel_v2(                (<masterinterface.MasterInterface>self.petmodel)            )
+    cpdef inline void calc_soilevapotranspiration_v1(self) noexcept nogil:
+        cdef double thresh
+        cdef numpy.int64_t k
+        for k in range(self.parameters.control.nmbhru):
+            if self.parameters.control.soil[k]:
+                self.sequences.fluxes.soilevapotranspiration[k] = self.sequences.fluxes.potentialsoilevapotranspiration[k]
+                if self.sequences.fluxes.soilevapotranspiration[k] > 0.0:
+                    if self.sequences.factors.soilwater[k] < 0.0:
+                        self.sequences.fluxes.soilevapotranspiration[k] = 0.0
+                    else:
+                        thresh = self.parameters.control.soilmoisturelimit[k] * self.parameters.control.maxsoilwater[k]
+                        if self.sequences.factors.soilwater[k] < thresh:
+                            self.sequences.fluxes.soilevapotranspiration[k] = self.sequences.fluxes.soilevapotranspiration[k] * (self.sequences.factors.soilwater[k] / thresh)
+            else:
+                self.sequences.fluxes.soilevapotranspiration[k] = 0.0
+    cpdef inline void update_soilevapotranspiration_v1(self) noexcept nogil:
+        cdef double excess
+        cdef double ei
+        cdef double ets
+        cdef double pet
+        cdef double pets
+        cdef double pei
+        cdef double r
+        cdef numpy.int64_t k
+        for k in range(self.parameters.control.nmbhru):
+            if self.parameters.control.soil[k]:
+                r = self.parameters.control.excessreduction[k]
+                pei = self.sequences.fluxes.potentialinterceptionevaporation[k]
+                pets = self.sequences.fluxes.potentialsoilevapotranspiration[k]
+                pet = r * pei + (1.0 - r) * (pei + pets) / 2.0
+                ets = self.sequences.fluxes.soilevapotranspiration[k]
+                ei = self.sequences.fluxes.interceptionevaporation[k]
+                excess = ets + ei - pet
+                if (pet >= 0.0) and (excess > 0.0):
+                    self.sequences.fluxes.soilevapotranspiration[k] = self.sequences.fluxes.soilevapotranspiration[k] - (r * excess)
+                elif (pet < 0.0) and (excess < 0.0):
+                    self.sequences.fluxes.soilevapotranspiration[k] = self.sequences.fluxes.soilevapotranspiration[k] - (r * excess)
+            else:
+                self.sequences.fluxes.soilevapotranspiration[k] = 0.0
+    cpdef inline void update_soilevapotranspiration_v2(self) noexcept nogil:
+        cdef numpy.int64_t k
+        for k in range(self.parameters.control.nmbhru):
+            if self.parameters.control.soil[k] and (self.sequences.factors.snowcover[k] < 1.0):
+                self.sequences.fluxes.soilevapotranspiration[k] = self.sequences.fluxes.soilevapotranspiration[k] * (1.0 - self.sequences.factors.snowcover[k])
+            else:
+                self.sequences.fluxes.soilevapotranspiration[k] = 0.0
+    cpdef inline void calc_airtemperature_tempmodel_v1(self, masterinterface.MasterInterface submodel) noexcept nogil:
+        cdef numpy.int64_t k
+        for k in range(self.parameters.control.nmbhru):
+            self.sequences.factors.airtemperature[k] = submodel.get_temperature(k)
+    cpdef inline void calc_airtemperature_tempmodel_v2(self, masterinterface.MasterInterface submodel) noexcept nogil:
+        cdef numpy.int64_t k
+        submodel.determine_temperature()
+        for k in range(self.parameters.control.nmbhru):
+            self.sequences.factors.airtemperature[k] = submodel.get_temperature(k)
+    cpdef inline void calc_interceptedwater_intercmodel_v1(self, masterinterface.MasterInterface submodel) noexcept nogil:
+        cdef numpy.int64_t k
+        for k in range(self.parameters.control.nmbhru):
+            self.sequences.factors.interceptedwater[k] = submodel.get_interceptedwater(k)
+    cpdef inline void calc_soilwater_soilwatermodel_v1(self, masterinterface.MasterInterface submodel) noexcept nogil:
+        cdef numpy.int64_t k
+        for k in range(self.parameters.control.nmbhru):
+            self.sequences.factors.soilwater[k] = submodel.get_soilwater(k)
+    cpdef inline void calc_snowcover_snowcovermodel_v1(self, masterinterface.MasterInterface submodel) noexcept nogil:
+        cdef numpy.int64_t k
+        for k in range(self.parameters.control.nmbhru):
+            self.sequences.factors.snowcover[k] = submodel.get_snowcover(k)
+    cpdef double get_waterevaporation(self, numpy.int64_t k) noexcept nogil:
+        return self.sequences.fluxes.waterevaporation[k]
+    cpdef double get_interceptionevaporation(self, numpy.int64_t k) noexcept nogil:
+        return self.sequences.fluxes.interceptionevaporation[k]
+    cpdef double get_soilevapotranspiration(self, numpy.int64_t k) noexcept nogil:
+        return self.sequences.fluxes.soilevapotranspiration[k]
+    cpdef inline void calc_airtemperature(self) noexcept nogil:
+        if self.tempmodel_typeid == 1:
+            self.calc_airtemperature_tempmodel_v1(                (<masterinterface.MasterInterface>self.tempmodel)            )
+        elif self.tempmodel_typeid == 2:
+            self.calc_airtemperature_tempmodel_v2(                (<masterinterface.MasterInterface>self.tempmodel)            )
+    cpdef inline void calc_potentialinterceptionevaporation(self) noexcept nogil:
+        if self.petmodel_typeid == 1:
+            self.calc_potentialinterceptionevaporation_petmodel_v1(                (<masterinterface.MasterInterface>self.petmodel)            )
+        elif self.petmodel_typeid == 2:
+            self.calc_potentialinterceptionevaporation_petmodel_v2(                (<masterinterface.MasterInterface>self.petmodel)            )
+    cpdef inline void calc_potentialwaterevaporation(self) noexcept nogil:
+        if self.petmodel_typeid == 1:
+            self.calc_potentialwaterevaporation_petmodel_v1(                (<masterinterface.MasterInterface>self.petmodel)            )
+        elif self.petmodel_typeid == 2:
+            self.calc_potentialwaterevaporation_petmodel_v2(                (<masterinterface.MasterInterface>self.petmodel)            )
+    cpdef inline void calc_waterevaporation(self) noexcept nogil:
+        cdef numpy.int64_t k
+        for k in range(self.parameters.control.nmbhru):
+            if self.parameters.control.water[k] and self.sequences.factors.airtemperature[k] > self.parameters.control.temperaturethresholdice[k]:
+                self.sequences.fluxes.waterevaporation[k] = self.sequences.fluxes.potentialwaterevaporation[k]
+            else:
+                self.sequences.fluxes.waterevaporation[k] = 0.0
+    cpdef inline void calc_interceptedwater(self) noexcept nogil:
+        if self.intercmodel_typeid == 1:
+            self.calc_interceptedwater_intercmodel_v1(                (<masterinterface.MasterInterface>self.intercmodel)            )
+    cpdef inline void calc_interceptionevaporation(self) noexcept nogil:
+        cdef numpy.int64_t k
+        for k in range(self.parameters.control.nmbhru):
+            if self.parameters.control.interception[k]:
+                self.sequences.fluxes.interceptionevaporation[k] = min(                    self.sequences.fluxes.potentialinterceptionevaporation[k], self.sequences.factors.interceptedwater[k]                )
+            else:
+                self.sequences.fluxes.interceptionevaporation[k] = 0.0
+    cpdef inline void calc_soilwater(self) noexcept nogil:
+        if self.soilwatermodel_typeid == 1:
+            self.calc_soilwater_soilwatermodel_v1(                (<masterinterface.MasterInterface>self.soilwatermodel)            )
+    cpdef inline void calc_snowcover(self) noexcept nogil:
+        if self.snowcovermodel_typeid == 1:
+            self.calc_snowcover_snowcovermodel_v1(                (<masterinterface.MasterInterface>self.snowcovermodel)            )
+    cpdef inline void calc_potentialsoilevapotranspiration(self) noexcept nogil:
+        if self.petmodel_typeid == 1:
+            self.calc_potentialsoilevapotranspiration_petmodel_v1(                (<masterinterface.MasterInterface>self.petmodel)            )
+        elif self.petmodel_typeid == 2:
+            self.calc_potentialsoilevapotranspiration_petmodel_v2(                (<masterinterface.MasterInterface>self.petmodel)            )
+    cpdef inline void calc_soilevapotranspiration(self) noexcept nogil:
+        cdef double thresh
+        cdef numpy.int64_t k
+        for k in range(self.parameters.control.nmbhru):
+            if self.parameters.control.soil[k]:
+                self.sequences.fluxes.soilevapotranspiration[k] = self.sequences.fluxes.potentialsoilevapotranspiration[k]
+                if self.sequences.fluxes.soilevapotranspiration[k] > 0.0:
+                    if self.sequences.factors.soilwater[k] < 0.0:
+                        self.sequences.fluxes.soilevapotranspiration[k] = 0.0
+                    else:
+                        thresh = self.parameters.control.soilmoisturelimit[k] * self.parameters.control.maxsoilwater[k]
+                        if self.sequences.factors.soilwater[k] < thresh:
+                            self.sequences.fluxes.soilevapotranspiration[k] = self.sequences.fluxes.soilevapotranspiration[k] * (self.sequences.factors.soilwater[k] / thresh)
+            else:
+                self.sequences.fluxes.soilevapotranspiration[k] = 0.0
+    cpdef inline void calc_interceptedwater_intercmodel(self, masterinterface.MasterInterface submodel) noexcept nogil:
+        cdef numpy.int64_t k
+        for k in range(self.parameters.control.nmbhru):
+            self.sequences.factors.interceptedwater[k] = submodel.get_interceptedwater(k)
+    cpdef inline void calc_soilwater_soilwatermodel(self, masterinterface.MasterInterface submodel) noexcept nogil:
+        cdef numpy.int64_t k
+        for k in range(self.parameters.control.nmbhru):
+            self.sequences.factors.soilwater[k] = submodel.get_soilwater(k)
+    cpdef inline void calc_snowcover_snowcovermodel(self, masterinterface.MasterInterface submodel) noexcept nogil:
+        cdef numpy.int64_t k
+        for k in range(self.parameters.control.nmbhru):
+            self.sequences.factors.snowcover[k] = submodel.get_snowcover(k)
+    cpdef void determine_interceptionevaporation_v1(self) noexcept nogil:
+        self.calc_potentialinterceptionevaporation_v3()
+        self.calc_interceptedwater_v1()
+        self.calc_interceptionevaporation_v1()
+    cpdef void determine_soilevapotranspiration_v1(self) noexcept nogil:
+        self.calc_soilwater_v1()
+        self.calc_snowcover_v1()
+        self.calc_potentialsoilevapotranspiration_v2()
+        self.calc_soilevapotranspiration_v1()
+        self.update_soilevapotranspiration_v1()
+        self.update_soilevapotranspiration_v2()
+    cpdef void determine_waterevaporation_v1(self) noexcept nogil:
+        self.calc_airtemperature_v1()
+        self.calc_potentialwaterevaporation_v1()
+        self.calc_waterevaporation_v1()
+    cpdef void determine_interceptionevaporation(self) noexcept nogil:
+        self.calc_potentialinterceptionevaporation_v3()
+        self.calc_interceptedwater_v1()
+        self.calc_interceptionevaporation_v1()
+    cpdef void determine_soilevapotranspiration(self) noexcept nogil:
+        self.calc_soilwater_v1()
+        self.calc_snowcover_v1()
+        self.calc_potentialsoilevapotranspiration_v2()
+        self.calc_soilevapotranspiration_v1()
+        self.update_soilevapotranspiration_v1()
+        self.update_soilevapotranspiration_v2()
+    cpdef void determine_waterevaporation(self) noexcept nogil:
+        self.calc_airtemperature_v1()
+        self.calc_potentialwaterevaporation_v1()
+        self.calc_waterevaporation_v1()