Beispiel #1
0
 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()
Beispiel #2
0
        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()
Beispiel #3
0
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()
Beispiel #4
0
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()
Beispiel #5
0
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)