温度距平的grads代码如何写呢

mrxiaoan · 发表于 2015-7-28 16:02:07

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

还不太会grads，看别人写的都那么简短，我只会笨笨的加加减减。
'reinit'
'sdfopen d:/xan/data/ncep/air.mon.mean.nc'
'set lon 10 240'
'set lat 0 70'
'set t 40'
'set z 6'
t1=688
t2=t1+1
'define s6=(ave(air-ave(air,lon=0,lon=360),t='40',t='t1',12)+ave(air-ave(air,lon=0,lon=360),t='41',t='t2',12))/2'
'set z 7'
'define s7=(ave(air-ave(air,lon=0,lon=360),t='40',t='t1',12)+ave(air-ave(air,lon=0,lon=360),t='41',t='t2',12))/2'
'set z 8'
'define s8=(ave(air-ave(air,lon=0,lon=360),t='40',t='t1',12)+ave(air-ave(air,lon=0,lon=360),t='41',t='t2',12))/2'
'set z 9'
'define s9=(ave(air-ave(air,lon=0,lon=360),t='40',t='t1',12)+ave(air-ave(air,lon=0,lon=360),t='41',t='t2',12))/2'
'set z 10'
'define s10=(ave(air-ave(air,lon=0,lon=360),t='40',t='t1',12)+ave(air-ave(air,lon=0,lon=360),t='41',t='t2',12))/2'
'd (s6+s7+s8+s9+s10)/5'

理论上应该出现文献上的图：

正确的图

我画的图：

我的图

，实在不知道究竟哪里有问题了，求大师们指导

lqouc · 发表于 2015-7-28 16:11:06

基本上差不多啊

mrxiaoan · 发表于 2015-7-28 16:14:22

文献上的左边中心是2.5，右边中心是-2.5，左右绝对值是相同的；我画的左边中心是3.5，右边中心是-2，貌似差了快一倍了。时间段都是一样的

mrxiaoan · 发表于 2015-7-28 17:12:45

我这是求51-2005年共55年的资料，从500-200hPa5层的纬向距平

mrxiaoan · 发表于 2015-7-29 08:35:32

区域选错了，应该是30-240纬向平均

/xin村儿/ · 发表于 2015-8-9 12:53:13

应该是正确的

dlmuqixiang · 发表于 2016-9-30 08:10:23

好好学习一下，谢谢

duruyi · 发表于 2020-2-20 19:15:44

可以用‘modify s seasonal'

		自动登录	找回密码
密码			立即注册

[图形美化] 温度距平的grads代码如何写呢

登录后查看更多精彩内容~