#reset gamma.ax = None gamma_th.ax = None #gamma.t = 0 #print gamma[0,:] #/elem['dt']] fig = plt.figure() gamma.ax = fig.add_subplot(223) gamma_th = Frame(np.array(soln['gammamax'][1:ny/3]),meta={'dx':dky,'x0':dky,'stationary':True,'nt':gamma.nt,'yscale':'symlog'}) gamma_th.ax =fig.add_subplot(223) gamma.shareax = True gamma.yscale = 'symlog' sigma = blob.raw_data.std(axis=2) frm_data1D = Frame(np.average(blob.raw_data,axis=2),meta={'sigma':sigma,'t_array':time}) frames= [frm_data,frm_fft,gamma_th,gamma,frm_amp] #FrameMovie(frames,fast=False,fig=fig) fig = plt.figure() # #let's examine the boundary # bc_left = Frame(blob.raw_data[:,0:20,:],meta={'mask':True}) # bc_right = Frame(blob.raw_data[:,-5:-1,:],meta={'mask':True})
lam_history = Frame(lam,meta={'dx':tchunk,'stationary':False,'fontsz':18,'ylabel':'', 'xlabel':r'$t$','ticksize':14,'title':r'$\lambda$','xlabel':r't'}) lam_history.ax = fig.add_subplot(111) lam_history.ax.set_ylim([0.0,100.*np.round((lam_history[-1]+100.)/100.)]) lam_history.render(fig,111) print lam_history.x fig.savefig(pp,format='pdf') plt.close(fig) pp.close() pp = PdfPages(save_path+'/'+key+'compare_lam.pdf') fig = plt.figure() lam_history.ax = None lam_history.yscale='linear' #lam_history = Frame(lam,meta={'stationary':False,'title':'','fontsz':18,'ylabel':'','xlabel':r'$t$', # 'ticksize':14,'yscale':'linear'}) lam_rough_history = Frame(lam_rough,meta={'stationary':False,'title':'','fontsz':18,'ylabel':'','xlabel':r'$t$', 'ticksize':14,'yscale':'linear','dx':tchunk}) lam_history.render(fig,111) lam_rough_history.render(fig,111) fig.savefig(pp,format='pdf') plt.close(fig) pp.close() pp = PdfPages(save_path+'/'+key+'loglam.pdf') fig = plt.figure()