Check Wavelength Solutions Against Sky OH Lines
We check the wavelength solutions derived from the arc lines against the sky OH lines.
The script is in /home/osrseng/syeh/wavelengthCalibration/skycubeCheck, titled check_[filterscale]sky[year].pro. E.g., check_Kc3100_sky2025.pro
Edit the lines in the script accordingly to specify which sky cube and OH line library to use:
pro check_Kc3100_sky2025, makecorr = makecorr, psfile = psfile
dir = '/home/osrseng/syeh/wavelengthCalibration/skycubeCheck/skycubes/'
skycubes = dir+['s250522_a011002_Kc3_100.fits']
print, skycubes
outdir = '/home/osrseng/syeh/wavelengthCalibration/skycubeCheck/'
ohfile = '/home/osrseng/syeh/wavelengthCalibration/skycubeCheck/OH-lines/ohlines_kn3.dat'Then edit the lines below:
pro kn3_ohlines_fit, psfile = psfile
; PURPOSE: load in a cube and do a fit to the OH lines with calibspec
file = '/home/osrseng/syeh/wavelengthCalibration/skycubeCheck/skycubes/s250522_a011002_Kc3_100.fits'
ohfile = '/home/osrseng/syeh/wavelengthCalibration/skycubeCheck/OH-lines/ohlines_kn3.dat'Compile and run the code:
IDL> .comp check_Kc3100_sky2025
% Compiled module: CHECK_KC3100_SKY2025.
% Compiled module: PLOT_FWHM_MAP.
% Compiled module: KN3_OHLINES_FIT.
IDL> check_Kc3100_sky2025, /makecorr, psfile='./Kc3100_sky2025.ps'The code will produce two plots. The left panel shows the wavelength solution shift in pixels in spatial directions. The right panel shows the wavelength solution shift in angstroms in spatial directions.
