rtc-tools 2.6.0b2__py3-none-any.whl → 2.7.0a1__py3-none-any.whl

{rtc_tools-2.6.0b2.dist-info → rtc_tools-2.7.0a1.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: rtc-tools
-Version: 2.6.0b2
+Version: 2.7.0a1
 Summary: Toolbox for control and optimization of water systems.
 Home-page: https://oss.deltares.nl/web/rtc-tools/home
 Author: Deltares
@@ -29,8 +29,9 @@ License-File: COPYING.LESSER
 Requires-Dist: casadi ==3.6.*,>=3.6.3
 Requires-Dist: numpy <1.26,>=1.16.0
 Requires-Dist: scipy <1.11,>=1.0.0
-Requires-Dist: pymoca ==0.9.1
+Requires-Dist: pymoca ==0.9.*,>=0.9.1
 Requires-Dist: rtc-tools-channel-flow >=1.1.0
+Requires-Dist: defusedxml >=0.7.0
 Provides-Extra: all
 Requires-Dist: netCDF4 ; extra == 'all'
 Provides-Extra: netcdf

{rtc_tools-2.6.0b2.dist-info → rtc_tools-2.7.0a1.dist-info}/RECORD

@@ -1,6 +1,6 @@
 rtctools/__init__.py,sha256=91hvS2-ryd2Pvw0COpsUzTwJwSnTZ035REiej-1hNI4,107
-rtctools/_version.py,sha256=dUPcNcf2xQ9Ynkv1d2c0aTxHuz_HBjht0yZ3i8NDzGc,499
-rtctools/rtctoolsapp.py,sha256=A6M0v6vymAZWdBPLzgRVyAUfahD0knC8B4IptJ7_sBM,4097
+rtctools/_version.py,sha256=My52FWqz9UB3_PwxQjfVCPUkHjej_S48gNc1rsNpaKc,499
+rtctools/rtctoolsapp.py,sha256=UnkuiJhv0crEEVs8H6PYvMuc2y_q6V_xLuyKEgXj9GM,4200
 rtctools/util.py,sha256=PaeKfDUA174ODZbY5fZjCTf-F-TdhW7yEuP189Ro190,9075
 rtctools/_internal/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 rtctools/_internal/alias_tools.py,sha256=XuQSAhhFuVtwn0yrAObZWIKPsSF4j2axXRtEmitIFPs,5310
@@ -10,8 +10,8 @@ rtctools/_internal/debug_check_helpers.py,sha256=UgQTEPw4PyR7MbYLewSSWQqTwQj7xr5
 rtctools/data/__init__.py,sha256=EllgSmCdrlvQZSd1VilvjPaeYJGhY9ErPiQtedmuFoA,157
 rtctools/data/csv.py,sha256=iYOEED3AaNxt4ews_aAkHfl9Tq9a-9vjxvYwjVR_lQE,5231
 rtctools/data/netcdf.py,sha256=xpk4Ggl7gItNG6lO7p3OJPR-elK8_CiCtxUI7cX0gwk,19109
-rtctools/data/pi.py,sha256=hBgK_dSLv_V9pWFGNzOEnG76ffLsjW3SBSWfgVJpXSA,45380
-rtctools/data/rtc.py,sha256=FBuUQ6aL4D3y6puKudyuLAeP_vmpDiliqzJg8f1kd7g,9043
+rtctools/data/pi.py,sha256=Ni1hBDdhQdcWYO-NUPhKA1WJdzbSXjZc5w5xauOBcJM,45437
+rtctools/data/rtc.py,sha256=1yGJZGq2Z36MYLiLuZaHnxupL4mgw-Wuu54PAG05kcM,9077
 rtctools/data/storage.py,sha256=67J4BRTl0AMEzlKNZ8Xdpy_4cGtwx8Lo_tL2n0G4S9w,13206
 rtctools/data/interpolation/__init__.py,sha256=GBubCIT5mFoSTV-lOk7cpwvZekNMEe5bvqSQJ9HE34M,73
 rtctools/data/interpolation/bspline.py,sha256=qevB842XWCH3fWlWMBqKMy1mw37ust-0YtSnb9PKCEc,948
@@ -20,7 +20,7 @@ rtctools/data/interpolation/bspline2d.py,sha256=ScmX0fPDxbUVtj3pbUE0L7UJocqroD_6
 rtctools/optimization/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 rtctools/optimization/collocated_integrated_optimization_problem.py,sha256=Bv1a1O4QTAdIpe9SOA9_oEDDph-oj-9q5yfN4SFMdW8,131343
 rtctools/optimization/control_tree_mixin.py,sha256=CC6TWp3kFQgMokx6213pRLx9iY17Fd0VcwG4Wpwa0Uo,8974
-rtctools/optimization/csv_lookup_table_mixin.py,sha256=xiyOKpPT3kQ6yuxSJRorGYW_QwnlBzW6vvrtvQNhHsE,17264
+rtctools/optimization/csv_lookup_table_mixin.py,sha256=h4WKuPR1BJbYuJfQ9lx4rXalB6yYK-zajINabRL7BrA,17298
 rtctools/optimization/csv_mixin.py,sha256=sRp5paHWfCw2bz-23Nw-HdFLS3CZTpVwaBdFo98DbvE,12252
 rtctools/optimization/goal_programming_mixin.py,sha256=GK25DrbAY_rMsra080pSCDZwzLQNN2Ppd-2d0_FEllg,32999
 rtctools/optimization/goal_programming_mixin_base.py,sha256=oh9CsEiyYTmthcfvRbX-9Z9bIo6SHv_DCiVt9kx0sjI,43781
@@ -30,21 +30,21 @@ rtctools/optimization/io_mixin.py,sha256=AsZQ7YOUcUbWoczmjTXaSje5MUEsPNbQyZBJ6qz
 rtctools/optimization/linearization_mixin.py,sha256=mG5S7uwvwDydw-eBPyQKnLyKoy08EBjQh25vu97afhY,1049
 rtctools/optimization/linearized_order_goal_programming_mixin.py,sha256=LQ2qpYt0YGLpEoerif4FJ5wwzq16q--27bsRjcqIU5A,9087
 rtctools/optimization/min_abs_goal_programming_mixin.py,sha256=WMOv9EF8cfDJgTunzXfI_cUmBSQK26u1HJB_9EAarfM,14031
-rtctools/optimization/modelica_mixin.py,sha256=ysVMayNA4sSFoHkSdhjWOxT6UzOVbN0ZeM4v-RpvZXE,17161
+rtctools/optimization/modelica_mixin.py,sha256=6rotTEge4l2oDib4tnG2ccXqG1ryEzNxbRkk_StFAHg,17769
 rtctools/optimization/netcdf_mixin.py,sha256=-zkXh3sMYE50c3kHsrmUVGWMSFm-0cXQpGrCm0yn-Tc,7563
 rtctools/optimization/optimization_problem.py,sha256=qzpc81NaZMeoXKuayFmBF15iXYuNAk5yxmaER_Gcz_A,44131
-rtctools/optimization/pi_mixin.py,sha256=63qda6i7hjtDuP3hL6RO29vCCP11aUpR9B4KoqlLFVI,11314
+rtctools/optimization/pi_mixin.py,sha256=n_ZyifuNPxcRoJXKR6Ks7DgK_b7-EToYJ23_khOwpzo,11786
 rtctools/optimization/planning_mixin.py,sha256=O_Y74X8xZmaNZR4iYOe7BR06s9hnmcapbuHYHQTBPPQ,724
 rtctools/optimization/single_pass_goal_programming_mixin.py,sha256=Zb9szg3PGT2o6gkGsXluSfEaAswkw3TKfPQDzUrj_Y4,25784
 rtctools/optimization/timeseries.py,sha256=nCrsGCJThBMh9lvngEpbBDa834_QvklVvkxJqwX4a1M,1734
 rtctools/simulation/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 rtctools/simulation/csv_mixin.py,sha256=rGDUFPsqGHmF0_dWdXeWzWzMpkPmwCNweTBVrwSh31g,6704
 rtctools/simulation/io_mixin.py,sha256=SJasNGI--OQ9Y-Z61oeeaGCxSrNddYz4AOVfJYbmf74,6209
-rtctools/simulation/pi_mixin.py,sha256=uwl61LYjb8dmMz910EB2-bC0KSuhLzsrJzk0hxWYEhk,9359
-rtctools/simulation/simulation_problem.py,sha256=gTAimG2MLw_TTkeHLkIMxpYgAmR-voqzvje7pcFnw4U,44556
-rtc_tools-2.6.0b2.dist-info/COPYING.LESSER,sha256=46mU2C5kSwOnkqkw9XQAJlhBL2JAf1_uCD8lVcXyMRg,7652
-rtc_tools-2.6.0b2.dist-info/METADATA,sha256=J_s90hgXnSTTIZg76jopXRUeL0V3VUJfOyvorAaE_1s,1452
-rtc_tools-2.6.0b2.dist-info/WHEEL,sha256=oiQVh_5PnQM0E3gPdiz09WCNmwiHDMaGer_elqB3coM,92
-rtc_tools-2.6.0b2.dist-info/entry_points.txt,sha256=-x622IB_l1duw2M6t6syfQ4yzOiQTp0IZxKGcYRgWgk,151
-rtc_tools-2.6.0b2.dist-info/top_level.txt,sha256=pnBrb58PFPd1kp1dqa-JHU7R55h3alDNJIJnF3Jf9Dw,9
-rtc_tools-2.6.0b2.dist-info/RECORD,,
+rtctools/simulation/pi_mixin.py,sha256=vcizArZsOn7tGlSFCsUmkjmRFA1A1kMNqUllLlkVc9Y,9831
+rtctools/simulation/simulation_problem.py,sha256=wGRIOm4AHAy75R7cbqeA13-Qy1ED3RY3Yc5Hvtmv6PY,48541
+rtc_tools-2.7.0a1.dist-info/COPYING.LESSER,sha256=46mU2C5kSwOnkqkw9XQAJlhBL2JAf1_uCD8lVcXyMRg,7652
+rtc_tools-2.7.0a1.dist-info/METADATA,sha256=p33wEk1F8t8jCHyV1bUe5R4JcedANOtp-9NsTH1mWBc,1494
+rtc_tools-2.7.0a1.dist-info/WHEEL,sha256=GJ7t_kWBFywbagK5eo9IoUwLW6oyOeTKmQ-9iHFVNxQ,92
+rtc_tools-2.7.0a1.dist-info/entry_points.txt,sha256=-x622IB_l1duw2M6t6syfQ4yzOiQTp0IZxKGcYRgWgk,151
+rtc_tools-2.7.0a1.dist-info/top_level.txt,sha256=pnBrb58PFPd1kp1dqa-JHU7R55h3alDNJIJnF3Jf9Dw,9
+rtc_tools-2.7.0a1.dist-info/RECORD,,

{rtc_tools-2.6.0b2.dist-info → rtc_tools-2.7.0a1.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: bdist_wheel (0.42.0)
+Generator: bdist_wheel (0.43.0)
 Root-Is-Purelib: true
 Tag: py3-none-any
 

rtctools/_version.py

@@ -8,11 +8,11 @@ import json
 
 version_json = '''
 {
- "date": "2024-01-18T16:11:20+0000",
+ "date": "2024-04-25T06:42:12+0000",
  "dirty": false,
  "error": null,
- "full-revisionid": "e868661a3b37bb27f61ffcf9cdf1d53d857d5108",
- "version": "2.6.0b2"
+ "full-revisionid": "e3d59506cabe8df7e05acd089ccf574e106a577f",
+ "version": "2.7.0a1"
 }
 '''  # END VERSION_JSON
 

rtctools/data/pi.py

@@ -5,6 +5,7 @@ import logging
 import os
 import xml.etree.ElementTree as ET
 
+import defusedxml.ElementTree as DefusedElementTree
 import numpy as np
 
 ns = {"fews": "http://www.wldelft.nl/fews", "pi": "http://www.wldelft.nl/fews/PI"}
@@ -30,7 +31,7 @@ class Diag:
         """
         self.__path_xml = os.path.join(folder, basename + ".xml")
 
-        self.__tree = ET.parse(self.__path_xml)
+        self.__tree = DefusedElementTree.parse(self.__path_xml)
         self.__xml_root = self.__tree.getroot()
 
     def get(self, level=ERROR_FATAL):
@@ -87,7 +88,7 @@ class DiagHandler(logging.Handler):
         self.__path_xml = os.path.join(folder, basename + ".xml")
 
         try:
-            self.__tree = ET.parse(self.__path_xml)
+            self.__tree = DefusedElementTree.parse(self.__path_xml)
             self.__xml_root = self.__tree.getroot()
         except Exception:
             self.__xml_root = ET.Element("{%s}Diag" % (ns["pi"],))
@@ -135,7 +136,7 @@ class ParameterConfig:
             basename = basename + ".xml"
         self.__path_xml = os.path.join(folder, basename)
 
-        self.__tree = ET.parse(self.__path_xml)
+        self.__tree = DefusedElementTree.parse(self.__path_xml)
         self.__xml_root = self.__tree.getroot()
 
     def get(self, group_id, parameter_id, location_id=None, model=None):
@@ -300,9 +301,8 @@ class ParameterConfig:
                 # get table contenstart_datetime
                 el_row = child.findall("pi:row", ns)
                 table = {
-                    columnId[key]: np.empty(len(el_row), columnType[key])  # initialize table
-                    for key in columnId
-                }
+                    columnId[key]: np.empty(len(el_row), columnType[key]) for key in columnId
+                }  # initialize table
 
                 i_row = 0
                 for row in el_row:
@@ -378,7 +378,7 @@ class Timeseries:
         if self.make_new_file:
             self.__reset_xml_tree()
         else:
-            self.__tree = ET.parse(self.__path_xml)
+            self.__tree = DefusedElementTree.parse(self.__path_xml)
             self.__xml_root = self.__tree.getroot()
 
         self.__values = [{}]
@@ -855,7 +855,7 @@ class Timeseries:
 
                 variable = self.__data_config.variable(header)
 
-                miss_val = float(header.find("pi:missVal", ns).text)
+                miss_val = header.find("pi:missVal", ns).text
                 values = self.__values[ensemble_member][variable]
 
                 # Update the header, which may have changed
@@ -868,11 +868,8 @@ class Timeseries:
                     self.__xml_root.remove(series)
                     continue
 
-                # Replace NaN with missing value
-                nans = np.isnan(values)
-                values[nans] = miss_val
-
                 # Write output
+                nans = np.isnan(values)
                 if self.__binary:
                     f.write(values.astype(self.__pi_dtype).tobytes())
                 else:
@@ -897,7 +894,10 @@ class Timeseries:
                         event = ET.Element("pi:event")
                         event.set("date", t.strftime("%Y-%m-%d"))
                         event.set("time", t.strftime("%H:%M:%S"))
-                        event.set("value", str(values[i]))
+                        if nans[i]:
+                            event.set("value", miss_val)
+                        else:
+                            event.set("value", str(values[i]))
                         series.append(event)
                         if self.dt:
                             t += self.dt
@@ -907,9 +907,6 @@ class Timeseries:
                         for i in range(len(values), len(events)):
                             series.remove(events[i])
 
-                # Restore NaN
-                values[nans] = np.nan
-
         if self.__binary:
             f.close()
 

rtctools/data/rtc.py

@@ -1,8 +1,9 @@
 import logging
 import os
-import xml.etree.ElementTree as ET
 from collections import namedtuple
 
+import defusedxml.ElementTree as DefusedElementTree
+
 ts_ids = namedtuple("ids", "location_id parameter_id qualifier_id")
 p_ids = namedtuple("ids", "model_id location_id parameter_id")
 
@@ -31,7 +32,7 @@ class DataConfig:
 
         path = os.path.join(folder, "rtcDataConfig.xml")
         try:
-            tree = ET.parse(path)
+            tree = DefusedElementTree.parse(path)
             root = tree.getroot()
 
             timeseriess1 = root.findall("./*/fews:timeSeries", ns)

rtctools/optimization/csv_lookup_table_mixin.py

@@ -2,7 +2,6 @@ import configparser
 import glob
 import logging
 import os
-import pickle
 from typing import Iterable, List, Tuple, Union
 
 import casadi as ca
@@ -323,7 +322,7 @@ class CSVLookupTableMixin(OptimizationProblem):
 
             # If tck file is newer than the csv file, first try to load the cached values from
             # the tck file
-            tck_filename = filename.replace(".csv", ".tck")
+            tck_filename = filename.replace(".csv", ".npz")
             valid_cache = False
             if os.path.exists(tck_filename):
                 if no_curvefit_options:
@@ -338,11 +337,13 @@ class CSVLookupTableMixin(OptimizationProblem):
                             output
                         )
                     )
-                    with open(tck_filename, "rb") as f:
-                        try:
-                            tck, function = pickle.load(f)
-                        except Exception:
-                            valid_cache = False
+                    try:
+                        with np.load(filename.replace(".csv", ".npz")) as data:
+                            tck = (data["arr_0"], data["arr_1"], int(data["arr_2"]))
+                        function = ca.Function.load(filename.replace(".csv", ".ca"))
+                    except Exception:
+                        valid_cache = False
+
             if not valid_cache:
                 logger.info("CSVLookupTableMixin: Recalculating tck values for {}".format(output))
 
@@ -446,11 +447,8 @@ class CSVLookupTableMixin(OptimizationProblem):
                 )
 
             if not valid_cache:
-                pickle.dump(
-                    (tck, function),
-                    open(filename.replace(".csv", ".tck"), "wb"),
-                    protocol=pickle.HIGHEST_PROTOCOL,
-                )
+                np.savez(filename.replace(".csv", ".npz"), *tck)
+                function.save(filename.replace(".csv", ".ca"))
 
     def lookup_tables(self, ensemble_member):
         # Call parent class first for default values.

rtctools/optimization/modelica_mixin.py

@@ -48,9 +48,21 @@ class ModelicaMixin(OptimizationProblem):
             else:
                 model_name = self.__class__.__name__
 
-        self.__pymoca_model = pymoca.backends.casadi.api.transfer_model(
-            kwargs["model_folder"], model_name, self.compiler_options()
-        )
+        compiler_options = self.compiler_options()
+        logger.info(f"Loading/compiling model {model_name}.")
+        try:
+            self.__pymoca_model = pymoca.backends.casadi.api.transfer_model(
+                kwargs["model_folder"], model_name, compiler_options
+            )
+        except RuntimeError as error:
+            if compiler_options.get("cache", False):
+                raise error
+            compiler_options["cache"] = False
+            logger.warning(f"Loading model {model_name} using a cache file failed: {error}.")
+            logger.info(f"Compiling model {model_name}.")
+            self.__pymoca_model = pymoca.backends.casadi.api.transfer_model(
+                kwargs["model_folder"], model_name, compiler_options
+            )
 
         # Extract the CasADi MX variables used in the model
         self.__mx = {}

rtctools/optimization/pi_mixin.py

@@ -277,3 +277,16 @@ class PIMixin(IOMixin):
         :class:`pi.Timeseries` object for holding the output data.
         """
         return self.__timeseries_export
+
+    def set_unit(self, variable: str, unit: str):
+        """
+        Set the unit of a time series.
+
+        :param variable:        Time series ID.
+        :param unit:            Unit.
+        """
+        assert hasattr(
+            self, "_PIMixin__timeseries_import"
+        ), "set_unit can only be called after read() in pre() has finished."
+        self.__timeseries_import.set_unit(variable, unit, 0)
+        self.__timeseries_export.set_unit(variable, unit, 0)

rtctools/rtctoolsapp.py

@@ -108,7 +108,8 @@ def download_examples(*args):
 
         opener = urllib.request.build_opener()
         urllib.request.install_opener(opener)
-        local_filename, _ = urllib.request.urlretrieve(url)
+        # The security warning can be dismissed as the url variable is hardcoded to a remote.
+        local_filename, _ = urllib.request.urlretrieve(url)  # nosec
     except HTTPError:
         sys.exit("Could not found examples for RTC-Tools version {}.".format(version))
 

rtctools/simulation/pi_mixin.py

@@ -240,3 +240,16 @@ class PIMixin(IOMixin):
     def get_timeseries(self, variable):
         _, values = self.io.get_timeseries(variable)
         return values
+
+    def set_unit(self, variable: str, unit: str):
+        """
+        Set the unit of a time series.
+
+        :param variable:        Time series ID.
+        :param unit:            Unit.
+        """
+        assert hasattr(
+            self, "_PIMixin__timeseries_import"
+        ), "set_unit can only be called after read() in pre() has finished."
+        self.__timeseries_import.set_unit(variable, unit, 0)
+        self.__timeseries_export.set_unit(variable, unit, 0)

rtctools/simulation/simulation_problem.py

@@ -82,9 +82,21 @@ class SimulationProblem(DataStoreAccessor):
                 model_name = self.__class__.__name__
 
         # Load model from pymoca backend
-        self.__pymoca_model = pymoca.backends.casadi.api.transfer_model(
-            kwargs["model_folder"], model_name, self.compiler_options()
-        )
+        compiler_options = self.compiler_options()
+        logger.info(f"Loading/compiling model {model_name}.")
+        try:
+            self.__pymoca_model = pymoca.backends.casadi.api.transfer_model(
+                kwargs["model_folder"], model_name, compiler_options
+            )
+        except RuntimeError as error:
+            if compiler_options.get("cache", False):
+                raise error
+            compiler_options["cache"] = False
+            logger.warning(f"Loading model {model_name} using a cache file failed: {error}.")
+            logger.info(f"Compiling model {model_name}.")
+            self.__pymoca_model = pymoca.backends.casadi.api.transfer_model(
+                kwargs["model_folder"], model_name, compiler_options
+            )
 
         # Extract the CasADi MX variables used in the model
         self.__mx = {}
@@ -319,6 +331,17 @@ class SimulationProblem(DataStoreAccessor):
         """
         Initialize state vector with default values
 
+        Initial values are first read from the given Modelica files.
+        If an initial value equals zero or is not provided by a Modelica file,
+        and the variable is not marked as fixed,
+        then the initial value is tried to be set with the initial_state method.
+        When using CSVMixin, this method by default looks for initial values
+        in an initial_state.csv file.
+        Furthermore, if a variable is not marked as fixed
+        and no initial value is given by the initial_state method,
+        the initial value can be overwritten using the seed method.
+        When a variable is marked as fixed, the initial value is only read from the Modelica file.
+
         :param config_file: Path to an initialization file.
         """
         if config_file:
@@ -393,29 +416,35 @@ class SimulationProblem(DataStoreAccessor):
         for var in itertools.chain(self.__pymoca_model.states, self.__pymoca_model.alg_states):
             var_name = var.symbol.name()
             var_nominal = self.get_variable_nominal(var_name)
+            start_values = {}
 
             # Attempt to cast var.start to python type
+            start_value_is_symbolic = False
             mx_start = ca.MX(var.start)
             if mx_start.is_constant():
                 # cast var.start to python type
-                start_val = var.python_type(mx_start.to_DM())
+                start_value_pymoca = var.python_type(mx_start.to_DM())
+                if start_value_pymoca is not None and start_value_pymoca != 0:
+                    start_values["modelica"] = start_value_pymoca
             else:
-                # var.start is a symbolic expression with unknown value
-                start_val = None
+                start_value_is_symbolic = True
+                start_values["modelica"] = mx_start
 
-            if start_val == 0.0 and not var.fixed:
+            if not var.fixed:
                 # To make initialization easier, we allow setting initial states by providing
                 # timeseries with names that match a symbol in the model. We only check for this
                 # matching if the start and fixed attributes were left as default
                 try:
-                    start_val = self.initial_state()[var_name]
+                    start_values["initial_state"] = self.initial_state()[var_name]
                 except KeyError:
                     pass
                 else:
                     # An initial state was found- add it to the constrained residuals
                     logger.debug(
                         "Initialize: Added {} = {} to initial equations "
-                        "(found matching timeseries).".format(var_name, start_val)
+                        "(found matching timeseries).".format(
+                            var_name, start_values["initial_state"]
+                        )
                     )
                     # Set var to be fixed
                     var.fixed = True
@@ -425,36 +454,75 @@ class SimulationProblem(DataStoreAccessor):
                 # timeseries with names that match a symbol in the model. We only check for this
                 # matching if the start and fixed attributes were left as default
                 try:
-                    start_val = self.seed()[var_name]
+                    start_values["seed"] = self.seed()[var_name]
                 except KeyError:
                     pass
                 else:
                     # An initial state was found- add it to the constrained residuals
                     logger.debug(
                         "Initialize: Added {} = {} as initial guess "
-                        "(found matching timeseries).".format(var_name, start_val)
+                        "(found matching timeseries).".format(var_name, start_values["seed"])
+                    )
+
+            # Set the start value based on the different inputs.
+            if "seed" in start_values:
+                input_source = "seed"
+                source_description = "seed method"
+            elif "modelica" in start_values:
+                input_source = "modelica"
+                source_description = "modelica file"
+            elif "initial_state" in start_values:
+                input_source = "initial_state"
+                source_description = "initial_state method (typically reads initial_state.csv)"
+            else:
+                start_values["modelica"] = start_value_pymoca
+                input_source = "modelica"
+                source_description = "modelica file or default value"
+            start_val = start_values.get(input_source, None)
+            is_numeric = start_val is not None and not start_value_is_symbolic
+            numeric_start_val = start_val if is_numeric else 0.0
+            if len(start_values) > 1:
+                logger.warning(
+                    "Initialize: Multiple initial values for {} are provided: {}.".format(
+                        var_name, start_values
                     )
+                    + " Value from {} will be used to continue.".format(source_description)
+                )
 
             # Attempt to set start_val in the state vector. Default to zero if unknown.
             try:
-                self.set_var(var_name, start_val if start_val is not None else 0.0)
+                self.set_var(var_name, numeric_start_val)
             except KeyError:
                 logger.warning(
                     "Initialize: {} not found in state vector. "
-                    "Initial value of {} not set.".format(var_name, start_val)
+                    "Initial value of {} not set.".format(var_name, numeric_start_val)
                 )
 
             # Add a residual for the difference between the state and its starting expression
-            start_expr = start_val if start_val is not None else var.start
+            start_expr = start_val
             if var.fixed:
                 # Set bounds to be equal to each other, such that IPOPT can
                 # turn the decision variable into a parameter.
+                if var.min != -np.inf or var.max != np.inf:
+                    logger.info(
+                        "Initialize: bounds of {} will be overwritten".format(var_name)
+                        + " by the start value given by {}.".format(source_description)
+                    )
                 var.min = start_expr
                 var.max = start_expr
             else:
                 # minimize residual
                 minimized_residuals.append((var.symbol - start_expr) / var_nominal)
 
+            # Check that the start_value is in between the variable bounds.
+            if start_val is not None and is_numeric:
+                if not (var.min <= start_val and start_val <= var.max):
+                    logger.warning(
+                        "Initialize: start value {} = {}".format(var_name, start_val)
+                        + " is not in between bounds {} and {}".format(var.min, var.max)
+                        + " and will be adjusted."
+                    )
+
         # Default start var for ders is zero
         for der_var in self.__mx["derivatives"]:
             self.set_var(der_var.name(), 0.0)
@@ -610,7 +678,15 @@ class SimulationProblem(DataStoreAccessor):
         # If unsuccessful, stop.
         return_status = solver.stats()["return_status"]
         if return_status not in {"Solve_Succeeded", "Solved_To_Acceptable_Level"}:
-            raise Exception('Initialization Failed with return status "{}"'.format(return_status))
+            if return_status == "Infeasible_Problem_Detected":
+                message = (
+                    "Initialization Failed with return status: {}. ".format(return_status)
+                    + "This means no initial state could be found "
+                    + "that satisfies all equations and constraints."
+                )
+            else:
+                message = "Initialization Failed with return status: {}. ".format(return_status)
+            raise Exception(message)
 
         # Update state vector with initial conditions
         self.__state_vector[: self.__n_states] = initial_state["x"][: self.__n_states].T