plot
#
	expand 0.75
	data 1614_na.plot
	read { i 1 j 2 k 7 m 8 p 9 q 10 }
	window 1 2 1 2
	ctype 2 ltype 0
#
        set x=5889.9512
        set y=(i/x)*0.984349
        set shift=y-1
        set v1=shift*300000
	set zk = -(LN(k))
	set zp = -(LN(p))
	set z1 = zk+zp
	erase
	limits -700 400 z1
	box 1 5 4 4 
	ctype 3 ltype 2
	connect v1 z1
#
	ctype 5 ltype 0
        set l=5895.9243
        set u=(i/l)*0.984349
        set shift=u-1
        set v2=shift*300000
	set zm = -(LN(m))
	set zq = -(LN(q))
	set z2 = zm+zq
	connect v2 z2
#
	relocate -650 1.45
	ctype 2
	label NGC1614 Na I
	relocate -650 1.35
	label Apparent Optical Depth
	relocate -650 0.3
	label \tau_a(v) = -ln \frac{I_{obs}(v)}{ I_o(v)}
	xlabel v (km/s)
	ylabel \tau_a
#
	data 1614legend.dat
	lines 1 5
	read { r 1 s 2 }
	ctype 2 ltype 0
	connect r s
	lines 6 7
	read { r 1 s 2 }
	ltype 0 ctype 5
	relocate 270 1.22
	label \lambda 5896
	connect r s
	lines 8 9 
	read { r 1 s 2 }
	ltype 2 ctype 3
	relocate 270 1.38
	label \lambda 5890
	connect r s
	ctype 2 ltype 0
#
	set nak = (3.768E14*z1)/(0.6311*i)
	set nam = (3.768E14*z2)/(0.3180*i)
	window 1 2 1 1
	limits -700 400 nam 
	box 1 5 4 4 
	ctype 5
	connect v2 nam
	ctype 3 ltype 2
	connect v1 nak
#
	ctype 2 ltype 0
	xlabel v (km/s)
	ylabel N_a (cm^{-2} [km/s]^{-1})
	relocate -650 1.4E11
	label NGC1614 Na I
	relocate -650 1.3E11
	label Apparent Column Density
	relocate -650 4E10
	label N_a(v)=\frac{m_ec}{\pi e^2}\frac{\tau_a(v)}{{\it f}\lambda(\AA)}