Notebook to examined 3.6 test runs
In [1]:
import netCDF4 as nc
import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline
In [2]:
d = nc.Dataset('/data/nsoontie/MEOPAR/SalishSea/results/2Ddomain/3.4-compare/SalishSea_1h_20020915_20020921_grid_T.nc')
sal = d.variables['vosaline'][:]
sal = np.ma.masked_values(sal,0)
deps = d.variables['deptht'][:]
temp = d.variables['votemper'][:]
temp = np.ma.masked_values(temp,0)
ssh = d.variables['sossheig'][:]
d = nc.Dataset('/data/nsoontie/MEOPAR/SalishSea/results/2Ddomain/3.4-compare/SalishSea_1h_20020915_20020921_grid_U.nc')
U = d.variables['vozocrtx'][:]
U = np.ma.masked_values(U,0)
depsU = d.variables['depthu'][:]
In [3]:
smin=22
smax=34
t=0
y=2
plt.pcolormesh(np.arange(0,sal.shape[3]),deps,sal[t,:,y,:],vmin=smin,vmax=smax)
plt.axis([0,1100,400,0])
plt.colorbar()
Out[3]:
<matplotlib.colorbar.Colorbar instance at 0x7fed8536ba28>
In [4]:
t=60
plt.pcolormesh(np.arange(0,sal.shape[3]),deps,sal[t,:,y,:],vmin=smin,vmax=smax)
plt.axis([0,1100,400,0])
plt.colorbar()
Out[4]:
<matplotlib.colorbar.Colorbar instance at 0x7fed854b5128>
In [5]:
t=120
plt.pcolormesh(np.arange(0,sal.shape[3]),deps,sal[t,:,y,:],vmin=smin,vmax=smax)
plt.axis([0,1100,400,0])
plt.colorbar()
Out[5]:
<matplotlib.colorbar.Colorbar instance at 0x7fed854c1c20>
In [6]:
tmin=3; tmax=10
t=0
plt.pcolormesh(np.arange(0,sal.shape[3]),deps,temp[t,:,y,:],vmin=tmin,vmax=tmax)
plt.axis([0,1100,400,0])
plt.colorbar()
Out[6]:
<matplotlib.colorbar.Colorbar instance at 0x7fed7ef7bea8>
In [7]:
t=120
plt.pcolormesh(np.arange(0,sal.shape[3]),deps,temp[t,:,y,:],vmin=tmin,vmax=tmax)
plt.axis([0,1100,400,0])
plt.colorbar()
Out[7]:
<matplotlib.colorbar.Colorbar instance at 0x7fed7ee4e128>
In [8]:
t=0
umin=-.5
umax=0.5
plt.pcolormesh(np.arange(0,U.shape[3]),depsU,U[t,:,y,:],vmin=umin,vmax=umax)
plt.axis([0,1100,400,0])
plt.colorbar()
Out[8]:
<matplotlib.colorbar.Colorbar instance at 0x7fed7e595830>
In [9]:
t=100
plt.pcolormesh(np.arange(0,U.shape[3]),depsU,U[t,:,y,:],vmin=umin,vmax=umax)
plt.axis([0,1100,400,0])
plt.colorbar()
Out[9]:
<matplotlib.colorbar.Colorbar instance at 0x7fed7e47cab8>
In [10]:
t=120
plt.pcolormesh(np.arange(0,U.shape[3]),depsU,U[t,:,y,:],vmin=umin,vmax=umax)
plt.axis([0,1100,400,0])
plt.colorbar()
Out[10]:
<matplotlib.colorbar.Colorbar instance at 0x7fed7e2f1128>
Tidal signal? Look at boundary
In [11]:
plt.plot(U[:,0,y,3])
Out[11]:
[<matplotlib.lines.Line2D at 0x7fed7e1c5710>]
I don't see a tidal signal at all. ssh?
In [12]:
plt.plot(ssh[:,y,10])
Out[12]:
[<matplotlib.lines.Line2D at 0x7fed7e102d10>]
The amplitude seemse reall small but this looks like a tidal response. Howevever, the frequency is weird. Shoudnt there be two mins a day? this looks like 3 or 4.
In [14]:
plt.pcolormesh(ssh[:,y,:],vmin=-3,vmax=3)
plt.colorbar()
Out[14]:
<matplotlib.colorbar.Colorbar instance at 0x7fed7de92cb0>
In [16]:
plt.pcolormesh(U[:,0,y,:],vmin=-3,vmax=3)
plt.colorbar()
Out[16]:
<matplotlib.colorbar.Colorbar instance at 0x7fed7c389e18>
Ok, I definitely have tides in 3.4 but they aren't there in 3.6... Check namelists.
In [ ]: