# 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.
tran_distances = tran_distances - coast_location
# Create basis vectors of new coordinate system
b_lon = trans_lon1[0] - lon0
b_lat = trans_lat1[0] - lat0
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=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.
tran_distances = tran_distances - coast_location
# Create basis vectors of new coordinate system
b_lon = trans_lon1[0] - lon0
b_lat = trans_lat1[0] - lat0
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 = 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.
tran_distances = tran_distances - coast_location
# Create basis vectors of new coordinate system
b_lon = trans_lon1[0] - lon0
b_lat = trans_lat1[0] - lat0
# Create transects
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.
tran_distances = tran_distances - coast_location
# Create basis vectors of new coordinate system
b_lon = trans_lon1[0] - lon0
b_lat = trans_lat1[0] - lat0
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=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.
tran_distances = tran_distances - coast_location
# Create basis vectors of new coordinate system
b_lon = trans_lon1[0] - lon0
b_lat = trans_lat1[0] - lat0