【新人求教】画站点图出现entire grids undefined问题

糖份

       RT。
       各位大神，本人新人一个，grads 是自学的。在计算EU与全国160站降水的相关系数过程中，相关系数已经计算完毕，且在FORTRAN里面已经完成向按GRADS要求排列的二进制资料文件，现在就差用GRADS画图了。
      相关系数的ctl文件：
dset e:/rcoreu2.grd
dtype station
stnmap e:/temperature.map
undef -999.0
title temperature
tdef 1  linear jun1951 1yr
vars 1
cor  0 99  correlation coefficient
endvars
;

站点的ctl文件
dset  e:\grid.dat
title  grid
undef 99999.9
xdef 71 linear 70.0 1.0
ydef 41 linear 15.0  1.0
zdef 1  levels 1000   
tdef 1  linear jun1951 1yr
vars 1
g  0 99 grid data1
endvars
;

运行的gs文件

'reinit'
'open e:/grid.ctl'
'open e:/cor_t.ctl'
'set grid off'
'set lon 73 135'
'set lat 15 55'
'set mpdset hires'
'd oacres(g,cor.2)'
;


不知为何最后一运行就出现




请各位大神指教一二，谢谢啦
这个是画图用的数据文件 1.rar (6.12 KB, 下载次数: 3)

2012-11-25 18:38 上传

点击文件名下载附件

做个霸气的木头

站点数据处理很多环节，你的肯定是数据处理那块出问题了，因为单画站点图，只有一个站，还是-0，看看这个帖子自己分析先http://bbs.06climate.com/forum.php?mod=viewthread&tid=7503，

river

估计是计算相关系数写成二进制的时候出问题了吧。把fortran代码也贴上来吧

糖份

做个霸气的木头 发表于 2012-11-25 19:29

                               
登录/注册后可看大图

站点数据处理很多环节，你的肯定是数据处理那块出问题了，因为单画站点图，只有一个站，还是-0，看看这个帖 ...

你好，站点数据处理的程序我自己检查过了，看不出来错误啊。fortran检查不出错，运行也没有问题
      IMPLICIT NONE
        INTEGER  NX,NT,it
        PARAMETER(NX=160,NT=58)
        real  lat(nx),lon(nx)
        real  aa(61),eu(nt),t(nx,nt),s(nx),ave(nx,nt),rr(nx),ry(nx)
        real  ave_t(nx)
        real  ave_eu,rx
        REAL  A,B,X,Y,tim
        INTEGER  I,J,K,nlev,nflag
        character*8  stid(nx)

      open(10,file='e:/t1601.txt',STATUS='OLD')
      read(10, *) ((t(i,it), i=1, nx), it=1,nt)

        do i=1,nx
      do it=1,nt
        ave_t(i)=ave_t(i)+t(i,it)
        enddo
        ave_t(i)=ave_t(i)/nt
        enddo
       
      
      open(14,file='e:/eu.grd',form='binary',status='old')
      do it=1,61
        read(14)aa(it)
      enddo
      write(*,*)'read data ok'

      do it=1,nt-3
        eu(it)=aa(it+3)
        enddo

      do it=1,nt
      ave_eu=ave_eu+eu(it)
      enddo
      ave_eu=ave_eu/nt



      do i=1,nx
      do it=1,nt
      s(i)=s(i)+(eu(it)-ave_eu)*(t(i,it)-ave_t(i))
      enddo;enddo

      do it=1,nt
      rx=rx+((eu(it)-ave_eu)*(eu(it)-ave_eu))
      enddo

      do i=1,nx
      do it=1,nt
      ry(i)=ry(i)+((t(i,it)-ave_t(i))*(t(i,it)-ave_t(i)))
      enddo;enddo

!     相关系数
      do i=1,nx
      rr(i)=s(i)/(sqrt(rx)*sqrt(ry(i)))
      write(30)rr(i)
        write(*,*)rr(i)
      enddo
      write(*,*)'write cor to file ok'



        open(1, file='e:/Lonlat.TXT',STATUS='OLD')
        do k = 1, nx
        read(1,*) stid(k),  lat(k), lon(k)
                write(*,*) stid(k),  lat(k), lon(k)
        enddo

        open(21,file='e:/rcoreu2.grd',form='binary')
        open(22,file='e:/rcoreu2.txt')

   
      tim   = 0.0
      nlev  = 1
      nflag = 1

      do i= 1,nx
      write(21) stid(i), lat(i), lon(i), tim, nlev, nflag, rr(i)
      write(*,*) stid(i), lat(i), lon(i), tim, nlev, nflag, rr(i)
  
      nlev = 0
      write(21) stid(i-1), lat(i-1), lon(i-1), tim, nlev, nflag   
      enddo
      
      end

糖份

这个是用到的数据文件

糖份

river 发表于 2012-11-25 19:39

                               
登录/注册后可看大图

估计是计算相关系数写成二进制的时候出问题了吧。把fortran代码也贴上来吧

嗯，虽然不知道是不是这个问题，不过谢谢你啊。我已经重新换了个方法算了，这回算出来了。

mofangbao

      do i= 1,nx
      write(21) stid(i), lat(i), lon(i), tim, nlev, nflag, rr(i)
      write(*,*) stid(i), lat(i), lon(i), tim, nlev, nflag, rr(i)
  
      nlev = 0
      write(21) stid(i-1), lat(i-1), lon(i-1), tim, nlev, nflag   
      enddo
这一段的enddo，放到nlev=0前面去，你是对站点进行循环的，这样写就成了nx个时次，一个站点的数据了，实际上应该是nx个站点，一个时次下次贴代码要记得挑重点贴，你这样一通代码全贴上来，谁闲着很多时间帮你一句句的看呢？

学业有成

请问你的问题解决了？我也遇到此类问题，望解决 qq；455528540