nc转dat (grads)后，作图是图形不能重合

胖胖

下载了一整年的温度数据格点数据（欧洲中心），具体格式如下

把此nc文件转换为grads可读的dat,使用gs如下：

'reinit'
'sdfopen J:\pang\1.nc'
'set gxout fwrite'
'set fwrite J:\pang\1.dat'
mt=1
while (mt<=365)
mz=1
while(mz<=37)
'set t 'mt''
'set z 'mz''
'set lon 69.75 140.25'
'set lat 15 55.5 '
'd t.1'
mz=mz+1
endwhile
mt=mt+1
endwhile
'disable fwrite'
'reinit'

得到dat 后描述文件为：

dset J:\pang\1.dat
undef -9.96921e+36
title temperature
xdef  95 linear 69.75 0.75
ydef   55 linear 15 0.75
zdef 37 levels 1000 975 950 925 900 875 850 825 800 775 750 700 650 600 550 500 450 400 350 300 250 225 200 175 150 125 100 70 50 30 20 10 7 5 3 2 1
tdef 365 linear JAN2001 1dy
vars   1
t 37 -999 tempetature
Endvars

sdfopen nc和open dat，作图得到的是两张完全不同的图。

这是哪里出了问题，以为是缺测值设置错误，试了好几个都不对。求解？

↘禹禹獨行…

'set lon 69.75 140.25'
'set lat 15 55.5 '
把这个放在第四行试试。
'd t.1'不应该是'd t'吗？

胖胖

↘禹禹獨行… 发表于 2013-9-26 22:39

                               
登录/注册后可看大图

'set lon 69.75 140.25'
'set lat 15 55.5 '
把这个放在第四行试试。

原来在外面呢，经纬度放在哪里，都是这个结果。t.1（温度）是因为时间的默认也是t，才这么写。我试试改成t 哦。

river

把这两句
'set lon 69.75 140.25'
'set lat 15 55.5 '
改成
'set x 1 95'
'set y 1 55'
试一下看看

胖胖

river 发表于 2013-9-26 22:44

                               
登录/注册后可看大图

把这两句
'set lon 69.75 140.25'
'set lat 15 55.5 '

原来用的就是这个，不到才改的，还是这样。这两个写法有什么区别吗，我觉得表达都是同一个意思。

胖胖

问题解决了，在d变量的时候，t.1改t，就可以。

冬日情愫+

脚本的确没有发现有什么问题。楼主可以试试输出数据看看。

胖胖

冬日情愫+ 发表于 2013-9-26 23:35

                               
登录/注册后可看大图

脚本的确没有发现有什么问题。楼主可以试试输出数据看看。

把d t.1改为 d t,就可以，我自己也觉得奇怪呢。t.1输出的比t小些，并不是时间。

river

胖胖 发表于 2013-9-26 22:49

                               
登录/注册后可看大图

原来用的就是这个，不到才改的，还是这样。这两个写法有什么区别吗，我觉得表达都是同一个意思。

表达的意思是同一个的，但是直接设置经纬度很容易出现错位的情况，是因为有的时候有的人想当然的设置经纬度，没有考虑到格点资料有没有这个点

虹hhh

你好 请问用你类似的程序运行之后，dat文件中为什么时间层上 只读取时间t=1这一层，t=2以后的都没提取