; batch file called from 	create_many_saturn_images.pro, calculate_incidence_angles.pro,
; calculate_ring_brightness_plus_planet.pro etc..

	;x-east in the imageplane,y-north,z-toward the spacecraft
	;help,subsolar_azimuth_deg,phase_angle_deg

	solar_direction=[cos((subsolar_azimuth_deg+90.)/180.*!pi)*sin(double(phase_angle_deg)/180.*!pi),$
					sin((subsolar_azimuth_deg+90.)/180.*!pi)*sin(double(phase_angle_deg)/180.*!pi),$
					cos(double(phase_angle_deg)/180.*!pi)]
	solar_direction_abs=(solar_direction##transpose(solar_direction))^0.5
	solar_direction=solar_direction/solar_direction_abs(0)