import scipy as sp import imp imp.load_source('transect_analysis', '/home/563/esh563/goulburn_NT/transect_analysis.py') import transect_analysis as ta # Specify start and end coords on coast. # Choose order so that (transect_axis, coastline_axis) forms a right hand coordinate system lon0 = 134.5293 lat0 = -12.4715 coast_lon1 = 133.3290 coast_lat1 = -12.1468 month = 12 trans_lon0, trans_lat0, trans_lon1, trans_lat1, n_points, n_trans, coast_distances, tran_distances = ta.define_transects( lon0, lat0, coast_lon1, coast_lat1, 453300, 25000) static_path = '/g/data/ua8/ARCCSS_Data/MCASClimate/v1-0/static/static.nc' static = xr.open_dataset(static_path).sel(latitude=slice(-12.75, -6), longitude=slice(130, 139)) static_tran = ta.calc_transects(static, trans_lon0, trans_lat0, trans_lon1, trans_lat1, n_points, n_trans) static_tran = static_tran.assign_coords(coastal_axis=coast_distances) static_tran = static_tran.assign_coords(transect_axis=tran_distances) # Calcualate distance where landmask drops below 0.5 coast_i = np.where( static_tran.mean('coastal_axis').LANDMASK.values < 0.5)[0][0] - 1 coast_location = tran_distances[coast_i] # Redefine tran_distances so that coastline occurs at 0.
import pyproj as pp import scipy as sp import imp imp.load_source('transect_analysis', '/home/563/esh563/goulburn_NT/transect_analysis.py') import transect_analysis as ta # Specify start and end coords on coast. # Choose order so that (transect_axis, coastline_axis) forms a right hand coordinate system lon0 = 134.5293 lat0 = -12.4715 coast_lon1 = 133.3290 coast_lat1 = -12.1468 trans_lon0, trans_lat0, trans_lon1, trans_lat1, n_points, n_trans, coast_distances, tran_distances = ta.define_transects( lon0, lat0, coast_lon1, coast_lat1, 453300, spacing=8 * 10**3) static_path = '/g/data/ua8/ARCCSS_Data/MCASClimate/v1-0/static/static.nc' static = xr.open_dataset(static_path).sel(latitude=slice(-12.75, -6), longitude=slice(130, 139)) static_tran = ta.calc_transects(static, trans_lon0, trans_lat0, trans_lon1, trans_lat1, n_points, n_trans) static_tran = static_tran.assign_coords(coastal_axis=coast_distances) static_tran = static_tran.assign_coords(transect_axis=tran_distances) # Calcualate distance where landmask drops below 0.5 coast_i = np.where( static_tran.mean('coastal_axis').LANDMASK.values < 0.5)[0][0] - 1 coast_location = tran_distances[coast_i] # Redefine tran_distances so that coastline occurs at 0.