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

冬日情愫+

挺好的，值得好好学习！

Forcast

当年弄过，怎么过的都记不起来了，感觉还是有点麻烦的，一直想再把它弄通做个笔记，看来还要学习

qxtlyf

GRADS绘图需要亲自操作，搞明白每一步的用处和容易出现的错误，才能学会；GRADS绘图主要还是数据的处理比较难，有数据就好绘图

sskkyy2008

网格的描述文件：
dset          d:\drought\grid.grd
undef        -999.0
title            Grid data
xdef          71          linear     70  1
ydef          41          linear     15   1
zdef           1           linear     1000    1
tdef            1          linear      jul1951     1mo
vars     1
g       0        99         grid data
endvars
脚本文件：
'reinit'
*'open d:\drought\precip.ctl'
'open d:\drought\grid.ctl'
'open d:\drought\index.ctl'
'enable print d:\drought\index.gmf'
*'set lon 73 135'
*'set lat 15 55'
*'set lat 20 36'
*'set lon 75 120'
*'set mpdset hires  sic'
*'province-basemap sic 1'
'define   a=oacres(g, ind .2, 1.5)'
'set gxout shaded'
'set gxout contour'
'd a'
'print'
'disable print'
;

1012893089

很有用

深度狙杀

zhaolanxiai 发表于 2011-6-27 16:19

                               
登录/注册后可看大图

我按照嘟嘟鱼的程序，怎么还是弄不出来呢？:Q

file:///C:/Documents%20and%20Settings/Administrator/My%20Documents/My%20Pictures/QQ截图未命名.png

gs不用写这句话的把

kittyhare

为什么我画出的图线没有小数点

时光⑩漫步

顶起guoyipeng！！！

mmpyb

第一个程序运行时出现了好多错误，

mmpyb

D:\wsh\gs\zhdian.for(13) : Warning: In the call to STNTOGRD, actual argument #1 does not match the type and kind of the corresponding dummy argument.
      call stntogrd(ind)