求助将高斯格点数据转常规dat数据

愉悦

最近在使用ncep的地面感热和潜热通量数据来求地面热通量，但是遇到了个问题，就是这个数据是gauss格点数据，是分辨率为1.875的数据，因为要求几十年的平均，所以想到把数据进行双线性插值成2.5*2.5的以便用fortran求这几十年的日平均，但是在数据格式转换上就出问题了。首先想到的是用grads的re()语句nc转成dat，同时也把数据转成了常规用的格距的数据，但是发现第一个时次还基本吻合，第二个时次就完全不同了呢。求有懂的人指教一下。反复翻了re()语句的规则，都觉得没有错了呢
用的gs如下：
'reinit'
'sdfopen L:\try\ghflux\gflux.sfc.gauss.2008.nc'
'set gxout fwrite'
'set fwrite L:\try\ghflux\25gflux.08.dat'
'set x 1 192'
'set y 1 94'
'set z 1'
it=1
while(it<=31)
'set t 'it''
'd re(gflux,144,linear,0,2.5,73,linear,-90,2.5,bl)'
it=it+1
endwhile
'disable fwrite'
;
还有就是我把这种分辨率的数据转换成dat数据，即便是不进行双线性插值我发现也有问题呢？是不是我的fwrite的程序有问题？不知道有没有谁用过高斯格点数据，对这类数据要求某年和30年平均的逐日距平该如何处理？

river

这个NC资料可以直接在grads里面用的吧，那就直接使用grads计算出结果接行了啊。grads求个平均什么的还是不在话下的

愉悦

river 发表于 2013-10-28 20:00

                               
登录/注册后可看大图

这个NC资料可以直接在grads里面用的吧，那就直接使用grads计算出结果接行了啊。grads求个平均什么的还是不在 ...

可是要用到三十年的数据 总共就有60个文件，grads不能打开这么多吧？

river

愉悦 发表于 2013-10-28 20:21

                               
登录/注册后可看大图

可是要用到三十年的数据 总共就有60个文件，grads不能打开这么多吧？

批量描述，听说过的吧？

愉悦

river 发表于 2013-10-28 21:56

                               
登录/注册后可看大图

批量描述，听说过的吧？

搞定了哈，谢谢。嘿嘿，不过那个re()还是没有搞懂呢~

deeli

敢问楼主怎么解决的啊~我也遇到问题了哇~~~

891971604

好多不懂啊，怎么办

spacer

891971604 发表于 2014-2-12 04:10
好多不懂啊，怎么办

可以直接写ctl。

HymntotheSea

楼主re的问题解决了么

愉悦

HymntotheSea 发表于 2014-9-20 19:30
楼主re的问题解决了么

没有 你成功了吗？