pro pdsplot1,source
  data = rdcol(source+'_current.pds',0L,20000L,[1,2])
  plot, data(0,*),data(1,*), $
    title=source,xtitle='Frequency (1/hr)', ytitle='Power/<Power>', $
    xrange=[0.0, 0.5]
  image = tvrd()
  write_gif,source+'_current_1.gif',image
  return
end

pro pdsplot2,source
  data = rdcol(source+'_current.pds',0L,20000L,[1,2])
  plot, data(0,*),data(1,*), $
    title=source,xtitle='Frequency (1/hr)', ytitle='Power/<Power>', $
    xrange=[0.0, 0.05]
  image = tvrd()
  write_gif,source+'_current_2.gif',image
  return
end

pro pdsweb
  on_error, 1
  set_plot, 'z'
  tek_color
  device, set_resolution=[640,240]

  pdsplot1, 'zand'
  pdsplot2, 'zand'

  return
end