IDL script to find the membrane mirror frequency response Data taken on August 3, 2003 at UH88. Reference data taken July 7, 2003 Freq Distance to focus (cm) 1740 47.70 1760 44.13 1780 42.55 1800 41.95 1820 41.79 1840 41.79 1860 42.86 1880 44.57 1900 45.84 1920 49.97 freq = [1740, 1760, 1780, 1800, 1820, 1840, 1860, 1880, 1900, 1920] dist = [47.70, 44.13, 42.55, 41.95, 41.79, 41.79, 42.86, 44.57, 45.84, 49.97] dist2 = dist + 3.175 freq2 = [1750, 1770, 1790, 1810, 1830, 1850, 1870, 1890, 1910, 1930] dist3 = [45.12, 42.70, 40.52, 40.16, 40.40, 41.12, 42.15, 43.10, 45.79, 49.33] dist4 = dist3+3.175 ; freq2 is a second set of data taken at the summit sum = dist2 + dist4 avg= sum/2.0 ; this averages the two data sets together freqhilo = [1800, 1810, 1820, 1830, 1840, 1850, 1860, 1870, 1880, 1890, 1900, 1910, 1920, 1930, 1940, 1950, 1960, 1970, 1980, 1990, 2000, 2020, 2040, 2060, 2080] disthilo = [36.23, 35.64, 34.73, 33.77, 33.14, 31.99, 30.96, 30.12, 29.45, 28.58, 27.74, 27.42, 26.63, 26.39, 26.11, 25.96, 26.00, 26.04, 26.31, 26.11, 26.55, 27.98, 29.85, 32.62, 35.60] set_plot, 'ps' Device, filename='membrane_summit_hilo.ps', Xoffset=0, Yoffset=2, Xsize=7.0, Ysize=7.0, /Inches sum = dist2 + dist4 avg= sum/2.0 plot, freq, avg, Xrange=[1720, 2100], Xstyle=1, Yrange=[20, 54], Ystyle=1, Xtitle = 'Frequency (Hz)',Ytitle = 'Focal length (cm)',Title='Focal length of IRTF membrane vs. Frequency at summit', Psym=-5 oplot, freqhilo, disthilo, Psym=-6 xvalues=[2000, 2030] yvalues=[50, 50] PlotS, xvalues, yvalues, Linestyle=0, Psym=-5 XYoutS, 0.85, 0.855, 'Summit', Alignment=0.0, /Normal, Charsize=1.0 xvalues=[2000, 2030] yvalues=[48, 48] PlotS, xvalues, yvalues, Linestyle=0, Psym=-6 XYoutS, 0.85, 0.80, 'Hilo', Alignment=0.0, /Normal, Charsize=1.0 device, /close ; closes device set_plot, 'x' ; makes it so that the plot goes to the screen again plot, freq, avg, Xrange=[1720, 1950], Xstyle=1, Yrange=[40, 54], Psym=-5 oplot, freqhilo-140, disthilo*1.68, Psym=-6