#
#txres = Ngl.Resources() # Text resources desired
#txres.txFontHeightF = 0.015
#txres.txFontColor = "OrangeRed"
istart = 356
jstart = 373
iend = 387
jend = 411
for file in lst_file:
print "Plotting "+file
nc = netCDF4.Dataset(file, "r")
temp = nc.variables["temp"][0,:,:,:]
time = nc.variables["ocean_time"][0]
myday = jday2date(time/86400.)
# date_tag = myday.strftime('%d %B %Y')
date_tag = myday.strftime('%Y_%m_%d')
print date_tag
plotout = date_tag + '.png'
# plot = Ngl.contour_map(wks, aice, res)
# Ngl.text_ndc(wks, date_tag, 0.85, 0.84, txres)
plt.clf()
pyroms_toolbox.transectview(temp, -1, istart, iend, jstart, \
jend, 'ARCTIC2', outfile=plotout, cmin=-1.8, cmax=10.25)
nc.close()
# Ngl.frame(wks)
#Ngl.end()
m.drawmapboundary(fill_color='#99ffff')
m.fillcontinents(color='0.3', lake_color='0.3')
# m.fillcontinents(color='coral',lake_color='aqua')
idiv = nc.variables["sh_d_hf"][it, :, :]
aice = nc.variables["aice"][it, :, :]
idiv = np.where(aice == 0.0, 0.0, idiv)
idiv *= 1.e5
cs = m.pcolormesh(x, y, idiv, vmin=-.6, vmax=.6, cmap=cmap)
# cs = m.contourf(x, y, idiv, levels=levels, cmap=cmap, extend='both')
parallels = np.arange(45., 90., 15.)
# labels = [left,right,top,bottom]
m.drawparallels(parallels)
meridians = np.arange(0., 360., 15.)
m.drawmeridians(meridians, labels=[1, 0, 0, 1])
# csa = m.contour(x, y, idiv, levels=levels, linewidths=(1,), colors='k')
# csa = m.contour(x, y, idiv, levels=levels, colors=('k',))
myday = jday2date(times[it] / 24.)
date_tag = myday.strftime('%d %B %Y')
plt.title(date_tag, fontsize=20)
fname = myday.strftime('movie_idiv/frame_%(number)04d.png' %
{'number': it})
print(date_tag)
cbar = plt.colorbar(cs, orientation='vertical')
cbar.ax.tick_params(labelsize=15)
plt.savefig(fname)
plt.close()
nc.close()
if plot_grid: plt.plot(x, z.T, 'k')
# if plot_grid: plt.plot(x, z.T, 'k', hold=True);
# plt.ylim(-700,1)
plt.ylim(-60, 3)
cbar = plt.colorbar(cs, orientation='horizontal')
cbar.ax.tick_params(labelsize=15)
#plt.xlim(400,600)
times = z_file.variables["time"][:]
ntimes = len(times)
cmap = plt.cm.get_cmap("seismic")
for record in range(ntimes - 10, ntimes):
#record = -1 # Last record
myday = jday2date(times[record] / 24.)
date_tag = myday.strftime('%H:%M, %d %B %Y')
plt.figure(figsize=(10, 6))
plt.suptitle(date_tag, fontsize=20)
# plt.subplot(1,2,1); plot_section(z_file, record, variable='u', vmin=-2., vmax=2.); plt.title('U');
# plt.subplot(1,2,2); plot_section(z_file, record, y=427, variable='v', vmin=-2., vmax=2.); plt.title('V');
plt.subplot(2, 2, 1)
plot_section(z_file, record, variable='u', vmin=-2., vmax=2.)
plt.title('U')
plt.subplot(2, 2, 2)
plot_section(z_file, record, y=359, variable='v', vmin=-2., vmax=2.)
plt.title('V')
plt.subplot(2, 2, 3)
plot_section(z_file, record, variable='salt', vmin=None, vmax=None)
plt.title('salt')
plt.subplot(2, 2, 4)