Detailed panel by panel look at coastaline of SSS150 and its river file
import cmocean.cm as cm
import copy
import matplotlib.pyplot as plt
import numpy as np
import xarray as xr
%matplotlib inline
cm_thermal = copy.copy(cm.thermal)
cm_thermal.set_bad('white')
cm_thermal.set_under('grey')
b_SSS150 = xr. open_dataset('../../../../grid/sss150/bathymetry_sss150_update202108210102AUM1.nc')
new_bathy = copy.copy(b_SSS150)
new_bathy
<xarray.Dataset> Dimensions: (x: 710, y: 826) Dimensions without coordinates: x, y Data variables: nav_lon (y, x) float32 ... nav_lat (y, x) float32 ... Bathymetry (y, x) float32 ...
[586460 values with dtype=float32]
[586460 values with dtype=float32]
[586460 values with dtype=float32]
r_SSS150 = xr.open_dataset('../../../../tools/I_ForcingFiles/Rivers/ncfiles/R_SSS150_Dailies_y2023m01d01.nc').rorunoff[0]
def make_the_plots(istart, iend, jstart, jend, bathy=b_SSS150):
fig, axs = plt.subplots(1, 3, figsize=(15, 5))
for ax in [axs[0], axs[2]]:
colours = ax.pcolormesh(np.arange(jstart, jend+1), np.arange(istart, iend+1),
b_SSS150.Bathymetry[istart:iend, jstart:jend],
cmap=cm_thermal, vmax=30, vmin=0.5);
axs[1].pcolormesh(np.arange(jstart, jend+1), np.arange(istart, iend+1),
b_SSS150.Bathymetry[istart:iend, jstart:jend]>0,
cmap='Greys');
axs[0].set_title('Bathymetry')
axs[1].set_title('Rivers')
axs[2].set_title('Both')
for i in range(istart, iend):
for j in range(jstart, jend):
if r_SSS150[i, j] > 0:
axs[1].plot(j+0.5, i+0.5, 'c+')
axs[2].plot(j+0.5, i+0.5, 'wx')
return fig, axs
istart, iend = 0, 826
jstart, jend = 0, 710
make_the_plots(istart, iend, jstart,jend)
istart, iend = 0, 50
jstart, jend = 0, 50
make_the_plots(istart, iend, jstart, jend)
istart, iend = 0, 50
jstart, jend = 150, 200
make_the_plots(istart, iend, jstart,jend)
istart, iend = 0, 50
jstart, jend = 200, 250
make_the_plots(istart, iend, jstart,jend)
istart, iend = 0, 50
jstart, jend = 230, 300
make_the_plots(istart, iend, jstart,jend)
istart, iend = 50, 100
jstart, jend = 200, 250
make_the_plots(istart, iend, jstart,jend)
istart, iend = 50, 100
jstart, jend = 250, 300
make_the_plots(istart, iend, jstart,jend)
istart, iend = 100, 150
jstart, jend = 199, 250
make_the_plots(istart, iend, jstart,jend)
istart, iend = 100, 150
jstart, jend = 250, 300
make_the_plots(istart, iend, jstart,jend)
istart, iend = 150, 200
jstart, jend = 150, 200
make_the_plots(istart, iend, jstart,jend)
istart, iend = 150, 200
jstart, jend = 200, 250
make_the_plots(istart, iend, jstart,jend)
istart, iend = 150, 200
jstart, jend = 250, 300
make_the_plots(istart, iend, jstart,jend)
istart, iend = 150, 200
jstart, jend = 300, 350
make_the_plots(istart, iend, jstart,jend)
istart, iend = 150, 200
jstart, jend = 350, 650
make_the_plots(istart, iend, jstart,jend)
istart, iend = 150, 200
jstart, jend = 650, 710
make_the_plots(istart, iend, jstart,jend)
istart, iend = 200, 250
jstart, jend = 150, 200
make_the_plots(istart, iend, jstart,jend)
istart, iend = 200, 250
jstart, jend = 200, 250
make_the_plots(istart, iend, jstart,jend)
istart, iend = 200, 250
jstart, jend = 250, 300
make_the_plots(istart, iend, jstart,jend)
istart, iend = 200, 250
jstart, jend = 300, 350
make_the_plots(istart, iend, jstart,jend)
istart, iend = 200, 250
jstart, jend = 350, 400
make_the_plots(istart, iend, jstart,jend)
istart, iend = 200, 250
jstart, jend = 550, 600
make_the_plots(istart, iend, jstart,jend)
istart, iend = 200, 250
jstart, jend = 600, 650
make_the_plots(istart, iend, jstart,jend)
istart, iend = 200, 250
jstart, jend = 650, 710
make_the_plots(istart, iend, jstart,jend)
istart, iend = 202, 212
jstart, jend = 700, 710
make_the_plots(istart, iend, jstart, jend)
istart, iend = 202, 212
jstart, jend = 700, 710
new_bathy.Bathymetry[203:211, -2:] = 0
make_the_plots(istart, iend, jstart, jend, bathy=new_bathy)
istart, iend = 250, 300
jstart, jend = 150, 220
make_the_plots(istart, iend, jstart,jend)
istart, iend = 250, 300
jstart, jend = 200, 250
make_the_plots(istart, iend, jstart,jend)
istart, iend = 250, 300
jstart, jend = 250, 300
make_the_plots(istart, iend, jstart,jend)
istart, iend = 230, 300
jstart, jend = 290, 350
make_the_plots(istart, iend, jstart,jend)
istart, iend = 250, 300
jstart, jend = 350, 400
make_the_plots(istart, iend, jstart,jend)
istart, iend = 250, 300
jstart, jend = 400, 450
make_the_plots(istart, iend, jstart,jend)
istart, iend = 250, 300
jstart, jend = 450, 500
make_the_plots(istart, iend, jstart,jend)
istart, iend = 250, 300
jstart, jend = 500, 550
make_the_plots(istart, iend, jstart,jend)
istart, iend = 250, 300
jstart, jend = 550, 600
make_the_plots(istart, iend, jstart,jend)
istart, iend = 250, 300
jstart, jend = 600, 650
make_the_plots(istart, iend, jstart,jend)
istart, iend = 300, 350
jstart, jend = 150, 200
make_the_plots(istart, iend, jstart,jend)
istart, iend = 300, 350
jstart, jend = 200, 250
make_the_plots(istart, iend, jstart,jend)
istart, iend = 300, 350
jstart, jend = 250, 300
make_the_plots(istart, iend, jstart,jend)
istart, iend = 300, 350
jstart, jend = 400, 450
make_the_plots(istart, iend, jstart,jend)
istart, iend = 300, 350
jstart, jend = 450, 500
make_the_plots(istart, iend, jstart,jend)
istart, iend = 350, 400
jstart, jend = 0, 50
make_the_plots(istart, iend, jstart,jend)
istart, iend = 350, 400
jstart, jend = 50, 100
make_the_plots(istart, iend, jstart,jend)
istart, iend = 350, 400
jstart, jend = 100, 150
make_the_plots(istart, iend, jstart,jend)
istart, iend = 350, 400
jstart, jend = 150, 200
make_the_plots(istart, iend, jstart,jend)
istart, iend = 350, 400
jstart, jend = 200, 250
make_the_plots(istart, iend, jstart,jend)
istart, iend = 350, 400
jstart, jend = 250, 300
make_the_plots(istart, iend, jstart,jend)
istart, iend = 350, 400
jstart, jend = 300, 370
make_the_plots(istart, iend, jstart,jend)
istart, iend = 350, 400
jstart, jend = 350, 400
make_the_plots(istart, iend, jstart,jend)
istart, iend = 350, 400
jstart, jend = 450, 500
make_the_plots(istart, iend, jstart,jend)
istart, iend = 395, 450
jstart, jend = 0, 50
make_the_plots(istart, iend, jstart,jend)
istart, iend = 400, 450
jstart, jend = 50, 100
make_the_plots(istart, iend, jstart,jend)
istart, iend = 400, 450
jstart, jend = 100, 150
make_the_plots(istart, iend, jstart,jend)
istart, iend = 400, 450
jstart, jend = 130, 200
make_the_plots(istart, iend, jstart,jend)
istart, iend = 400, 450
jstart, jend = 300, 350
make_the_plots(istart, iend, jstart,jend)
istart, iend = 400, 455
jstart, jend = 350, 400
make_the_plots(istart, iend, jstart,jend)
istart, iend = 400, 450
jstart, jend = 450, 500
make_the_plots(istart, iend, jstart,jend)
istart, iend = 400, 450
jstart, jend = 500, 550
make_the_plots(istart, iend, jstart,jend)
istart, iend = 450, 500
jstart, jend = 0, 50
make_the_plots(istart, iend, jstart,jend)
istart, iend = 450, 500
jstart, jend = 50, 100
make_the_plots(istart, iend, jstart,jend)
istart, iend = 450, 500
jstart, jend = 100, 150
make_the_plots(istart, iend, jstart,jend)
istart, iend = 450, 500
jstart, jend = 150, 200
make_the_plots(istart, iend, jstart,jend)
istart, iend = 450, 500
jstart, jend = 350, 400
make_the_plots(istart, iend, jstart,jend)
istart, iend = 450, 500
jstart, jend = 500, 550
make_the_plots(istart, iend, jstart,jend)
istart, iend = 450, 500
jstart, jend = 550, 600
make_the_plots(istart, iend, jstart,jend)
istart, iend = 500, 550
jstart, jend = 0, 50
make_the_plots(istart, iend, jstart,jend)
istart, iend = 500, 550
jstart, jend = 50, 100
make_the_plots(istart, iend, jstart,jend)
istart, iend = 500, 550
jstart, jend = 95, 150
make_the_plots(istart, iend, jstart,jend)
istart, iend = 500, 550
jstart, jend = 150, 200
make_the_plots(istart, iend, jstart,jend)
istart, iend = 500, 550
jstart, jend = 350, 400
make_the_plots(istart, iend, jstart,jend)
istart, iend = 500, 550
jstart, jend = 500, 550
make_the_plots(istart, iend, jstart,jend)
istart, iend = 500, 550
jstart, jend = 550, 600
make_the_plots(istart, iend, jstart,jend)
istart, iend = 550, 600
jstart, jend = 0, 50
make_the_plots(istart, iend, jstart,jend)
istart, iend = 550, 600
jstart, jend = 50, 100
make_the_plots(istart, iend, jstart,jend)
istart, iend = 550, 602
jstart, jend = 100, 150
make_the_plots(istart, iend, jstart,jend)
istart, iend = 550, 600
jstart, jend = 150, 200
make_the_plots(istart, iend, jstart,jend)
istart, iend = 550, 600
jstart, jend = 350, 400
make_the_plots(istart, iend, jstart,jend)
istart, iend = 550, 600
jstart, jend = 500, 550
make_the_plots(istart, iend, jstart,jend)
## Tile Twelve
* Pitt
* no rivers
istart, iend = 550, 600
jstart, jend = 550, 600
make_the_plots(istart, iend, jstart,jend)
istart, iend = 600, 650
jstart, jend = 0, 50
make_the_plots(istart, iend, jstart,jend)
istart, iend = 600, 650
jstart, jend = 50, 100
make_the_plots(istart, iend, jstart,jend)
istart, iend = 600, 650
jstart, jend = 95, 150
make_the_plots(istart, iend, jstart,jend)
istart, iend = 600, 650
jstart, jend = 145, 200
make_the_plots(istart, iend, jstart,jend)
istart, iend = 600, 650
jstart, jend = 500, 550
make_the_plots(istart, iend, jstart,jend)
istart, iend = 650, 700
jstart, jend = 100, 150
make_the_plots(istart, iend, jstart,jend)
istart, iend = 650, 700
jstart, jend = 150, 200
make_the_plots(istart, iend, jstart,jend)
istart, iend = 650, 700
jstart, jend = 500, 550
make_the_plots(istart, iend, jstart,jend)
istart, iend = 700, 750
jstart, jend = 100, 150
make_the_plots(istart, iend, jstart,jend)
istart, iend = 700, 750
jstart, jend = 150, 200
make_the_plots(istart, iend, jstart,jend)
## Tile Three
* Howe Sound
* river okay
istart, iend = 750, 800
jstart, jend = 100, 150
make_the_plots(istart, iend, jstart,jend)
## Tile Four
* Howe Sound
* rivers okay
istart, iend = 750, 800
jstart, jend = 150, 200
make_the_plots(istart, iend, jstart,jend)
## Tile Five-Fourteen
* Land
# Row Seventeen
## Tile One-Three
* Land
## Tile Four
* Squamish
istart, iend = 800, 825
jstart, jend = 170, 200
fig, axs = make_the_plots(istart, iend, jstart,jend)
axs[0].plot(176, 807, 'r*');
axs[0].plot(177, 807, 'co');
axs[0].plot(176, 806, 'co');
print (b_SSS150.Bathymetry[806, 176].values, b_SSS150.Bathymetry[807, 177].values,
b_SSS150.Bathymetry[807, 176].values)
4.0 4.0 0.0
istart, iend = 790, 825
jstart, jend = 170, 200
new_bathy.Bathymetry[807, 176] = 4
make_the_plots(istart, iend, jstart, jend, bathy=new_bathy);
new_bathy['history'] = 'From bathymetry_sss150_update202108210102AUM1.nc with Squamish River corrected and Fraser River closed on East'
new_bathy.to_netcdf('bathymetry_sss150_20240822.nc')