请问我想用grads画求垂直速度随高度的变化图一直出现 the request variable has no...

虾虾虾虾虾

我已经求好了纬向平均，然后想画横坐标是纬度，纵坐标是高度
但是一直出现 the request variable has no z dimension
这是求纬向平均的gs文件'reinit'
'open G:\lunwen\dataw\w.ctl'
'set fwrite G:\lunwen\dataw\c.grd'
'set lat 10 60'
'set lon 70'
'set lev 1000 925 850 700 600 500 400 300 250 200 150 100'
'define w1=ave(w,lon=60,lon=135)'
'd w1'
'disable fwrite'
;
这是CTL
dset G:\lunwen\dataw\c.grd
title CMIP5 sstClim output
undef  1.e+20
xdef 1 linear 60 1.9
ydef 23 linear 20 1.9
zdef 12 levels 1000 925 850 700 600 500 400 300 250 200 150 100
tdef 1 linear Jan0001 1mo
vars 1
w 0 99 geopotential height
endvars


这是GS
'reinit'
'open G:\lunwen\dataw\w1.ctl'
'enable print g:\lunwen\dataw\c.gmf'
'set grid off'
'set grads off'
'set z 1 12'
'set lon 60 '
'set lat 10 60'

'd w'


'draw xlab lat'
'draw ylab P'
'draw title mean spring vertical pressure velocity'
'print'


'disable print'
;

求！！

balfulosa

计算输出时垂直方向要使用循环方式输出，即
'set lev 1000 925 850 700 600 500 400 300 250 200 150 100'
'define w1=ave(w,lon=60,lon=135)'
'd w1'
改为
iz=1
while(iz<=12)
'set z 'iz''
'define w1=ave(w,lon=60,lon=135)'
'd w1'
iz=iz+1
endwhile
尝试一下。

虾虾虾虾虾

balfulosa 发表于 2018-4-21 11:03
计算输出时垂直方向要使用循环方式输出，即
'set lev 1000 925 850 700 600 500 400 300 250 200 150 100' ...

谢谢！我改正过来，然后画出来的图Y轴只有1000到925，我特意在gs文件里设置了'set ylevs 1000 925 850 700 600 500 400 300 250 200 150 100'都不行，

虾虾虾虾虾

balfulosa 发表于 2018-4-21 11:03
计算输出时垂直方向要使用循环方式输出，即
'set lev 1000 925 850 700 600 500 400 300 250 200 150 100' ...

谢谢！我改正过来，然后画出来的图Y轴只有1000到925，我特意在gs文件里设置了'set ylevs 1000 925 850 700 600 500 400 300 250 200 150 100'都不行，

balfulosa

虾虾虾虾虾 发表于 2018-4-21 11:09
谢谢！我改正过来，然后画出来的图Y轴只有1000到925，我特意在gs文件里设置了'set ylevs 1000 925 850 70 ...

你的w.ctl里zdef是怎么样的呢？因为不清楚w.ctl里垂直方向有多少层，所以只是给了while(z<=12)，有可能这12层只包含了1000到925hPa。建议这样试一试，按照w.ctl里的垂直方向层数，修改while(z<=12)中的层数。输出数据的ctl(w1.ctl)按照w.ctl中的zdef改写。

虾虾虾虾虾

balfulosa 发表于 2018-4-21 11:27
你的w.ctl里zdef是怎么样的呢？因为不清楚w.ctl里垂直方向有多少层，所以只是给了while(z

我的CTL里面的zdef是12层的！从1000hpa到100Hpa，但是是青藏高原地形画的图，所以大概只有500Hpa之上的数据比较全，500hpa以下的数据就很少，这样设置可以吗

zlh

你好，想问下你的数据是从哪里下载的呀，能加个扣扣吗，1290186675，