Updating Hatteras scripts, developed pipelines for road and island offset

Author: HannahAline

cascade/beach_dune_manager.py

@@ -209,11 +209,11 @@ def filter_overwash(
         Percent overwash removed from barrier interior to dunes
         [%, overwash_filter+overwash_to_dune <=100]
     post_storm_xyz_interior_grid: grid,
-        Interior barrier island topography [for Barrier3d, decameters MHW]
+        Interior barrier island topography_dunes [for Barrier3d, decameters MHW]
     pre_storm_xyz_interior_grid: grid,
-        Interior barrier island topography [for Barrier3d, decameters MHW]
+        Interior barrier island topography_dunes [for Barrier3d, decameters MHW]
     post_storm_yxz_dune_grid: grid,
-        Dune topography [for Barrier3d, decameters above the berm elevation]
+        Dune topography_dunes [for Barrier3d, decameters above the berm elevation]
     artificial_maximum_dune_height: float,
         The maximum dune height than can be created by bulldozer after a storm
         [for Barrier3d, dam above berm elevation]

cascade/roadway_manager.py

@@ -61,10 +61,10 @@ def bulldoze(
     time_index: int,
         Time index for drowning error message
     xyz_interior_grid: array
-        Interior barrier island topography [z units specified by dz; for Barrier3d,
+        Interior barrier island topography_dunes [z units specified by dz; for Barrier3d,
         dz=10, decameters MHW]
     yxz_dune_grid:
-        Dune topography [z units specified by dz; for Barrier3d, dz=10, decameters
+        Dune topography_dunes [z units specified by dz; for Barrier3d, dz=10, decameters
         above the berm elevation]
     road_ele: float
         Road elevation [m; needs to be in same reference frame as xyz; for Barrier3d,
@@ -193,7 +193,7 @@ def rebuild_dunes(
     Parameters
     ----------
     yxz_dune_grid: ndarray, shape (ny, nx)
-        Dune topography [z units specified by dz; for Barrier3D, dz=10, decameters
+        Dune topography_dunes [z units specified by dz; for Barrier3D, dz=10, decameters
         above the berm elevation]
     max_dune_height: float, optional
         Maximum dune height for dune rebuilding [m]
@@ -258,7 +258,7 @@ def set_growth_parameters(
     Parameters
     ----------
     yxz_dune_grid:
-        Dune topography [units must be the same as Dmax]
+        Dune topography_dunes [units must be the same as Dmax]
     Dmax: float
         Maximum natural equilibrium dune height [default in Barrier3D is decameters]
     growthparam:
@@ -318,7 +318,7 @@ def get_road_relocation_elevation(
     time_index: int,
         Time index for drowning error message
     xyz_interior_grid: array
-        Interior barrier island topography [z units specified by dz; for Barrier3d,
+        Interior barrier island topography_dunes [z units specified by dz; for Barrier3d,
         dz=10, decameters MHW]
     road_width: int
         Width of roadway [m]

cascade/tools/plotters.py

@@ -78,7 +78,7 @@ def plot_XShoreTransects(barrier3d, TMAX):
     plt.hlines(barrier3d._SL, -1, len(CrossElev + 1), colors="dodgerblue")
     plt.xlabel("Cross-Shore Distance (dam)")
     plt.ylabel("Elevation (m)")
-    plt.title("Cross-shore Topo Transects")
+    plt.title("Cross-shore PEA_2011 Transects")
     plt.show()
 
     return fig
@@ -1177,7 +1177,7 @@ def plot_nonlinear_stats_RoadwayManager(
     plt.subplot(2, 2, 2)
     if dunes_rebuilt is not None:
         plt.plot(dunes_rebuilt[0:tmax_sim], "k")
-        plt.ylabel("Dunes Rebuilt")
+        plt.ylabel("PEA_2011 Rebuilt")
         plt.xlabel("Time (yr)")
 
     # road relocated
@@ -1517,7 +1517,7 @@ def plot_nonlinear_stats_BeachDuneManager(
     plt.subplot(2, 2, 2)
     if dunes_rebuilt is not None:
         plt.plot(dunes_rebuilt, "k")
-        plt.ylabel("Dunes Rebuilt")
+        plt.ylabel("PEA_2011 Rebuilt")
         plt.xlabel("Time (yr)")
 
     # nourishment volumes