nc转grd出错了吗？？

chaseechoes

nc转grd的gs文件：

'reinit'
'sdfopen f:\HadISST_ICE\HadISST_ice.nc'
'set gxout fwrite'
'set fwrite f:\HadISST_ICE\area.grd'
i=1
while(i<=1701)
'set x 1 360'
'set y 91 180'
'set t 'i
'set z 1'
'd sea_ice_area_fr'
i=i+1
endwhile
'close 1'
'disable fwrite'

数据是海冰密集度，数值在0~1，直接用nc文件画图是对的，但是用ctl画的不对
DSET f:\HadISST_ICE\area.grd
TITLE MONTHLY HadISST_ICE
UNDEF 32767
XDEF 360 LINEAR -179.5 1.0
YDEF 90 LINEAR 0.5 1.0
ZDEF 1 LEVELS 1000
TDEF 10 LINEAR JAN1870 1MO
VARS 1
FR 0 99 sea_ice_area_fr
ENDVARS

用grads    d fr能出图，但是数据不对，contouring: -1e+30 to 0 interval 1e+29
数据为什么会出问题啊，怎么会变得这么大？谢谢各位了~~file:///C:\Users\ZWX\Documents\Tencent Files\470721876\Image\6VW0IBJOTCVNIG6C$(6GLYM.jpg
file:///C:\Users\ZWX\Documents\Tencent Files\470721876\Image\6VW0IBJOTCVNIG6C$(6GLYM.jpg
file:///C:\Users\ZWX\Documents\Tencent Files\470721876\Image\6VW0IBJOTCVNIG6C$(6GLYM.jpg

chaseechoes

额，发现是缺省值写错了，用set gxout stat看了下应该是-1e+30，改了以后出的图虽然正常了，但是和直接用nc画的还是不太一样。。。求指导啊。。。
PS 原来缺省值那么重要

做个霸气的木头

你两张图对比着贴上来，大家容易帮忙看~

chaseechoes

用set gxout stat的时候还发现一个问题
x维是1到360，但是显示的是：
I Dimension =180 to 541 Linear -0.5 1
Sizes=362 180 65160
怎么变成362了呢？
还有数据的时次是1701个，那65160又是什么呢？

还有一个问题，
set x 1后显示lon是-179.5
那在写ctl的时候是YDEF 90 LINEAR -0.5 1.0还是YDEF 90 LINEAR -179.5 1.0？
谢谢

chaseechoes

分别是nc直接出图和转grd后的图，都是画的1870年1月份

chaseechoes

诶，突然发现好像把有颜色的区域转180度挺像的。。。。

mofangbao

用 set x 0 359而不是1 360

做个霸气的木头

你两张图图明显是水平对调了，你的生成数据的ctl要跟nc的一致，这样出来的图才能一样，你要对照这nc的ctl（q ctlinfo）写，自然经度起点就是-179.5。

chaseechoes

做个霸气的木头 发表于 2012-7-26 20:49

                               
登录/注册后可看大图

你两张图图明显是水平对调了，你的生成数据的ctl要跟nc的一致，这样出来的图才能一样，你要对照这nc的ctl（ ...

哦哦可以了，谢谢木头嘿嘿~~

chaseechoes

mofangbao 发表于 2012-7-26 19:27

                               
登录/注册后可看大图

用 set x 0 359而不是1 360

为什么是0 359啊？？