extract sun zenith and azimuth angle

[rese_support]

I was wondering if we could extract sun zenith and azimuth angle for all images at once when we could provide image acquisition year and time as well as latitude and longitude for all image at the same time.

You'd need to write a small IDL script to do so; the function 'sun' from within atcor. Usage is as follows:


 angles = sun(month=month,day=day,hour=hour,minute=minute, lat=lat,lon=lon)
 azi = angles[1]
 zen = angles[0]
 dayofyear=angles[2]