forrest@0: program to_read_int
forrest@0: integer, parameter :: nlat = 180 * 20, nlon = 360 * 20, &
forrest@0:        recl = 2*nlat*nlon
forrest@0: integer(2), dimension(nlat,nlon) :: data
forrest@0: open(11,file='Npp_0.05deg_mean.int16',form='unformatted',&
forrest@0:     access='direct',recl=recl)
forrest@0: read(11,rec=1)data
forrest@0: close(11)
forrest@0: do i = 1, nlon
forrest@0:   do j = 1, nlat
forrest@0:      write(*,*)data(j,i)
forrest@0:   end do
forrest@0: enddo
forrest@0: end