Up to date version for meeting with KA. Modified location of sandbag …

…functions
UNC-CECL · Apr 26, 2024 · 3c853d4 · 3c853d4
1 parent 59d4ea4
commit 3c853d4
Show file tree

Hide file tree

Showing 4 changed files with 18 additions and 16 deletions.
diff --git a/cascade/cascade.py b/cascade/cascade.py
@@ -610,6 +610,21 @@ def update(self):
         # human dynamics modules
         ###############################################################################
 
+        # ~~~~~~~~~~~~~~ SandbagManager ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+        # If sandbag dynamics are enabled check if conditions are met for sandbag
+        # emplacement. When sandbag conditions are met, SandbagManager will rebuild
+        # dunes if they fall below a user defined threshold.
+        for iB3D in range(self._ny):
+            if self._sandbag_management_on[iB3D] == True:
+                sandbag_emplacement = check_sandbag_need(dune_road_distance=self._roadways[iB3D]._road_setback,
+                                                              design_elevation = self._sandbag_elevation,
+                                                              barrier3d = self._barrier3d[iB3D],
+                                                              )
+                self._sandbag_Need_TS[iB3D] = np.append(self._sandbag_Need_TS[iB3D],sandbag_emplacement)
+
+            # Check sandbag need
+
+
         # ~~~~~~~~~~~~~~ RoadwayManager ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
         # Remove overwash from roadway after each model year, place on the dune,
         # rebuild dunes if fall below height threshold, and check if dunes should
@@ -696,19 +711,6 @@ def update(self):
                     + (self._initial_beach_width[iB3D] / 10)  # dam
                 )
 
-        # ~~~~~~~~~~~~~~ SandbagManager ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-        # If sandbag dynamics are enabled check if conditions are met for sandbag
-        # emplacement. When sandbag conditions are met, SandbagManager will rebuild
-        # dunes if they fall below a user defined threshold.
-        for iB3D in range(self._ny):
-            if self._sandbag_management_on[iB3D] == True:
-                sandbag_emplacement = check_sandbag_need(dune_road_distance=self._roadways[iB3D]._road_setback,
-                                                              design_elevation = self._sandbag_elevation,
-                                                              barrier3d = self._barrier3d[iB3D],
-                                                              )
-                self._sandbag_Need_TS[iB3D] = np.append(self._sandbag_Need_TS[iB3D],sandbag_emplacement)
-
-            # Check sandbag need
 
         # ~~~ CHOM coupler (in development) ~~~
         # Provide agents in the Coastal Home Ownership Model (CHOM) with variables

diff --git a/cascade/roadway_manager.py b/cascade/roadway_manager.py
@@ -437,7 +437,7 @@ def check_sandbag_need(
         dune_road_distance,
         design_elevation,
         barrier3d,
-        threshold_elevation = 0.25,
+        threshold_elevation = 0.05,
 ):
     time_index = barrier3d.time_index -1
     print('Time index is '+str(time_index))

diff --git a/scripts/ocracoke_ms/plot_script.py b/scripts/ocracoke_ms/plot_script.py
@@ -326,7 +326,7 @@ def plot_ElevAnimation_CASCADE(
 
 os.chdir('C:\\Users\\frank\\PycharmProjects\\CASCADE\\Run_output')
 # run_name='Wreck_ACC_RSLR3_S3' # 5 Length
-run_name = "Sandbag_Change_80"  # 4 length
+run_name = "Sandbag_80"  # 4 length
 # run_name='Metompkin_Marsh_S10_3'
 # run_name='Smith_S10_3' # 5
 

diff --git a/scripts/ocracoke_ms/run_script.py b/scripts/ocracoke_ms/run_script.py
@@ -14,7 +14,7 @@
 s_file = 'C:\\Users\\frank\\PycharmProjects\\CASCADE\\data\\Ocracoke_init_data\\S1.npy'
 #d_file = 'C:\\Users\\frank\\PycharmProjects\\CASCADE\\data\\Ocracoke_init_data\\b3d_high-elevations.csv'
 #s_file = 'StormList_0_baseline.npy'
-run_name = 'Sandbag_Change_80'
+run_name = 'Sandbag_80'
 
 #offsets = np.load('C:\\Users\\frank\\PycharmProjects\\CASCADE\\data\\Ocracoke_init_data\\mod_offset.npy')
 #offsets = offsets[2:]