def update(val):
ratio=so2.val
constants['ratio_copy']=ratio
o3,o2,o,J_o2,J_o3,o3_running=steady(constants,False)
m.set_ydata(J_o2*1E10)
l.set_ydata(o3/1E12)
n.set_ydata(J_o3*1E3)
du=o3_running/2.69E16
#ax0.annotate(str(int(round(du))), xy=(0.99,0.01),xycoords='axes fraction',horizontalalignment='right', verticalalignment='bottom')
fig.canvas.draw_idle()
def update(val):
ratio=so2.val
constants['M_surf']=2.5E19*(ratio+0.79)
constants['M']=constants['M_surf']*np.exp(-constants['heights']/constants['H'])
o3,o2,o,J_o2,J_o3,o3_running=steady(constants,False)
m.set_ydata(J_o2*1E10)
l.set_ydata(o3/1E12)
n.set_ydata(J_o3*1E3)
du=o3_running/2.69E16
#ax0.annotate(str(int(round(du))), xy=(0.99,0.01),xycoords='axes fraction',horizontalalignment='right', verticalalignment='bottom')
fig.canvas.draw_idle()
def logoxyoz(constants):
ratios=np.logspace(-6,0,50)
fig, ax=plt.subplots()
for i in range(len(ratios)):
constants['ratio']=ratios[i]
constants['M_surf']=2.5E19*(constants['ratio']+0.79)
constants['M']=constants['M_surf']*np.exp(-constants['heights']/constants['H'])
o3,o2,o,J_o2,J_o3,o3_running=steady(constants,False)
du=o3_running
# du=o3_running/2.69E16
ratio=constants['ratio']
plt.scatter(ratio,du,marker="x",color='r')
ax.set_xscale('log')
#ax.set_xlim([0,1])
ax.set_yscale('log')
plt.xlabel(r'Oxygen Content (PAL)')
plt.ylabel(r'Ozone Column / cm-2')
#plt.show()
plt.savefig('dobson_log.png')
plt.close()
def linoxyoz(constants):
ratios=np.linspace(0.01,1.00,100)
fig, ax=plt.subplots()
for i in range(len(ratios)):
constants['ratio']=ratios[i]
constants['M_surf']=2.5E19*(constants['ratio']+0.79)
constants['M']=constants['M_surf']*np.exp(-constants['heights']/constants['H'])
o3,o2,o,J_o2,J_o3,o3_running=steady(constants,False)
du=o3_running
du=o3_running/2.69E16
rat=100*constants['ratio']
plt.scatter(rat,du,marker=".",color='navy')
x,y=(21,21),(0,1100)
plt.plot(x,y,color='k')
ax.set_xlim([0,80])
ax.set_ylim([550,1100])
plt.xlabel(r'Oxygen Content / %')
plt.ylabel(r'Ozone Column / DU')
plt.savefig('dobson.png')
plt.close()
def steady_state_chapman(constants):
print "CALCULATING STEADY STATE OZONE COLUMN (CHAPMAN)"
o3,o2,o,J_o2,J_o3,o3_running=steady(constants,True)
altconc(constants,o3,o3_running)