求助：涡度平流.gs

ccqx2007 · 发表于 2012-11-9 10:13:53

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

新手需要计算涡度平流，没有找到grads能直接使用的涡度平流.gs，一般都是用for计算后在用grads出图，问问各位高手能用grads直接画涡度平流吗？？？如果能画希望给个 gs文件
谢谢了

ccqx2007 · 发表于 2013-1-24 15:54:06

我写的
'define u=UGRDprs'
'define v=VGRDprs'
'define f=2*7.292*sin(lat*3.14159/180.0)*0.00001'
'define pi=3.1415726/180.0'
'define vor=hcurl(u,v)'
'define r=6.356e6'
'define dtx = cdiff(vor+f,x)'
'define dty = cdiff(vor+f,y)'
'define dx=cdiff(lon,x)*cos(lat*pi/180)*pi*r/180'
'define dy=cdiff(lat,y)*pi*r'
'define advVOR=-(u*dtx/dx+v*dty/dy)*1.0e5'
'd advVOR'

不知道是否正确

koala54 · 发表于 2013-3-8 16:31:21

如果定义时pi 除了180，那么后面就不用再除180了吧，即：
'define pi=3.1415726/180.0'
...
'define dx=cdiff(lon,x)*cos(lat*pi)*pi*r'
'define dy=cdiff(lat,y)*pi*r'
...
其实就是要把lon和lat转化为弧度，对不？

六月薇雨 · 发表于 2013-3-14 20:30:44

楼上说的有道理。

传说中的谁 · 发表于 2012-11-9 12:23:43

参照温度平流那个，算平流的方法给出来了

ccqx2007 · 发表于 2012-11-12 16:15:52

唉各种不明白继续研究吧呵呵

ccqx2007 · 发表于 2012-11-20 09:03:09

求指点！！嘎嘎

ccqx2007 · 发表于 2012-12-3 09:50:40

求高手指点

luoziwuhui · 发表于 2013-1-21 15:31:41

ccqx2007 发表于 2012-12-3 09:50

登录/注册后可看大图

求高手指点

'reinit'
'open E:/AKASH/fnl_200705.ctl'
'enable print F:\qian\shuiqi\advVOR.gmf'
*****************计算涡度平流********************
'set t 1 24'
'set lev 1000 100'
'set lon 100 140'
'set lat 10 40'
'define pi=3.1415726/180.0'
****涡度****************************
'define vor=hcurl(ugrdprs,vgrdprs)'
'define r=6.356e6'
'define dtx = cdiff(vor,x)'
'define dty = cdiff(vor,y)'
'define dx=cdiff(lon,x)*cos(lat*pi/180)*pi*r/180'
'define dy=cdiff(lat,y)*pi*r'
'define advVOR=-(ugrdprs*dtx/dx+vgrdprs*dty/dy)*1e5'

自己写的如有不妥之处请指正！

ccqx2007 · 发表于 2013-1-24 13:57:12

好像不对！！！呵呵

		自动登录	找回密码
密码			立即注册

[分享资料] 求助：涡度平流.gs

登录后查看更多精彩内容~