Skip to content

Commit

Permalink
New beach nourishment functions and storms
Browse files Browse the repository at this point in the history
  • Loading branch information
franklin1895 committed Nov 11, 2024
1 parent 8eea9ef commit 0c8781c
Show file tree
Hide file tree
Showing 108 changed files with 1,527 additions and 72 deletions.
12 changes: 6 additions & 6 deletions data/Ocracoke_init_data/Ocracoke-CASCADE-parameters.yaml
Original file line number Diff line number Diff line change
Expand Up @@ -20,7 +20,7 @@ GrowthParamStart: false
HdDiffu: 0.75
Ki: 7.5e-06
Kr: 7.5e-05
LShoreface: 523.9734307609483
LShoreface: 524.8348346135758
MHW: 0.46
MaxShrubHeight: 5.3
MaxUpSlope: 0.25
Expand All @@ -31,15 +31,15 @@ Qs_min: 1.0
Qshrub_max: 0.15
RSLR_Constant: true
RSLR_const: 0.0056
Rat: 0
Rat: -2.0
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: 2099.9999999999995
ShrubEl_max: 2.3
ShrubEl_min: 1.2
Shrub_ON: 0
Expand All @@ -53,8 +53,8 @@ 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
elevation_file: C:\Users\frank\PycharmProjects\CASCADE\data\Ocracoke_init_data\elevations\Topography_2019\Domain_49_topography_2019.npy
dune_file: C:\Users\frank\PycharmProjects\CASCADE\data\Ocracoke_init_data\dunes\Dunes_2019\Domain_41_dune_2019.npy
elevation_file: C:\Users\frank\PycharmProjects\CASCADE\data\Ocracoke_init_data\elevations\Topography_2019\Domain_41_topography_2019.npy
enable_sandbags: false
growth_param_file: growthparam_1000dam.npy
k_sf: 19324.753176079226
Expand All @@ -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_9_baseline.npy
storm_file: C:\Users\frank\PycharmProjects\CASCADE\data\Ocracoke_init_data\storms\Synthetic_Storms\OCR_Future_StormList_2_baseline.npy
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
4 changes: 2 additions & 2 deletions scripts/ocracoke_ms/Forward_Simulations_Mac.py
Original file line number Diff line number Diff line change
Expand Up @@ -11,7 +11,7 @@
os.chdir('/Users/ceclmac/PycharmProjects/CASCADE')

# Set the number of years to simulate
run_years = 100
run_years = 126

# Set the start year
start_year = 2024
Expand Down Expand Up @@ -41,7 +41,7 @@
# Specify the number of B3D domains and control the number of domains
buffer_enabled = True
island_grid_number = 39
Total_B3D_Number = 70
Total_B3D_Number = 69

# Set whether inlet is accretional or erosional:
source_sink_load_name = '/Users/ceclmac/PycharmProjects/CASCADE/data/Ocracoke_init_data/Source_Sink/S32.csv'
Expand Down
73 changes: 23 additions & 50 deletions scripts/ocracoke_ms/Hindcast_Test_Script.py
Original file line number Diff line number Diff line change
Expand Up @@ -31,47 +31,18 @@
os.chdir('C:\\Users\\frank\\PycharmProjects\\CASCADE\\Run_output')

run_name_batch = []
'''
# Changed Value
run_name_batch.append('NCB_S32_SR_NR_1974_F')
run_name_batch.append('OCR_S32_SR_NR_1974_F')
run_name_batch.append('Duck_S32_SR_NR_1974_F')
# Original Value
run_name_batch.append('NCB_S32_MR_OR_1974')
run_name_batch.append('OCR_1_0_S32_MR_OR_1974')
run_name_batch.append('Duck_S32_SR_OR_1974')
run_name_batch.append('OCR_30_S38_SR_NR_NS_1974_Final')
run_name_batch.append('OCR_30_S38_SR_NR_NS_1997_Final')
run_name_batch.append('OCR_S32_SR_NR_NS_1974_Final')
run_name_batch.append('OCR_S32_SR_NR_NS_1997_Final')

run_name_batch.append('OCR_40_S38_SR_NR_NS_1974_Final')
run_name_batch.append('OCR_40_S38_SR_NR_NS_1997_Final')
run_name_batch.append('OCR_I_No_Management_ST0_Sink2')
run_name_batch.append('OCR_I_Status_Quo_ST0_Sink2')
#run_name_batch.append('OCR_I_Nourishment_S0_Accretional_Sink_Test_N3')

run_name_batch.append('OCR_40_S39_SR_NR_NS_1974_Final')
run_name_batch.append('OCR_40_S39_SR_NR_NS_1997_Final')
'''

run_name_batch.append('OCR_S32_SR_NR_NS_1974_Threshold_1_75_1_0')
run_name_batch.append('OCR_S32_SR_NR_NS_1974_Threshold_1_75_1_1')
run_name_batch.append('OCR_S32_SR_NR_NS_1974_Threshold_1_75_1_2')
run_name_batch.append('OCR_S32_SR_NR_NS_1974_Threshold_1_75_1_3')
run_name_batch.append('OCR_S32_SR_NR_NS_1974_Threshold_1_75_1_4')
run_name_batch.append('OCR_S32_SR_NR_NS_1974_Threshold_1_75_1_5')

run_name_batch.append('OCR_S32_SR_NR_NS_1997_Threshold_1_75_1_0')
run_name_batch.append('OCR_S32_SR_NR_NS_1997_Threshold_1_75_1_1')
run_name_batch.append('OCR_S32_SR_NR_NS_1997_Threshold_1_75_1_2')
run_name_batch.append('OCR_S32_SR_NR_NS_1997_Threshold_1_75_1_3')
run_name_batch.append('OCR_S32_SR_NR_NS_1997_Threshold_1_75_1_4')
run_name_batch.append('OCR_S32_SR_NR_NS_1997_Threshold_1_75_1_5')
#run_name_batch.append('OCR_I_Nourishment_S0_Erosional_Sink_Test_N1')
#run_name_batch.append('OCR_I_Nourishment_S0_Erosional_Sink_Test_N2')
#run_name_batch.append('OCR_I_Nourishment_S0_Erosional_Sink_Test_N3')

#nt_run = 46
#nt_run = 32
nt_run = 23
nt_run = 95
number_barrier3d_models = 70
buffer_length = 15
All_EP_Change = []
Expand Down Expand Up @@ -100,7 +71,6 @@
ny = np.size(b3d)

b3d_list.append(copy.deepcopy(b3d))
directory = "C:\\Users\\frank\\PycharmProjects\\CASCADE\\"
# TMax_MGMT = Needed 0
# TMAX_Sim = Last simulation year of the model 99
TMax_Sim = nt_run # Give length of simulation
Expand All @@ -115,7 +85,7 @@
all_shoreline_change = cascade._brie_coupler.brie.x_s_save

All_Year_1_Shoreline_Position = all_shoreline_change[:,1]
All_Year_44_Shoreline_Position = all_shoreline_change[:,-1]
All_Year_44_Shoreline_Position = all_shoreline_change[:,nt_run]

Year_1_Shoreline_Positions = All_Year_1_Shoreline_Position[buffer_length:-buffer_length]
Year_1_Shoreline_Positions[0] = 1624
Expand Down Expand Up @@ -244,36 +214,39 @@

# Plot all the different OCR storm intensities

plt.plot(domain_nums, LRR_74_97, label = 'Historic Change',color='grey')
#plt.plot(domain_nums, LRR_74_97, label = 'Historic Change',color='grey')
#plt.plot(domain_nums, LRR_97_20, label = 'Historic Change',color='grey')

plt.axhline(y = 0, color = 'k', linestyle = '--')
#plt.plot(domain_nums, All_EP_Change[0], label= '0%',color='#1f77b4')
#plt.plot(domain_nums, All_EP_Change[1], label= '10%',color='#ff7f0e')
plt.plot(domain_nums,All_EP_Change[2], label = '20%', color = '#2ca02c')
plt.plot(domain_nums, All_EP_Change[0], label= 'Natural',color='#1f77b4')
plt.plot(domain_nums, All_EP_Change[1], label= 'Status Quo',color='#ff7f0e')
#plt.plot(domain_nums, All_EP_Change[2], label = '3', color = '#2ca02c')
#plt.plot(domain_nums,All_EP_Change[3], label = '30%', color = '#d62728')
#plt.plot(domain_nums,All_EP_Change[4], label = '40%', color = '#9467bd')
#plt.plot(domain_nums,All_EP_Change[5], label = '50%', color = '#8c564b')

plt.legend()
plt.title('Historic vs Modeled Change: 1974-1997')
plt.tight_layout()
plt.title('Beach Nourishment Values: I')
plt.ylabel('Shoreline Change Rate (m/yr)')
plt.xlabel('B3D Domain')
plt.show()

plt.plot(domain_nums, LRR_97_20, label = 'Historic Change',color='grey')

print('Hello Purr')
'''
#plt.plot(domain_nums, LRR_97_20, label = 'Historic Change',color='grey')
plt.axhline(y = 0, color = 'k', linestyle = '--')
#plt.plot(domain_nums, All_EP_Change[6], label= '0%',color='#1f77b4')
#plt.plot(domain_nums, All_EP_Change[7], label= '10%',color='#ff7f0e')
plt.plot(domain_nums,All_EP_Change[8], label = '20%', color = '#2ca02c')
plt.plot(domain_nums, All_EP_Change[3], label= '1',color='#1f77b4')
plt.plot(domain_nums, All_EP_Change[4], label= '2',color='#ff7f0e')
plt.plot(domain_nums, All_EP_Change[5], label = '3', color = '#2ca02c')
#plt.plot(domain_nums,All_EP_Change[9], label = '30%', color = '#d62728')
#plt.plot(domain_nums,All_EP_Change[10], label = '40%', color = '#9467bd')
#plt.plot(domain_nums,All_EP_Change[11], label = '50%', color = '#8c564b')
#plt.plot(domain_nums,All_EP_Change[13], label = '2.0')#, color = '#1f77b4')
plt.legend()
plt.title('Historic vs Modeled Change: 1997-2020')
plt.tight_layout()
plt.title('Beach Nourishment Values: I')
plt.ylabel('Shoreline Change Rate (m/yr)')
plt.xlabel('B3D Domain')
plt.show()
Expand Down Expand Up @@ -316,7 +289,7 @@
plt.xlabel('B3D Domain')
plt.show()
'''
# Plot Ocracoke Storms
#plt.plot(domain_nums, LRR_74_97, label = 'Historic Change',color='grey')
plt.plot(domain_nums, LRR_97_20, label = 'Historic Change',color='grey')
Expand Down
Loading

0 comments on commit 0c8781c

Please sign in to comment.