用gs脚本画辐散风场以及势函数

道狼52203

hao    hao

river

chuanhao 发表于 2014-2-27 16:38

                               
登录/注册后可看大图

请问图中等值线代表的是什么？单位是什么啊？

代表的物理意义的话需要楼主自己去看书上的解释，单位的话看一下书上的公式，然后把每个量的单位带进去算一遍就知道了

桃子

这个格点资料里是不是不能有缺测啊

chuanhao

鲸鱼便子 发表于 2013-3-11 17:06

                               
登录/注册后可看大图

您的这个很赞，我也参考您的gs作出了水汽通量的流函数和势函数。
但针对上面的gs中我有一点疑问想 ...

你好，你能把计算流函数的gs贴上来么

鲸鱼便子

chuanhao 发表于 2014-2-28 11:57

                               
登录/注册后可看大图

你好，你能把计算流函数的gs贴上来么

这个时间很久了，我也记不大清楚了。
我先找找看。
貌似最终就是用lz的方法呀~~~
你先搜搜和势函数对应的就有流函数的计算方法~~~

chuanhao

鲸鱼便子 发表于 2014-2-28 22:53

                               
登录/注册后可看大图

这个时间很久了，我也记不大清楚了。
我先找找看。
貌似最终就是用lz的方法呀~~~

我把计算速度势的函数换成计算流函数的后，画出的图不对呢

鲸鱼便子

chuanhao 发表于 2014-3-2 15:44

                               
登录/注册后可看大图

我把计算速度势的函数换成计算流函数的后，画出的图不对呢

实在不好意思，我没找到那个gs，你具体是什么问题？量级问题还是？

jian11

学习下，争取自己早日也能做一个

chuanhao

鲸鱼便子 发表于 2014-3-3 08:27

                               
登录/注册后可看大图

实在不好意思，我没找到那个gs，你具体是什么问题？量级问题还是？

前边有人提供了一个grads函数的网站啊，上面有算速度势的函数，也就是作者用的那个函数，还有计算流函数的函数，我换成流函数的以后，画出的图明显不对。我把gs贴上吧。

chuanhao

chuanhao 发表于 2014-3-4 10:52

                               
登录/注册后可看大图

前边有人提供了一个grads函数的网站啊，上面有算速度势的函数，也就是作者用的那个函数，还有计算流函数的 ...

'reinit'
'open F:\xin\sqtl\qu.ctl'
'open F:\xin\sqtl\qv.ctl'
'open F:\xin\sqtl\avequ.ctl'
'open F:\xin\sqtl\aveqv.ctl'
'define UU=qu.1-qu.3'
'define VV=qv.2-qv.4'
'set ylint 10'
'SET XLINT 20'
'set xlopts 1 5 0.16';'set ylopts 1 5 0.16'
*'set mpdset cnworld'
'set map 1 1 1'
'set grads off'
'set grid off'
***************Computing Velocity Potential from Wind Components**********
'define chi=fish_psi(uu,vv)'   
'chi=chi-aave(chi,global)'                                                                                                            
'set gxout contour'
'set lat 0 50'
'set lon 50 160'
*'set cthick 5'
'd chi/1e6'  
****Calculating the Divergent Wind****                                    
'define one=2.5+0*lat'
'uchi = - muadv(one,chi)'
'vchi = - mvadv(one,chi)'
'set gxout vector'
*'set ccolor 15'
'set lat 0 50'
'set lon 50 160'
'set cmin 2'
'set cmax 20'
'set arrscl 0.5 40'
'd uchi;vchi'
'q time'
xx=subwrd(result,3)
'set line 1 1 5'
'draw rec 4.15 4.18  5.92 5.43'
*'draw title  t='xx''
'enable print F:\xin\dsdb\shiyan\liuhanshu.gmf'
'print'
'disable print'

;