GrADS如何输出text，csv文件； 和如果快速获取每日平均值

rouqiudog

新手GrADS，已经逐渐感觉GrADS的方便了。

1.  今天设置了单个像元的长时间序列数据，想保存出来Text，费了半天劲没有搞定。 在网上查了很多，没有解决这个问题。
也试了那个网上有的fprintf.gs 但是保存出来的只有2行，保存不全。

请问哪位大侠帮忙解答一下，指点一下如何把数据输出为 TXT文件。 谢谢。


2. 另外GrADS的平均值ave计算的时候，我看到是设置t=1到end来限定需要平均的时段。 由于我的数据是hourly数据，而我只需要计算出日平均，那么需要设置t=1, t=24 (代表第一天的)。这样不能每一天，就写一个吧。

请问有没有简单的方面 直接得出 365天每天的平均值啊，和 每月的平均值啊？

谢谢！

river

关于第一个问题，楼主可以打开gs文件自己去看用法



关于第二个问题，楼主可以利用while循环来每次增加24就可以实现

mofangbao

我能说这两个问题论坛都有现成的解决方案么，不过得你自己找了

rouqiudog

river 发表于 2014-11-25 08:19
关于第一个问题，楼主可以打开gs文件自己去看用法

非常感谢 River的回复。 问题已经解决，我发现我之前在网上粘的fprintf有问题，这次试用了 整合版本里面自带的fprintf，搞定了。

第二个问题，我试试。 我还想问一下，如果文件有缺失是否影响GrADS对于t的对应和设置？ 这个问题应该会影响到平均值的计算。

具体情况是这样的：
我原则上我的数据是1年的每天24小时的文件 （24 * 365 = 8760个），但是，在随机的不连续的若干天上，文件有缺失，有的是20个，有的是15个之类的。

那么我为了方便读取设置CTL为:
DSET ^data_2011%j3_%h2.dat
TITLE data
options little_endian yrev
options template
UNDEF -9999.0
XDEF     10  LINEAR   115.  0.125
YDEF     10  LINEAR    35.75  0.125
ZDEF 1 LINEAR 1 1
TDEF 8760 LINEAR 00:00Z01jan2011 1hr
VARS 1
inso  1  99  solar radiation            
ENDVARS

请问GrADS会自动跳过没有的数据吗？ 也就是说我这里设置的是8760个tdef，但是实际上的文件数为8000个。请问GrADS会识别正确对应的t和文件名吗？
（例如 t=1 应该对应的是2011年1月1日的0小时数据，即使这个数据是缺失的话；真正有的数据从2011年1月1日的1小时数据， 请问这个时候GrADS对应的是t=2 还是t=1）.

谢谢！

rouqiudog

mofangbao 发表于 2014-11-25 08:46
我能说这两个问题论坛都有现成的解决方案么，不过得你自己找了

谢谢管理员。
我也是先在论坛里搜索查找，没有找到才提问的。 可能是答案包含在其他的问题了。

问题已经解决一部分。

river

rouqiudog 发表于 2014-11-25 18:26
非常感谢 River的回复。 问题已经解决，我发现我之前在网上粘的fprintf有问题，这次试用了 整合版本里面 ...

很不幸，GrADS没有那么智能，你设置了8760个tdef，但是实际上的文件数为8000个，那么GrADS在计算过程中会报错的吧·······

rouqiudog

river 发表于 2014-11-25 21:29
很不幸，GrADS没有那么智能，你设置了8760个tdef，但是实际上的文件数为8000个，那么GrADS在计算过程中会 ...

谢谢您的回复。

我自己试验了一下，GrADS可以的，确实是智能了啊。

我真实的数据开头几天都是完整的，中间和最后几天都是没有文件，发现使用哪个tdef 8760，提取的结果是正确的，在没有数据的那几天提取的直接就是undef了。

确实GrADS惊艳了一把。

river

rouqiudog 发表于 2014-11-25 23:44
谢谢您的回复。

我自己试验了一下，GrADS可以的，确实是智能了啊。

我倒是没有试过，还以为没有数据的那几天会提示出错呢······
看来grads还是有好多高级的功能我们都不知道呢。
我觉得你可以把这个当成一个GrADS的一个小技巧，发个帖子让大家学习一下。