pro l2land_wrapper

;  This example routine demonstrates use of the HDF-EOS read utilites with a
;  MISR Level 2 land surface 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.

filename = dialog_pickfile(TITLE="Select MISR Level 2 Land Surface File")
gridname = 'SubregParamsLnd'
fieldname = 'LandBRF'

;  get the start and end block numbers for the file

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

;  read the first five blocks of the Land BRF data field

status = read_grid_field(filename,gridname,fieldname,landbrf_data,BLOCKRANGE=[start_block,start_block+4])

;  get the resolution (in meters) from grid attributes along with
;  the specific grid attributes available for Land BRF

status = read_grid_attr(filename,gridname,"Block_size.resolution_x",xres)
status = read_grid_attr(filename,gridname,"Block_size.resolution_y",yres)

status = read_grid_attr(filename,gridname,"Min LandBRF",landbrf_min)
status = read_grid_attr(filename,gridname,"Max LandBRF",landbrf_max)
status = read_grid_attr(filename,gridname,"Scale LandBRF",landbrf_scale)
status = read_grid_attr(filename,gridname,"Offset LandBRF",landbrf_offset)
status = read_grid_attr(filename,gridname,"Fill LandBRF",landbrf_fillvalue)
status = read_grid_attr(filename,gridname,"Underflow LandBRF",landbrf_underflowvalue)
status = read_grid_attr(filename,gridname,"Overflow LandBRF",landbrf_overflowvalue)

stop

end