program to_read_int

 integer, parameter :: nlat = 180 * 20, nlon = 360 * 20, &

        recl = 2*nlat*nlon

 integer(2), dimension(nlat,nlon) :: data

 open(11,file='Npp_0.05deg_mean.int16',form='unformatted',&

     access='direct',recl=recl)

 read(11,rec=1)data

 close(11)

 do i = 1, nlon

   do j = 1, nlat

      write(*,*)data(j,i)

   end do

 enddo

 end