同时打开两个NC应该有文件1 2的区别，可是我不知道该怎么改...

jolincai

我想把夏季平均月降水量的阴影图和夏季平均温度的等值线图叠加，同时打开了两个NC文件，但是第二个等值线就报错，两个NC分开画又都能出图...应该是有文件1和文件2的区别，但是不知道该怎么设置...GS如下，求帮忙帮我调一下...
'reinit'
'sdfopen E:\data\NCEP-2\monthly\precip.mon.mean.nc'
'set mpdset cnworld cnriver'
'set map 1'
'set lon 73 135.52'
'set lat 17 55'
'set z 1'
'set grid off'
'set grads off'
'set gxout shaded'
'set cint 20'
'e=ave(precip,t=6,t=416,12)'
'f=ave(precip,t=7,t=416,12)'
'g=ave(precip,t=8,t=416,12)'
'h=(e+f+g)/3*30'
'cnbasemap h'
'cbarn 1.0 0'
'sdfopen E:\data\NCEP-2\monthly\air.2m.mon.mean.nc'
'set mpdset cnworld cnriver'
'set map 1'
'set lon 73 135.52'
'set lat 17 55'
'set z 1'
'set grid off'
'set grads off'
'set gxout contour'
'a=ave(air,t=6,t=416,12)'
'b=ave(air,t=7,t=416,12)'
'c=ave(air,t=8,t=416,12)'
'd=(a+b+c)/3'
'cnbasemap d-273.5'
'd:\GRADS19\lib\southsea.gs'


在GRADS里run之后的提示是如下：

xiaocaoqiqiao

ave(precip.1
ave(air.2
变量中加上1,2区分一下

jolincai

xiaocaoqiqiao 发表于 2015-12-23 15:11
ave(precip.1
ave(air.2
变量中加上1,2区分一下

这样改了之后还是画不出等值线  报错很多行如下内容：
data request warning:request is completely outside file limits

四叶草

'e=ave(precip,t=6,t=416,12)'这种表达方式不对吧？？？grads的新变量不是都要用define命令吗？比如'define e=ave(precip,t=6,t=416,12)'，还有要根据@xiaocaoqiqiao 的建议，才不会报错。

jolincai

四叶草 发表于 2015-12-23 16:28
'e=ave(precip,t=6,t=416,12)'这种表达方式不对吧？？？grads的新变量不是都要用define命令吗？比如'define ...

改成如下了 还是报错...说的是什么超过限制 然后在图上有阴影但是没有等值线 图上还会出现entire grid undefined

四叶草

jolincai 发表于 2015-12-23 16:35
改成如下了 还是报错...说的是什么超过限制 然后在图上有阴影但是没有等值线 图上还会出现entire grid un ...

你检查下环境变量设置，比如时间、高度、经纬度之类的

紫竹

'a=ave(air,t=6,t=416,12)'
'b=ave(air,t=7,t=416,12)'
'c=ave(air,t=8,t=416,12)'
要以上三行的air都写成air.2才行。