def _dumpic(self): """ Dump initial condition plots For each variable (uc,vc,T,S,eta): - surface plot - seabed plot """ varnames = ['uc','vc','temp','salt','eta'] sun = Spatial(self.suntanspath+'/'+self.icfile,klayer=[0]) # Plot the surface variables for vv in varnames: h=plt.figure() sun.variable=vv sun.loadData() sun.clim = [sun.data.min(),sun.data.max()] sun.plot() outfile = '%s/IC_%s_surface.png'%(self.plotdir,vv) sun.savefig(outfile) del h # Plot the seabed variables sun.klayer=[-1] for vv in varnames: h=plt.figure() sun.variable=vv sun.loadData() sun.clim = [sun.data.min(),sun.data.max()] sun.plot() outfile = '%s/IC_%s_seabed.png'%(self.plotdir,vv) sun.savefig(outfile) del h
def _dumpic(self): """ Dump initial condition plots For each variable (uc,vc,T,S,eta): - surface plot - seabed plot """ varnames = ['uc', 'vc', 'temp', 'salt', 'eta'] sun = Spatial(self.suntanspath + '/' + self.icfile, klayer=[0]) # Plot the surface variables for vv in varnames: h = plt.figure() sun.variable = vv sun.loadData() sun.clim = [sun.data.min(), sun.data.max()] sun.plot() outfile = '%s/IC_%s_surface.png' % (self.plotdir, vv) sun.savefig(outfile) del h # Plot the seabed variables sun.klayer = [-1] for vv in varnames: h = plt.figure() sun.variable = vv sun.loadData() sun.clim = [sun.data.min(), sun.data.max()] sun.plot() outfile = '%s/IC_%s_seabed.png' % (self.plotdir, vv) sun.savefig(outfile) del h
Generates a series of animations """ from sunpy import Spatial import numpy as np import matplotlib.pyplot as plt runname='SFBay3D' ncfile = '%s/%s_0*.nc'%('data',runname) k = 1 # depth layer plt.figure() sun = Spatial(ncfile,klayer=[k],variable='salt') sun.tstep=np.arange(0,len(sun.time)) sun.loadData() sun.clim = [28.0,32.0] sun.animate(vector_overlay=False) sun.saveanim('plots/%s_salt.mov'%runname) sun.clim=None plt.figure() sun = Spatial(ncfile,klayer=[k],variable='uc') sun.tstep=np.arange(0,len(sun.time)) sun.loadData() sun.animate(vector_overlay=False) sun.saveanim('plots/%s_uc.mov'%runname) sun.clim=None plt.figure() sun = Spatial(ncfile,klayer=[k],variable='eta') sun.tstep=np.arange(0,len(sun.time))