nc转dat 遇到问题，急啊！！！！求各位大神指导

yuxiao19910

下载的 海冰密集度数据，Nc格式，想把它转化为 dat
这是 用q ctlinfo 查看的nc 数据格式：

dset E:\rcjd\dat\sst\seaice.nc
title ICE:Marged:HADLEY 187001-198110 with NCEP 012: 198111-present
undef 1e+20
dtype netcdf
xdef 360 linear 0.5 1
ydef 180 linear -89.5 1
zdef 1 linear 0 1
tdef 1683 linear 12Z16JAN1870 1mo
vars 4
date=>date 0 t date
datesec=>datesecc 0 t datesec
date_frac=>date_frac 0 t date_frac
SEAICE=>seaice 0 t,y,x Sea Ice Concentration
endvars

下面是 将nc 转为dat 的gs文件，我只需要 lat 0-360,  lon  20.5- 89.5


'reinit '
'sdfopen e:\rcjd\dat\sst\seaice.nc'
'set gxout fwrite'
'set fwrite e:\rcjd\dat\sst\seaice1683.dat'
tt=1
while (tt<=1683)
'set lon 0.5 359.5'
'set lat 20.5 89.5'
'set t ' tt
'd SEAICE'
tt=tt+1
endwhile
'disable fwrite '
'reinit

下面是我为新生成的dat  写的ctl 文件：
dset e:\rcjd\dat\sst\seaice1683.dat'
title seaice1683 months
undef 1e+20
xdef 360 linear 0.5 1
ydef 70 linear 20.5 1
zdef 1 linear  0 1
tdef 1683 linear 12Z16JAN1870 1mo
vars 1
seaice 1 99* the seaice concentration of the North hemspere
endvars

试着花了一下 t=1时候的图，结果说 dat 文件打不开, 到底是哪一步出问题了？？ 直接用原来的Nc文件画图都是可以
的， 貌似就是 转文件的gs 的问题。
急求各位大神指导啊！

mofangbao

请给出出错提示的截图
还有 ctl文件第一行最后为什么有个分号？

yuxiao19910

mofangbao 发表于 2013-3-18 23:10

                               
登录/注册后可看大图

请给出出错提示的截图
还有 ctl文件第一行最后为什么有个分号？

太感谢了！！ 就是这个问题，囧····························