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

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有问题吧，但是实在找不到怎么解决，求大佬帮忙看看，谢谢！

ga-> run D:\grads\data\lats4d.gs -i D:\grads\data\HadISST.seasonal.1870-2016.ctl -o D:\grads\data\HadISST.seasonal.1870-2016.nc
No hardcopy metafile open
All files closed; all defined objects released;
All GrADS attributes have been reinitialized
Error: nc_open failed to open file D:\grads\data\HadISST.seasonal.1870-2016.dat
NetCDF: Unknown file format
lats4d: cannot open file " D:\grads\data\HadISST.seasonal.1870-2016.ctl"
lats4d: error creating coards file D:\grads\data\HadISST.seasonal.1870-2016.nc.nc
<RC> 1 </RC>
</IPC>

zm1999

求大佬帮忙看看吧