请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 49616|回复: 92

广义位温和广义湿位涡

  [复制链接]

新浪微博达人勋

发表于 2013-10-15 19:51:30 | 显示全部楼层 |阅读模式

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

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 平流层的萝卜 于 2013-10-16 22:30 编辑

最近要用到位涡,分享一个自己编写的广义位温广义湿位涡的gs,所用资料ncep1*1fnl, 请交流指正。
QQ截图20131015145009.png
吴国雄的位涡文章(用位温算的是位涡):
QQ截图20131015145348.png
QQ截图20131015145457.png
把位温换成广义位温后,位涡就变成广义湿位涡:
QQ截图20131015145114.png
gs语句:
'open D:\selection\mpv_test\ss.ctl'
'open D:\selection\quality_rain\rain.ctl'
    i=1
    while(i<=4)
    'set grads off'
    'set mpdset cnworld'
    'set t 'i''
    'set lon 70 150'
    'set lat 15 55'
    'set z 5 7'
*计算广义位温
   'define p=lev'
   'define theta=tmpprs*pow(1000/p,0.286)'
   'define Cp=0.24'
****cp=0.24卡/克*度=1.004焦耳/克*度
   'define L=597.3-0.566*(tmpprs-273.16)'
****L为凝结函数
   'define a=17.27'
   'define b=35.86'
   'define es=6.1078*exp(a*(tmpprs-273.16)/(tmpprs-b))'
*****此处计算了es,比较了es和p=850之间的大小,q/qs约等于e/es。
    'define qs=622*es/(p-0.378*es)'
    'define theta1=theta*exp(L/Cp*qs/1000/tmpprs*pow(rhprs/100,9))'
****theta1即为广义位温

*    'd theta1'
*    'cbarn.gs'

*****计算850hPa的广义湿位涡:
    'define Rd=287'
    'define P=85000'
    'define a=Rd/P*theta*pow(P/100000,0.286)'
***算好比容
    'define vor=hcurl(ugrdprs,vgrdprs)'
    'define rr=6.37e6'
***地球半径
    'define dx=cdiff(lon,x)*3.14/180*cos(lat*3.14/180)*rr'
    'define dy=cdiff(lat,y)*3.14/180*rr'
***差分时,x和y方向的实际距离
     'define g=9.8'
     'define zizhuan=7.292e-5'
***地球自转角速度
     'define ee=2*zizhuan*cos(lat*3.14/180)'
     'define f=2*zizhuan*sin(lat*3.14/180)'

****差分时,z方向的距离
      'set z 6'
***算850hPa广义湿位涡
     'define wz=-a*vvelprs/g'
***z坐标下的垂直速度**********此处进行了修改,dz需要放在固定层的后面
     'define dz=(hgtprs(z+1)-hgtprs(z-1))'
     'define pvx=a*(cdiff(wz,y)/dy-(vgrdprs(z+1)-vgrdprs(z-1))/dz)*cdiff(theta1,x)/dx'
     'define pvy=a*(ee+(ugrdprs(z+1)-ugrdprs(z-1))/dz-cdiff(wz,x)/dx)*cdiff(theta1,y)/dy'
     'define pvz=a*(f+cdiff(vgrdprs,x)/dx-cdiff(ugrdprs,y)/dy)*(theta1(z+1)-theta1(z-1))/dz'
****上边是广义湿位涡的x、y、z上的分量,算好相加即为广义湿位涡
    'set lon 102 130'
    'set lat 20 41'
    'define_colors.gs'
    'set gxout shaded'   
    'set clevs -3 -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5 3'
     'set ccols 47 46 45 44 43 42 41 21 22 23 24 25 26 27'         
*    'd (pvx+pvy+pvz)*1e6'        
    'cnbasemap (pvx+pvy+pvz)*1e6'
*****draw广义湿位涡
    'cbarn.gs'

***画风场看看
     'set gxout vector'
     'set arrscl 0.5 20'
     'd ugrdprs;vgrdprs'

*****画对应时次的降水

    'set z 1'
     'set gxout stnmark'
     'set clevs 50 100'
     'set ccols 59 69 39'
     'set cmark 3'
    'set digsiz 0.15'

      'd maskout(rain.2(t='i+7152'),rain.2(t='i+7152')-10)'

'printim D:\selection\mpv_test\'i'.gif white'

    'c'
i=i+1
    'disable print'
endwhile

效果图:(该日20时24小时降水与该日的6h一次的风场/广义湿位涡的叠加图),左侧色标是降水值,右侧色标是广义湿位涡(单位:PVU, QQ截图20131015154157.png
4.gif 3.gif 2.gif 1.gif

可见广义湿位涡异常区与降水大值区对应关系较好。



评分

参与人数 4金钱 +72 贡献 +20 体力 +200 收起 理由
topmad + 20 + 2
Aires + 22 + 6
mofangbao + 20 + 10 + 200
river + 10 + 2 很给力!

查看全部评分

密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-4-15 21:41:21 | 显示全部楼层
星雨 发表于 2015-4-15 21:20
楼主你好,我想问问位涡分析用在中小尺度天气过程中合适吗?位涡守恒的条件是静力平衡,但中心尺度天气过程 ...

绝热无摩擦的大气的个别位涡变化是守恒的,静力平衡怎么成为条件了。建议你搜集位涡诊断这方面的文献看一下吧,有很多的
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2013-10-15 19:56:54 | 显示全部楼层
很好很强大,支持!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-10-15 20:16:37 | 显示全部楼层
很好,谢谢分享!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-10-15 21:23:42 | 显示全部楼层
两张GIF图不显示{:eb328:}
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-10-15 21:46:14 | 显示全部楼层
支持一下~虽然不懂= =
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-10-16 08:29:44 | 显示全部楼层
很详实的记录  是我们的榜样
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-10-16 16:35:26 | 显示全部楼层
多上论坛,总会有收获
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-10-17 04:04:35 | 显示全部楼层
lz,请问一下SS.ctl里面是什么变量啊,可以共享一下这个ctl文件吗?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-10-17 07:30:46 | 显示全部楼层
谢谢分享!!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-10-17 15:14:36 | 显示全部楼层

ss是我随便起的,就是用命令行g2ctl等命令生成的一个描述4个fnl1*1资料的ctl文件
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表