;read Rayleigh Lidar netcdf files

pro readRayleighLidarNetcdf, file, title, numt, numz, month, day, year, hour, z, gvlat, gvlon, gvalt, temperature, temperr


id=NCDF_OPEN(file)
NCDF_ATTGET, id, /GLOBAL, 'Title', title
NCDF_ATTGET, id, /GLOBAL, 'Month', month
NCDF_ATTGET, id, /GLOBAL, 'DayOfMonth', day
NCDF_ATTGET, id, /GLOBAL, 'Year', year
NCDF_ATTPUT, id, /GLOBAL, 'NumTimes', numt
NCDF_ATTPUT, id, /GLOBAL, 'NumAlt', numz

title=string(title)
month=string(month)
day=string(day)
year=string(year)
print,'Global Attributes:'
for i=0,5 do begin
  result = NCDF_ATTNAME(id, i,/global)
  print,i,' ',result
endfor
print,'Variable Names:'
for i=0,6 do begin
  result = NCDF_VARINQ(id, i)
  NCDF_ATTGET, id, i, 'units',units

  print,i, '. ',result.name, ', ',string(units)
  
endfor

tempid = NCDF_VARID(id, 'Temperature')
NCDF_VARGET, id, tempid, temperature
temperrid = NCDF_VARID(id, 'TempErr')
NCDF_VARGET, id, temperrid, temperr
timeid = NCDF_VARID(id, 'Hour')
NCDF_VARGET, id, timeid, hour
zid = NCDF_VARID(id, 'Altitude')
NCDF_VARGET, id, zid, z
gvaltid = NCDF_VARID(id, 'GValt')
NCDF_VARGET, id, gvaltid, gvalt
latid = NCDF_VARID(id, 'Latitude')
NCDF_VARGET, id, latid, gvlat
lonid = NCDF_VARID(id, 'Longitude')
NCDF_VARGET, id, lonid, gvlon

NCDF_CLOSE, id

return

end