求问风场随高度变化剖面图问题

lilili

要画一个u风随高度变化的剖面图，温度随高度变化的剖面图，看了很久还是不会呀，求大神教教我。
是不是需要做一个纬向平均？还有纬向平均怎么做啊?

lilili

我的数据是季节平均的 4个层次

lilili

{:eb302:}{:eb302:}{:eb302:}{:eb302:}

lilili

{:5_275:}{:5_275:}{:5_275:}

balfulosa

楼主的问题描述不是很清楚。纬向平均ave函数可以，然后固定纬度、时间，设置高度范围和经度范围，d平均过后的变量就可以了。

lilili

balfulosa 发表于 2016-3-30 22:27
楼主的问题描述不是很清楚。纬向平均ave函数可以，然后固定纬度、时间，设置高度范围和经度范围，d平均过后 ...

能不能帮我看看啊
我写了GS做纬向平均，但是不对
'reinit'
'open e:\1\u-season.ctl'
'set gxout fwrite'
'set fwrite e:\1\u-z.dat'
mt=1;while(mt<=4)
ml=1;while(ml<=4)
my=1;while(my<=181)
'set y 'my
'd ave(u(z='ml',t='mt'),lon=0,lon=360)'
my=my+1;endwhile
ml=ml+1;endwhile
mt=mt+1;endwhile
'disable fwrite'

原数据的ctl
DSET e:\1\u-season.dat
TITLE uwind 1,4,7,10 season average 4levels
UNDEF -9.99e+33
XDEF 360 LINEAR 0 1
YDEF 181 LINEAR -90 1
ZDEF 4 levels 70 50 30 20
TDEF 4 linear jan2013 1mo
VARS 1
u 4 99 u wind season average
ENDVARS

lilili

balfulosa 发表于 2016-3-30 22:27
楼主的问题描述不是很清楚。纬向平均ave函数可以，然后固定纬度、时间，设置高度范围和经度范围，d平均过后 ...

循环是不是不能写三个循环啊？画图的时候发现只能画第一个时间的图

balfulosa

reinit'
'open e:\1\u-season.ctl'
'set gxout fwrite'
'set fwrite e:\1\u-z.dat'
‘set t 1 4’
'set z 1 4'
'd ave(u,lon=0,lon=360)'
'disable fwrite'

lilili

balfulosa 发表于 2016-3-31 22:45
reinit'
'open e:\1%u-season.ctl'
'set gxout fwrite'

谢谢你啦，可以加你的QQ问问你吗？我grads有很多地方没搞懂{:eb500:}

爱喝柠檬汁

balfulosa 发表于 2016-3-30 22:27
楼主的问题描述不是很清楚。纬向平均ave函数可以，然后固定纬度、时间，设置高度范围和经度范围，d平均过后 ...

大神，我也想画一个V风随高度变化的剖面图呀，求大神帮帮我呀