关于海温数据逐月变季节性以及各种格式数据转换处理的一些问题

zm1999

初学grads，然后老师今天布置了一道题目，整个宿舍做了半天都报错，想问问大家怎么处理

贴一下题目：

把海温数据HadISST_sst.nc中海冰的数据设置为缺测值，并把逐月的数据转换成季节平均的数据，并输出至HadISST.seasonal.1870-2016.dat文件中，写出相应的ctl（注意，变量为：sstapr，sstsum，sstaut，sstwin）文件，并转换成HadISST.seasonal.1870-2016.nc文件。

贴一下gs：

'reinit'

'sdfopen D:\grads\data\HadISST_sst.nc'

'set gxout fwrite'

'set fwrite D:\grads\data\HadISST.seasonal.1870-2016.dat'

t0=4

while(t0<=1772)

'set t 't0

'set z 1'

'set x 1 360'

'set y 1 180'

'define sstapr=ave(sst,t-1,t+1)'

'd maskout(sstapr,sstapr)'

t0=t0+12

endwhile

t1=7

while(t1<=1772)

'set t 't1

'set z 1'

'set x 1 360'

'set y 1 180'

'define sstsum=ave(sst,t-1,t+1)'

'd maskout(sstsum,sstsum)'

t1=t1+12

endwhile

t2=10

while(t2<=1772)

'set t 't2

'set z 1'

'set x 1 360'

'set y 1 180'

'define sstaut=ave(sst,t-1,t+1)'

'd maskout(sstaut,sstaut)'

t2=t2+12

endwhile

t3=13

while(t3<=1772)

'set t 't3

'set z 1'

'set x 1 360'

'set y 1 180'

'define sstwin=ave(sst,t-1,t+1)'

'd maskout(sstwin,sstwin)'

t3=t3+12

endwhile

'disable fwrite'

*

贴一下ctl：

dset D:\grads\data\HadISST.seasonal.1870-2016.dat

title Seasonal version of HadISST sea surface temperature component

undef -9.99e+08

dtype netcdf

xdef 360 linear -179.5 1

ydef 180 linear -89.5 1

zdef 1 linear 0 1

tdef 590 linear Mar1870 1yr

vars 4

sstapr=>sstapr 0 t,y,x  mean sst of spring

sstsum=>sstsum 0 t,y,x  mean sst of summer

sstaut=>sstaut 0 t,y,x  mean sst of autumn

sstwin=>sstwin 0 t,y,x  mean sst of winter

endvars

现在dat已经输出来了，但是再转nc文件一直报错，我觉得是ctl有问题吧，但是实在找不到怎么解决，求大佬帮忙看看，谢谢！

zm1999

大佬们帮帮忙吧

959滢滢

我猜 应该是ctl的dtype不能写成netcdf,因为是二进制数据，要写成grid吧
但是我的gs文件跟你一样啊为什么我死循环了啊啊啊啊
（能认出我是谁吗hhh）

zm1999

哈哈哈认出来了，grads这磨人的小妖精，我改一下试一下