PARAMETER(II=304, JJ=448, NREL=II*JJ*2) integer*2 data(ii*jj) real seaice(ii,jj) OPEN(UNIT,FILE='northern/monthly/bt_20060901_f13_v02_n.bin',status='old' & ,recl=nrel, access='DIRECT',form='unformatted') print*,'reading the data' read(unit,rec=1) data close(unit) do 100 i=1,ii do 100 j=1,jj L=ii*(j-1)+i seaice(i,j)=float(data(L)/10.) c print*,seaice(i,j) 100 continue print*,seaice(151,201) STOP END