points_good = 1
level_elev, thresh, levpercent = ro.calc_level_ice(
elevation_km, pint, pwidth, min_ridge_height)
elevation_km = elevation_km - level_elev
thresh = thresh - level_elev
if (azi1 > 300):
mask_km = where((azi_km > azi1) | (azi_km < azi2))
else:
mask_km = where((azi_km > azi1) | (azi_km < azi2))
if (size(mask_km) > 100):
#MASK IF WEIRDLY DOESN'T HAVE ATM POINTS ON THE SIDE. THINK ABOUT THIS NUMBER MORE!
xatm1D1, yatm1D1, elevation1D1, xatm2m, yatm2m, odist2m, elevS2m, tptsN, tpts, dist_diffs = ro.atm_analysis_1D(
xatm_km,
yatm_km,
azi_km,
elevation_km,
min_ridge_height,
mask_km,
atrkres=atrkres)
peak_locs = odist2m[tptsN]
elevS2mN = elevS2m[tptsN]
xatm2mN = xatm2m[tptsN]
yatm2mN = yatm2m[tptsN]
print mean(dist_diffs)
#may need to check upper percentile distance?
if (size(tptsN) > 0):
found_ridges = 1
if (size(tptsN) > 1):
