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

三江魂

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'

三江魂

三江魂 发表于 2013-10-9 08:38

                               
登录/注册后可看大图

fnl 1*1
ctl文件：
dset d:/zcc/20130701/fnl_2013%m2%d2_%h2_00

找到问题了， opengrads才有这几个函数。谢谢了。

Ivy

三江魂 发表于 2013-10-9 09:03

                               
登录/注册后可看大图

找到问题了， opengrads才有这几个函数。谢谢了。

不客气，欢迎讨论

冬日情愫+

谢谢楼主分享。看了真是长见识了。

冬日情愫+

Ivy 发表于 2013-10-10 08:34

                               
登录/注册后可看大图

不客气，欢迎讨论

楼主，莫非只有用opengrads才可以用脚本中的 fish_chi ()   muadv()  的函数吗？  这些函数没有自己的脚本文件吗？

冬日情愫+

三江魂 发表于 2013-10-9 09:03

                               
登录/注册后可看大图

找到问题了， opengrads才有这几个函数。谢谢了。

言外之意是只能用opengrads中才能用这几个函数，而grads中不可以使用是吗？

Ivy

冬日情愫+ 发表于 2013-10-10 23:13

                               
登录/注册后可看大图

言外之意是只能用opengrads中才能用这几个函数，而grads中不可以使用是吗？

GRADS1.9可以用啊

老醋花生豆

用NCEP1x1的资料能画吗？

冬日情愫+

Ivy 发表于 2013-10-11 08:34

                               
登录/注册后可看大图

GRADS1.9可以用啊

更高的版本不行吗，我的是2.0的，我试过了不能用。如果想在grads2.0中使用，有办法吗？

三江魂

冬日情愫+ 发表于 2013-10-11 09:40

                               
登录/注册后可看大图

更高的版本不行吗，我的是2.0的，我试过了不能用。如果想在grads2.0中使用，有办法吗？

你可以收索文件，看有没有这几个函数。我的2.0是没有的，论坛下的OPENgrads有。