pro grainger,nopr,iout,indir,outdir,lstrn,iset,lset,$
  nlines,wcm,wavelength,rndat,ridat,iwave

; *****
; Will work with Roy Grainger's volcanic ash indices

; *****
; The input netCDF file
   fil=strarr(1)
   fil='grainger_volcanic_ash.nc'
   filnm=strcompress(indir+fil)

; For graphs
   lset=lstrn

; *****
; Read in from the netCDF file
   rdgrainger,nopr,iout,filnm,$
   ntext,text,$
   nlines,wcm,wavelength,rndat,ridat,rnerr,rierr

; *****
  return
  end