Removing print screen and updating plot scripts to allow for specific 20% plot functions

franklin1895 · franklin1895 · commit 03b3d9ac34cc · 2025-06-11T18:10:22.000-04:00
diff --git a/cascade/roadway_manager.py b/cascade/roadway_manager.py
@@ -159,10 +159,10 @@ def bulldoze(
     if ((seaside_water_cells > percent_water_cells_touching_road) or (
         bayside_water_cells > percent_water_cells_touching_road)) and  allow_causeway==False:
         roadway_drown = True
-        print(
-            f"Roadway width drowned at {time_index - 1} years, "
-            f"{percent_water_cells_touching_road * 100.0}% of road borders water"
-        )
+        # print(
+        #     f"Roadway width drowned at {time_index - 1} years, "
+        #     f"{percent_water_cells_touching_road * 100.0}% of road borders water"
+        # )
     else:
         roadway_drown = False
 
@@ -236,7 +236,7 @@ def rebuild_dunes(
     )
 
     rebuild_dune_volume = np.sum(new_dune_domain - yxz_dune_grid)
-
+    z = 20
     return new_dune_domain, rebuild_dune_volume
 
 def build_interior_dunes(
@@ -430,10 +430,10 @@ def get_road_relocation_elevation(
     # or relocated
     if road_ele <= 0:
         roadway_drown = 1
-        print(
-            f"Roadway cannot be relocated at {time_index - 1} years b/c the road "
-            "would be at or below MSL"
-        )
+        # print(
+        #     f"Roadway cannot be relocated at {time_index - 1} years b/c the road "
+        #     "would be at or below MSL"
+        # )
 
     return road_ele, roadway_drown
 
@@ -486,7 +486,7 @@ def road_relocation_checks(
 
         if road_setback < 0:
             road_relocated = 1
-            print('Roadway relocated')
+            #print('Roadway relocated')
 
             # relocate the road only if the width of the island allows it
             if (
@@ -498,10 +498,10 @@ def road_relocation_checks(
             ):
                 relocation_break = 1
                 # time_index - 1 because B3D advances time step at end of dune_update
-                print(
-                    "Island is too narrow for roadway to be relocated. Roadway "
-                    f"eaten up by dunes at {time_index - 1} years"
-                )
+                # print(
+                #     "Island is too narrow for roadway to be relocated. Roadway "
+                #     f"eaten up by dunes at {time_index - 1} years"
+                # )
             else:
                 road_setback = road_relocation_setback
 
@@ -524,6 +524,7 @@ def check_sandbag_need(
                 for cell in range(0, len(barrier3d.DuneDomain[time_index, :, width])):
                     if barrier3d.DuneDomain[time_index, cell, width] < threshold_elevation:
                         barrier3d._DuneRestart[width][cell] = design_elevation / 10
+                        c = 10
 
         if exceeds_min_dune_threshold == True:
             sandbag_need = True
@@ -771,10 +772,10 @@ def update(self, barrier3d, trigger_dune_knockdown):
         # road cannot be below 0 m MHW (sea level); stop managing!
         if self._road_ele < 0:
             self._drown_break = 1
-            print(
-                f"Roadway drowned in place at {self._time_index - 1} years due to "
-                "SLR - road cannot be below 0 m MHW"
-            )
+            # print(
+            #     f"Roadway drowned in place at {self._time_index - 1} years due to "
+            #     "SLR - road cannot be below 0 m MHW"
+            # )
 
             # an adaptation solution may be to knock down the dunes so that they
             # are small and can easily be overwashed
diff --git a/data/Ocracoke_init_data/Ocracoke-CASCADE-parameters.yaml b/data/Ocracoke_init_data/Ocracoke-CASCADE-parameters.yaml
@@ -20,7 +20,7 @@ GrowthParamStart: false
 HdDiffu: 0.75
 Ki: 7.5e-06
 Kr: 7.5e-05
-LShoreface: 523.9734307609483
+LShoreface: 524.5623739370453
 MHW: 0.46
 MaxShrubHeight: 5.3
 MaxUpSlope: 0.25
@@ -31,29 +31,29 @@ Qs_min: 1.0
 Qshrub_max: 0.15
 RSLR_Constant: true
 RSLR_const: 0.0056
-Rat: 0
+Rat: 20
 Rin_i: 0.1
 Rin_r: 2.0
 SalineLimit: 5
 Sandbag_elevation: 1.9
 SeededRNG: true
 Seedmax: 1000.0
 Seedmin: 100.0
-ShorefaceToe: 1560.3448275999995
+ShorefaceToe: 1919.9999999999995
 ShrubEl_max: 2.3
 ShrubEl_min: 1.2
 Shrub_ON: 0
 SprayDist: 170
 StormSeries: []
 StormStart: 2
-TMAX: 126
+TMAX: 100
 TideAmp: 1.2
 TimeFruit: 5.0
 UprootLimit: -0.2
 beta: 0.04
 disp_mu: -0.721891
 disp_sigma: 1.5
-dune_file: C:\Users\frank\PycharmProjects\CASCADE\data\Ocracoke_init_data\dunes\Sample_1_dune.npy
+dune_file: C:\Users\frank\PycharmProjects\CASCADE\data\Ocracoke_init_data\dunes\Dunes_2019\Domain_49_dune_2019.npy
 elevation_file: C:\Users\frank\PycharmProjects\CASCADE\data\Ocracoke_init_data\elevations\Topography_2019\Domain_49_topography_2019.npy
 enable_sandbags: false
 growth_param_file: growthparam_1000dam.npy
@@ -65,4 +65,4 @@ rmin: 0.55
 s_sf_eq: 0.01698590034494604
 sandbag_elevation: 1.8
 sandbag_need: false
-storm_file: C:\Users\frank\PycharmProjects\CASCADE\data\Ocracoke_init_data\storms\Synthetic_Storms\OCR_Future_StormList_99_baseline.npy
+storm_file: C:\Users\frank\PycharmProjects\CASCADE\data\Ocracoke_init_data\storms\Synthetic_Storms\OCR_Future_StormList_1_baseline.npy
diff --git a/scripts/ocracoke_ms/Future_Runs_Test.py b/scripts/ocracoke_ms/Future_Runs_Test.py
@@ -15,7 +15,7 @@
 Save_Path = 'C:\\Users\\frank\\OneDrive - University of North Carolina at Chapel Hill\\Chapter 3\\Model Runs\\Summary_Values\\'
 
 #Management_Name = ['Status_Quo','Natural', 'Nourishment']
-Management_Name = ['Preemptive_Relocation']
+Management_Name = ['Status_Quo']
 #RSLR_Rate = ['IL','I','IH']
 RSLR_Rate = ['I']
 Sink_Name = ['Erosional_Sink','Accretional_Sink']
diff --git a/scripts/ocracoke_ms/Test_Bay_Interior_Location.py b/scripts/ocracoke_ms/Test_Bay_Interior_Location.py
@@ -9,7 +9,7 @@
 
 os.chdir('E:\\Model_Runs')
 
-run_name ='OCR_IL_Status_Quo_S2_Erosional_Sink'
+run_name ='OCR_IL_Nourishment_S2_Erosional_Sink'
 
 output = np.load(run_name + ".npz", allow_pickle=True)
 cascade = output["cascade"]
diff --git a/scripts/ocracoke_ms/run_forward_batch_forward_simulations.py b/scripts/ocracoke_ms/run_forward_batch_forward_simulations.py
@@ -11,7 +11,7 @@
 os.chdir('C:\\Users\\frank\\PycharmProjects\\CASCADE')
 
 # Set the number of years to simulate
-run_years = 126
+run_years = 100
 
 # Set the start year
 start_year = 2024
@@ -22,8 +22,8 @@
 Total_B3D_Number = 69
 
 # Set scenario type
-preemptive_relocation = True
-status_quo = False
+preemptive_relocation = False
+status_quo = True
 nourishment = False
 
 if preemptive_relocation == True:
@@ -59,7 +59,7 @@
 Storms = 'Baseline'
 
 # RSLR Data
-RSLR_Type = 'I'
+RSLR_Type = 'IH'
 
 if RSLR_Type == 'IL':
     RSLR_Data = np.load('C:\\Users\\frank\\PycharmProjects\\CASCADE\\data\\Ocracoke_init_data\\RSLR\\Int_Low_SLR.npy')
@@ -85,7 +85,7 @@
 run_name = []
 for snames in range(start_num,end_num):
     if Storms == 'Baseline':
-        name_base = 'OCR_'+str(RSLR_Type)+str(Management_name)+'S'+str(snames)
+        name_base = 'OCR_'+str(RSLR_Type)+str(Management_name)+'S'+str(snames)+'_T'
     else:
         name_base = 'OCR_'+str(RSLR_Type)+str(Management_name)+'S'+str(snames)+'_10'
     Temp_Name = []
@@ -371,7 +371,7 @@ def alongshore_uniform(run_name,
     background_erosion = background_erosion_list #list(background_erosion_list) # m/yr, background shoreline erosion
     rebuild_dune_threshold = 1  # m
     sandbag_management_on = sandbag_cells
-    sandbag_elevation = 4 # m
+    sandbag_elevation = 2 # m
 
     # baseline models for comparison -- all roadways ----------------------------------------
     roads_on = road_cells
@@ -422,7 +422,7 @@ def alongshore_uniform(run_name,
         allow_causeway = allow_causeway
     )
 
-for k in range(68,len(run_name)):
+for k in range(0,len(run_name)):
     for l in range(0,len(background_threshold_list)):
         alongshore_uniform(run_name=run_name[k][l],
                            s_file=s_file[k],
