def run(self):
self.status = 'in progress'
if (self.each != 0):
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt
fig = plt.figure()
for i in range(0,self.maxiter):
self.hydrostep()
self.time += self.dt
if (self.each != 0):
if (i % self.each == 0):
I1 = (self.rho/self.m)**2*AIA.t171(self.T)*1e8
I2 = (self.rho/self.m)**2*AIA.t193(self.T)*1e8
I3 = (self.rho/self.m)**2*AIA.t211(self.T)*1e8
self.I171 += [np.sum((I1*self.dx)[np.where(np.abs(self.s-self.L*0.5) < self.L*0.5-self.L*0.15)])/self.L]
self.I193 += [np.sum((I2*self.dx)[np.where(np.abs(self.s-self.L*0.5) < self.L*0.5-self.L*0.15)])/self.L]
self.I211 += [np.sum((I3*self.dx)[np.where(np.abs(self.s-self.L*0.5) < self.L*0.5-self.L*0.15)])/self.L]
plt.subplot(211)
plt.title('Density 10^4 / Temperature', size = 18)
plt.plot(self.s, np.log10(self.rho/self.m)-4)
plt.plot(self.s, np.log10(self.T))
plt.axis([0,np.max(self.s),4,7])
#plt.axis([0,1e9,4,7])
plt.subplot(212)
#plt.plot(self.s, self.rhou)
#plt.axis([0,1e9,-1e-8,1e-8])
plt.title('AIA 171/193/211 intensity', size = 18)
plt.plot(self.s,np.log10(I1))
plt.plot(self.s,np.log10(I2))
plt.plot(self.s,np.log10(I3))
plt.axis([0,np.max(self.s),-2,4])
fname = (img_dir + str(np.int(i/self.each)).zfill(np.ceil(np.log10(self.maxiter/self.each)).astype(np.int)) + '.png')
fig.savefig(fname)
fig.clf()
self.status = 'done'
