...
| Code Block | ||||
|---|---|---|---|---|
| ||||
class SkyModel(object) def __init__(self, params, kwargs): Set#Set observatory parameters, filter throughput curves, which sky components to compute use#use the kwargs to set any "magic numbers" that are needed to tune the results to the site. load#load up template spectra def _loadMoonSED(self): load#load up a grid of scattered moon light SEDs def interpolateMoon(self, moonPhase, moonAlt, moon2targetDist): Interpolate#Interpolate the moon SEDs to the requested paramter values. #similar load/interpolate methods for the other sky components--zodiacal light, lower atmosphere, upper atmosphere, airglow continuum, twilight, def setRaDecMjd(self, ra, dec, mjd, kwargs): set the MJD of the observations and the locations to calculate sky brightnesses at. additional kwargs for things like PWV, solar activity, etc. Calculate things like airmass for each sky point, distance to moon for each sky point def calcSky(self, mags=True): use the template spectra to interpolate a best guess for the sky SED at each RA,Dec point. Return SEDs or magnitudes |
...
- We will need to add a twilight sky model from scratch since there is not one in the ESO model
- The Canon data does not extend into the IR, so we will need to rely on the photodiode data (pointed at zenith), to quantify the sky brightness variation in y-band.
- Does anyone have a favourite way to interpolate template spectra? I'd obviously like to avoid python loops.