program getdata c c program to retrieve Jones and Briffa (1992) CRU c land/sst monthly anomalies 1854-1993 c c INTEGER i, j, k, DATA(72), month,year real lon(72),lat(36) open (unit=4,file='anomalies-new', $ status='old') do i=1,72 lon(i)=2.5+(i-1)*5.0 end do do j=1,36 lat(j)=87.5-(j-1)*5.0 end do imonth = 12 iyear = 140 imax = imonth*iyear do n=1,imax read (4,*,END=800) month,year do k=1,36 read (4,550) (DATA(j), j=1,18) read (4,550) (DATA(j), j=19,36) read (4,550) (DATA(j), j=37,54) read (4,550) (DATA(j), j=55,72) end do end do 550 FORMAT (18I5) 560 format (2i7,2f8.1) 800 nmonth = n-1 write (6,*) nmonth, ' months of data read in.' close (unit=4) stop end