Plots from abort and init files from changing 1_namelist.agrif test of Haro Strait 1
import numpy as np
import netCDF4 as nc
import xarray as xr
import matplotlib.pyplot as plt
import numpy.ma as ma
%matplotlib inline
agrif_abort = nc.Dataset('/data/vdo/MEOPAR/completed-runs/changing_agrif/1_output.abort.nc')
agrif_init = nc.Dataset('/data/vdo/MEOPAR/completed-runs/changing_agrif/1_output.init.nc')
agrif_mask = nc.Dataset('/data/vdo/MEOPAR/completed-runs/changing_agrif/1_mesh_mask.nc')
fig,(ax1,ax2) = plt.subplots(1,2,figsize=((10,8)))
s = ma.masked_array(agrif_abort.variables['vosaline'][0,27,:,:], mask = 1 - agrif_mask.variables['tmask'][0,27,:,:])
mesh = ax1.pcolormesh(s)
fig.colorbar(mesh,ax=ax1)
ax1.plot(121,532,'ro')
ax1.title('')
s = ma.masked_array(agrif_abort.variables['vosaline'][0,:,532,:], mask = 1 - agrif_mask.variables['tmask'][0,:,532,:])
mesh = ax2.pcolormesh(s)
ax2.set_ylim((32,0))
fig.colorbar(mesh,ax=ax2)
ax2.plot(121,27,'ro')
[<matplotlib.lines.Line2D at 0x7f21a3e900b8>]
fig,(ax1,ax2) = plt.subplots(1,2,figsize=((10,8)))
s = ma.masked_array(agrif_abort.variables['vosaline'][0,0,:,:], mask = 1 - agrif_mask.variables['tmask'][0,0,:,:])
mesh = ax1.pcolormesh(s,#vmin = 15, vmax = 33
)
fig.colorbar(mesh,ax=ax1)
ax1.plot(121,532,'ro')
s = ma.masked_array(agrif_abort.variables['vosaline'][0,:,532,:], mask = 1 - agrif_mask.variables['tmask'][0,:,532,:])
mesh = ax2.pcolormesh(s)
ax2.set_ylim((32,0))
fig.colorbar(mesh,ax=ax2)
ax2.plot(121,27,'ro')
[<matplotlib.lines.Line2D at 0x7f21a3ba72e8>]
fig,(ax1,ax2) = plt.subplots(1,2,figsize=((10,8)))
s = ma.masked_array(agrif_abort.variables['votemper'][0,27,:,:], mask = 1 - agrif_mask.variables['tmask'][0,27,:,:])
mesh = ax1.pcolormesh(s,#vmin = 15, vmax = 33
)
fig.colorbar(mesh,ax=ax1)
ax1.plot(121,532,'ro')
s = ma.masked_array(agrif_abort.variables['votemper'][0,:,532,:], mask = 1 - agrif_mask.variables['tmask'][0,:,532,:])
mesh = ax2.pcolormesh(s)
ax2.set_ylim((32,0))
fig.colorbar(mesh,ax=ax2)
ax2.plot(121,27,'ro')
[<matplotlib.lines.Line2D at 0x7f21ab3b9cf8>]
fig,(ax1,ax2) = plt.subplots(1,2,figsize=((10,8)))
s = ma.masked_array(agrif_abort.variables['vozocrtx'][0,27,:,:], mask = 1 - agrif_mask.variables['umask'][0,27,:,:])
mesh = ax1.pcolormesh(s,#vmin = 15, vmax = 33
)
fig.colorbar(mesh,ax=ax1)
ax1.plot(121,532,'ro')
s = ma.masked_array(agrif_abort.variables['vozocrtx'][0,:,532,:], mask = 1 - agrif_mask.variables['umask'][0,:,532,:])
mesh = ax2.pcolormesh(s)
ax2.set_ylim((32,0))
fig.colorbar(mesh,ax=ax2)
ax2.plot(121,27,'ro')
[<matplotlib.lines.Line2D at 0x7f21a3874630>]
fig,(ax1,ax2) = plt.subplots(1,2,figsize=((10,8)))
s = ma.masked_array(agrif_abort.variables['vomecrty'][0,27,:,:], mask = 1 - agrif_mask.variables['vmask'][0,27,:,:])
mesh = ax1.pcolormesh(s,#vmin = 15, vmax = 33
)
fig.colorbar(mesh,ax=ax1)
ax1.plot(121,532,'ro')
s = ma.masked_array(agrif_abort.variables['vomecrty'][0,:,532,:], mask = 1 - agrif_mask.variables['vmask'][0,:,532,:])
mesh = ax2.pcolormesh(s)
ax2.set_ylim((32,0))
fig.colorbar(mesh,ax=ax2)
ax2.plot(121,27,'ro')
[<matplotlib.lines.Line2D at 0x7f21a35fb978>]
fig,ax = plt.subplots(1,1,figsize=((10,8)))
ssh = ma.masked_array(agrif_abort.variables['sossheig'][0,:,:],mask = 1-agrif_mask.variables['tmask'][0,0,:,:])
mesh = ax.pcolormesh(ssh,vmin = 0, vmax = 1.4)
fig.colorbar(mesh)
<matplotlib.colorbar.Colorbar at 0x7f21a2f3df98>
expanded = nc.Dataset('/data/vdo/MEOPAR/completed-runs/expanded_again/1_output.abort.nc')
expanded_mask = nc.Dataset('/data/vdo/MEOPAR/completed-runs/expanded_again/1_mesh_mask.nc')
fig,ax = plt.subplots(1,1,figsize=((10,8)))
ssh = ma.masked_array(expanded.variables['sossheig'][0,:,:],mask = 1-expanded_mask.variables['tmask'][0,0,:,:])
mesh = ax.pcolormesh(ssh,vmin = 0, vmax = 1.4)
fig.colorbar(mesh)
<matplotlib.colorbar.Colorbar at 0x7f21a2e88438>
fig,(ax1,ax2) = plt.subplots(1,2,figsize=((10,8)))
s = ma.masked_array(expanded.variables['votemper'][0,0,:,:], mask = 1 - expanded_mask.variables['tmask'][0,0,:,:])
mesh = ax1.pcolormesh(s)
fig.colorbar(mesh,ax=ax1)
ax1.plot(121,532,'ro')
#ax1.title('')
s = ma.masked_array(expanded.variables['votemper'][0,:,532,:], mask = 1 - expanded_mask.variables['tmask'][0,:,532,:])
mesh = ax2.pcolormesh(s)
ax2.set_ylim((32,0))
fig.colorbar(mesh,ax=ax2)
ax2.plot(121,27,'ro')
[<matplotlib.lines.Line2D at 0x7f21a286fb70>]
expanded2 = nc.Dataset('/data/vdo/MEOPAR/completed-runs/expanded_again2/1_output.abort.nc')
expanded2_mask = nc.Dataset('/data/vdo/MEOPAR/completed-runs/expanded_again2/1_mesh_mask.nc')
fig,ax = plt.subplots(1,1,figsize=((10,8)))
ssh = ma.masked_array(expanded2.variables['sossheig'][0,:,:],mask = 1-expanded2_mask.variables['tmask'][0,0,:,:])
mesh = ax.pcolormesh(ssh,#vmin = 0, vmax = 1.4
)
fig.colorbar(mesh)
<matplotlib.colorbar.Colorbar at 0x7f219bc9a240>
fig,(ax1,ax2) = plt.subplots(1,2,figsize=((14,12)))
s = ma.masked_array(expanded2.variables['votemper'][0,27,:,:], mask = 1 - expanded2_mask.variables['tmask'][0,27,:,:])
mesh = ax1.pcolormesh(s)
fig.colorbar(mesh,ax=ax1)
ax1.plot(121,532,'ro')
#ax1.title('')
s = ma.masked_array(expanded2.variables['votemper'][0,:,532,:], mask = 1 - expanded2_mask.variables['tmask'][0,:,532,:])
mesh = ax2.pcolormesh(s)
ax2.set_ylim((32,0))
fig.colorbar(mesh,ax=ax2)
ax2.plot(121,27,'ro')
[<matplotlib.lines.Line2D at 0x7f21a0ca7390>]
fig,(ax1,ax2) = plt.subplots(1,2,figsize=((14,12)))
s = ma.masked_array(expanded2.variables['vosaline'][0,27,:,:], mask = 1 - expanded2_mask.variables['tmask'][0,27,:,:])
mesh = ax1.pcolormesh(s)
fig.colorbar(mesh,ax=ax1)
ax1.plot(121,532,'r*')
#ax1.title('')
s = ma.masked_array(expanded2.variables['vosaline'][0,:,532,:], mask = 1 - expanded2_mask.variables['tmask'][0,:,532,:])
mesh = ax2.pcolormesh(s)
ax2.set_ylim((32,0))
fig.colorbar(mesh,ax=ax2)
ax2.plot(121,27,'r*')
[<matplotlib.lines.Line2D at 0x7f21a0ab7828>]
fig,(ax1,ax2) = plt.subplots(1,2,figsize=((14,12)))
s = ma.masked_array(expanded2.variables['vozocrtx'][0,27,:,:], mask = 1 - expanded2_mask.variables['umask'][0,27,:,:])
mesh = ax1.pcolormesh(s,vmin = -1, vmax = 1
)
fig.colorbar(mesh,ax=ax1)
ax1.plot(121,532,'r*')
s = ma.masked_array(expanded2.variables['vozocrtx'][0,:,532,:], mask = 1 - expanded2_mask.variables['umask'][0,:,532,:])
mesh = ax2.pcolormesh(s, #vmin = -0.1, vmax = 0.07
)
ax2.set_ylim((32,0))
fig.colorbar(mesh,ax=ax2)
ax2.plot(121,27,'r*')
[<matplotlib.lines.Line2D at 0x7f21a0ce0d30>]
fig,(ax1,ax2) = plt.subplots(1,2,figsize=((14,12)))
s = ma.masked_array(expanded2.variables['vomecrty'][0,27,:,:], mask = 1 - expanded2_mask.variables['vmask'][0,27,:,:])
mesh = ax1.pcolormesh(s,#vmin = -1, vmax = 1.25
)
fig.colorbar(mesh,ax=ax1)
ax1.plot(121,532,'r*')
#ax1.title('')
s = ma.masked_array(expanded2.variables['vomecrty'][0,:,532,:], mask = 1 - expanded2_mask.variables['vmask'][0,:,532,:])
mesh = ax2.pcolormesh(s)
ax2.set_ylim((32,0))
fig.colorbar(mesh,ax=ax2)
ax2.plot(121,27,'r*')
[<matplotlib.lines.Line2D at 0x7f21a0284ba8>]
print(expanded2.variables['vomecrty'][0,27,532,121],
expanded2.variables['vozocrtx'][0,27,532,121],
expanded2.variables['vovecrtz'][0,27,532,121])
0.0 0.0 0.0
print(expanded.variables['vozocrtx'][0,27,532,121],
expanded.variables['vomecrty'][0,27,532,121],
expanded.variables['vovecrtz'][0,27,532,121])
0.0 0.0 0.0
print(agrif_abort.variables['vozocrtx'][0,27,532,121],
agrif_abort.variables['vomecrty'][0,27,532,121],
agrif_abort.variables['vovecrtz'][0,27,532,121])
0.0 0.0 0.0