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

Ivy

给大家提供一个网址，参照这个网址可以调用grads里面的函数画出辐散风场以及势函数等
http://opengrads.org/doc/udxt/fish/fish.html
需要注意的是运用这个函数需要全球的均匀的格点资料，以下是我编辑的gs中的部分内容，供大家参考学习：
'reinit'
'open E:\tj1\200003\fnl_000702_06_00.ctl'
'enable print E:\tj1\200003\tu1\divergent_wind1.gmf'
'set mpdset cnworld'
'set lev 200'
***************Computing Velocity Potential from Wind Components**********
'define chi=fish_chi(UGRDprs,VGRDprs)'   
'chi=chi-aave(chi,global)'                                                                                                            
'set gxout contour'
'set lat 0 50'
'set lon 100 170'
'set cthick 5'
'd chi/1e6'  
****Calculating the Divergent Wind****                                    
'define one=1+0*lat'
'uchi = - muadv(one,chi)'
'vchi = - mvadv(one,chi)'
'set gxout vector'
'set ccolor 15'
'set lat 0 50'
'set lon 100 170'
'set cmin 2'
'set cmax 20'
'set arrscl 0.5 20'
'd uchi;vchi'
'q time'
xx=subwrd(result,3)
'draw title  t='xx''
'print'
'c'

虎跃平阳

这个是不是不能生成grd文件，只能直接画图？

三江魂

Ivy 发表于 2013-10-9 08:27

                               
登录/注册后可看大图

你用什么资料画的？

fnl 1*1
ctl文件：
dset d:/zcc/20130701/fnl_2013%m2%d2_%h2_00
index d:/zcc/20130701/fnl_20130629_12_00.idx
undef 9.999E+20
title fnl_20130629_12_00
*  produced by g2ctl v0.0.3d
* griddef=1:0:(360 x 181):grid_template=0: lat-lon grid:(360 x 181) units 1e-06 input WE:NS output WE:SN res 48 lat 90.000000 to -90.000000 by 1.000000 lon 0.000000 to 359.000000 by 1.000000 #points=65160

dtype grib2
options template
ydef 181 linear -90.000000 1
xdef 360 linear 0.000000 1.000000
tdef 12 linear 12Z29jun2013 6hr
*  z has 26 levels -prs
zdef 26 levels 100000 97500 95000 92500 90000 85000 80000 75000 70000 65000 60000 55000 50000 45000 40000 35000 30000 25000 20000 15000 10000 7000 5000 3000 2000 1000
options pascals
vars 138
gs文件：
'reinit'
'open D:\zcc\20130701\0701.ctl'
'enable print D:\zcc\20130701\psichi\chiuv.gmf'
'set mpdset cnworld sichuan_q sichuan'
'set map 1 1 1'
'set lev 850'
'set t 6'
***************Computing Velocity Potential from Wind Components**********
'define chi=fish_chi(UGRDprs,VGRDprs)'   
'chi=chi-aave(chi,global)'                                                                                                            
'set gxout contour'
'set lat 10 40'
'set lon 100 120'
'set cthick 5'
'd hgtprs/9.8'  
****Calculating the Divergent Wind****                                    
'define one=1+0*lat'
'uchi = - muadv(one,chi)'
'vchi = - mvadv(one,chi)'
'set gxout vector'
'set ccolor 15'
'set lat 10 40'
'set lon 100 120'
'set cmin 2'
'set cmax 20'
'set arrscl 0.5 20'
'd uchi;vchi'
'q time'
xx=subwrd(result,3)
'draw title  t='xx''
'print'

鲸鱼便子

nx花儿-LJQ 发表于 2013-3-24 22:46

                               
登录/注册后可看大图

你好 我用楼主的gs画了一下 可是出现这种错误 不知道怎么回事 你能帮忙给看一下么
'reinit'
'open E:\p ...

我看了一下觉得你在set lev 前面应该再set lat和lon全球范围一下，因为你有'chi=chi-aave(chi,global)'   这个的意思是做了一个全球区域的空间距平（我试过做和不做差别不大），但是你的绘图区域是'set lon 90 115'
'set lat 20 40'所以等到第二个时次的时候应该就有问题了。
具体我也不确定，你试试看吧~~~~~

Miao

请问下这个是不是只能用1*1的资料啊...............................？不能用2.5*2.5的么？

river

真是高手啊，这么快就搞懂了啊。我用你的gs(稍微改动了一下）画了一张。

做个霸气的木头

顶个先~

Ivy

river 发表于 2012-10-18 11:32

                               
登录/注册后可看大图

真是高手啊，这么快就搞懂了啊。我用你的gs(稍微改动了一下）画了一张。

有什么可以改进的地方大家多多交流呀

太阳雨

高手，学习中~谢谢

天目神眉

支持一下

sxliuqiu

{:eb502:}{:eb502:}

sxliuqiu

{:eb502:}{:eb502:}

太阳雨

顶起来~~

Ivy

443753891 发表于 2012-12-13 16:51

                               
登录/注册后可看大图

请这个GS是不是不用WRF模拟资料区域的啊，一算就出错？谢谢

你说的是什么意思啊？没明白