read into the .pro file when scaling day to night with "solar_contin_lower" and "solar_continu_upper" defining franhofer line so the synthetic solar spectrum can be shifted ==========5577============= openr, lun, 'solar_spectrum_6563.dat', /GET_LUN synth_solar=fltarr(2,72679) solar_contin_lower = 5575.1 solar_contin_upper = 5576.1 franhofer_lower = (5576.095 - 0.2) + (radial_vel[k] * 5575.9 / 3.0e5) franhofer_upper = (5576.095 + 0.2) + (radial_vel[k] * 5576.3 / 3.0e5) ==========6563============= openr, lun, 'solar_spectrum_6563.dat', /GET_LUN synth_solar=fltarr(2,2760) solar_contin_lower = 6565 solar_contin_upper = 6567 franhofer_lower = (6568.88 - 0.2) + (radial_vel[k] * 6568.4 / 3.0e5) franhofer_upper = (6568.88 + 0.2) + (radial_vel[k] * 6569.3 / 3.0e5) ------solar spectrum plot------- oplot, [6553, 6575], [0, 0], linestyle = 0, color =0 ----calculating scatter for each plot----- oplot, [6563, 6575], [sigma_continuum_upper, sigma_continuum_upper], THICK = 4, LINESTYLE = 1, COLOR = 10 oplot, [6563, 6575], [sigma_continuum_lower, sigma_continuum_lower], THICK = 4, LINESTYLE = 1, COLOR = 10 oplot, [6563, 6575], [mean_continuum, mean_continuum], THICK = 4, LINESTYLE = 0, COLOR = 10 ----plotting in normalized---- plot, lambda, Venus_night_final[*,i], xrange=[6558,6569], /xstyle, $ ==========3914============= openr, lun, 'solar_spectrum_3914.dat', /GET_LUN fluxcal=fltarr(2,2720) IF feature_true_lambda[k] EQ 3914.4 THEN BEGIN solar_contin_lower = 3917.2 solar_contin_upper = 3918 franhofer_lower = (3913.3 - 0.2) + (radial_vel[k] * 3913.8 / 3.0e5) franhofer_upper = (3913.3 + 0.2) + (radial_vel[k] * 3912.8 / 3.0e5)