;!!!!!!!!!!!!!! NOTE:
; this is taken out of saturn_brightness.pro to save computations
; and not to iterate it over pixels
; Tables from Donese 93
; NOTE !!!! I put the phase angle 180 there myself, copied the values for 150
; !!!! THIS IS NOT PHYSICALLY REALISTIC, JUST A PLACEHOLDER
		ab_belt_blue=[[0.    ,30.   , 60   ,90    ,120   ,150   ,180] $
					 ,[0.6736,0.6197,0.5726,0.5948,0.7869,2.1323,2.1323] $
					 ,[1.1339,1.1254,1.1552,1.1970,1.2385,1.4539,1.4539]]
		ab_zone_blue=[[0.    ,30.   , 60   ,90    ,120   ,150   ,180] $
					 ,[0.5950,0.5659,0.5406,0.5290,0.6015,1.5886,1.5886] $
					 ,[1.0948,1.0938,1.1372,1.1606,1.1636,1.3727,1.3727]]
		ab_belt_red= [[0.    ,30.   , 60   ,90    ,120   ,150   ,180] $
					 ,[1.6854,1.5946,1.4529,1.3386,1.3653,2.2539,2.2539] $
					 ,[1.4812,1.4779,1.4586,1.4224,1.3637,1.3385,1.3385]]
		ab_zone_red= [[0.    ,30.   , 60   ,90    ,120   ,150   ,180] $
					 ,[1.6911,1.5995,1.4567,1.3436,1.3680,2.2021,2.2021] $
					 ,[1.4825,1.4794,1.4592,1.4221,1.3618,1.3341,1.3341]]

		ab_barkstorm_versus_phase_angle=$
			(ab_belt_blue+ab_zone_blue+ab_belt_red+ab_zone_red)/4.

		; CONVERT ANGLES TO RADIANS
		ab_barkstorm_versus_phase_angle(*,0)= $
			(180.-ab_barkstorm_versus_phase_angle(*,0))*!pi/180.
		;plot,ab_barkstorm_versus_phase_angle(*,0),ab_barkstorm_versus_phase_angle(*,1)
		;oplot,ab_barkstorm_versus_phase_angle(*,0),ab_barkstorm_versus_phase_angle(*,2),color=150