program test parameter(NX=360,NY=181,NT=72) real water(NX,NY) open(1,file='water_ss_cdas_9398.dat',form='formatted') do IT=1,NT read(1,10)((water(i,j),i=1,nx),j=1,ny) enddo 10 format(10f7.2) stop end