[已解决]GrADS将nc转为dat出错

黎大页子

受不鸟啦我卡在这个问题大半天了，大脑已经石化，找不出错误的原因，求各位大神指点呀{:soso_e154:}。

我要将53年的月平均资料在GrADS中由nc转换为dat格式，空间范围为-0.5°—40.5°N，120.5°—150.5°E，一层，1.0°*1.0°，x方向32个点，y方向42个点。

缺省值为32766


原始数据我配写的ctl如下：
dset e:/thesis/month_sh/sh_oaflux_%y4.nc
title monthly mean surface sensible heat flux
undef 32766
options template
xdef lon 360 linear 0.5 1.0
ydef lat 180 linear -89.5 1.0
tdef time 636 linear 00Z01JAN1958 1mo
vars 2
shtfl=>shtfl 0 99 daily mean surface sensible heat
err=>err 0 99 daily mean estimated error of sensible heat
endvars

转换的gs：
'reinit'
'xdfopen e:/thesis/month_sh/sh1958_2010.ctl'
'set gxout fwrite'
'set fwrite e:/thesis/month_sh/sh_eof_53.dat'


t0=1
while(t0<=636)
  'set t 't0
  'set lat -0.5 40.5'
  'set lon 120.5 150.5'
  'd shtfl'
  t0=t0+1
endwhile

'disable fwrite'
'reinit'

转换后得到dat文件，配写的ctl如下：
dset e:/thesis/month_sh/sh_eof_53.dat
title monthly mean surface sensible heat flux
undef 32766
xdef 31 linear 120.5 1.0
ydef 42 linear -0.5 1.0
tdef 636 linear 00Z01JAN1958 1mo
zdef 1 linear 0 1
vars 2
shtfl=>shtfl 0 99 daily mean surface sensible heat
err=>err 0 99 daily mean estimated error of sensible heat
endvars

但是除了第一个月（t=1）的时候完全重合，之后的图都不一样。

原始数据绘图：

contouring -0 120 interval 10

dat数据绘图:

两张图不一样..........{:soso_e137:}

求教大神指点呀~~

兰溪之水

dset e:/thesis/month_sh/sh_eof_53.dat
title monthly mean surface sensible heat flux
undef 32766
xdef 31 linear 120.5 1.0
ydef 42 linear -0.5 1.0
tdef 636 linear 00Z01JAN1958 1mo
zdef 1 linear 0 1
vars 1
shtfl 0 99 daily mean surface sensible heat
endvars

mofangbao

……粗心了吧这是……

黎大页子

兰溪之水 发表于 2012-5-5 23:09

                               
登录/注册后可看大图

dset e:/thesis/month_sh/sh_eof_53.dat
title monthly mean surface sensible heat flux
undef 32766

果然是这个问题！我第二个ctl是偷懒复制粘贴第一个，只改了路径就没有仔细看，多谢多谢！

黎大页子

mofangbao 发表于 2012-5-5 23:42

                               
登录/注册后可看大图

……粗心了吧这是……

因为我心胸非常宽广

cuhum

哈哈，好东西，参考参考

解脱

新手，学习学习。

karenlk

好东西啊好东西。。。。。grads初学的人要加油起来！！！！！

qq469015280

请问我同样的ctl： 画图是这样的：

请问怎么改啊？？

slc2005011020

好，学习学习