st-bender /
sciapy
| @@ 70-79 (lines=10) @@ | ||
| 67 | iers.conf.iers_auto_url = "https://astroconda.org/aux/astropy_mirror/iers_a_1/finals2000A.all" |
|
| 68 | ||
| 69 | ||
| 70 | def solar_zenith_angle(alt, lat, lon, time): |
|
| 71 | atime = Time(time) |
|
| 72 | loc = coord.EarthLocation.from_geodetic( |
|
| 73 | height=alt * units.km, |
|
| 74 | lat=lat * units.deg, |
|
| 75 | lon=lon * units.deg, |
|
| 76 | ) |
|
| 77 | altaz = coord.AltAz(location=loc, obstime=atime) |
|
| 78 | sun = coord.get_sun(atime) |
|
| 79 | return sun.transform_to(altaz).zen.value |
|
| 80 | ||
| 81 | ||
| 82 | def read_spectra(year, orbit, spec_base=None, skip_upleg=True): |
|
| @@ 69-78 (lines=10) @@ | ||
| 66 | iers.conf.iers_auto_url = "https://astroconda.org/aux/astropy_mirror/iers_a_1/finals2000A.all" |
|
| 67 | ||
| 68 | ||
| 69 | def solar_zenith_angle(alt, lat, lon, time): |
|
| 70 | atime = Time(time) |
|
| 71 | loc = coord.EarthLocation.from_geodetic( |
|
| 72 | height=alt * units.km, |
|
| 73 | lat=lat * units.deg, |
|
| 74 | lon=lon * units.deg, |
|
| 75 | ) |
|
| 76 | altaz = coord.AltAz(location=loc, obstime=atime) |
|
| 77 | sun = coord.get_sun(atime) |
|
| 78 | return sun.transform_to(altaz).zen.value |
|
| 79 | ||
| 80 | ||
| 81 | def read_spectra(year, orbit, spec_base=None, skip_upleg=True): |
|
