|
|
|
|
@ -159,7 +159,11 @@ def find_beach_slope(chain_list, z_ave_list, ref_index, check_value):
|
|
|
|
|
while i>0 and z_ave_list[i]<check_value+1:
|
|
|
|
|
i=i-1
|
|
|
|
|
|
|
|
|
|
# Hide numpy divide-by-zero warnings
|
|
|
|
|
np.seterr(all='ignore')
|
|
|
|
|
slope=(z_ave_list[i]-z_ave_list[ref_index])/(chain_list[i]-chain_list[ref_index])
|
|
|
|
|
np.seterr(all=None)
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
beach_ave=interpolate.interp1d([min(chain_list),max(chain_list)], [(min(chain_list)-chain_list[ref_index])*slope+z_ave_list[ref_index], (z_ave_list[ref_index]-(chain_list[ref_index]-max(chain_list))*slope)])
|
|
|
|
|
|
|
|
|
|
|
