关于多时次插值后保存问题咨询大家

男紫汗

这两天在做站点插值，幸亏家园的诸多教程，目前已经到了插值阶段。我的资料是气候中心的64年台站降水资料，想要插值后保存一下各时次（不知道先计算再插值画图关系大不大），论坛上有很多分享如清风大神http://bbs.06climate.com/forum.php?mod=viewthread&tid=7357（这里貌似是单时次），和http://bbs.06climate.com/forum.php?mod=viewthread&tid=10006&extra=&page=1（这里貌似没看到保存）等等。技术渣实在无能为力了，恳请大家指点！‘’‘’‘’‘’‘’‘’‘’‘’‘’‘’
之前一直在Fortran下，忘了去grads那里看看了，看了下ave.grd生成了。。。昨晚上生成了一个ave.grd,今天发现前面步骤有点问题，改了之后运行gs就生成不了。要么就是这种问题(图片），后来又版搜了下，改上面第二个链接中'open grid.ctl'
'open station.ctl'
'set lat 15 55'
'set lon 70 140'
i=1
while (i<=62)
'set t 'i''
'set gxout fwrite'cccccccccccccccccccccccccccccc一般set fwrite 和set gxout fwrite 是反过来写的，在这里在颠倒顺序，并且在多次里需要将set fwrite写到时次的循环里，才能得到准确的62个时次的数据，否则只有同样的值或者一个时次的值。
'set fwrite ave.grd'

'define a=oacres(g,t.2,1.5)'
'define a1=maskout(a,g-0.5)'
'define aa=smth9(a1)'
'd aa*0.1'
i=i+1
endwhile
'disable fwrite'
;
格式如下
'reinit'
'open D:\fortran\station\station\stationgrid.ctl'
'open D:\fortran\station\station\station.ctl'
'set lon 70 140'
'set lat 15 55'
'set z 1'
'set fwrite D:\fortran\station\station\ave.grd'
'set gxout fwrite'
i=1
while(i<=768)
'set t 'i
'define a=oacres(g,rain.2)'
'define a1=maskout(a,g-0.5)'
'define aa=smth9(a1)'
'd aa'
i=i+1
endwhile
'disable fwrite'
;

就生成了grd，应了第二链接中Sylvia所说的可能和电脑有关

四叶草

http://bbs.06climate.com/forum.p ... p;extra=&page=1
这个帖子里已经讲了站点插值为格点数据，并保存多个时次格点数据的方法，不是楼主想要的？

男紫汗

四叶草 发表于 2015-12-15 21:11
http://bbs.06climate.com/forum.php?mod=viewthread&tid=10006&extra=&page=1
这个帖子里已经讲了站点插 ...

一直在Fortran下处理，忘了去grads文件夹看了

四叶草

男紫汗 发表于 2015-12-16 14:09
一直在Fortran下处理，忘了去grads文件夹看了