HydPy 6.2.dev1__cp313-cp313-win_amd64.whl

hydpy/cythons/autogen/c_meteo_sun_morsim.pyx

@@ -0,0 +1,519 @@
+#!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 DerivedParameters:
+    pass
+@cython.final
+cdef class FixedParameters:
+    pass
+@cython.final
+cdef class Sequences:
+    pass
+@cython.final
+cdef class InputSequences:
+    cpdef inline void load_data(self, numpy.int64_t idx)  noexcept nogil:
+        cdef numpy.int64_t k
+        if self._globalradiation_inputflag:
+            self.globalradiation = self._globalradiation_inputpointer[0]
+        elif self._globalradiation_diskflag_reading:
+            self.globalradiation = self._globalradiation_ncarray[0]
+        elif self._globalradiation_ramflag:
+            self.globalradiation = self._globalradiation_array[idx]
+    cpdef inline void save_data(self, numpy.int64_t idx)  noexcept nogil:
+        cdef numpy.int64_t k
+        if self._globalradiation_diskflag_writing:
+            self._globalradiation_ncarray[0] = self.globalradiation
+        if self._globalradiation_ramflag:
+            self._globalradiation_array[idx] = self.globalradiation
+    cpdef inline set_pointerinput(self, str name, pointerutils.PDouble value):
+        if name == "globalradiation":
+            self._globalradiation_inputpointer = value.p_value
+@cython.final
+cdef class FactorSequences:
+    cpdef inline void load_data(self, numpy.int64_t idx)  noexcept nogil:
+        cdef numpy.int64_t k
+        if self._earthsundistance_diskflag_reading:
+            self.earthsundistance = self._earthsundistance_ncarray[0]
+        elif self._earthsundistance_ramflag:
+            self.earthsundistance = self._earthsundistance_array[idx]
+        if self._solardeclination_diskflag_reading:
+            self.solardeclination = self._solardeclination_ncarray[0]
+        elif self._solardeclination_ramflag:
+            self.solardeclination = self._solardeclination_array[idx]
+        if self._timeofsunrise_diskflag_reading:
+            self.timeofsunrise = self._timeofsunrise_ncarray[0]
+        elif self._timeofsunrise_ramflag:
+            self.timeofsunrise = self._timeofsunrise_array[idx]
+        if self._timeofsunset_diskflag_reading:
+            self.timeofsunset = self._timeofsunset_ncarray[0]
+        elif self._timeofsunset_ramflag:
+            self.timeofsunset = self._timeofsunset_array[idx]
+        if self._possiblesunshineduration_diskflag_reading:
+            self.possiblesunshineduration = self._possiblesunshineduration_ncarray[0]
+        elif self._possiblesunshineduration_ramflag:
+            self.possiblesunshineduration = self._possiblesunshineduration_array[idx]
+        if self._dailypossiblesunshineduration_diskflag_reading:
+            self.dailypossiblesunshineduration = self._dailypossiblesunshineduration_ncarray[0]
+        elif self._dailypossiblesunshineduration_ramflag:
+            self.dailypossiblesunshineduration = self._dailypossiblesunshineduration_array[idx]
+        if self._unadjustedsunshineduration_diskflag_reading:
+            self.unadjustedsunshineduration = self._unadjustedsunshineduration_ncarray[0]
+        elif self._unadjustedsunshineduration_ramflag:
+            self.unadjustedsunshineduration = self._unadjustedsunshineduration_array[idx]
+        if self._sunshineduration_diskflag_reading:
+            self.sunshineduration = self._sunshineduration_ncarray[0]
+        elif self._sunshineduration_ramflag:
+            self.sunshineduration = self._sunshineduration_array[idx]
+        if self._dailysunshineduration_diskflag_reading:
+            self.dailysunshineduration = self._dailysunshineduration_ncarray[0]
+        elif self._dailysunshineduration_ramflag:
+            self.dailysunshineduration = self._dailysunshineduration_array[idx]
+        if self._portiondailyradiation_diskflag_reading:
+            self.portiondailyradiation = self._portiondailyradiation_ncarray[0]
+        elif self._portiondailyradiation_ramflag:
+            self.portiondailyradiation = self._portiondailyradiation_array[idx]
+    cpdef inline void save_data(self, numpy.int64_t idx)  noexcept nogil:
+        cdef numpy.int64_t k
+        if self._earthsundistance_diskflag_writing:
+            self._earthsundistance_ncarray[0] = self.earthsundistance
+        if self._earthsundistance_ramflag:
+            self._earthsundistance_array[idx] = self.earthsundistance
+        if self._solardeclination_diskflag_writing:
+            self._solardeclination_ncarray[0] = self.solardeclination
+        if self._solardeclination_ramflag:
+            self._solardeclination_array[idx] = self.solardeclination
+        if self._timeofsunrise_diskflag_writing:
+            self._timeofsunrise_ncarray[0] = self.timeofsunrise
+        if self._timeofsunrise_ramflag:
+            self._timeofsunrise_array[idx] = self.timeofsunrise
+        if self._timeofsunset_diskflag_writing:
+            self._timeofsunset_ncarray[0] = self.timeofsunset
+        if self._timeofsunset_ramflag:
+            self._timeofsunset_array[idx] = self.timeofsunset
+        if self._possiblesunshineduration_diskflag_writing:
+            self._possiblesunshineduration_ncarray[0] = self.possiblesunshineduration
+        if self._possiblesunshineduration_ramflag:
+            self._possiblesunshineduration_array[idx] = self.possiblesunshineduration
+        if self._dailypossiblesunshineduration_diskflag_writing:
+            self._dailypossiblesunshineduration_ncarray[0] = self.dailypossiblesunshineduration
+        if self._dailypossiblesunshineduration_ramflag:
+            self._dailypossiblesunshineduration_array[idx] = self.dailypossiblesunshineduration
+        if self._unadjustedsunshineduration_diskflag_writing:
+            self._unadjustedsunshineduration_ncarray[0] = self.unadjustedsunshineduration
+        if self._unadjustedsunshineduration_ramflag:
+            self._unadjustedsunshineduration_array[idx] = self.unadjustedsunshineduration
+        if self._sunshineduration_diskflag_writing:
+            self._sunshineduration_ncarray[0] = self.sunshineduration
+        if self._sunshineduration_ramflag:
+            self._sunshineduration_array[idx] = self.sunshineduration
+        if self._dailysunshineduration_diskflag_writing:
+            self._dailysunshineduration_ncarray[0] = self.dailysunshineduration
+        if self._dailysunshineduration_ramflag:
+            self._dailysunshineduration_array[idx] = self.dailysunshineduration
+        if self._portiondailyradiation_diskflag_writing:
+            self._portiondailyradiation_ncarray[0] = self.portiondailyradiation
+        if self._portiondailyradiation_ramflag:
+            self._portiondailyradiation_array[idx] = self.portiondailyradiation
+    cpdef inline set_pointeroutput(self, str name, pointerutils.PDouble value):
+        if name == "earthsundistance":
+            self._earthsundistance_outputpointer = value.p_value
+        if name == "solardeclination":
+            self._solardeclination_outputpointer = value.p_value
+        if name == "timeofsunrise":
+            self._timeofsunrise_outputpointer = value.p_value
+        if name == "timeofsunset":
+            self._timeofsunset_outputpointer = value.p_value
+        if name == "possiblesunshineduration":
+            self._possiblesunshineduration_outputpointer = value.p_value
+        if name == "dailypossiblesunshineduration":
+            self._dailypossiblesunshineduration_outputpointer = value.p_value
+        if name == "unadjustedsunshineduration":
+            self._unadjustedsunshineduration_outputpointer = value.p_value
+        if name == "sunshineduration":
+            self._sunshineduration_outputpointer = value.p_value
+        if name == "dailysunshineduration":
+            self._dailysunshineduration_outputpointer = value.p_value
+        if name == "portiondailyradiation":
+            self._portiondailyradiation_outputpointer = value.p_value
+    cpdef inline void update_outputs(self) noexcept nogil:
+        if self._earthsundistance_outputflag:
+            self._earthsundistance_outputpointer[0] = self.earthsundistance
+        if self._solardeclination_outputflag:
+            self._solardeclination_outputpointer[0] = self.solardeclination
+        if self._timeofsunrise_outputflag:
+            self._timeofsunrise_outputpointer[0] = self.timeofsunrise
+        if self._timeofsunset_outputflag:
+            self._timeofsunset_outputpointer[0] = self.timeofsunset
+        if self._possiblesunshineduration_outputflag:
+            self._possiblesunshineduration_outputpointer[0] = self.possiblesunshineduration
+        if self._dailypossiblesunshineduration_outputflag:
+            self._dailypossiblesunshineduration_outputpointer[0] = self.dailypossiblesunshineduration
+        if self._unadjustedsunshineduration_outputflag:
+            self._unadjustedsunshineduration_outputpointer[0] = self.unadjustedsunshineduration
+        if self._sunshineduration_outputflag:
+            self._sunshineduration_outputpointer[0] = self.sunshineduration
+        if self._dailysunshineduration_outputflag:
+            self._dailysunshineduration_outputpointer[0] = self.dailysunshineduration
+        if self._portiondailyradiation_outputflag:
+            self._portiondailyradiation_outputpointer[0] = self.portiondailyradiation
+@cython.final
+cdef class FluxSequences:
+    cpdef inline void load_data(self, numpy.int64_t idx)  noexcept nogil:
+        cdef numpy.int64_t k
+        if self._extraterrestrialradiation_diskflag_reading:
+            self.extraterrestrialradiation = self._extraterrestrialradiation_ncarray[0]
+        elif self._extraterrestrialradiation_ramflag:
+            self.extraterrestrialradiation = self._extraterrestrialradiation_array[idx]
+        if self._clearskysolarradiation_diskflag_reading:
+            self.clearskysolarradiation = self._clearskysolarradiation_ncarray[0]
+        elif self._clearskysolarradiation_ramflag:
+            self.clearskysolarradiation = self._clearskysolarradiation_array[idx]
+        if self._dailyglobalradiation_diskflag_reading:
+            self.dailyglobalradiation = self._dailyglobalradiation_ncarray[0]
+        elif self._dailyglobalradiation_ramflag:
+            self.dailyglobalradiation = self._dailyglobalradiation_array[idx]
+    cpdef inline void save_data(self, numpy.int64_t idx)  noexcept nogil:
+        cdef numpy.int64_t k
+        if self._extraterrestrialradiation_diskflag_writing:
+            self._extraterrestrialradiation_ncarray[0] = self.extraterrestrialradiation
+        if self._extraterrestrialradiation_ramflag:
+            self._extraterrestrialradiation_array[idx] = self.extraterrestrialradiation
+        if self._clearskysolarradiation_diskflag_writing:
+            self._clearskysolarradiation_ncarray[0] = self.clearskysolarradiation
+        if self._clearskysolarradiation_ramflag:
+            self._clearskysolarradiation_array[idx] = self.clearskysolarradiation
+        if self._dailyglobalradiation_diskflag_writing:
+            self._dailyglobalradiation_ncarray[0] = self.dailyglobalradiation
+        if self._dailyglobalradiation_ramflag:
+            self._dailyglobalradiation_array[idx] = self.dailyglobalradiation
+    cpdef inline set_pointeroutput(self, str name, pointerutils.PDouble value):
+        if name == "extraterrestrialradiation":
+            self._extraterrestrialradiation_outputpointer = value.p_value
+        if name == "clearskysolarradiation":
+            self._clearskysolarradiation_outputpointer = value.p_value
+        if name == "dailyglobalradiation":
+            self._dailyglobalradiation_outputpointer = value.p_value
+    cpdef inline void update_outputs(self) noexcept nogil:
+        if self._extraterrestrialradiation_outputflag:
+            self._extraterrestrialradiation_outputpointer[0] = self.extraterrestrialradiation
+        if self._clearskysolarradiation_outputflag:
+            self._clearskysolarradiation_outputpointer[0] = self.clearskysolarradiation
+        if self._dailyglobalradiation_outputflag:
+            self._dailyglobalradiation_outputpointer[0] = self.dailyglobalradiation
+@cython.final
+cdef class LogSequences:
+    pass
+@cython.final
+cdef class Model(masterinterface.MasterInterface):
+    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:
+        self.__hydpy_reuse_process_radiation_v1__ = False
+    cpdef void load_data(self, numpy.int64_t idx) noexcept nogil:
+        self.idx_sim = idx
+        self.sequences.inputs.load_data(idx)
+    cpdef void save_data(self, numpy.int64_t idx) noexcept nogil:
+        self.idx_sim = idx
+        self.sequences.inputs.save_data(idx)
+        self.sequences.factors.save_data(idx)
+        self.sequences.fluxes.save_data(idx)
+    cpdef inline void run(self) noexcept nogil:
+        self.calc_earthsundistance_v1()
+        self.calc_solardeclination_v2()
+        self.calc_timeofsunrise_timeofsunset_v1()
+        self.calc_dailypossiblesunshineduration_v1()
+        self.calc_possiblesunshineduration_v2()
+        self.update_loggedglobalradiation_v1()
+        self.calc_dailyglobalradiation_v2()
+        self.calc_extraterrestrialradiation_v2()
+        self.calc_clearskysolarradiation_v1()
+        self.calc_dailysunshineduration_v2()
+        self.calc_portiondailyradiation_v1()
+        self.adjust_clearskysolarradiation_v1()
+        self.calc_unadjustedsunshineduration_v1()
+        self.update_loggedunadjustedsunshineduration_v1()
+        self.calc_sunshineduration_v2()
+    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:
+        self.sequences.factors.update_outputs()
+        self.sequences.fluxes.update_outputs()
+    cpdef inline void calc_earthsundistance_v1(self) noexcept nogil:
+        self.sequences.factors.earthsundistance = 1.0 + 0.033 * cos(            2 * self.parameters.fixed.pi / 366.0 * (self.parameters.derived.doy[self.idx_sim] + 1)        )
+    cpdef inline void calc_solardeclination_v2(self) noexcept nogil:
+        self.sequences.factors.solardeclination = 0.41 * cos(            2.0 * self.parameters.fixed.pi * (self.parameters.derived.doy[self.idx_sim] - 171.0) / 365.0        )
+    cpdef inline void calc_timeofsunrise_timeofsunset_v1(self) noexcept nogil:
+        cdef double d_dt
+        self.sequences.factors.timeofsunrise = (12.0 / self.parameters.fixed.pi) * acos(            tan(self.sequences.factors.solardeclination) * tan(self.parameters.derived.latituderad)            + 0.0145            / cos(self.sequences.factors.solardeclination)            / cos(self.parameters.derived.latituderad)        )
+        self.sequences.factors.timeofsunset = 24.0 - self.sequences.factors.timeofsunrise
+        d_dt = (self.parameters.derived.utclongitude - self.parameters.control.longitude) * 4.0 / 60.0
+        self.sequences.factors.timeofsunrise = self.sequences.factors.timeofsunrise + (d_dt)
+        self.sequences.factors.timeofsunset = self.sequences.factors.timeofsunset + (d_dt)
+    cpdef inline void calc_dailypossiblesunshineduration_v1(self) noexcept nogil:
+        self.sequences.factors.dailypossiblesunshineduration = self.sequences.factors.timeofsunset - self.sequences.factors.timeofsunrise
+    cpdef inline void calc_possiblesunshineduration_v2(self) noexcept nogil:
+        cdef double d_t1
+        cdef double d_t0
+        cdef double d_stc
+        d_stc = self.parameters.derived.sct[self.idx_sim]
+        d_t0 = max((d_stc - self.parameters.derived.hours / 2.0), self.sequences.factors.timeofsunrise)
+        d_t1 = min((d_stc + self.parameters.derived.hours / 2.0), self.sequences.factors.timeofsunset)
+        self.sequences.factors.possiblesunshineduration = max(d_t1 - d_t0, 0.0)
+    cpdef inline void update_loggedglobalradiation_v1(self) noexcept nogil:
+        cdef numpy.int64_t idx
+        for idx in range(self.parameters.derived.nmblogentries - 1, 0, -1):
+            self.sequences.logs.loggedglobalradiation[idx] = self.sequences.logs.loggedglobalradiation[idx - 1]
+        self.sequences.logs.loggedglobalradiation[0] = self.sequences.inputs.globalradiation
+    cpdef inline void calc_dailyglobalradiation_v2(self) noexcept nogil:
+        cdef numpy.int64_t idx
+        self.sequences.fluxes.dailyglobalradiation = 0.0
+        for idx in range(self.parameters.derived.nmblogentries):
+            self.sequences.fluxes.dailyglobalradiation = self.sequences.fluxes.dailyglobalradiation + (self.sequences.logs.loggedglobalradiation[idx])
+        self.sequences.fluxes.dailyglobalradiation = self.sequences.fluxes.dailyglobalradiation / (self.parameters.derived.nmblogentries)
+    cpdef inline void calc_extraterrestrialradiation_v2(self) noexcept nogil:
+        cdef double d_sunsethourangle
+        d_sunsethourangle = (self.sequences.factors.timeofsunset - self.sequences.factors.timeofsunrise) * self.parameters.fixed.pi / 24.0
+        self.sequences.fluxes.extraterrestrialradiation = (            self.parameters.fixed.solarconstant * self.sequences.factors.earthsundistance / self.parameters.fixed.pi        ) * (            d_sunsethourangle            * sin(self.sequences.factors.solardeclination)            * sin(self.parameters.derived.latituderad)            + cos(self.sequences.factors.solardeclination)            * cos(self.parameters.derived.latituderad)            * sin(d_sunsethourangle)        )
+    cpdef inline void calc_clearskysolarradiation_v1(self) noexcept nogil:
+        cdef numpy.int64_t idx
+        idx = self.parameters.derived.moy[self.idx_sim]
+        self.sequences.fluxes.clearskysolarradiation = self.sequences.fluxes.extraterrestrialradiation * (            self.parameters.control.angstromconstant[idx] + self.parameters.control.angstromfactor[idx]        )
+    cpdef inline void calc_dailysunshineduration_v2(self) noexcept nogil:
+        self.sequences.factors.dailysunshineduration = self.return_sunshineduration_v1(            self.sequences.fluxes.dailyglobalradiation,            self.sequences.fluxes.extraterrestrialradiation,            self.sequences.factors.dailypossiblesunshineduration,        )
+    cpdef inline void calc_portiondailyradiation_v1(self) noexcept nogil:
+        cdef double d_temp
+        cdef double d_p
+        cdef double d_tlp
+        cdef double d_dt
+        cdef numpy.int64_t i
+        cdef double d_fac
+        d_fac = 2.0 * self.parameters.fixed.pi / 360.0
+        self.sequences.factors.portiondailyradiation = 0.0
+        for i in range(2):
+            if i:
+                d_dt = self.parameters.derived.hours / 2.0
+            else:
+                d_dt = -self.parameters.derived.hours / 2.0
+            d_tlp = (100.0 * d_fac) * (                (self.parameters.derived.sct[self.idx_sim] + d_dt - self.sequences.factors.timeofsunrise)                / (self.sequences.factors.timeofsunset - self.sequences.factors.timeofsunrise)            )
+            if d_tlp <= 0.0:
+                d_p = 0.0
+            elif d_tlp < 100.0 * d_fac:
+                d_p = 50.0 - 50.0 * cos(1.8 * d_tlp)
+                d_temp = 3.4 * sin(3.6 * d_tlp) ** 2
+                if d_tlp <= 50.0 * d_fac:
+                    d_p = d_p - (d_temp)
+                else:
+                    d_p = d_p + (d_temp)
+            else:
+                d_p = 100.0
+            if i:
+                self.sequences.factors.portiondailyradiation = self.sequences.factors.portiondailyradiation + (d_p)
+            else:
+                self.sequences.factors.portiondailyradiation = self.sequences.factors.portiondailyradiation - (d_p)
+    cpdef inline void adjust_clearskysolarradiation_v1(self) noexcept nogil:
+        self.sequences.fluxes.clearskysolarradiation = self.sequences.fluxes.clearskysolarradiation * ((            self.parameters.derived.nmblogentries * self.sequences.factors.portiondailyradiation / 100.0        ))
+    cpdef inline void calc_unadjustedsunshineduration_v1(self) noexcept nogil:
+        self.sequences.factors.unadjustedsunshineduration = self.return_sunshineduration_v1(            self.sequences.inputs.globalradiation,            self.sequences.fluxes.extraterrestrialradiation            * self.parameters.derived.nmblogentries            * self.sequences.factors.portiondailyradiation            / 100.0,            self.sequences.factors.possiblesunshineduration,        )
+    cpdef inline void update_loggedunadjustedsunshineduration_v1(self) noexcept nogil:
+        cdef numpy.int64_t idx
+        for idx in range(self.parameters.derived.nmblogentries - 1, 0, -1):
+            self.sequences.logs.loggedunadjustedsunshineduration[idx] = (                self.sequences.logs.loggedunadjustedsunshineduration[idx - 1]            )
+        self.sequences.logs.loggedunadjustedsunshineduration[0] = self.sequences.factors.unadjustedsunshineduration
+    cpdef inline void calc_sunshineduration_v2(self) noexcept nogil:
+        cdef numpy.int64_t idx
+        cdef double d_denom
+        cdef double d_nom
+        d_nom = self.sequences.factors.unadjustedsunshineduration * self.sequences.factors.dailysunshineduration
+        if d_nom == 0.0:
+            self.sequences.factors.sunshineduration = 0.0
+        else:
+            d_denom = 0.0
+            for idx in range(self.parameters.derived.nmblogentries):
+                d_denom = d_denom + (self.sequences.logs.loggedunadjustedsunshineduration[idx])
+            self.sequences.factors.sunshineduration = min(d_nom / d_denom, self.sequences.factors.possiblesunshineduration)
+    cpdef void process_radiation_v1(self) noexcept nogil:
+        if not self.__hydpy_reuse_process_radiation_v1__:
+            self.run()
+            self.__hydpy_reuse_process_radiation_v1__ = True
+    cpdef double get_possiblesunshineduration_v1(self) noexcept nogil:
+        return self.sequences.factors.possiblesunshineduration
+    cpdef double get_sunshineduration_v1(self) noexcept nogil:
+        return self.sequences.factors.sunshineduration
+    cpdef double get_clearskysolarradiation_v1(self) noexcept nogil:
+        return self.sequences.fluxes.clearskysolarradiation
+    cpdef double get_globalradiation_v2(self) noexcept nogil:
+        return self.sequences.inputs.globalradiation
+    cpdef inline double return_sunshineduration_v1(self, double globalradiation, double extraterrestrialradiation, double possiblesunshineduration) noexcept nogil:
+        cdef double d_sd
+        cdef numpy.int64_t idx
+        if extraterrestrialradiation <= 0.0:
+            return possiblesunshineduration
+        idx = self.parameters.derived.moy[self.idx_sim]
+        d_sd = (possiblesunshineduration / self.parameters.control.angstromfactor[idx]) * (            globalradiation / extraterrestrialradiation - self.parameters.control.angstromconstant[idx]        )
+        return min(max(d_sd, 0.0), possiblesunshineduration)
+    cpdef inline void calc_earthsundistance(self) noexcept nogil:
+        self.sequences.factors.earthsundistance = 1.0 + 0.033 * cos(            2 * self.parameters.fixed.pi / 366.0 * (self.parameters.derived.doy[self.idx_sim] + 1)        )
+    cpdef inline void calc_solardeclination(self) noexcept nogil:
+        self.sequences.factors.solardeclination = 0.41 * cos(            2.0 * self.parameters.fixed.pi * (self.parameters.derived.doy[self.idx_sim] - 171.0) / 365.0        )
+    cpdef inline void calc_timeofsunrise_timeofsunset(self) noexcept nogil:
+        cdef double d_dt
+        self.sequences.factors.timeofsunrise = (12.0 / self.parameters.fixed.pi) * acos(            tan(self.sequences.factors.solardeclination) * tan(self.parameters.derived.latituderad)            + 0.0145            / cos(self.sequences.factors.solardeclination)            / cos(self.parameters.derived.latituderad)        )
+        self.sequences.factors.timeofsunset = 24.0 - self.sequences.factors.timeofsunrise
+        d_dt = (self.parameters.derived.utclongitude - self.parameters.control.longitude) * 4.0 / 60.0
+        self.sequences.factors.timeofsunrise = self.sequences.factors.timeofsunrise + (d_dt)
+        self.sequences.factors.timeofsunset = self.sequences.factors.timeofsunset + (d_dt)
+    cpdef inline void calc_dailypossiblesunshineduration(self) noexcept nogil:
+        self.sequences.factors.dailypossiblesunshineduration = self.sequences.factors.timeofsunset - self.sequences.factors.timeofsunrise
+    cpdef inline void calc_possiblesunshineduration(self) noexcept nogil:
+        cdef double d_t1
+        cdef double d_t0
+        cdef double d_stc
+        d_stc = self.parameters.derived.sct[self.idx_sim]
+        d_t0 = max((d_stc - self.parameters.derived.hours / 2.0), self.sequences.factors.timeofsunrise)
+        d_t1 = min((d_stc + self.parameters.derived.hours / 2.0), self.sequences.factors.timeofsunset)
+        self.sequences.factors.possiblesunshineduration = max(d_t1 - d_t0, 0.0)
+    cpdef inline void update_loggedglobalradiation(self) noexcept nogil:
+        cdef numpy.int64_t idx
+        for idx in range(self.parameters.derived.nmblogentries - 1, 0, -1):
+            self.sequences.logs.loggedglobalradiation[idx] = self.sequences.logs.loggedglobalradiation[idx - 1]
+        self.sequences.logs.loggedglobalradiation[0] = self.sequences.inputs.globalradiation
+    cpdef inline void calc_dailyglobalradiation(self) noexcept nogil:
+        cdef numpy.int64_t idx
+        self.sequences.fluxes.dailyglobalradiation = 0.0
+        for idx in range(self.parameters.derived.nmblogentries):
+            self.sequences.fluxes.dailyglobalradiation = self.sequences.fluxes.dailyglobalradiation + (self.sequences.logs.loggedglobalradiation[idx])
+        self.sequences.fluxes.dailyglobalradiation = self.sequences.fluxes.dailyglobalradiation / (self.parameters.derived.nmblogentries)
+    cpdef inline void calc_extraterrestrialradiation(self) noexcept nogil:
+        cdef double d_sunsethourangle
+        d_sunsethourangle = (self.sequences.factors.timeofsunset - self.sequences.factors.timeofsunrise) * self.parameters.fixed.pi / 24.0
+        self.sequences.fluxes.extraterrestrialradiation = (            self.parameters.fixed.solarconstant * self.sequences.factors.earthsundistance / self.parameters.fixed.pi        ) * (            d_sunsethourangle            * sin(self.sequences.factors.solardeclination)            * sin(self.parameters.derived.latituderad)            + cos(self.sequences.factors.solardeclination)            * cos(self.parameters.derived.latituderad)            * sin(d_sunsethourangle)        )
+    cpdef inline void calc_clearskysolarradiation(self) noexcept nogil:
+        cdef numpy.int64_t idx
+        idx = self.parameters.derived.moy[self.idx_sim]
+        self.sequences.fluxes.clearskysolarradiation = self.sequences.fluxes.extraterrestrialradiation * (            self.parameters.control.angstromconstant[idx] + self.parameters.control.angstromfactor[idx]        )
+    cpdef inline void calc_dailysunshineduration(self) noexcept nogil:
+        self.sequences.factors.dailysunshineduration = self.return_sunshineduration_v1(            self.sequences.fluxes.dailyglobalradiation,            self.sequences.fluxes.extraterrestrialradiation,            self.sequences.factors.dailypossiblesunshineduration,        )
+    cpdef inline void calc_portiondailyradiation(self) noexcept nogil:
+        cdef double d_temp
+        cdef double d_p
+        cdef double d_tlp
+        cdef double d_dt
+        cdef numpy.int64_t i
+        cdef double d_fac
+        d_fac = 2.0 * self.parameters.fixed.pi / 360.0
+        self.sequences.factors.portiondailyradiation = 0.0
+        for i in range(2):
+            if i:
+                d_dt = self.parameters.derived.hours / 2.0
+            else:
+                d_dt = -self.parameters.derived.hours / 2.0
+            d_tlp = (100.0 * d_fac) * (                (self.parameters.derived.sct[self.idx_sim] + d_dt - self.sequences.factors.timeofsunrise)                / (self.sequences.factors.timeofsunset - self.sequences.factors.timeofsunrise)            )
+            if d_tlp <= 0.0:
+                d_p = 0.0
+            elif d_tlp < 100.0 * d_fac:
+                d_p = 50.0 - 50.0 * cos(1.8 * d_tlp)
+                d_temp = 3.4 * sin(3.6 * d_tlp) ** 2
+                if d_tlp <= 50.0 * d_fac:
+                    d_p = d_p - (d_temp)
+                else:
+                    d_p = d_p + (d_temp)
+            else:
+                d_p = 100.0
+            if i:
+                self.sequences.factors.portiondailyradiation = self.sequences.factors.portiondailyradiation + (d_p)
+            else:
+                self.sequences.factors.portiondailyradiation = self.sequences.factors.portiondailyradiation - (d_p)
+    cpdef inline void adjust_clearskysolarradiation(self) noexcept nogil:
+        self.sequences.fluxes.clearskysolarradiation = self.sequences.fluxes.clearskysolarradiation * ((            self.parameters.derived.nmblogentries * self.sequences.factors.portiondailyradiation / 100.0        ))
+    cpdef inline void calc_unadjustedsunshineduration(self) noexcept nogil:
+        self.sequences.factors.unadjustedsunshineduration = self.return_sunshineduration_v1(            self.sequences.inputs.globalradiation,            self.sequences.fluxes.extraterrestrialradiation            * self.parameters.derived.nmblogentries            * self.sequences.factors.portiondailyradiation            / 100.0,            self.sequences.factors.possiblesunshineduration,        )
+    cpdef inline void update_loggedunadjustedsunshineduration(self) noexcept nogil:
+        cdef numpy.int64_t idx
+        for idx in range(self.parameters.derived.nmblogentries - 1, 0, -1):
+            self.sequences.logs.loggedunadjustedsunshineduration[idx] = (                self.sequences.logs.loggedunadjustedsunshineduration[idx - 1]            )
+        self.sequences.logs.loggedunadjustedsunshineduration[0] = self.sequences.factors.unadjustedsunshineduration
+    cpdef inline void calc_sunshineduration(self) noexcept nogil:
+        cdef numpy.int64_t idx
+        cdef double d_denom
+        cdef double d_nom
+        d_nom = self.sequences.factors.unadjustedsunshineduration * self.sequences.factors.dailysunshineduration
+        if d_nom == 0.0:
+            self.sequences.factors.sunshineduration = 0.0
+        else:
+            d_denom = 0.0
+            for idx in range(self.parameters.derived.nmblogentries):
+                d_denom = d_denom + (self.sequences.logs.loggedunadjustedsunshineduration[idx])
+            self.sequences.factors.sunshineduration = min(d_nom / d_denom, self.sequences.factors.possiblesunshineduration)
+    cpdef void process_radiation(self) noexcept nogil:
+        if not self.__hydpy_reuse_process_radiation_v1__:
+            self.run()
+            self.__hydpy_reuse_process_radiation_v1__ = True
+    cpdef double get_possiblesunshineduration(self) noexcept nogil:
+        return self.sequences.factors.possiblesunshineduration
+    cpdef double get_sunshineduration(self) noexcept nogil:
+        return self.sequences.factors.sunshineduration
+    cpdef double get_clearskysolarradiation(self) noexcept nogil:
+        return self.sequences.fluxes.clearskysolarradiation
+    cpdef double get_globalradiation(self) noexcept nogil:
+        return self.sequences.inputs.globalradiation
+    cpdef inline double return_sunshineduration(self, double globalradiation, double extraterrestrialradiation, double possiblesunshineduration) noexcept nogil:
+        cdef double d_sd
+        cdef numpy.int64_t idx
+        if extraterrestrialradiation <= 0.0:
+            return possiblesunshineduration
+        idx = self.parameters.derived.moy[self.idx_sim]
+        d_sd = (possiblesunshineduration / self.parameters.control.angstromfactor[idx]) * (            globalradiation / extraterrestrialradiation - self.parameters.control.angstromconstant[idx]        )
+        return min(max(d_sd, 0.0), possiblesunshineduration)

hydpy/cythons/autogen/c_meteo_temp_io.pxd

@@ -0,0 +1,112 @@
+#!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
+ctypedef void (*CallbackType) (Model)  noexcept nogil
+cdef class CallbackWrapper:
+    cdef CallbackType callback
+@cython.final
+cdef class Parameters:
+    cdef public ControlParameters control
+    cdef public DerivedParameters derived
+@cython.final
+cdef class ControlParameters:
+    cdef public numpy.int64_t nmbhru
+    cdef public double[:] hruarea
+    cdef public double[:] temperatureaddend
+@cython.final
+cdef class DerivedParameters:
+    cdef public double[:] hruareafraction
+@cython.final
+cdef class Sequences:
+    cdef public InputSequences inputs
+    cdef public FactorSequences factors
+@cython.final
+cdef class InputSequences:
+    cdef public double temperature
+    cdef public numpy.int64_t _temperature_ndim
+    cdef public numpy.int64_t _temperature_length
+    cdef public bint _temperature_ramflag
+    cdef public double[:] _temperature_array
+    cdef public bint _temperature_diskflag_reading
+    cdef public bint _temperature_diskflag_writing
+    cdef public double[:] _temperature_ncarray
+    cdef public bint _temperature_inputflag
+    cdef double *_temperature_inputpointer
+    cpdef inline void load_data(self, numpy.int64_t idx)  noexcept nogil
+    cpdef inline void save_data(self, numpy.int64_t idx)  noexcept nogil
+    cpdef inline set_pointerinput(self, str name, pointerutils.PDouble value)
+@cython.final
+cdef class FactorSequences:
+    cdef public double[:] temperature
+    cdef public numpy.int64_t _temperature_ndim
+    cdef public numpy.int64_t _temperature_length
+    cdef public numpy.int64_t _temperature_length_0
+    cdef public bint _temperature_ramflag
+    cdef public double[:,:] _temperature_array
+    cdef public bint _temperature_diskflag_reading
+    cdef public bint _temperature_diskflag_writing
+    cdef public double[:] _temperature_ncarray
+    cdef public double meantemperature
+    cdef public numpy.int64_t _meantemperature_ndim
+    cdef public numpy.int64_t _meantemperature_length
+    cdef public bint _meantemperature_ramflag
+    cdef public double[:] _meantemperature_array
+    cdef public bint _meantemperature_diskflag_reading
+    cdef public bint _meantemperature_diskflag_writing
+    cdef public double[:] _meantemperature_ncarray
+    cdef public bint _meantemperature_outputflag
+    cdef double *_meantemperature_outputpointer
+    cpdef inline void load_data(self, numpy.int64_t idx)  noexcept nogil
+    cpdef inline void save_data(self, numpy.int64_t idx)  noexcept nogil
+    cpdef inline set_pointeroutput(self, str name, pointerutils.PDouble value)
+    cpdef inline void update_outputs(self) noexcept nogil
+@cython.final
+cdef class Model(masterinterface.MasterInterface):
+    cdef public Parameters parameters
+    cdef public Sequences sequences
+    cpdef inline void simulate(self, numpy.int64_t idx)  noexcept nogil
+    cpdef void reset_reuseflags(self) noexcept nogil
+    cpdef void load_data(self, numpy.int64_t idx) noexcept nogil
+    cpdef void save_data(self, numpy.int64_t idx) noexcept nogil
+    cpdef inline void run(self) noexcept nogil
+    cpdef inline void update_inlets(self) noexcept nogil
+    cpdef inline void update_outlets(self) noexcept nogil
+    cpdef inline void update_receivers(self, numpy.int64_t idx) noexcept nogil
+    cpdef inline void update_senders(self, numpy.int64_t idx) noexcept nogil
+    cpdef inline void update_outputs(self) noexcept nogil
+    cpdef inline void calc_temperature_v1(self) noexcept nogil
+    cpdef inline void adjust_temperature_v1(self) noexcept nogil
+    cpdef inline void calc_meantemperature_v1(self) noexcept nogil
+    cpdef void determine_temperature_v1(self) noexcept nogil
+    cpdef double get_temperature_v1(self, numpy.int64_t s) noexcept nogil
+    cpdef double get_meantemperature_v1(self) noexcept nogil
+    cpdef inline void calc_temperature(self) noexcept nogil
+    cpdef inline void adjust_temperature(self) noexcept nogil
+    cpdef inline void calc_meantemperature(self) noexcept nogil
+    cpdef void determine_temperature(self) noexcept nogil
+    cpdef double get_temperature(self, numpy.int64_t s) noexcept nogil
+    cpdef double get_meantemperature(self) noexcept nogil