grads格点资料插值站点时为什么不能输出数据？

勇敢的悲伤

用grads将格点资料插值到站点的时候遇到一些问题，希望各位高手能给解析一下。
r.ctl为站点背景资料，14071914gedian.ctl为格点资料。
1、程序虽然能运行，但是运行之后并没有创建的14071914zhandian.dat文件出现
2、格点插值站点的时候不能用‘set z ’ zz 这种格式的循环吗？
3、用grads2ascii 能把不同高度的资料放在一个文件里面吗？
4、输出的站点资料除了站点号，经度、维度和变量的值外，还有一个值，这个值是什么呢，如果是高度的话，为什么输出所有文件里这个值都是固定的，没有变呢？
刚刚接触格点插值站点，所以问题有点儿多，包涵~~
'reinit'
'open d:\test\r.ctl'
'open d:\test\14071914gedian.ctl'

'set lon 110 116'
'set lat 30 37'
'set gxout fwrite'
'set fwrite d:\test\14071914zhandian.dat'

zz=1
while(zz<=11)
'd gr2stn(h.2(z='zz'),g1)'
'd gr2stn(t.2(z='zz'),g1)'
'd gr2stn(ttd.2(z='zz'),g1)'
'd gr2stn(windx.2(z='zz'),g1)'
'd gr2stn(wind.2(z='zz'),g1)'

*讲转成站点的数据输入到文档中
'grads2ascii gr2stn(h.2(z='zz'),g1) d:\test\h2'%zz%'.txt %8.0f'
'grads2ascii gr2stn(t.2(z='zz'),g1) d:\test\t2'%zz%'.txt %8.0f'
'grads2ascii gr2stn(ttd.2(z='zz'),g1) d:\test\ttd2'%zz%'.txt %8.0f'
'grads2ascii gr2stn(windx.2(z='zz'),g1) d:\test\windx2'%zz%'.txt %8.0f '
'grads2ascii gr2stn(wind.2(z='zz'),g1) d:\test\wind2'%zz%'.txt %8.0f'

zz=zz+1
endwhile
'disable fwrite'
;

mofangbao

楼主的步子跨的有点大
你为什么不觉得是你上面 d的时候就出了问题？

胖胖

我自己做过类似的工作，确实不能输出。就这我写过解决方法，很麻烦。当时提问过，还是没弄明白！

勇敢的悲伤

mofangbao 发表于 2014-7-27 09:42
楼主的步子跨的有点大
你为什么不觉得是你上面 d的时候就出了问题？

你是说'd gr2stn(h.2(z='zz'),g1)' ？这个有什么问题吗？

勇敢的悲伤

胖胖 发表于 2014-7-27 10:12
我自己做过类似的工作，确实不能输出。就这我写过解决方法，很麻烦。当时提问过，还是没弄明白！

哎，弄这个头都大了，难道只能一个一个输出，然后再用Fortran重新建一个文件吗。。。

lqouc

勇敢的悲伤 发表于 2014-7-27 11:18
哎，弄这个头都大了，难道只能一个一个输出，然后再用Fortran重新建一个文件吗。。。

关于兰溪的那个脚本，你的意思是不能循环往里面写数据？新版里增加了一个参数，可以不进行覆盖而是追加。具体用法看帖吧。
http://bbs.06climate.com/forum.php?mod=viewthread&tid=11515

勇敢的悲伤

lqouc 发表于 2014-7-27 11:29
关于兰溪的那个脚本，你的意思是不能循环往里面写数据？新版里增加了一个参数，可以不进行覆盖而是追加。 ...

恩恩。谢谢~用了这个之后可以循环输入了。
但是很奇怪高度场为什么一直都是500，而且一行只有4个数据就换行了。我已经定义过每行放的数据个数了，可还是没有成功。

四叶草

楼主可以试着用set gxout print 看看。

qx2864973225

回复一下 混个金钱

HAHAHAH

同样的困扰