Added calculation of @kaltevuus from @korkeus as a fallback

Author: s-hiitola

Scripts/assignment/emme_assignment.py

@@ -109,6 +109,7 @@ def prepare_network(self,
                 pass
             self.emme_project.modeller.emmebank.create_function(
                 idx, param.volume_delay_funcs[idx])
+        self._calculate_gradients()
         self.emme_project.create_extra_function_parameters(el1="@kaltevuus")
 
     def init_assign(self, 
@@ -319,6 +320,11 @@ def _add_bus_stops(self):
         modified_network: Network = mnw.add_bus_stops(network)
         self.mod_scenario.publish_network(modified_network)
 
+    def _calculate_gradients(self):
+        network: Network = self.mod_scenario.get_network()
+        modified_network: Network = mnw.calculate_gradients(network)
+        self.mod_scenario.publish_network(modified_network)
+
     def _create_matrices(self, time_period, id_hundred, id_ten):
         """Create EMME matrices for storing demand and impedance.
 

Scripts/utils/modify_network.py

@@ -19,6 +19,18 @@ def warn(self, message):
 
 # Utility functions that modify a network. Functions take a network as input and return the modified network
 
+def calculate_gradients(network):
+    for link in network.links():
+        if link['@kaltevuus'] == 0 and not link.i_node.is_centroid and not link.j_node.is_centroid:
+            try:
+                if link.i_node['@korkeus'] != 0 and link.j_node['@korkeus'] != 0:
+                    gradient = ((link.i_node['@korkeus'] - link.j_node['@korkeus']) / link.length)*0.1
+                    link['@kaltevuus'] = gradient
+                    log.debug(f"Calculated @kaltevuus for link {link.id}: {gradient}")
+            except KeyError:
+                log.info("@korkeus extra_attribute has not been defined. Skipping adjustment of @kaltevuus values")
+    return network
+
 def add_bus_stops(network):
     # Initialize an empty dictionary to store line IDs and maximum stop distances
     data = {"line_id": [], "maximum_stop_distance": [], "is_motorway": [], "loops": []}

Scripts/utils/validate_loaded_network.py

@@ -2,7 +2,7 @@
 
 def validate_loaded(network, fares=None):
     for line in network.transit_lines():
-        if str(line.vehicle.mode) in "mr":
+        if str(line.vehicle.mode) in "mrj":
             segments = list(line.segments())
             for seg1,seg2 in zip(segments[:-1],segments[1:]):
                 if seg1.data1 == 0 and (seg2.allow_boardings == 1 or seg2.allow_alightings ==1):