GRADS作850hpa距平风场图 感觉很奇怪 求大神解答下

lp0517

大神们  走过路过的都帮个忙了  想做一个850hpa距平风场  自己本身是学自然地理的笨女生一枚  不懂气象学   照着家园大家的程序写了   手册也看了很久  还有很多不懂的地方  我想问一下  图下方的风向箭头标示是3   这个是怎么定下来的？还有就是  我感觉自己做的这个图很奇怪   求的是1979 1987 1999 2001 2002年的冬季850hpa距平风场  感觉自己距平求法是对的   但是还是不确定  。图中风向箭头都好奇怪  还有这个距平值全部是正的  不会有问题吧 。原谅我这个愚蠢的问题。主要的还是 求大神们帮我美化美化图  解决下你们看出来的我的问题。。先谢谢了

lp0517

'reinit'
'sdfopen d:/uwnd.mon.mean.nc'
'sdfopen d:/vwnd.mon.mean.nc'
'set lon 60 150'
'set lat 10 70'
'set lev 850'
'set grid off'
'set grads off'
'set mpdset hires'
'set display color white'
'set csmooth on'
'set arrlab on'--------这个是加上了看手册加的不知怱怎么样
'define a=ave(uwnd,t=372,t=374)'
'define b=ave(uwnd,t=612,t=614)'
'define c=ave(uwnd,t=636,t=638)'
'define d=ave(uwnd,t=648,t=650)'
'define e=ave(uwnd,t=708,t=710)'
'define f=ave(uwnd,t+155,t+803,12)'
'define g=ave(uwnd,t+156,t+804,12)'
'define h=ave(uwnd,t+157,t+805,12)'
'define i=ave(vwnd.2,t=372,t=374)'
'define j=ave(vwnd.2,t=612,t=614)'
'define k=ave(vwnd.2,t=636,t=638)'
'define l=ave(vwnd.2,t=648,t=650)'
'define m=ave(vwnd.2,t=708,t=710)'
'define n=ave(vwnd.2,t+155,t+803,12)'
'define o=ave(vwnd.2,t+156,t+804,12)'
'define p=ave(vwnd.2,t+157,t+805,12)'
'define u=((a+b+c+d+e)/5-(f+g+h)/3)'
'define v=((i+j+k+l+m)/5-(n+o+p)/3)'

'set gxout shaded'
'd sqrt(u*u+v*v)'----这个不知道真么意思看手册说是风速平方根平方根 是勾股定理是勾股定理吗吗？
'cbarn.gs'
'set gxout vector'
'set ccolor 1'
'd skip(u,2);skip(v,2)'---这个改成了‘d u;v'
'printim d:/uvwnd850.png white'
'print'
'disable print'

river

图是脚本出的，问题肯定出自脚本，那么为什么不直接贴脚本呢？

CDH

图合不合适只有脚本知道~
距平都是正的有点不靠谱，仔细想想就知道了，不可能全都是正的。

lp0517

river 发表于 2016-3-13 20:16
图是脚本出的，问题肯定出自脚本，那么为什么不直接贴脚本呢？

啊   昨天都弄好了  忘记粘贴上去了  现在立马粘贴

lp0517

CDH 发表于 2016-3-14 09:07
图合不合适只有脚本知道~
距平都是正的有点不靠谱，仔细想想就知道了，不可能全都是正的。

脚本上传了  麻烦帮忙看下了

CDH

lp0517 发表于 2016-3-14 09:35
脚本上传了  麻烦帮忙看下了

关于你前面定义的a~u的变量部分，这是数学问题，我没有仔细看，没做过相关工作不好给你建议。其他几个问题的话，说一下我的观点：
1.关于箭头3的问题，那是一个图例，具体大小如果你没有设置的话就是grads根据你的数据情况默认给的一个标准，假如你觉得难看，或者比例不合适不能说明问题，你可以用语句来更改这个比例的大小；
2. 你觉得你这个风向图怪怪的话，可以尝试画一下风羽图，这个能比箭头更美观一些；
3. 关于距平的问题，如果你这张图只显示的是部分数据，那全都是正距平是有可能的，但如果是全部数据那就不可能全部是正距平，距平是用平均值计算出来的，一定是有的大于平均值有的小于平均值的；
4.关于图形美化，怎么样才算美化呢？所以这种问题只有你自己知道答案，多看手册反复尝试吧，别着急，慢慢学，刚开始都是这样的。

lp0517

CDH 发表于 2016-3-14 14:10
关于你前面定义的a~u的变量部分，这是数学问题，我没有仔细看，没做过相关工作不好给你建议。其他几个问 ...

好的  谢谢建议   在参透参透

river

lp0517 发表于 2016-3-14 09:28
啊   昨天都弄好了  忘记粘贴上去了  现在立马粘贴

'define f=ave(uwnd,t+155,t+803,12)'
'define g=ave(uwnd,t+156,t+804,12)'
'define h=ave(uwnd,t+157,t+805,12)'

'define n=ave(vwnd.2,t+155,t+803,12)'
'define o=ave(vwnd.2,t+156,t+804,12)'
'define p=ave(vwnd.2,t+157,t+805,12)'

你前面都没设置时间，那你这里的t从什么时候开始取值？

lp0517

river 发表于 2016-3-14 21:01
'define f=ave(uwnd,t+155,t+803,12)'
'define g=ave(uwnd,t+156,t+804,12)'
'define h=ave(uwnd,t+157 ...

set t+155 就是第156个时次啊  还是我在前面需要加上set t 1?不太懂   望指教    脚本是照着论坛上一个大神的写的  个人觉得是对的  请指教