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): (no_sulphur+biogenic) -(no_sulphur)'
ind_title = 'Indirect effect (W/m2):(no_sulphur+biogenic) -(no_sulphur)'
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_3(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 -- 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) - month = ' + str(j)
ind_title = 'Indirect effect (W/m2) - month = ' + str(j)
#dir_title = 'Direct effect (W/m2) - (no_sulphur+biogenic)-(no_sulphur)'
#ind_title = 'Indirect effect (W/m2) - (no_sulphur+biogenic)-(no_sulphur)'
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_3(di_eff, dir_title, -1, 1, '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))
save_title_1 = 'direct_month_' + str(j) + '.eps'
save_title_2 = 'indirect_month_' + str(j) + '.eps'
plt.savefig(save_title_1, dp1=500)
#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, -4, 4, 'seismic')
plt.savefig(save_title_2, dp1=500)
#plt.savefig('indir_effect_nosulphur_2.eps',dp1 = 500)
print('Direct radiative forcing -- SHORTWAVE , ',
di_eff_2 = di_eff_2[6,:,:,:]
#di_eff.units = 'W m-2'
di_eff = di_eff*-1
in_eff = in_eff*-1
dir_title = 'Direct forcing (W/m2): (PD) -(PI)'
ind_title = 'Indirect forcing (W/m2):(PD)-(PI)'
di_eff.units = 'W m-2'
in_eff.units = 'W m-2'
image_path = '/home/users/eeara/no_sulphur_new/images/'
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_3(di_eff,dir_title,-7,7,'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(image_path+'dir_effect_nosulphur.png',dpi = 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,-7,7,'coolwarm')
plt.savefig(image_path + 'indir_effect_nosulphur_pdpi.png',dpi = 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): (no_sulphur+biogenic) -(no_sulphur)'
ind_title = 'Indirect effect (W/m2):(no_sulphur+biogenic) -(no_sulphur)'
di_eff.units = 'W m-2'
in_eff.units = 'W m-2'
image_path = '/home/users/eeara/no_sulphur_new/images/'
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_3(di_eff, dir_title, -2, 2, '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(image_path + 'dir_effect_nosulphur.png', dpi=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, -5, 5, 'coolwarm')
plt.savefig(image_path + 'indir_effect_nosulphur_2.png', dpi=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)