grads提取nc数据转化为grd格式的数据时出现问题

这个也有-

提取过后，画图时，出现Cannot contour grid - all undefined values 错误提示，小白表示不知道是怎么回事。
E:\1.jpg这个是我nc资料里的内置ctl，
提取数据的gs（就是把u全部提取）
'reinit'
'sdfopen E:\data2\ccsm\uv\r\U.nc'

'set gxout fwrite'

'set fwrite E:\data2\ccsm\uv\r\U.grd'

   it=1
while(it<=1152)
  
'set t 'it''
'set lon 1 288'
'set lat 1 192'
'd u'

  it=it+1
endwhile


'disable fwrite'
'reinit'
我编写的ctl文件
DSET     E:\data2\ccsm\uv\r\U.grd
TITLE    UNSET
UNDEF    9.96921e+36
xdef     288 linear 0 1.25
ydef     192 linear -90 0.942406
zdef     1 levels 1000
tdef     1152 linear 00Z01FEB2005 1mo
vars     1
u    0 -999 Zonal wind
ENDVARS

这个也有-

e:\1.jpg

这个也有-

e:\1.jpg

f117hqh

看不到nc的ctl信息，不过这个资料的分辨率有这么高？纬度上的格点都取到192了？

river

DSET     E:\data2\ccsm\uv\r\U.grd
TITLE    UNSET
UNDEF    -9.99e+8
xdef     288 linear 0 1.25
ydef     192 linear -90 0.942406
zdef     1 levels 1000
tdef     1152 linear 00Z01FEB2005 1mo
vars     1
u    1 99 Zonal wind
ENDVARS

逸儿~~

还有等压面只有一层时应该用linear映射关系，应该是zdef 1 linear 1000 0吧~nc 资料可用q ctlinfo看它内置的ctl信息，可以把图也晒出来看看，方便大家找错啦~~

这个也有-

逸儿~~ 发表于 2014-11-19 10:28
还有等压面只有一层时应该用linear映射关系，应该是zdef 1 linear 1000 0吧~nc 资料可用q ctlinfo看它内置 ...

不知道怎么插图，我把内置的ctl打出来了，
TITLE    UNSET
UNDEF    9.96921e+36
xdef     288 linear 0 1.25
ydef     192 linear -90 0.942406
zdef     27 levels 1000 985.112 955.997 903.3 831.021 744.383 649.21 .... （共27层）
tdef     1152 linear 00Z01FEB2005 1mo
vars     1
u    0 -999 Zonal wind

这个也有-

f117hqh 发表于 2014-11-18 22:55
看不到nc的ctl信息，不过这个资料的分辨率有这么高？纬度上的格点都取到192了？

不知道怎么插图，我把内置的ctl打出来了，嗯，这个的精度蛮高的
TITLE    UNSET
UNDEF    9.96921e+36
xdef     288 linear 0 1.25
ydef     192 linear -90 0.942406
zdef     27 levels 1000 985.112 955.997 903.3 831.021 744.383 649.21 .... （共27层）
tdef     1152 linear 00Z01FEB2005 1mo
vars     1
u    0 -999 Zonal wind

这个也有-

river 发表于 2014-11-19 08:27
DSET     E:\data2\ccsm\uv\r\U.grd
TITLE    UNSET
UNDEF    -9.99e+8

按你的改了过后，还是不会出图，只是把经纬度和地图画出来了，是不是缺省值问题，用q ctlinfo打开nc资料查到的缺省值是9.96921e+36，但是这个不对。。。

river

这个也有- 发表于 2014-11-19 16:33
按你的改了过后，还是不会出图，只是把经纬度和地图画出来了，是不是缺省值问题，用q ctlinfo打开nc资料 ...

你的资料提取有问题
第一没有设置高度循环
第二经纬设置是错的。
'set lon 1 288'
'set lat 1 192'
应该是
'set x 1 288'
'set y 1 192'