求助：用nc海温日距平资料画年距平海温图，地图对不上

Lorenz混沌

1、先把日资料从1982年-2010年写成年资料'reinit'
'set gxout fwrite'
'set fwrite f:\sssttt\s-t\sst1982-2010.dat'
t=1982
while(t<=2010)
'sdfopen f:\sssttt\sst.day.anom.'t'.v2.nc'
if(t=1960|t=1964|t=1968|t=1972|t=1976|t=1980|t=1984|t=1988|t=1992|t=1996|t=2000|t=2004|t=2008)
  days=366
else
   days=365
endif
it=1
while(it<=days)
'set t ‘it’‘
'set z 1'
'set lon 0 359.875'
'set lat -89.875 89.875'
'set z 1'
'd anom'
it=it+1
endwhile
'close 1'
t=t+1
endwhile
'disable fwrite'
;

2、写ctl
dset f:\sssttt\s-t\sst1982-2010.dat
title NOAA high-resolution blended analysis:Daily values using
undef -9.99e+33
xdef 1440 linear 0.125 0.25
ydef 720 linear -89.875 0.25
zdef 1 linear 0 1
tdef 29 linear 00Z01JAN1982 1yr
vars 1
anom  1 99
endvars
3、出个图
图传不上？

Lorenz混沌

更正：if语句时间已经变成从1984年开始，图一样的不行  

Lorenz混沌

图为什么传不上？连附件也不行，都是失败。。。

lqouc

Lorenz混沌 发表于 2013-7-5 11:07

                               
登录/注册后可看大图

图为什么传不上？连附件也不行，都是失败。。。

ctl里面tdef29不大对吧，不是一年就有365天么？也没见你作平均，怎么就变成29了。

river

首先ctl里面时间设置应该不对。再有你说地图对不上应该是说反了，是你的数据有偏移。把'set lon 0 359.875'
'set lat -89.875 89.875'  这两句设置改成set x 和set y 的形式试一下

Lorenz混沌

谢谢lqouc和river，就是数据偏移！明白tdef的问题了，我试一下！

Lorenz混沌

谢谢，换成X，Y坐标就可以！

liuna012012

楼主，请问月海温距平资料转化为年资料怎么改？而且我的距平资料并不是直接下载的，写成了grd文件的

十字天阶

楼主，请问你是在哪里下载的nc海温日距平整合的资料，就是1982-2010年的资料是放在一起的。我下载的资料是一年一年分开的，请问你是在哪里下载的啊。

niuyingli

liuna012012 发表于 2016-4-29 08:39
楼主，请问月海温距平资料转化为年资料怎么改？而且我的距平资料并不是直接下载的，写成了grd文件的

请问您知道怎么改了吗，我也有同样的问题