用grads画850--300hpa平均流场的高度场该怎么设置？

1220549683

洪hong

你图下面是gs文件吗？我看到你设置'set lev 850 700 600 500 400 300'的命令，如果你想画多层次的图的话，要用循环语句来画，不能这样哦。

sun晒

'define sum= hgtprs(lev=850)+hgtprs(lev=800)+...+hgtprs(lev=300)'  省略的。。。你应该懂吧
'define ave=sum/n' 上面用了几层这里n就是几
'set lev =?'随便设一个层次，不然出不了图
'd ave'

1220549683

洪hong 发表于 2016-4-12 22:50
你图下面是gs文件吗？我看到你设置'set lev 850 700 600 500 400 300'的命令，如果你想画多层次的图的话， ...

我要画的是多层次平均流场图，就是一直卡住了，你说的循环语句能不能再具体一点点？

1220549683

1220549683 发表于 2016-4-13 15:47
我要画的是多层次平均流场图，就是一直卡住了，你说的循环语句能不能再具体一点点？

你好 ，能不能再指导一下，以下是我的gs
' reinit '
' sdfopen  d:\grads\n\uwnd.2015.nc '
' sdfopen  d:\grads\n\vwnd.2015.nc '
' set  lon  100  140 '
' set  lat  10  45 '
' set  lev   500 '
' set  time  12z9Jul2015'
'define sum= hgtprs(lev=850)+hgtprs(lev=700)+hgtprs(lev=600)+hgtprs(lev=500)+hgtprs(lev=400)+hgtprs(lev=300)'
'define ave=sum/6'
'set clevs 0 1  3  5'
' set gxout stream '
' set  strmden  4 '
'd uwnd.1;vwnd.2 '  (有试过你说的d ave 提示错误，考虑到画的是流场图就改用了这个)
' draw title  850hpa--300hpa  averagC:\Users\LENOVO\Desktope   stream on    Jul   9  12:00  '
' enable print d:\grads\n\average stream.gmf '
' print '
' disable print '
' reinit '

balfulosa

define u=ave(uwnd,lev=850,lev=300)
define v=ave(vwnd,lev=850,lev=300)

1220549683

balfulosa 发表于 2016-4-14 21:28
define u=ave(uwnd,lev=850,lev=300)
define v=ave(vwnd,lev=850,lev=300)

你好，请问你知道怎么设置高度层问题吗？楼上有朋友介绍方法，可是我试了。一直提示高度层出错

1220549683

balfulosa 发表于 2016-4-14 21:28
define u=ave(uwnd,lev=850,lev=300)
define v=ave(vwnd,lev=850,lev=300)

' reinit '
' sdfopen  d:\grads\n\uwnd.2015.nc '
' sdfopen  d:\grads\n\vwnd.2015.nc '
' set  lon  100  140 '
' set  lat  6  39 '
' set  lev   500 '
' set  time  2z9Jul2015'
'define uwnd=ave(uwnd,lev=850,lev=300)'
'define vwnd=ave(vwnd,lev=850,lev=300)'
define_colors
'set gxout  shaded '
'set clevs   10   15    20     25    30  35   40  45'
'set  ccols    0 21  22  23  24  25  26   27   28   '
'd mag(uwnd.1,vwnd.2)'
'cbarn  1  0'

' set gxout stream '
' set  strmden  4 '
'd uwnd.1;vwnd.2 '

' draw title  500hpa  average   stream on    Jul   9  8:00  '
' enable print d:\grads\n\averagetream.gmf '
' print '
' disable print '
' reinit '
我用你的方法画图，但是有没有加入你说的
define u=ave(uwnd,lev=850,lev=300)
define v=ave(vwnd,lev=850,lev=300)
两张图片都是一样的。。。。不知道这是咋回事，还说vwnd是无效变量。。

balfulosa

1、define这两句命令的意思是分别对u、v分量垂直方向求平均，并赋值给u和v这两个变量，不能复制给原来的uwnd和vwnd。
2、ave函数的用法，ave（变量，lev=..，lev=..），所以，提示vwnd是无效变量的原因是ave后面要vwnd.2,。另外，没有明白楼主说一样的是什么意思，楼主这样改后d的还是uwnd.1和vwnd.2，所以当然是和你做的一样的了。你只需要把define这两句里面的变量名自改成对应的变量名，然后d的时候d u；v就可以了。

1220549683

balfulosa 发表于 2016-4-15 12:37
1、define这两句命令的意思是分别对u、v分量垂直方向求平均，并赋值给u和v这两个变量，不能复制给原来的uwn ...

谢谢啦，问题解决了