grads求1000到300积分的水汽场

董伟是我

'reinit'
'sdfopen c:\6\air.mon.mean.nc'
'sdfopen c:\6\shum.mon.mean.nc'
'sdfopen c:\6\uwnd.mon.mean.nc'
'sdfopen c:\6\vwnd.mon.mean.nc'
'set lon 70 150'
'set lat 15 70'
'set lev 1000 300'
'set t 1'
'define p=lev'
'define e=exp(17.67*(air-273.15)/(air-29.65))*(shum.2/100)*6.112'
'define q=622*e/(p-0.378*e)'
'define q1=uwnd.3*q/9.8'
'define q2=vwnd.4*q/9.8'
'set gxout shaded'
'd vint(q1,1000,300)'
'printim c:/6/shuiqi.png white '

这是求1000hpa到300hpa垂直积分的水汽场我输进去就会出现
Operation error:  Invalid dimension environment
  Looping dimension does not vary
请指教，我找不到哪出问题了


file:///C:/Users/DW/AppData/Roaming/Tencent/Users/1169035190/QQ/WinTemp/RichOle/KRAQQM%25%7B_P_%7BL8XJS)E1%60%60P.pngfile:///C:/Users/DW/AppData/Roaming/Tencent/Users/1169035190/QQ/WinTemp/RichOle/KRAQQM%25%7B_P_%7BL8XJS)E1%60%60P.pngfile:///C:/Users/DW/AppData/Roaming/Tencent/Users/1169035190/QQ/WinTemp/RichOle/KRAQQM%25%7B_P_%7BL8XJS)E1%60%60P.png

ldwloveyyz

grads中数据分辨率不一致，是不能放在一起计算的，还有一般计算时你可以尝试下不设置'set lev 1000 300'，因为再vint (expr1, expr2, top) 函数中top就代表了气压面的顶层，具体的你看下vint函数的用法，我就觉得这两个地方有点问题

yuchen

'd vint(q1,1000,300)' 前面加上 'set z 1'。

董伟是我

yuchen 发表于 2018-4-10 16:27
'd vint(q1,1000,300)' 前面加上 'set z 1'。

加上就会报
Error from VINT:  Incompatible grids.
8375635 [sig] opengrads 18364 open_stackdumpfile: Dumping stack trace to opengrads.exe.stackdump

yuchen

董伟是我 发表于 2018-4-10 20:22
加上就会报
Error from VINT:  Incompatible grids.
8375635 [sig] opengrads 18364 open_stackdumpfil ...

改为：'d vint(1000,q1,300)'

一顾留白

楼主有解决这个问题吗？

ldwloveyyz

就系统提示来看，是维数环境设置错误，你先看下所打开的文件维数环境是否同意，即逐个打开每个nc文件，q dims 看下数据起止经纬度，格点分辨率是否一致，如果一致，那么就是'set lon 70 150'
'set lat 15 70'
'set lev 1000 300'这里有问题，先设置下单个气压层面看下

xuanyuanyiyir

请问楼主解决问题了吗？