pro plot_vline, x, linestyle=linestyle
;
;
;
xarr = [1,1]*x
yarr = [!cymin, !cymax]
cmd = 'oplot, xarr, yarr'
if (keyword_set(linestyle)) then cmd = cmd + ', linestyle=linestyle'
stat = execute(cmd)
;
end