#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()
