多变量单时次观测数据grads绘图求助

liuzf406

第7步, 编写 gs 脚本进行插值后作出图形, 脚本如下:
'reinit'
'set grads off'
'open c:/isotope_grads/grid.ctl'
'open c:/isotope_grads/sta.ctl'
'set lon 70 140'
'set lat 15 55'
'set mpdset hires'
'define grid1=oacres(d18o,dd18o.2)'
'define grid2=oacres(d2h,d2h.2)'
'define grid3=oacres(dex,dex.2)'
'define grid4=oacres(annp,annp.2)'
'define grid5=oacres(annt,annt.2)'
'define grid6=oacres(dd18o,dd18o.2)'
'set gxout shaded'
'page q1'
'd grid1'
'draw title d18o'
'set gxout contour'
'cnbasemap rgrid'
'cbarn'

'set gxout shaded'
'page q2'
'd grid2'
'draw title d2h'
'set gxout contour'
'cnbasemap rgrid'
'cbarn'

'set gxout shaded'
'page q3'
'd grid3'
'draw title dex'
'set gxout contour'
'cnbasemap rgrid'
'cbarn'

'set gxout shaded'
'page q4'
'd grid6'
'draw title dd18o'
'set gxout contour'
'cnbasemap rgrid'
'cbarn'
;
运行后竟然出图了, 啊, 真不敢相信,本以是求大家帮忙看问题的. 前两天一直没弄出来, 感谢楼上两位的指导啊, 正如版主"做个霸气的木头"所说, 可能是我的以前你grid.ctl文件出错了. 再次对"清风"和"做个霸气的木头"表示深深的感谢.
我不懂fortran, 只是照搬清风的. 对grads,也只是初学者.
上图让大家看看吧: