import xarray as xr
import matplotlib.pyplot as plt
%matplotlib inline
import numpy as np
VAT25 = xr.open_dataset('/data/vdo/MEOPAR/completed-runs/VAT25/SalishSea_1h_20150507_20150509_grid_T_20150508-20150508.nc')
VAT27 = xr.open_dataset('/data/vdo/MEOPAR/completed-runs/VAT27/SalishSea_1h_20150507_20150509_grid_T_20150508-20150508.nc')
fig,(ax1,ax2) = plt.subplots(1,2,figsize=((10,8)))
VAT25.sossheig.isel(time_counter=23).plot(ax=ax1,vmin = -3, vmax=3, cmap = 'bwr')
VAT27.sossheig.isel(time_counter=23).plot(ax=ax2,vmin = -3, vmax=3, cmap = 'bwr')
<matplotlib.collections.QuadMesh at 0x7f39a1ff3588>
(VAT25.sossheig.isel(time_counter=23) - VAT27.sossheig.isel(time_counter=23)).plot(figsize=((10,8)),
vmin = -0.02,
vmax = 0.02,
cmap = 'bwr')
<matplotlib.collections.QuadMesh at 0x7f399ff2e7f0>
VAT26 = xr.open_dataset('/data/vdo/MEOPAR/completed-runs/VAT26/SalishSea_1h_20150507_20150509_grid_T_20150508-20150508.nc')
VAT28 = xr.open_dataset('/data/vdo/MEOPAR/completed-runs/VAT28/SalishSea_1h_20150507_20150509_grid_T_20150508-20150508.nc')
fig,(ax1,ax2) = plt.subplots(1,2,figsize=((10,8)))
VAT26.sossheig.isel(time_counter=23).plot(ax=ax1,vmin = -3, vmax=3, cmap = 'bwr')
VAT28.sossheig.isel(time_counter=23).plot(ax=ax2,vmin = -3, vmax=3, cmap = 'bwr')
<matplotlib.collections.QuadMesh at 0x7f67534cc3c8>
(VAT26.sossheig.isel(time_counter=23) - VAT28.sossheig.isel(time_counter=23)).plot(figsize=((10,8)),
vmin = -0.02,
vmax = 0.02,
cmap = 'bwr')
<matplotlib.collections.QuadMesh at 0x7f675137ccf8>