import numpy as np
import netCDF4 as nc
import xarray as xr
%matplotlib inline
import matplotlib.pyplot as plt
control = xr.open_dataset('/data/vdo/MEOPAR/completed-runs/VAT21/SalishSea_1h_20150507_20150514_grid_T_20150507-20150507.nc')
VAT20 = xr.open_dataset('/data/vdo/MEOPAR/completed-runs/VAT20/SalishSea_1h_20150507_20150514_grid_T_20150507-20150507.nc')
VAT24 = xr.open_dataset('/data/vdo/MEOPAR/completed-runs/VAT24/SalishSea_1h_20150507_20150514_grid_T_20150507-20150507.nc')
VAT25 = xr.open_dataset('/data/vdo/MEOPAR/completed-runs/VAT25/SalishSea_1h_20150507_20150509_grid_T_20150507-20150507.nc')
control2 = xr.open_dataset('/data/vdo/MEOPAR/completed-runs/VAT25/SalishSea_1h_20150507_20150509_grid_V_20150508-20150508.nc')
VAT252 = xr.open_dataset('/data/vdo/MEOPAR/completed-runs/VAT25/SalishSea_1h_20150507_20150509_grid_V_20150508-20150508.nc')
figs, (ax1,ax2) = plt.subplots(1,2,sharey=True,figsize = ((20,8)))
control.vosaline.isel(deptht = 0, time_counter = 10).plot(ax=ax1)
VAT20.vosaline.isel(deptht = 0, time_counter = 10 ).plot(ax=ax2)
#VAT24.vosaline.isel(deptht = 0,time_counter = 10).plot(ax=ax3)
<matplotlib.collections.QuadMesh at 0x7fcd0af82cc0>
(VAT20.vosaline.isel(deptht = 0, time_counter = 10 ) - control.vosaline.isel(deptht = 0, time_counter = 10)).plot(figsize=(10,8))
<matplotlib.collections.QuadMesh at 0x7fcd1b0a7048>
figs, (ax1,ax2) = plt.subplots(1,2,sharey=True,figsize = ((20,8)))
control.vosaline.isel(deptht = 0, time_counter = 10).plot(ax=ax1)
VAT24.vosaline.isel(deptht = 0,time_counter = 10).plot(ax=ax2)
<matplotlib.collections.QuadMesh at 0x7fcd08c08940>
(VAT24.vosaline.isel(deptht = 0, time_counter = 10 ) - control.vosaline.isel(deptht = 0, time_counter = 10)).plot(figsize=(10,8))
<matplotlib.collections.QuadMesh at 0x7fcd0748f198>
figs, (ax1,ax2,ax3) = plt.subplots(1,3,sharey=True,figsize = ((10,8)))
control.vosaline.isel(deptht = 0, time_counter = 23).plot(ax=ax1)
VAT20.vosaline.isel(deptht = 0, time_counter = 23).plot(ax=ax2)
VAT24.vosaline.isel(deptht = 0,time_counter = 23).plot(ax=ax3)
<matplotlib.collections.QuadMesh at 0x7fcd34a06860>
difference1 = VAT20.vosaline.isel(deptht = 0, time_counter = 23) - control.vosaline.isel(deptht = 0, time_counter = 23)
difference1.plot()
<matplotlib.collections.QuadMesh at 0x7fcd2d443a58>
difference2 = VAT24.vosaline.isel(deptht = 0, time_counter = 23) - control.vosaline.isel(deptht = 0, time_counter = 23)
difference2.plot()
<matplotlib.collections.QuadMesh at 0x7fcd2c7f4f28>
figs, (ax1,ax2,ax3) = plt.subplots(1,3,sharey=True,figsize = ((10,8)))
control.votemper.isel(deptht = 0, time_counter = 23).plot(ax=ax1)
VAT20.votemper.isel(deptht = 0, time_counter = 23).plot(ax=ax2)
VAT24.votemper.isel(deptht = 0,time_counter = 23).plot(ax=ax3)
<matplotlib.collections.QuadMesh at 0x7fcd2ba142e8>
difference3 = VAT20.votemper.isel(deptht = 0, time_counter = 23) - control.votemper.isel(deptht = 0, time_counter = 23)
difference3.plot(figsize=(10,8))
<matplotlib.collections.QuadMesh at 0x7fcd05c7d630>
difference4 = VAT24.votemper.isel(deptht = 0, time_counter = 23) - control.votemper.isel(deptht = 0, time_counter = 23)
difference4.plot(figsize=(10,8))
<matplotlib.collections.QuadMesh at 0x7fcd050377f0>
fig,(ax1,ax2) = plt.subplots(1,2,figsize=(10,8))
control.vosaline.isel(deptht = 0, time_counter = 23).plot(ax=ax1)
VAT25.vosaline.isel(deptht = 0, time_counter = 23).plot(ax=ax2)
<matplotlib.collections.QuadMesh at 0x7fcd04368cc0>
difference5 = VAT25.vosaline.isel(deptht = 0, time_counter = 23) - control.vosaline.isel(deptht = 0, time_counter = 23)
difference5.plot(figsize=(10,8))
<matplotlib.collections.QuadMesh at 0x7fcd0ecce588>
figs,(ax1,ax2)=plt.subplots(1,2,figsize=(10,8))
control.sossheig.isel(time_counter = 23).plot(ax=ax1)
VAT25.sossheig.isel(time_counter = 23).plot(ax=ax2)
<matplotlib.collections.QuadMesh at 0x7fcd02b4a470>
difference6 = VAT25.sossheig.isel(time_counter = 23) - control.sossheig.isel(time_counter = 23)
difference6.plot(figsize=(10,8))
<matplotlib.collections.QuadMesh at 0x7fcd0f91b048>
figs,(ax1,ax2)=plt.subplots(1,2,figsize=(10,8))
control2.vomecrty.isel(time_counter = 23,depthv =0).plot(ax=ax1)
VAT252.vomecrty.isel(time_counter = 23,depthv = 0).plot(ax=ax2)
<matplotlib.collections.QuadMesh at 0x7fcd013134e0>
difference7 = VAT252.vomecrty.isel(time_counter = 23,depthv = 0) - control2.vomecrty.isel(time_counter = 23,depthv =0)
difference7.plot()
<matplotlib.collections.QuadMesh at 0x7fcd235b1cc0>
control3 = xr.open_dataset('/data/vdo/MEOPAR/completed-runs/VAT21/SalishSea_1h_20150507_20150514_grid_T_20150508-20150508.nc')
VAT253 = xr.open_dataset('/data/vdo/MEOPAR/completed-runs/VAT25/SalishSea_1h_20150507_20150509_grid_T_20150508-20150508.nc')
figs,(ax1,ax2) = plt.subplots(1,2,figsize=(10,8))
control3.sossheig.isel(time_counter=23).plot(ax=ax1)
VAT253.sossheig.isel(time_counter=23).plot(ax=ax2)
<matplotlib.collections.QuadMesh at 0x7fccffb70d30>
difference8 = VAT253.sossheig.isel(time_counter=23) - control3.sossheig.isel(time_counter=23)
difference8.plot(figsize=(10,8))
<matplotlib.collections.QuadMesh at 0x7fcd104dce48>
np.abs(difference8).max()
<xarray.DataArray 'sossheig' ()> array(0.035715579986572266) Coordinates: time_centered datetime64[ns] 2015-05-08T23:30:00 time_counter datetime64[ns] 2015-05-08T23:30:00
figs,(ax1,ax2)=plt.subplots(1,2,figsize=(10,8))
control3.vosaline.isel(time_counter=23,deptht=0).plot(ax=ax1)
VAT253.vosaline.isel(time_counter=23,deptht=0).plot(ax=ax2)
<matplotlib.collections.QuadMesh at 0x7fccfe367c18>
difference9 = VAT253.vosaline.isel(time_counter=23,deptht=0)-control3.vosaline.isel(time_counter=23,deptht=0)
difference9.plot(figsize=(10,8))
<matplotlib.collections.QuadMesh at 0x7fcd11122da0>
fig,(ax1,ax2)=plt.subplots(1,2,figsize=(10,8))
control3.votemper.isel(time_counter=23,deptht=0).plot(ax=ax1)
VAT253.votemper.isel(time_counter=23,deptht=0).plot(ax=ax2)
<matplotlib.collections.QuadMesh at 0x7fcd0429d7f0>
difference10 = VAT253.votemper.isel(time_counter=23,deptht=0)-control3.votemper.isel(time_counter=23,deptht=0)
difference10.plot()
<matplotlib.collections.QuadMesh at 0x7fcd210f2438>
np.abs(VAT253.votemper-control3.votemper).max()
<xarray.DataArray 'votemper' ()> array(3.087312698364258)
np.abs(VAT253.vosaline-control3.vosaline).max()
<xarray.DataArray 'vosaline' ()> array(13.474448204040527)
np.abs(VAT253.sossheig-control3.sossheig).max()
<xarray.DataArray 'sossheig' ()> array(0.07220423221588135)
VAT20a = xr.open_dataset('/data/vdo/MEOPAR/completed-runs/VAT20/SalishSea_1h_20150507_20150514_grid_T_20150513-20150513.nc')
controla = xr.open_dataset('/data/vdo/MEOPAR/completed-runs/VAT21/SalishSea_1h_20150507_20150514_grid_T_20150513-20150513.nc')
fig,(ax1,ax2)= plt.subplots(1,2,figsize=((10,8)))
controla.vosaline.isel(deptht=0,time_counter=23).plot(ax=ax1)
VAT20a.vosaline.isel(deptht = 0,time_counter = 23).plot(ax=ax2)
<matplotlib.collections.QuadMesh at 0x7fcd0d3c8f28>
difference01 = VAT20a.vosaline.isel(deptht = 0,time_counter = 23) - controla.vosaline.isel(deptht=0,time_counter=23)
difference01.plot(figsize=(10,8))
<matplotlib.collections.QuadMesh at 0x7fcd0bc520b8>