如何把中国区域的降水格点资料 cmap 2.5x2.5的txt数据转换成grads画图用的dat数据

云大小子

我的降水资料 cmap资料 是全球的2.5x2.5资料，请教各位，有什么程序可以将格点的降水资料 转换成grads画图用的dat数据

言深深

http://bbs.06climate.com/forum.p ... &extra=page%3D1看看这里面有没有你想要的

言深深

http://bbs.06climate.com/forum.p ... &extra=page%3D1看看这里面有没有你想要的

mofangbao

如果能把数据用fortran读进去就好转了 转grads的方法是通用的，格点和站点的两种方法都有分享过

云大小子

mofangbao 发表于 2011-12-2 21:54

                               
登录/注册后可看大图

如果能把数据用fortran读进去就好转了 转grads的方法是通用的，格点和站点的两种方法都有分享过

我用那种把站点资料转成dat数据的方法 读取网格点得txt数据 得到了一个dat数据，但是我写的一个ctl文件访问不了 说是missing xdef record  missing ydef record  等等

云大小子

言深深 发表于 2011-12-2 21:54

                               
登录/注册后可看大图

http://bbs.06climate.com/forum.php?mod=viewthread&tid=623&extra=page%3D1看看这里面有没有你想要的

我看了觉得读取站点 格点方法应该是一样的 但是我得到了数据结果后写了一个ctl文件打不开dat数据

云大小子

dset d:\dong.dat
title monthly jiangshui
undef -999999.
xdef 144 linear 1.25 2.5
ydef 72 -88.75 2.5
zdef 1 linear 1000 1
tdef 2 linear 00Z01jan2005 1mo
vars 1
pred        0   99    rainfall
endvars
这就是我的ctl文件

云大小子

dong.txt (525.48 KB, 下载次数: 31)

2011-12-2 22:04 上传

点击文件名下载附件 以下是我自己写的读取数据的fortran程序
   program main
  implicit none
  real :: iyear,imonth,rlat,rlon,rval,tim,iyrold,imnold
  
  integer  nflag ,nlev,iflag
   
  open (1,file='d:\xia.txt',status='old')
  open (2,file='d:\xia.dat',status='new',form='binary')
           
          nflag=0
   10    read(1,*,end=90)iyear,imonth,rlat,rlon,rval
           if(iflag.eq.0) then
          iflag=1
          iyrold=iyear
          imnold=imonth
           endif
       if (iyrold.ne.iyear.or.imnold.ne.imonth) then
          nlev=0
       write(2) rlat,rlon,tim,nflag
    endif
             iyrold=iyear
        imnold=imonth
              tim=0.0
      nlev=1
      iflag=1
      write(2)rlat,rlon,tim,nlev,nflag
      write(2)rval
  goto 10
  90    continue
          nlev=0
     write(2) rlat,rlon,tim,nlev,nflag
    stop
    close(1)
    close(2)
   end
降水格点资料 也附上  ，直接运行程序能得到结果 ，但是我的ctl文件却打不开dat数据 错在哪里呢

mofangbao

云大小子 发表于 2011-12-2 22:00

                               
登录/注册后可看大图

dset d:\dong.dat
title monthly jiangshui
undef -999999.

写格点的方法和站点不一样，比站点简单多了，你这是格点资料，用些格点的方法就行了
比如
open(1,file='test.grd',status='replace',form='binary')
do i=1,t
  do j=1,lat
    do k=1,lon
      write(1)val(i,j,k)
     enddo
  enddo
enddo
close(1)
大概是这个样子，顺序永远是 x y z var t  越后面的在循环中越外层，所以如果有高度层的话也是按照这个规则的，这个在中科院的一本 实用绘图里面有讲到类似的
   

mofangbao

云大小子 发表于 2011-12-2 22:05

                               
登录/注册后可看大图

以下是我自己写的读取数据的fortran程序
   program main
  implicit none

因为你按照站点去写格点数据了，改成我刚才贴的那个样子试试吧，你的ctl没啥问题，除了ydef少了一个 linear