IDL script for plotting Membrane focal length, units, and volts

Units     F.L.(")   Volts (V)
6000      36.15      0.40
7000      30.15      0.45
8000      27.14      0.55
9000      24.12      0.60
10000     21.56      0.65
12000     17.65      0.80
15000     14.07      1.00
20000     10.57      1.35
30000      7.10      2.00
40000      5.38      2.70
50000      2.80      5.50


units = [6000, 7000, 8000, 9000, 10000, 12000, 15000, 20000, 30000, 40000, 50000]

length = [36.15, 30.15, 27.14, 24.12, 21.56, 17.65, 14.07, 10.57, 7.10, 5.38, 2.80]
length = length * 2.54

volts = [0.40, 0.45, 0.55, 0.60, 0.65, 0.80, 1.00, 1.35, 2.00, 2.70, 5.50] 



set_plot, 'ps'
Device, filename='units_volts.ps', Xoffset=0, Yoffset=2, Xsize=7.0, Ysize=7.0, /Inches

plot, units, volts, Xrange=[0, 55000], Xstyle=1, Yrange=[0., 6.], Psym=-5, Xtitle = 'Control Units',Ytitle = 'Voltage (V)',Title='Applied Voltage to Membrane vs. Control Units'

device, /close  ; closes device
set_plot, 'x'   ; makes it so that the plot goes to the screen again





set_plot, 'ps'
Device, filename='units_stroke.ps', Xoffset=0, Yoffset=2, Xsize=7.0, Ysize=7.0, /Inches

logunits = alog10(units)
loglength = alog10(length)
plot, logunits, loglength, Xrange=[3.6, 4.8], Xstyle=1,Ystyle=1, Yrange=[0.5, 2.3], Psym=-5, Xtitle = 'Log Control Units',Ytitle = 'Log Focal length (cm)',Title='Focal length of IRTF membrane vs. Control Units'

device, /close  ; closes device
set_plot, 'x'   ; makes it so that the plot goes to the screen again