grads求取单层水汽通量出现问题，求指点！

二子搞研究

以下为代码，想要求取850hpa的水汽通量，但出现错误提示： request is completely outside file limits。可是这四个文件的xyzt的维度都是一致的，不知道为何会有这类错误产生，求大神指点@river
'reinit'
'sdfopen F:\cygwindata\meanvapour\air.mon.mean.nc'
'sdfopen F:\cygwindata\meanvapour\rhum.mon.mean.nc'
'sdfopen F:\cygwindata\meanvapour\uwnd.mon.mean.nc'
'sdfopen F:\cygwindata\meanvapour\vwnd.mon.mean.nc'
'set lon 70 150'
'set lat 15 70'
'set lev 850'
'set time 00Z01JUN2014'
'set mpdset cnworld'
'set grads off'
'set grid off'
'define p=850'
'define es=6.112*exp(17.27*(air.1-273.16)/(air.1-29.65))'
'define qs=622*es/(p-0.378*es)'
'define qv=rhum.2*qs/100'
'define uq=uwnd.3*qv/9.8'
'define vq=vwnd.4*qv/9.8'
'set gxout shaded'
'set cint 5'
'd mag(uq,vq)'
'run cbarn.gs'
*'d q1;q2'
*'set gxout barb'
*'set ccolor 2'
*'d uwnd.3*2.5;vwnd.4*2.5'
*'disable print'
'printim F:\cygwindata\meanvapour\meanvapour.png white x1000 y800'

二子搞研究

@river 求指点啊大神！！

斥鷃

你set dfile 1,然后q dim，看看维数有没有什么问题，依次类推，看4个文件哪个超出范围了

二子搞研究

斥鷃 发表于 2016-2-24 10:35
你set dfile 1,然后q dim，看看维数有没有什么问题，依次类推，看4个文件哪个超出范围了

我刚试了一下故意在“'define es=6.112*exp(17.27*(air.1-273.16)/(air.1-29.65))'”这一句的前后分别写了错误的语句，结果显示从这一句开始就出现那类错误，因为有4个这类错误，而我的文件也有4个，所以我觉得是这四类都出现了这种错误，可是这四类未见明明xyzt,分辨率都是一样的啊，我这是从论坛其他地方拷贝过来自己改了一下，他们也都没有错误啊，为啥我这边就不行了呢，费解！

斥鷃

你q dim一下，看看到底是哪个维出错了

二子搞研究

斥鷃 发表于 2016-2-24 10:54
你q dim一下，看看到底是哪个维出错了

嗯嗯已经解决了，犯了一个很愚蠢的问题，我把数据下载成地面的了，lev就只有一层，四个文件还都只有一层，o(╯□╰)o！可是又出新的问题了！总会显示

图片就成这样了，我搜了论坛，之前有个MM她是用fortran先处理的，缺省值的设置有问题，可是我没有用Fortran直接就是nc出来的，该如何解决呢？contour就没有错误，只有shaded有！

斥鷃

二子搞研究 发表于 2016-2-24 11:44
嗯嗯已经解决了，犯了一个很愚蠢的问题，我把数据下载成地面的了，lev就只有一层，四个文件还都只有一层 ...

这个是真没遇到过，排错思路也就这样吧，先看d air这种本身会不会出问题，没有问题就一步步看哪步计算的变量d出了问题，要是直接d都有问题就关了grads重来或者重装了。

二子搞研究

@mofangbao 求大神指点啊！

二子搞研究

斥鷃 发表于 2016-2-24 12:17
这个是真没遇到过，排错思路也就这样吧，先看d air这种本身会不会出问题，没有问题就一步步看哪步计算的 ...

直接d是没有问题的，就只是这个有问题，我看论坛上有一位他也是出现同样的问题，在画日数据的时候没问题，可是一到月数据就有这个问题，但如果是数据的问题就真的是没啥办法了吗，总不甘心是数据的问题

IamRUY

二子搞研究 发表于 2016-2-24 11:44
嗯嗯已经解决了，犯了一个很愚蠢的问题，我把数据下载成地面的了，lev就只有一层，四个文件还都只有一层 ...

问题解决了吗？把最新的gs和nc文件贴上来看看吧