PRO GRID_PROJECT_EXAMPLE, MAPGRID

;- Load the 25 km equal area gridded data (814880 32-bit signed integers)
;- (could use sds_read to read this from HDF)
openr, lun, 'nobs.dat', /get_lun
nobs = lonarr(814880)
readu, lun, nobs
free_lun, lun

;- Convert the data to float
nobs = float(nobs)

;- Create the array required by grid_resample
datagrid = fltarr(4, 814880)

;- Copy the number of observations into the array
datagrid[0, *, *] = nobs

;- Call grid_resample (only need to do this once for each datagrid)
grid_resample, datagrid, mapgrid

;- Call grid_project (can do this many times for each mapgrid)
grid_project, mapgrid

END