pro l2aerosol_wrapper

;  This example routine demonstrates use of the HDF-EOS read utilites with a
;  MISR Level 2 aerosol file to access a data field, global attributes, and
;  grid attributes.  Once the stop statement is reached, the variables
;  containing the values read can be examined in the idlde variable watch
;  window or with the IDL HELP command.

l2aerfile = DIALOG_PICKFILE(TITLE="Select MISR Level 2 aerosol data file:")
grid = "RegParamsAer"
field = "RegBestEstimateSpectralOptDepth"

;  get the start and end block numbers for the file

status = read_global_attr(l2aerfile,"Start_block",start_block)
status = read_global_attr(l2aerfile,"End block",end_block)
start_block = start_block[0]
end_block = end_block[0]

;  read the aerosol optical depth

status = read_grid_field(l2aerfile,grid,field,optdepth_data)

;  get the data resolution (in meters) for this grid (all fields
;  in the grid are at the same resolution)

status = read_grid_attr(l2aerfile,grid,"Block_size.resolution_x",xres)
status = read_grid_attr(l2aerfile,grid,"Block_size.resolution_y",yres)

stop
end