; HDF5 file testing
; Author: Alex Schaefer
; Last Modified: 29 June 2016


begin

  setfileoption("hdf5", "FileStructure", "Advanced")

  f = addfile("/Users/schaefer/Desktop/NASA/CATS/CATS-ISS_L1B_D-M7.2-V2-07.2015-08-26T16-27-44T17-14-51UTC.hdf5","r")

  ; Lats and Lons load correctly with 56476 points each

  lat = f->/geolocation/CATS_Fore_FOV_Latitude 
  lon = f->/geolocation/CATS_Fore_FOV_Longitude
  nlats = dimsizes(lat)
  print(nlats)

  ; Check how many lats and lons are non-zero
  print("The number of non-zero lats: " + tostring(num(lat.ne.0)))
  print("The number of non-zero lons: " + tostring(num(lon.ne.0)))
  

  ; Load the variable -- Should be 56476 x 533

  backscatter_1064 = f->/channel_1064/FFOV/Total_Attenuated_Backscatter1064_Fore_FOV 
  backscatter_dims = dimsizes(backscatter_1064)
  print(backscatter_dims)

  ; Checking to see how many of the points are non-zero values
  ; This takes awhile when checking, uncomment following line to check values

  ;print(num(backscatter_1064.ne.0))

end