Updating Hatteras scripts, developed pipelines for road and island offset

Author: HannahAline

data/hatteras_init/Hatteras-CASCADE-parameters.yaml

@@ -20,7 +20,7 @@ GrowthParamStart: false
 HdDiffu: 0.75
 Ki: 7.5e-06
 Kr: 7.5e-05
-LShoreface: 524.1959986328857
+LShoreface: 524.6532522894681
 MHW: 0.46
 MaxShrubHeight: 5.3
 MaxUpSlope: 0.25
@@ -39,14 +39,14 @@ Sandbag_elevation: 1.9
 SeededRNG: true
 Seedmax: 1000.0
 Seedmin: 100.0
-ShorefaceToe: 1099.9999999999995
+ShorefaceToe: 1452.7454838709673
 ShrubEl_max: 2.3
 ShrubEl_min: 1.2
 Shrub_ON: 0
 SprayDist: 170
 StormSeries: []
 StormStart: 2
-TMAX: 20
+TMAX: 19
 TideAmp: 1.2
 TimeFruit: 5.0
 UprootLimit: -0.2

data/hatteras_init/Hatteras-CASCADE-parameters_BACKUP.yaml

@@ -0,0 +1,69 @@
+BarrierLength: 500
+BayDepth: 3.0
+BermEl: 1.9
+BurialLimit: 0.75
+C1: 8.8
+C2: 4.6
+Cbb_i: 0.85
+Cbb_r: 0.7
+Cx: 10.0
+DShoreface: 8.9
+Dmaxel: 3.4
+Dshrub: 2.75
+Dstart: 0.5
+DuneParamStart: true
+DuneRestart: 0.075
+DuneWidth: 20.0
+Female: 0.5
+GermRate: 0.6
+GrowthParamStart: false
+HdDiffu: 0.75
+Ki: 7.5e-06
+Kr: 7.5e-05
+LShoreface: 524.3969226015213
+MHW: 0.46
+MaxShrubHeight: 5.3
+MaxUpSlope: 0.25
+OWss_i: 2
+OWss_r: 1
+Qs_bb_min: 1
+Qs_min: 1.0
+Qshrub_max: 0.15
+RSLR_Constant: true
+RSLR_const: 0.0056
+Rat: 0.0
+Rin_i: 0.1
+Rin_r: 2.0
+SalineLimit: 5
+Sandbag_elevation: 1.9
+SeededRNG: true
+Seedmax: 1000.0
+Seedmin: 100.0
+ShorefaceToe: 1699.3199999999995
+ShrubEl_max: 2.3
+ShrubEl_min: 1.2
+Shrub_ON: 0
+SprayDist: 170
+StormSeries: []
+StormStart: 2
+TMAX: 20
+TideAmp: 1.2
+TimeFruit: 5.0
+UprootLimit: -0.2
+beta: 0.04
+disp_mu: -0.721891
+disp_sigma: 1.5
+dune_file: C:\Users\hanna\PycharmProjects\CASCADE\data\hatteras_init\buffer\sample_1_dune.npy
+elevation_file: C:\Users\hanna\PycharmProjects\CASCADE\data\hatteras_init\buffer\sample_1_topography.npy
+enable_sandbags: false
+growth_param_file: C:\Users\hanna\PycharmProjects\CASCADE\data\default_cascade_variables\growthparam_1000dam.npy
+k_sf: 19324.753176079226
+mm: 2.0
+nn: 0.5
+rmax: 0.95
+rmin: 0.55
+s_sf_eq: 0.01698590034494604
+sandbag_elevation: 1.8
+sandbag_need: false
+storm_file: C:\Users\hanna\PycharmProjects\CASCADE\data\hatteras_init\storms\hindcast_storms\storms_1978_1997.npy
+threshold_in: 0.25

data/hatteras_init/duneline_change/duneline_change.py

@@ -0,0 +1,46 @@
+import pandas as pd
+
+# Input dune-line CSVs (raw, per-transect)
+dune1978 = pd.read_csv(
+    r"C:\Users\hanna\PycharmProjects\CASCADE\data\hatteras_init\island_offset\1978\1978_duneline_offset_raw.csv"
+)
+dune1997 = pd.read_csv(
+    r"C:\Users\hanna\PycharmProjects\CASCADE\data\hatteras_init\island_offset\1997\1997_duneline_offset_raw.csv"
+)
+
+# ---- 1. Aggregate each year to one value per domain (mean ORIG_LEN) ----
+d78 = (
+    dune1978
+    .groupby("domain_id")["ORIG_LEN"]
+    .mean()
+    .reset_index()
+    .rename(columns={"ORIG_LEN": "dune_1978_m"})
+)
+
+d97 = (
+    dune1997
+    .groupby("domain_id")["ORIG_LEN"]
+    .mean()
+    .reset_index()
+    .rename(columns={"ORIG_LEN": "dune_1997_m"})
+)
+
+# ---- 2. Merge per-domain means on domain_id ----
+df = d78.merge(d97, on="domain_id", how="inner")
+
+# ---- 3. Compute dune-line change (1997 - 1978) ----
+# NEW: positive = seaward / accretion, negative = landward / retreat
+df['dune_change_m'] = df['dune_1978_m'] - df['dune_1997_m']
+
+
+# ---- 4. Sort and save ----
+df = df.sort_values("domain_id")
+
+outpath = (
+    r"C:\Users\hanna\PycharmProjects\CASCADE\data\hatteras_init\duneline_change"
+    r"\dune_change_1978_1997_by_domain.csv"
+)
+df.to_csv(outpath, index=False)
+
+print("Saved dune-line change table to:", outpath)
+print(df.head())