ind_title = 'Indirect effect (W/m2): (sulphurless + biogenic) - (sulphurless) '
dir_title = 'Direct effect (W/m2) (sulphurless - present day) - NEW'
ind_title = 'Indirect effect (W/m2)(sulphurless - present day) - NEW'
dir_title = 'Direct effect (W/m2) (PD+BIOGEN - PD) - NEW'
ind_title = 'Indirect effect (W/m2)(PD+BIOGEN - PD) - NEW'
di_eff.units = 'W m-2'
in_eff.units = 'W m-2'
print('max value dir effect = ', np.max(di_eff.data))
print('min value dir effect = ', np.min(di_eff.data))
print('mean value dir effect = ', np.mean(di_eff.data))
#pltfunc.plot_diff(di_eff,dir_title,-2,2,'coolwarm')
pltfunc.plot_diff(di_eff,dir_title,-2,2,'coolwarm')
print('max value indir effect = ', np.max(in_eff.data))
print('min value indir effect = ', np.min(in_eff.data))
print('mean value indir effect = ', np.mean(in_eff.data))
plt.savefig('dir_effect_nosulphur_2.eps',dp1 = 500)
#pltfunc.plot_diff(in_eff,ind_title,-10,10,'coolwarm')
pltfunc.plot_diff(in_eff,ind_title,-2,2,'coolwarm')
plt.savefig('indir_effect_nosulphur_2.eps',dp1 = 500)
print('mean direct effect ',np.mean(di_eff.data))
print('mean indirect effect ',np.mean(in_eff.data))
#print 'Surface albedo forcing -- SHORTWAVE , ',np.mean(surf_alb_sw)
#print 'Surface albedo forcing -- LONGWAVE , ',np.mean(surf_alb_lw)
plt.show()
if i % 3 == 0 and i < 49:
dat1 = get_model_data(p_day1, data_level)
dat2 = get_model_data(p_day2, data_level)
a = dat1
b = dat2
#a = dat1*(29.0/150.0)/1e-12 # convert it to ppt
#b = dat2*(29.0/150.0)/1e-12
c = (b - a) * 100.0 / a
#c = b
#c = a
#c = b-a
#d=dat3
#e=(d-a)*100/a
print('the highest is = ', np.max(c.data))
print('the lowest is = ', np.min(c.data))
print('mean = ', np.mean(c.data))
temp_str = 'Percentage increase in Ozone concentration (1km) - present day'
#c.data=np.where(c.data>=0,c.data,1e-20) #IMPORTANT LINE TO MASK UNNCESSARY DATA
#plot_diff(c,temp_str,1e-20,1e0,20)
pltfunc.plot_diff(c, temp_str, -30, 30, 'seismic_r')
plt.savefig('ozone_nosulphur_2.eps', dpi=500)
plt.show()
# ------------------------------------------------------------------
di_eff_2 = di_eff_2.collapsed('time', iris.analysis.MEAN)
#di_eff.units = 'W m-2'
di_eff = di_eff * -1
in_eff = in_eff * -1
dir_title = 'Direct effect (W/m2): PI_sulphurless - PI'
ind_title = 'Indirect effect (W/m2): PI_sulphurless - PI'
di_eff.units = 'W m-2'
in_eff.units = 'W m-2'
print('max value dir effect = ', np.max(di_eff.data))
print('min value dir effect = ', np.min(di_eff.data))
print('mean value dir effect = ', np.mean(di_eff.data))
#plot_diff(di_eff,dir_title,-120,121,10)
pltfunc.plot_diff(di_eff, dir_title, -14, 14, 'coolwarm')
#pltfunc.plot_diff(di_eff,dir_title,np.min(di_eff.data),np.max(di_eff.data),'seismic')
print('max value indir effect = ', np.max(in_eff.data))
print('min value indir effect = ', np.min(in_eff.data))
print('mean value indir effect = ', np.mean(in_eff.data))
plt.savefig('dir_effect_nosulphur_2.eps', dp1=500)
#plot_diff(in_eff,ind_title,-550,551,100)
pltfunc.plot_diff(in_eff, ind_title, -45, 45, 'coolwarm')
#pltfunc.plot_diff(in_eff,ind_title,np.min(in_eff.data),np.max(in_eff.data),'seismic')
plt.savefig('indir_effect_nosulphur_2.eps', dp1=500)
print('Direct radiative forcing -- SHORTWAVE , ', np.mean(di_eff.data))
print('Direct radiative forcing -- LONGWAVE , ', np.mean(in_eff.data))
#print 'Surface albedo forcing -- SHORTWAVE , ',np.mean(surf_alb_sw)
#print 'Surface albedo forcing -- LONGWAVE , ',np.mean(surf_alb_lw)
#di_eff.units = 'W m-2'
di_eff = di_eff*-1
in_eff = in_eff*-1
dir_title = 'Direct effect (W/m2) '
ind_title = 'Indirect effect (W/m2)'
di_eff.units = 'W m-2'
in_eff.units = 'W m-2'
print('max value dir effect = ', np.max(di_eff.data))
print('min value dir effect = ', np.min(di_eff.data))
print('mean value dir effect = ', np.mean(di_eff.data))
#plot_diff(di_eff,dir_title,-120,121,10)
pltfunc.plot_diff(di_eff,dir_title,-2,2,'seismic')
#pltfunc.plot_diff(di_eff,dir_title,np.min(di_eff.data),np.max(di_eff.data),'seismic')
print('max value indir effect = ', np.max(in_eff.data))
print('min value indir effect = ', np.min(in_eff.data))
print('mean value indir effect = ', np.mean(in_eff.data))
plt.savefig('dir_effect_nosulphur_2.eps',dp1 = 500)
#plot_diff(in_eff,ind_title,-550,551,100)
#pltfunc.plot_diff(in_eff,ind_title,np.min(in_eff.data),np.max(in_eff.data),'seismic')
pltfunc.plot_diff_3(in_eff,ind_title,-2,2,'seismic')
plt.savefig('indir_effect_nosulphur_2.eps',dp1 = 500)
print('Direct radiative forcing = ',np.mean(di_eff.data))
print('Indirect radiative forcing = ',np.mean(in_eff.data))
#print 'Surface albedo forcing -- SHORTWAVE , ',np.mean(surf_alb_sw)
#print 'Surface albedo forcing -- LONGWAVE , ',np.mean(surf_alb_lw)
di_eff_2 = di_eff_2.collapsed('time', iris.analysis.MEAN)
#di_eff.units = 'W m-2'
di_eff = di_eff * -1
in_eff = in_eff * -1
dir_title = 'Direct effect (W/m2) - present day'
ind_title = 'Indirect effect (W/m2) - present day'
di_eff.units = 'W m-2'
in_eff.units = 'W m-2'
print('max value dir effect = ', np.max(di_eff.data))
print('min value dir effect = ', np.min(di_eff.data))
print('mean value dir effect = ', np.mean(di_eff.data))
#plot_diff(di_eff,dir_title,-120,121,10)
pltfunc.plot_diff(di_eff, dir_title, -14, 14, 'seismic')
#pltfunc.plot_diff(di_eff,dir_title,np.min(di_eff.data),np.max(di_eff.data),'seismic')
print('max value indir effect = ', np.max(in_eff.data))
print('min value indir effect = ', np.min(in_eff.data))
print('mean value indir effect = ', np.mean(in_eff.data))
plt.savefig('dir_effect_nosulphur_2.eps', dp1=500)
#plot_diff(in_eff,ind_title,-550,551,100)
#pltfunc.plot_diff(in_eff,ind_title,np.min(in_eff.data),np.max(in_eff.data),'seismic')
pltfunc.plot_diff(in_eff, ind_title, -45, 45, 'seismic')
plt.savefig('indir_effect_nosulphur_2.eps', dp1=500)
print('Direct radiative forcing -- SHORTWAVE , ', np.mean(di_eff.data))
print('Direct radiative forcing -- LONGWAVE , ', np.mean(in_eff.data))
#print 'Surface albedo forcing -- SHORTWAVE , ',np.mean(surf_alb_sw)
#print 'Surface albedo forcing -- LONGWAVE , ',np.mean(surf_alb_lw)
di_eff_2 = di_eff_2.collapsed('time', iris.analysis.MEAN)
#di_eff.units = 'W m-2'
di_eff = di_eff * -1
in_eff = in_eff * -1
dir_title = 'Direct effect (W/m2)'
ind_title = 'Indirect effect (W/m2)'
di_eff.units = 'W m-2'
in_eff.units = 'W m-2'
print('max value dir effect = ', np.max(di_eff.data))
print('min value dir effect = ', np.min(di_eff.data))
print('mean value dir effect = ', np.mean(di_eff.data))
#plot_diff(di_eff,dir_title,-120,121,10)
pltfunc.plot_diff(di_eff, dir_title, np.min(di_eff.data),
np.max(di_eff.data), 'seismic')
print('max value indir effect = ', np.max(in_eff.data))
print('min value indir effect = ', np.min(in_eff.data))
print('mean value indir effect = ', np.mean(in_eff.data))
plt.savefig('dir_effect_nosulphur_2.eps', dp1=500)
#plot_diff(in_eff,ind_title,-550,551,100)
pltfunc.plot_diff(in_eff, ind_title, np.min(in_eff.data),
np.max(in_eff.data), 'seismic')
plt.savefig('indir_effect_nosulphur_2.eps', dp1=500)
print('Direct radiative forcing -- SHORTWAVE , ', np.mean(di_eff.data))
print('Direct radiative forcing -- LONGWAVE , ', np.mean(in_eff.data))
#print 'Surface albedo forcing -- SHORTWAVE , ',np.mean(surf_alb_sw)
#print 'Surface albedo forcing -- LONGWAVE , ',np.mean(surf_alb_lw)
plt.show()