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

liwei871204

算出的t检验值，是txt格式的， ttongji.txt (174.52 KB, 下载次数: 53)

2012-4-14 11:58 上传

点击文件名下载附件 。编了一个fortran程序，转化为grd格式。 tzhuan.f90 (290 Bytes, 下载次数: 53)

2012-4-14 11:59 上传

点击文件名下载附件 。之后为grd格式的写了一个ctl文件，dset e:\slp\ttongji.grdtitle tjiangyan tzhi
undef -999999.
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 1 linear 1000 1
tdef 1 linear 00Z01jan2005 1mo
vars 1  
t        0   99    ttongjizhi
endvars
为什么在grads中打开ctl文件，图并没有显示t的值。只有底图。
求各位大神帮忙。。。。谢谢。。。。。。

river

liwei871204 发表于 2012-4-14 13:38

                               
登录/注册后可看大图

txt中存放的 先是纬度的144个，再从经度1开始，到73，纬度应该写在里面，经度应该写在外面的啊，我是这么理 ...

呵呵，楼主啊，经度是从0°到360°，纬唯独是从-90°到90°，看也知道是哪个范围大了吧。怎么可能是纬度的144个，再从经度1开始，到73呢。还有就是grads能识别的二进制资料的排列必须是从西往东，从南到北，从底层到高层，这是死的。按照这个顺序fortran中的循环应该先是时间，再来层次，然后是纬度，最后才是经度

mofangbao

你的fortran写的很有问题，那个 144 和73的数组只给最后一个数字复制了，其余全部没有数值，而且楼主对GrADS中的数据存放方式是否了解？应该纬度在外，经度在内，在不了解这些基础东西之前就按照自己的想法写当然会出错了。。。  

liwei871204

txt中存放的 先是纬度的144个，再从经度1开始，到73，纬度应该写在里面，经度应该写在外面的啊，我是这么理解的。grads是初学，理解不了必须纬度在外，经度在内，是这么严格要求的么？

liwei871204

结果出来了，果然是前面的程序，只给最后一个数值赋值，其余的没有。我从新设置了一下循环。

Crystal_PP

麻烦楼主能不能分享一下求t检验值的程序？

liwei871204

是合成差值场的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)=(dat(i,j,11)+dat(i,j,12)+dat(i,j,17)+dat(i,j,18)+dat(i,j,19)+dat(i,j,20)+dat(i,j,22)+dat(i,j,45)+dat(i,j,51)+dat(i,j,58))/n2
   enddo
enddo
!!!!!!!!!!!!!!!!!!!!!       

!偏早年和偏晚年的差值
do j=1,y
   do i=1,x
      dif(i,j)=earlyave(i,j)-lateave(i,j)
   enddo
enddo
!!!!!!!!!!!!!!!!!!!!!

!偏早年和偏晚年的方值
do j=1,y
   do i=1,x
      earlyvar(i,j)=((dat(i,j,9)-earlyave(i,j))**2+(dat(i,j,14)-earlyave(i,j))**2+(dat(i,j,24)-earlyave(i,j))**2+(dat(i,j,27)-earlyave(i,j))**2+(dat(i,j,33)-earlyave(i,j))**2+(dat(i,j,35)-earlyave(i,j))**2+(dat(i,j,37)-earlyave(i,j))**2+(dat(i,j,42)-earlyave(i,j))**2+(dat(i,j,44)-earlyave(i,j))**2+(dat(i,j,52)-earlyave(i,j))**2)/n1
             latevar(i,j)=((dat(i,j,11)-lateave(i,j))**2+(dat(i,j,12)-lateave(i,j))**2+(dat(i,j,17)-lateave(i,j))**2+(dat(i,j,18)-lateave(i,j))**2+(dat(i,j,19)-lateave(i,j))**2+(dat(i,j,20)-lateave(i,j))**2+(dat(i,j,22)-lateave(i,j))**2+(dat(i,j,45)-lateave(i,j))**2+(dat(i,j,51)-lateave(i,j))**2+(dat(i,j,58)-lateave(i,j))**2)/n2
   enddo
enddo
!!!!!!!!!!!!!!!!!!!!       

!计算t统计量
do j=1,y
   do i=1,x
      ts(i,j)=(earlyave(i,j)-lateave(i,j))/(sqrt(((n1-1)*earlyvar(i,j)+(n2-1)*latevar(i,j))/(n1+n2-2))*sqrt(1/n1+1/n2))
   enddo
enddo
!!!!!!!!!!!!!!!!!!!!       
       
!资料输出       
open(200,file='d:\u200junedif.grd',form='binary')
do j=1,y
   do i=1,x
      write(200) dif(i,j)
   enddo
enddo
close(200)

open(300,file='d:\u200junet-test.grd',form='binary')
do j=1,y
   do i=1,x
      write(300) ts(i,j)
   enddo
enddo
close(300)
!资料输出完毕

stop
end        program main

天黑了又亮。

顶

东篱

顶

又是那隻貓

学习一下下