; This procdure reads the TEMPO data and returns the variables 
; stored in a data structure. The inputs are file path

PRO read_TEMPO_hourly, files, output

  id = ncdf_open(files, /nowrite)
  parentId = NCDF_GROUPSINQ(id)                ;get the parent directory ID's

  ; verifies that ParentID[0] is the 'PRODUCT' subdirectory
  ;for i=0, n_elements(parentId)-1 do begin
  ;  PRODUCT_NOMENCLATURE=ncdf_groupname(ParentID[i])
  ;  print, i, PRODUCT_NOMENCLATURE
  ;endfor
 
   ; read PRODUCT/ variables 
   SUPPORT_ID = ncdf_groupsinq(ParentID[0])

   ncdf_varget,ParentID[0],'latitude', lat
   ncdf_varget,ParentID[0],'longitude', lon
   ncdf_varget,ParentID[0],'time', time
   ncdf_varget,ParentID[0],'longitude_bounds', lon_bnd
   ncdf_varget,ParentID[0],'latitude_bounds', lat_bnd

   ntimes = n_elements(time)
   nlon = n_elements(lat[*,0])

   ncdf_varget,ParentID[0],'solar_zenith_angle', sza
   ncdf_varget,ParentID[0],'terrain_height', height
   ncdf_varget,ParentID[0],'tropopause_pressure',p_trop

   year   = '2014' ;strmid(files[nf], strpos(files[nf], 'V01')+4,4)
   month  = '08'; strmid(files[nf], strpos(files[nf], '2014')+4,2)
   day    = strmid(file_basename(files), strpos(file_basename(files),'2014')+6,2)
   hour   = strmid(file_basename(files), strpos(file_basename(files),'2014')+9,2)
   minute = '00';strmid(files[nf], strpos(files[nf], 'V01')+15,2)
   seconds= '00';strmid(files[nf], strpos(files[nf], 'V01')+17,2)
   swath  = '0';strmid(files[nf], strpos(files[nf], 'V01')+22,3) & swath = float(swath)
   ;print, year, month, day, ' Time: ', hour, ':', minute, ':',seconds, ' Swath: ', swath

   SUPPORT_ID1 = ncdf_groupsinq(ParentID[1])
   ncdf_varget,ParentID[1],'column_amount', column
   ncdf_varget,ParentID[1],'tropospheric_column_amount', trop_column
   ncdf_varget,ParentID[1],'stratospheric_column_amount', strat_column
   ncdf_varget,ParentID[1],'main_data_quality_flag',qa_flag


   SUPPORT_ID3 = ncdf_groupsinq(ParentID[3])
   ncdf_varget,ParentID[3],'amf_molecule_tropospheric', amf_molec
   ncdf_varget,ParentID[3],'amf_molecule_specific', amf_molec_spec

   SUPPORT_ID4 = ncdf_groupsinq(ParentID[4])
   ncdf_varget,ParentID[4],'column_amount', true_column
   ncdf_varget,ParentID[4],'gas_profile', true_gas_profile
   ncdf_varget,ParentID[4],'tropospheric_column_amount', true_trop_column

  ncdf_close,id


 output =  { lat:lat, lon:lon, time:time, ntimes:ntimes, nlon:nlon, lat_bnd:lat_bnd, lon_bnd:lon_bnd, $
	     year:year, month:month, day:day, hour:hour, minute:minute, seconds:seconds, swath:swath, $
	     column:column, qa_flag:qa_flag, trop_column:trop_column, strat_column:strat_column, $
	     amf_molec:amf_molec, amf_molec_spec:amf_molec_spec, $
	     true_column:true_column, true_gas_profile:true_gas_profile, true_trop_column:true_trop_column}

END