python读取grd数据_一个关于grads显示grd格式数据的问题。

是合成差值场的t检验程序。

program main

implicit none

integer i,j,t,irec

integer,parameter :: x=145

integer,parameter :: y=73

integer,parameter :: n=58 !总体样本数

real,parameter :: n1=10.0 !偏早年份

real,parameter :: n2=10.0 !偏晚年份

real :: dat(x,y,n),earlyave(x,y),lateave(x,y),earlyvar(x,y),latevar(x,y),ts(x,y),dif(x,y)

!资料读取

open(100,file='d:\48-05u200june.grd',form='unformatted',access='direct',recl=x*y*4)

irec=1

do t=1,n

read(100,rec=irec) ((dat(i,j,t),i=1,x),j=1,y)

irec=irec+1

enddo

close(100)

!资料读取完毕

!偏早年和偏晚年的均值

do j=1,y

do i=1,x

earlyave(i,j)=(dat(i,j,9)+dat(i,j,14)+dat(i,j,24)+dat(i,j,27)+dat(i,j,33)+dat(i,j,35)+dat(i,j,37)+dat(i,j,42)+dat(i,j,44)+dat(i,j,52))/n1

lateave(i,j)=(da