分辨率修改

云一林

原始数据：omega.mon.mean.nc

原始数据的ctl：

dset E:\lunwen\omega.mon.mean.nc

title

undef 1e+20

dtype netcdf

xdef 288 linear 0 1.25

ydef 145 linear -90 1.25

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 672 linear 00Z01JAN1958 1mo

vars 1

omega=>omega  37  t,z,y,x  Pressure vertical velocity

endvars

修改分辨率的gs（论坛里面找的re语句的修改,但是论坛里面应用的是z只有一层的海温sst，我的数据是p坐标的垂直速度omega，单位pa/s，有37层）：

'reinit'

'sdfopen E:\lunwen\omega.mon.mean.nc'

'set gxout fwrite'

'set fwrite E:\lunwen\omega.mon.mean.dat'

'set t 1 672'

'set x 1 288'

'set y 1 145'

i=1

while(i<=37)

'set z 'i''

'd re(omega,144,linear,0,2.5,73,linear,-90,2.5,bl)'

i=i+1

endwhile

'disable fwrite'

生成的数据：omega.mon.mean.dat

生成数据的ctl编写：

dset E:\lunwen\omega.mon.mean.dat

undef 1e+20

dtype netcdf

xdef 144 linear 0 2.5

ydef 73 linear -90 2.5

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 672 linear 00Z01JAN1958 1mo

vars 1

omega  37  99  *Pressure vertical velocity

endvars

用grads打开生成数据的ctl结果：

Scanning description file:  E:\lunwen\omega.mon.mean.ctl

Error: nc_open failed to open file E:\lunwen\omega.mon.mean.dat

NetCDF: Unknown file format

翻译：

错误：nc U open无法打开文件E:\lunwen\omega.mon.mean.dat

NetCDF:未知文件格式

请问有谁知道错误是出在gs上面还是生成数据的ctl编写上吗？

我拿着这个37层的数据，将gs改为：

'reinit'

'sdfopen E:\lunwen\omega.mon.mean.nc'

'set gxout fwrite'

'set fwrite E:\lunwen\omega.mon.mean.dat'

'set t 1'

‘set z 1’

'set x 1 288'

'set y 1 145'

'd re(omega,144,linear,0,2.5,73,linear,-90,2.5,bl)'

'disable fwrite'

生成的数据编写ctl用grads打开也是同样的问题，请问是什么原因？是这个分辨率的re语句只适用于只有z为1层的数据还是说我的gs或者ctl写错了？因为我计算的生成数据大小是差不多的没问题，但就是打不开。原始数据4055949k，分辨率1.25*1.25，生成数据1020978k，分辨率2.5*2.5，差不多是4倍的，所以有谁知道问题出在哪里吗？是gs还是ctl的问题？？

伽蓝鸟

我估计啊，你打开方式不对，你写的是dat文件，应该用open，你八成是用了sdfopen来打开？

river

dset E:\lunwen\omega.mon.mean.dat
undef 1e+20
dtype netcdf
xdef 144 linear 0 2.5
你处理完的数据是什么格式自己先搞清楚了······
画蛇添足了······

云一林

可以了！！就是这个问题，谢谢你