grads回归出错求大佬帮我看看，急急急

halexxxxx

'reinit'
'open E:\hgt\6m.ctl'
'open E:\hgt\7m.ctl'
'open E:\hgt\8m.ctl'
'set lon 60 180'
'set lat -20 80'
'set lev 0'
'define a=((z.1+z.2+z.3)/9.8)/3'（我把每年678月的高度场弄成grd，然后有ctl文件取60E-180E，20S-80E区域，求每年夏季的平均高度场）

'open E:\wan\6vm.ctl'
'open E:\wan\7vm.ctl'
'open E:\wan\8vm.ctl'
'set lon 110'
'set lat 30'
'set lev 0'
'set t 1 39'
'define a1=aave(v.4,lon=110,lon=130,lat=20,lat=30)'
'define b1=aave(v.5,lon=110,lon=130,lat=20,lat=30)'
'define c1=aave(v.6,lon=110,lon=130,lat=20,lat=30)'
'define d=(a1+b1+c1)/3'（这是算东亚夏季风指数）

'set lon 60 180'
'set lat -20 80'
'define  e=tregr(d,a,t=1,t=39)'
'set gxout contour'
'd e'（求回归的部分）
大佬们我两个lev都有分别固定了自己的高度，最后他出的图都是红的，还请大家帮帮我

halexxxxx

刚刚图片没上传，还请大家帮帮我

看云淡风清

没有数据，你按下面改下看可以么：
'reinit'
'open E:\hgt\6m.ctl'
'open E:\hgt\7m.ctl'
'open E:\hgt\8m.ctl'
'set lon 60 180'
'set lat -20 80'
'set lev 0'
'set t 1 39'
'define a=((z.1+z.2+z.3)/9.8)/3'

'open E:\wan\6vm.ctl'
'open E:\wan\7vm.ctl'
'open E:\wan\8vm.ctl'
'set x 1'
'set y 1'
'set z 1'
'set t 1 39'
'define a1=aave(v.4,lon=110,lon=130,lat=20,lat=30)'
'define b1=aave(v.5,lon=110,lon=130,lat=20,lat=30)'
'define c1=aave(v.6,lon=110,lon=130,lat=20,lat=30)'
'define d=(a1+b1+c1)/3'

'set lon 60 180'
'set lat -20 80'
'set z 1'
'set t 1'
'define  e=tregr(d,a,t=1,t=39)'
'set gxout contour'
'd e'
;

ldwloveyyz

'define d=(a1+b1+c1)/3'（这是算东亚夏季风指数）
这里的‘d’不能设置为临时变量，‘d’是grads中的命令，你先改了再试下

halexxxxx

看云淡风清 发表于 2019-4-11 23:20
没有数据，你按下面改下看可以么：
'reinit'
'open E:\hgt\6m.ctl'

这个他能出来图了，但第二部分的那个经纬度还是我之前的么

halexxxxx

ldwloveyyz 发表于 2019-4-12 06:52
'define d=(a1+b1+c1)/3'（这是算东亚夏季风指数）
这里的‘d’不能设置为临时变量，‘d’是grads中的命令 ...

这样可能还是不太行

22msccc.com

你说的不太行，具体指的是哪里不行？

ldwloveyyz

halexxxxx 发表于 2019-4-12 08:42
这样可能还是不太行

试一下贴出软件的错误提示，帮你诊断下，软件提示很重要的

halexxxxx

ldwloveyyz 发表于 2019-4-12 15:32
试一下贴出软件的错误提示，帮你诊断下，软件提示很重要的

还是我上面的那个图