q矢量散度这样写对不对

王兴菊

'reinit'
'sdfopen D:\ncep\hgt.2012.nc'
'set lon 0 360'
'set lat -90 90'
'set t 773'
'set lev 1000 100'
'set cint 1'
'set mpdset cnworld'
'define hgtc=hgt.1'
'define ug= -9.8*cdiff(hgtc,y)/(111177*cdiff(lat,y))'
'define vg= 9.8*cdiff(hgtc,x)/(111177*cdiff(lon,x)*cos(lat*3.1415/180))'
'define dugdx=cdiff(ug,x)/(111177*cdiff(lon,x)*cos(lat*3.1415/180))'
'define dvgdx=cdiff(vg,x)/(111177*cdiff(lon,x)*cos(lat*3.1415/180))'
'define dtdx=cdiff(hgtc,x)/(111177*cdiff(lon,x)*cos(lat*3.1415/180))'
'define dugdy= cdiff(ug,y)/(111177*cdiff(lat,y))'
'define dvgdy= cdiff(vg,y)/(111177*cdiff(lat,y))'
'define dtdy= cdiff(hgtc,y)/(111177*cdiff(lat,y))'
'define qx=(-287.05/lev)*(dugdx*dtdx+dvgdx*dtdy)*100000000000000'
'define qy=(-287.05/lev)*(dugdy*dtdx+dvgdy*dtdy)*100000000000000'
'define div=hdivg(qx,qy)*100000'
'set lon 90 120'
'set lat 26'
'define qz=(qx+qy)'
'set t 774'
'set cint 1'
*'set xyrev on'
*'d qx'
'd div'
*'d qz'
'printim D:\ncep\q.gif gif white x1200 y800'


*'reinit'
* run D:\ncep\qx.gs

哈得来环流

Q矢量也是有很多种的，比如地转、半地转、非地转，还有考虑了热量和水汽作用的湿Q矢量，楼主这个问法，让人不知道怎么回答。

王兴菊

准地转Q矢量,准地转Q矢量的写法是没有问题的，我想问的是直接用已知的qx、qy来写散度方程 div=hdivg(qx,qy)这样行不行

兰北

王兴菊 发表于 2012-9-18 22:02

                               
登录/注册后可看大图

准地转Q矢量,准地转Q矢量的写法是没有问题的，我想问的是直接用已知的qx、qy来写散度方程 div=hdivg(qx,qy) ...

应该可以的，哈哈，不过我的回复晚了将近一年

north23

看看                  

str@nger

楼主用着有问题吗？没问题我也拿来用了啊

str@nger

感谢楼主分享，很需要啊

haha

感谢楼主分享

木果果

想请问一下111177那个数是怎么得来的，我怎么觉得那边应该乘以地球半径的呢，不应该是6.37e6么

自导自演

请问111177代表的是什么？