pro find_singularity
  restore,"temp/level_1a_trial_0160_orbit_000000_jeppesen_cam_py.sav"
  for i=0,n_elements(data)-1 do begin
    ssa=upsample(*(data[i].scattering_angle),10)
    sva=upsample(*(data[i].viewing_angle),10)
    sza=upsample(*(data[i].zenith_angle),10)
    rp=rayleigh_phase(ssa)
    albedo=calc_albedo(39,0.608,rp,sva,sza)
    plot,sza,chap(49.11248,sza,8.37*0.608),charsize=3,title=string(i),psym=3,/ylog
    oplot,sza,abs(1.0/cos(sza*!const.dtor)),psym=3,color=254
    wait,3
  end
end