利用站点资料画出等值线图的GrADS示例

qhmemory

查了好多文献，在有些过程上说的很模糊，让我云里雾里！感谢楼主，讲的很细！

空想夜

artmunich 发表于 2011-7-5 00:47

                               
登录/注册后可看大图

发现了文章中的一些小疏漏，嘿嘿
1，在写grid.grd的ctl文件时，vars后面有个number没写，如下
vars

不明白版主所说的：vars后面有个number没写，他后面不是有个1么？然后我按着他说的：《接着我们需要给该二进制数据写一个描述性文件如下（新建一个写字板，取名为"index.ctl"）》：这个就不懂了，“新建一个写字板”？？？ctl文件不是个程序文件么？可以用写字板新建？

空想夜

嘟嘟鱼 发表于 2011-7-5 09:13

                               
登录/注册后可看大图

很好，谢谢你指出文中的错误！当时写的太粗心了！

不好意思，后来我明白问题所在了

xu850202

说得好详细好，必须顶个

耐宁.

嘟嘟鱼 发表于 2011-6-27 17:28

                               
登录/注册后可看大图

插值半径不能取的太大，否则即使插出来，也是不真实的值，这个需要你不断的尝试然后再跟实际情况对比，找 ...

楼主  我是初学者  很感谢你能提供这么一个学习的个例
有个问题 在十进制转二进制的过程中
open(3,file='indexs.grd',form='binary')
tim=0.0
nlev=1
nflag=1
do i=1,160
write(3) stid(i),lat(i),lon(i),tim,nlev,nflag,x(i)
end do
nlev=0
write(3) stid(i-1),lat(i-1),lon(i-1),tim,nlev,nflag
close(3)

问题：
1.这个里面stid是什么啊？（是站点号？）
2.多出了变量tim nlev nflag，不明白为什么要把的这几个变量写进去
3.为什么还要分nlev=1和nlev=0两种情况 而且在nlev=0的情况下为什么stid(i-1),lat(i-1),lon(i-1),都减了1

初学者 尚有不明之处
希望楼主看到此贴能给予回复
非常感谢

耐宁.

我按照楼主的步骤全部做了一遍
出现了与上面同样的问题
就是GS文件不出图


于是我把GS文件中的指令一条一条的在GRADS编辑器上实验
发现grid出现了问题， 在打开grid.ctl的时候 系统报错 说无法打开该文件

是不是grid.ctl有问题？

我的2个CTL分别为




GS内容为
'reinit'
'open    d:/drought/grid'
'open    d:/drough/index'
'enable print d:/drought/index.gmf'
  'set lon 73 135'
'set lat 15 55'
'set mpdset hires  cnworld'
'set gxout shaded'
'set gxout contour'
  'd oacres(g,ind.2,1.5)'
  'print'
   'disable print'
   
    就是这样了
出不了图 图示空白的

耐宁.

坐等楼主  和斑竹大神们回复

dcy_23

{:eb513:}{:eb513:}{:eb513:}{:eb513:}

lishidi1988

非常好用。多多支持

Abrahamhong

我也很想学习GrADS画图，不知从何学起，有没有这方面的书？