;read land mask for blended winds, 0.5 x 0.5 degree grid, 
;from 0.5E, 88.0S to 360.0E, 88.0N:  land=1,ocean=0;
;file contains dimensions and grid info

;version 1.0 (Feb 5, 2004)
;version 2.0 (Feb 22, 2008)

fname = 'mask_blended_winds.v02.dat'

openr,lunit,fname,/get_lun

readf,lunit,nlon
readf,lunit,nlat

xlon = fltarr(nlon)
ylat = fltarr(nlat)
mask = intarr(nlon,nlat)   

readf,lunit,xlon
readf,lunit,ylat
readf,lunit,mask,format='(80i1)'
free_lun,lunit

print,'x lon=',xlon
print,'y lat=',ylat

end