pro create_many_saturn_images,phase_angles_deg, $
		ring_facing_angles_deg,$
		incidences_rings_deg , $
		planet_size_pix,field_versus_planet_size,out_dir $
		,phase_curve,inc_emiss_azim_angles_deg,$
		rho_albedo_phase_index_ring_versus_r_sat,ring_phase_lookup_table $
		,label_lookup,tau_lookup,albedo_lookup,solar_elev_lookup,obs_elev_lookup,obs_azim_lookup,ss_phase_index_lookup

	for i_inc=0,n_elements(incidences_rings_deg)-1 do begin
	for i_ring=0,n_elements(ring_facing_angles_deg)-1 do begin
	for i_ph=0,n_elements(phase_angles_deg)-1 do begin
		;print,'ss'
		phase_angle_deg=phase_angles_deg(i_ph)
		ring_facing_angle_deg=ring_facing_angles_deg(i_ring)
		incidence_rings_deg=incidences_rings_deg(i_inc)
		print,' phase_angle_deg',phase_angle_deg,' ring_facing_angle_deg',ring_facing_angle_deg,'incidence_rings_deg',incidence_rings_deg;,'subsolar_azimuths_deg', subsolar_azimuths_deg,'range_incident_rings',range_incident_rings

		if(abs(ring_facing_angle_deg+incidence_rings_deg) lt phase_angle_deg $
			and abs(ring_facing_angle_deg-incidence_rings_deg) lt phase_angle_deg $
			and abs(180.+(ring_facing_angle_deg-incidence_rings_deg)) gt phase_angle_deg $
			and abs(180.-(ring_facing_angle_deg-incidence_rings_deg)) gt phase_angle_deg) then begin
			create_single_image_inc_emis_phase,phase_angle_deg, $
				ring_facing_angle_deg,$
				incidence_rings_deg , $
				planet_size_pix,field_versus_planet_size,out_dir, $
				phase_curve,inc_emiss_azim_angles_deg,$
				rho_albedo_phase_index_ring_versus_r_sat,ring_phase_lookup_table $
				,label_lookup,tau_lookup,albedo_lookup,solar_elev_lookup,obs_elev_lookup,obs_azim_lookup,ss_phase_index_lookup

		endif
	endfor
	endfor
	endfor
end