Here is a fortran snippet to read the formatted ascii data: dimension antarc(360,51) do iyear=1950,2002 do imonth=1,12 read (11,2099) nmonth,nyear read (11,2100) ((antarc(i,j),i=1,360),j=1,51) end do end do 2099 format(i4,i6) 2100 format(8f15.3) The gridpoint antarc (1,1) should be at 0E, 90S.