program readmask

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

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

      parameter (nlon=720, nlat=353)

      real xlon(nlon), ylat(nlat)
      integer mask(nlon,nlat)

      character dfname*80

      dfname = 'mask_blended_winds.v02.dat'
      write(6,1010) dfname
 1010 format("input= ",a80)
      open(11,file=dfname,form='formatted')

      read(11,'(i8)') nlonin      
      read(11,'(i8)') nlatin
      write(6,*) 'Input: nlon=',nlonin,' nlat= ',nlatin
      if (nlonin.ne.nlon .or. nlatin.ne.nlat) then
         write(6,*) 'Error in input dimensions: nlon,nlat=',nlon,nlat
         stop
      endif      
      read(11,'(6f13.5)') xlon
      read(11,'(6f13.5)') ylat
      write(6,*) 'lon=',xlon
      write(6,*) 'lat=',ylat

      read(11,'(80i1)') mask


c     print mediterranean black sea, and caspian sea (0-55E)

c     write(6,*) 'write from ylat=',ylat(237),' to ',ylat(271)
c     do j=271,237,-1 
c        write(6,'(115i1)') (mask(i,j),i=1,110)
c     enddo

      close(11)
      stop
      end