tav=readx('av')

pscl,file='lnu.ps'
pset
jscolors
efacl=1000
efach=10000
efacf=100000
plot,tav(0,*),tav(2,*)/1000,psym=3,xtitle='!8l!X',ytitle='Frequency (mHz)',title='15 year MDI Medium-!8l!X mode coverage',yrange=[0.8,4.5],ystyle=1,xrange=[-5,305],xstyle=1,symsize=0.5,xmargin=[8,2]
wl=where(tav(0,*) lt 75)
wh=where(tav(0,*) ge 75)
wl=where(tav(1,*) ge 12)
wh=where((tav(1,*) lt 12) and (tav(1,*) ne 0))
wf=where(tav(1,*) eq 0)
errplot,tav(0,wl),(tav(2,wl)-efacl*tav(7,wl))/1000,(tav(2,wl)+efacl*tav(7,wl))/1000,color=1
errplot,tav(0,wh),(tav(2,wh)-efach*tav(7,wh))/1000,(tav(2,wh)+efach*tav(7,wh))/1000,color=2
errplot,tav(0,wf),(tav(2,wf)-efacf*tav(7,wf))/1000,(tav(2,wf)+efacf*tav(7,wf))/1000,color=3

errplot,[200],[4.15],[4.25],color=1
xyouts,205,4.15,strcompress(string(efacl),/remove_all)+'!ms'
errplot,[200],[3.95],[4.05],color=2
xyouts,205,3.95,strcompress(string(efach),/remove_all)+'!ms'
errplot,[200],[3.75],[3.85],color=3
xyouts,205,3.75,strcompress(string(efacf),/remove_all)+'!ms'

device,color=0
loadct,0
psend1
xwin

tav75=readx('av75')

pscl,file='lnu75.ps'
pset
jscolors
efacl=1000
efach=10000
efacf=100000
plot,tav75(0,*),tav75(2,*)/1000,psym=3,xtitle='!8l!X',ytitle='Frequency (mHz)',title='15 year MDI Medium-!8l!X mode coverage',yrange=[0.8,4.5],ystyle=1,xrange=[-5,305],xstyle=1,symsize=0.5,xmargin=[8,2]
wl=where(tav75(0,*) lt 75)
wh=where(tav75(0,*) ge 75)
wl=where(tav75(1,*) ge 12)
wh=where((tav75(1,*) lt 12) and (tav75(1,*) ne 0))
wf=where(tav75(1,*) eq 0)
errplot,tav75(0,wl),(tav75(2,wl)-efacl*tav75(7,wl))/1000,(tav75(2,wl)+efacl*tav75(7,wl))/1000,color=1
errplot,tav75(0,wh),(tav75(2,wh)-efach*tav75(7,wh))/1000,(tav75(2,wh)+efach*tav75(7,wh))/1000,color=2
errplot,tav75(0,wf),(tav75(2,wf)-efacf*tav75(7,wf))/1000,(tav75(2,wf)+efacf*tav75(7,wf))/1000,color=3

errplot,[200],[4.15],[4.25],color=1
xyouts,205,4.15,strcompress(string(efacl),/remove_all)+'!ms'
errplot,[200],[3.95],[4.05],color=2
xyouts,205,3.95,strcompress(string(efach),/remove_all)+'!ms'
errplot,[200],[3.75],[3.85],color=3
xyouts,205,3.75,strcompress(string(efacf),/remove_all)+'!ms'

device,color=0
loadct,0
psend1
xwin

end