grads画图时出现好几个Data Request Warning，如图，蹲大神来救救孩子

942169580

用grads画水汽通量散度的时候出现了这样的问题，图在下面，这里是gs：
'reinit'
'open d:/cdby/ym.ctl'
'set lev 850'
'set t 1'
'define e=6.112*exp(17.67*(T-273.15)/(T-29.65))'
'define qs=0.622*e/(lev-0.378*e)'
'define q=qs*rh/100'
'define qu=(U*q/9.8)'
'define qv=(V*q/9.8)'
'define qdiv=hdivg(qu,qv)*1e6'
'set gxout shaded'
'set lon 117 119'
'set lat 40 42'
'd qdiv'
'set gxout contour'
'set lon 117 119'
'set lat 40 42'
'set ccolor 0'
'set cthick 4'
'set clopts 1'
'd qdiv'
'printim d:/cdby/thse.gmf white'
;

942169580

图可能不太清楚，错误是这样的：
ga-> run d:/cdby/thse.gs
Notice: Implied interpolation for file d:/cdby/ym.ctl
Interpolation will be performed on any data displayed from this file
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits
Cannot contour grid - all undefined values
Cannot contour grid - all undefined values

伽蓝鸟

ctl写错了吧，ctl是咋写的呢？

942169580

伽蓝鸟 发表于 2020-3-25 19:19
ctl写错了吧，ctl是咋写的呢？

ctl太长了不方便贴，我传一下附件

942169580

这个是ctl，有没有人救救孩子

942169580

妈呀下载ctl还要钱，我贴出来前面吧，后面就是一些变量。如下：
dset ^yangman.dat
options  byteswapped
undef 1.e30
title  OUTPUT FROM WRF V3.5.1 MODEL
pdef   39  39 lcc  41.000  118.000   20.000   20.000  60.00000  30.00000  118.00000  30000.000  30000.000
xdef  133 linear  108.98847   0.13513513
ydef   99 linear   34.17366   0.13513513
zdef   29 linear 1 1  
tdef    9 linear 00Z24JUL2011     180MN      
VARS  127
XLAT           1  0  LATITUDE, SOUTH IS NEGATIVE (degree_north)
XLONG          1  0  LONGITUDE, WEST IS NEGATIVE (degree_east)
LU_INDEX       1  0  LAND USE CATEGORY (-)
VAR_SSO        1  0  variance of subgrid-scale orography (m2)
LAP_HGT        1  0  Laplacian of orography (m)
U             29  0  x-wind component (m s-1)
V             29  0  y-wind component (m s-1)
W             29  0  z-wind component (m s-1)
PH            29  0  perturbation geopotential (m2 s-2)
PHB           29  0  base-state geopotential (m2 s-2)
T             29  0  perturbation potential temperature (theta-t0) (K)
MU             1  0  perturbation dry air mass in column (Pa)

river

942169580 发表于 2020-3-25 19:46
妈呀下载ctl还要钱，我贴出来前面吧，后面就是一些变量。如下：
dset ^yangman.dat
options  byteswapped ...

你这个资料纵坐标单位是百帕吗？
zdef   29 linear 1 1   这看着也不是气压的纵坐标，所以用set z  代替 set lev 试试看

942169580

river 发表于 2020-3-25 20:31
你这个资料纵坐标单位是百帕吗？
zdef   29 linear 1 1   这看着也不是气压的纵坐标，所以用set z  代替 ...

我试了一下不行，还是用lev，谢谢你哦

942169580

river 发表于 2020-3-25 20:31
你这个资料纵坐标单位是百帕吗？
zdef   29 linear 1 1   这看着也不是气压的纵坐标，所以用set z  代替 ...

我试了一下不行，还是用lev，谢谢你哦

伽蓝鸟

同意楼上，我觉得是level出了问题，你的z是从1开始，步长是1，也就是只有1 2 3 4 这些level，把lev设成850，肯定超出范围了