求救，nc转gs出现错误

wssummer

各位大神，额，这个问题从昨晚一直纠结着，怎么尝试都不对啊，论坛里貌似也没出现过这样的错误。请大家帮忙看看，指点一下

我把nc文件转成dat文件，这是我的nc文件的q ctlinfo结果



然后是我的转化的gs文件


但是，我把生成的dat文件，写了一个ctl，然后画出来，结果是不正确的


这是我的ctl


这个是我设t=1时，画出来的全球降水距平分布图




这个是set t=100时刻的，正常情况下应该跟t=1时差不多的，但是经度却差了很多，而且等值线形状是一模一样的


所以我猜是不是因为我x设错了，请大神们看看，是不是我的转化的gs文件出了错误
因为我把转化好的dat文件用fortran读写了一下，发现每个时刻的值都是一模一样的，哭死了

herenso

我觉得gs文件写错了吧，这样写的话岂不是经纬度成了外循环，时间成了内循环了。应该把t循环放到最外面，set x/y/z放到里面，这样写出来的dat才是按先时间后空间的顺序排列的，符合ctl的描述。

herenso

我觉得gs文件写错了吧，这样写的话岂不是经纬度成了外循环，时间成了内循环了。应该把t循环放到最外面，set x/y/z放到里面，这样写出来的dat才是按先时间后空间的顺序排列的，符合ctl的描述。

river

gs没看出什么问题啊，楼主把循环稍微改一下试试吧比如改成
i=1
while(i<=765)
'set t 'i''
'd precip'
i=i+1
还有就是楼主尽量不要用简写吧'set t' t，有的时候多一个空格少一个空格的就不对了

做个霸气的木头

你这个问题好生奇怪？话说为什么1和100的局平会一样？

mofangbao

应该就是3楼说的那个设置少了个空格

wssummer

mofangbao 发表于 2013-3-22 08:43

                               
登录/注册后可看大图

应该就是3楼说的那个设置少了个空格

搞定了，我又二了，gs里不能用循环，set t 1 765就好了，TT

kongfeng0824

请问楼主你的降水数据是在哪里下载的呢？是日降水量还是月值数据呢？

wssummer

kongfeng0824 发表于 2013-3-25 08:35

                               
登录/注册后可看大图

请问楼主你的降水数据是在哪里下载的呢？是日降水量还是月值数据呢？

necp/ncar网站上可以下载的，是月平均场

kongfeng0824

wssummer 发表于 2013-4-10 15:15

                               
登录/注册后可看大图

necp/ncar网站上可以下载的，是月平均场

哦。其实我是想问问日值的数据的下载地址的呵呵