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

云大小子

mofangbao 发表于 2011-12-2 22:10

                               
登录/注册后可看大图

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

恩 好的 我试试

云大小子

mofangbao 发表于 2011-12-2 22:10

                               
登录/注册后可看大图

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

program main
       implicit none
      real :: lat,lon
       integer t
      real :: val，var
          open (4,file='d:\dong.txt',status='old')
        open(2,file='test.dat',status='new',form='binary')
       read(4,*) t,lat,lon,var   
                                                    我想令var是变量 ，但是val也是指降水量而且是数组，这个该怎么解决
          do t=1,2
            do lat=-88.75,88.75 ,2.5
              do lon=1.25  ,358.75 , 2.5
              write(2)val(t,lat,lon)
              enddo
              enddo
           enddo
            close(2)
            close(4)
             end

mofangbao

自己稍微思考下就知道了，通过我给的例子应该能举一反三的

云大小子

mofangbao 发表于 2011-12-2 22:49

                               
登录/注册后可看大图

自己稍微思考下就知道了，通过我给的例子应该能举一反三的

program main
implicit none
real :: lat,lon
integer t
real :: var
open (4,file='d:\dong.txt',status='old')
open(2,file='test.dat',status='new',form='binary')
do t=1,2
  do lat=-88.75,88.75 ,2.5
    do lon=1.25  ,358.75 , 2.5
      write(2)var
     enddo
  enddo
enddo
close(2)
close(4)
end
我这样行不 ，我试了一下，得到了一个数据但是感觉不对

luoziwuhui

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

                               
登录/注册后可看大图

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

因为你按照站点去写格点数据了，mofangbao 正解！
你的 ctl zdef 1 linear 1000 1
应该是 zdef 1 linear 1000 1000

冬日情愫+

对了，我没有查手册，貌似高度只有一层的时候用得是levels吧？ zdef 1 levels 1000 1000   是吗？

chuanhao

请问你这个资料在哪里下载的，能把网址发给我吗

chuanhao

你这个资料从哪里下载的？能告诉我怎么下载吗