xax=300.+findgen(21)*20
v=xax^(-4)*10.^10

set_plot,'ps'
device,filename='rayleighxsecgraph.ps'
device,/landscape
!p.font=0
device,/palatino
!p.charsize=2.
plot,xax,v/v(20),xtitle='wavelength, nm', ytitle='relative Rayleigh cross section', $
thick=3
oplot,[400,400],[0,40],linestyle=2
xyouts,500,25,charsize=2.,'Visible'  &  xyouts,340,25,charsize=2,'UV'
device,/close

set_plot,'x'
SPAWN,'/app/scivis/netpbm/bin/pstopnm -land -xsize 400 rayleighxsecgraph.ps'
SPAWN,'/app/scivis/netpbm/bin/ppmtogif rayleighxsecgraph001.ppm > '+ $
    'rayleighxsecgraph.gif'


END