program readdump PARAMETER ( NI = 144, NJ = 73, nd=592) real a(ni,nj) real prate(ni,nj,nd) open(11,file='dump_prate.m02.may.1981.cfs.data', 1 form='unformatted',status='old') do iday=1,nd read(11) a do i=1,144 do j=1,73 prate(i,j,iday)=a(i,j) c from S Pole to N Pole end do end do end do stop end