use SFCTA's bpr curve for collectors and locals

BayAreaMetro · Aug 2, 2024 · d34898c · d34898c
1 parent 172a183
commit d34898c
Showing 1 changed file with 8 additions and 2 deletions.
diff --git a/tm2py/components/network/create_tod_scenarios.py b/tm2py/components/network/create_tod_scenarios.py
@@ -163,18 +163,24 @@ def _create_highway_scenarios(self):
             # "(el1 + 60 * (0.25 *(put(put((volau + volad)/el2) - 1) + "
             # "(((get(2)*get(2) + (16 * el3 * get(1)^0.5))))"
         )
+        # SFCTA bpr curve https://github.com/BayAreaMetro/tm2py/issues/147#issuecomment-2261058755
+        bpr_local_tmplt="el1 * 1.8 * (1 + 0.60 * ((volau + volad)/el2)^8.5)"
         for f_id in ["fd1", "fd2"]:
             if emmebank.function(f_id):
                 emmebank.delete_function(f_id)
             emmebank.create_function(
                 f_id, bpr_tmplt + reliability_tmplt.format(**parameters["freeway"])
             )
+        for f_id in ["fd6", "fd7"]:
+            if emmebank.function(f_id):
+                emmebank.delete_function(f_id)
+            emmebank.create_function(
+                f_id, bpr_local_tmplt + reliability_tmplt.format(**parameters["road"])
+            )
         for f_id in [
             "fd3",
             "fd4",
             "fd5",
-            "fd6",
-            "fd7",
             "fd9",
             "fd10",
             "fd11",