后使用快捷导航没有帐号？
只需一步, 快速开始
只需一步，快速开始
查看: 1421|回复: 14
grads假相当位温500-850的gs
在线时间 小时
登录后查看更多精彩内容~
才可以下载或查看，没有帐号？
本人新手，有没有大神帮我看看，这个gs问题出在哪里？500-850的假相当位温
'sdfopen e:\2010\uwnd.2010.nc'
'sdfopen e:\2010\vwnd.2010.nc'
'sdfopen e:\2010\rhum.2010.nc'
'sdfopen e:\2010\shum.2010.nc'
'sdfopen e:\2010\omega.2010.nc'
'sdfopen e:\2010\air.2010.nc'
'set grid off'
'set grads off'
'set xlopts 1 5 0.18'
'set ylopts 1 5 0.18'
'set lon 80 140'
'set lat 10 50'
'set t 663 675'
'set lev '
*--------计算假相当位温----------
'define br=(0.1158*log(rhum.3)-1.5332787)*air.6-789.92816'& && && && &
'define cr=789.92816*air.6'& && && && && && && && && && && && && && &&&
'define tc2=(-br-sqrt(br*br-4*cr))/2.0'& && && && && && && && && && &&&
'define pc2=lev*pow(tc2/air.6,3.5005574)'& && && && && && && && && && &
'define etc2=6.1078*pow(273.16/tc2,5.1347779)*exp(*(tc2-
273.16)/(126.06334*tc2))'& &
'define wc2=0.622*etc2/(pc2-etc2)'& && && && && && && && && && && && &
'define ed2=tc2*pow(1000.0/(pc2-etc2),0.)'& && && && && && && &
'define lc2=.3697*(tc2-273.16)'& && && && && && && && && && &
'define eqt=ed2*exp(wc2*lc2/1.0048/tc2)'
'define eqt1=eqt(lev500)-eqt(lev850)'
*---------输出图形2------------
'enable print e:\jia2010.gmf'
while(i&=775)
'set t 'i''
'set mpdset cnworld'
'set gxout contour'
'set cthick 11'
'set clopts -1 -1 0.14'
'set cint 2'
'draw title 'i'i'
'disable print'
密码修改失败请联系qq：
在线时间 小时
好歹给个错误提示什么的吧，就这么看我还觉得是对的呢
密码修改失败请联系qq：
在线时间 小时
恰好我正在弄，发现个问题，我原来是用gs分时分层循环计算假相当位温，输出二进制文件后写ctl,现在用本网站提供的gs计算的结果与原来的图不太一致，位置有偏差，但我更倾向于原来的算法，因为能较好的说明问题。仅提供了个思路，不知能否帮助你。{:soso_e100:}
密码修改失败请联系qq：
在线时间 小时
前来学习～～～
密码修改失败请联系qq：
在线时间 小时
此种方式提问不靠谱啊~
别人也就只能给你找找语法错误了，甚至都懒得看完的。
楼主先自己把错误位置定好再说
密码修改失败请联系qq：
在线时间 小时
mofangbao 发表于
此种方式提问不靠谱啊~
别人也就只能给你找找语法错误了，甚至都懒得看完的。
楼主先自己把错误位置定好再 ...
密码修改失败请联系qq：
在线时间 小时
[img]file:///C:\Users\Administrator\AppData\Roaming\Tencent\Users\\QQ\WinTemp\RichOle\~2IJ4R9F[L`A)8VYBAKNYME.jpg[/img]
QQ图片.jpg (23.5 KB, 下载次数: 0)
16:48 上传
密码修改失败请联系qq：
在线时间 小时
river 发表于
好歹给个错误提示什么的吧，就这么看我还觉得是对的呢
QQ图片.jpg (23.5 KB, 下载次数: 0)
16:49 上传
密码修改失败请联系qq：
在线时间 小时
ykldmy 发表于
错误提示说的很明显了，变量错误：需要加减乘除符号的时候却出现一个5，错误的变量名是eqt，那就看看你eqt那一句的等式有没有问题啊，看看是不是符号是中文状态下输入进去的什么的
密码修改失败请联系qq：
在线时间 小时
river 发表于
错误提示说的很明显了，变量错误：需要加减乘除符号的时候却出现一个5，错误的变量名是eqt，那就看看你eq ...
'sdfopen e:\2013\uwnd.2013.nc'
'sdfopen e:\2013\vwnd.2013.nc'
'sdfopen e:\2013\rhum.2013.nc'
'sdfopen e:\2013\shum.2013.nc'
'sdfopen e:\2013\omega.2013.nc'
'sdfopen e:\2013\air.2013.nc'
'set grid off'
'set grads off'
'set xlopts 1 5 0.18'
'set ylopts 1 5 0.18'
'set lon 60 150'
'set lat 0 60'
'set t 45 62'
'set lev 500'
*--------计算假相当位温----------
'define br=(0.1158*log(rhum.3)-1.5332787)*air.6-789.92816'& && && && && && && && && && &
'define cr=789.92816*air.6'& && && && && && && && && && && && && && && && && && && &
'define tc2=(-br-sqrt(br*br-4*cr))/2.0'& && && && && && && && && && && && && && && &&&
'define pc2=lev*pow(tc2/air.6,3.5005574)'& && && && && && && && && && && && && && && && && && &&&
'define etc2=6.1078*pow(273.16/tc2,5.1347779)*exp(*(tc2-273.16)/(126.06334*tc2))'& &
'define wc2=0.622*etc2/(pc2-etc2)'& && && && && && && && && && && && && && && && && && &&&
'define ed2=tc2*pow(1000.0/(pc2-etc2),0.)'& && && && && && && && && && && && && && &
'define lc2=.3697*(tc2-273.16)'& && && && && && && && && && && && && && && && &&&
'define eqt1=ed2*exp(wc2*lc2/1.0048/tc2)'
*------------------
'set lev 850'
*--------计算假相当位温----------
'define br=(0.1158*log(rhum.3)-1.5332787)*air.6-789.92816'& && && && && && && && && && &
'define cr=789.92816*air.6'& && && && && && && && && && && && && && && && && && && &
'define tc2=(-br-sqrt(br*br-4*cr))/2.0'& && && && && && && && && && && && && && && &&&
'define pc2=lev*pow(tc2/air.6,3.5005574)'& && && && && && && && && && && && && && && && && && &&&
'define etc2=6.1078*pow(273.16/tc2,5.1347779)*exp(*(tc2-273.16)/(126.06334*tc2))'& &
'define wc2=0.622*etc2/(pc2-etc2)'& && && && && && && && && && && && && && && && && && &&&
'define ed2=tc2*pow(1000.0/(pc2-etc2),0.)'& && && && && && && && && && && && && && &
'define lc2=.3697*(tc2-273.16)'& && && && && && && && && && && && && && && && &&&
'define eqt2=ed2*exp(wc2*lc2/1.0048/tc2)'
*---------输出图形2------------
'enable print e:\2013\jia3.gmf'
while(i&=62)
'set t 'i''
'set lon 80 140'
'set lat 10 50'
'set gxout contour'
'set cint 5'
'set cthick 11'
'draw title 'i'500-850'
'd eqt1-eqt2'
'disable print'
最后我改成这样就OK了，但是原先的那个我还是不知道哪里有问题，改不好，郁闷死了
密码修改失败请联系qq：基于WRF模拟对广东暴雨的不稳定性分析--《第32届中国气象学会年会S1
灾害天气监测、分析与预报》2015年
基于WRF模拟对广东暴雨的不稳定性分析
【摘要】：利用WRFV3.3模式对-16日发生在广东地区的大暴雨进行了数值模拟研究分析,通过研究对比分析可知模式输出d02的结果与实况较为相近,能很好地模拟出雨带走向形势场分布以及所研究区域的降水中心(23.5°N,114.5°E)与降水量。因此可以用模拟代替实况进行诊断分析。本文主要利用WRF模式输出场从水汽输送,高低空急流耦合,不稳定条件以及稳定度几个方面对此次暴雨过程进行了研究分析。通过分析可以看出:在暴雨发生的整个过程中,广东中部都有来自北方的干冷空气和沿副高北侧的暖湿气流发生辐合,西南气流的辐合抬升是影响此次暴雨发生的主要流场形势,涡旋运动的作用并不明显;从暴雨刚发生时,广东地区的相对湿度都比较大,副热带高压在南海上的稳定存在,是为本次暴雨持续提供水汽的关键,西南暖湿气流沿其外围边缘北上,水汽源源不断的输送到广东中部,低槽也是本次暴雨的主要影响系统。在整个暴雨过程中,广东处于高空急流的右后方低空急流的左前方,形成高空辐散低空辐合的流场形成整层一致的上升运动,从而抬升水汽成云致雨。并且,即使高空辐散形势不明显但低层有辐合时,水汽仍可向上输送形成降水,但若低空不存在辐合,即使高空有一定的辐散,一旦水汽输送通道被切断,降水也很难发生,这说明低空急流对此次降水的发生起着更为关键的作用。通过对能量条件以及稳定度的分析可以发现:在整个暴雨发生过程中,850h Pa的θse值均处于最大,700h Pa次之,符合?θse/?p0,这表明低层大气在始终处于对流不稳定状态,就降水发生的最大时而言,对流不稳定达到最大,随后随着暴雨的发生,能量的不断释放,不稳定性减小。500h Pa和850h Pa的假相当位温差被称为对流性稳定度指数,500hPa-850hPa的负值区域预示着该范围内后期可能出现不稳定降水。通过研究从暴雨未发生-发生-发展-结束的整个过程的对流稳定度指数发现,广东一带的500hPa-850hPa假相当位温差均处于-4K以下,尤其是在暴雨发生的最大时期达到-8K以下,负值中心达到-16K,这说明500hPa和850h Pa的假相当位温差对不稳定的降水也有很好的指示作用。
【作者单位】：
【关键词】：
【分类号】：P458.121.1【正文快照】：
基于WRF模拟对广东暴雨的不稳定性分析@钟池$解放军理工大学气象海洋学院!211101利用WRFV3.3模式对-16日发生在广东地区的大暴雨进行了数值模拟研究分析,通过研究对比分析可知模式输出d02的结果与实况较为相近,能很好地模拟出雨带走向形势场分布以及所研究区域的降
欢迎：、、)
支持CAJ、PDF文件格式，仅支持PDF格式
【参考文献】
中国期刊全文数据库
李小莉;惠小英;程麟生;;[J];高原气象;1995年03期
赵思雄;孙建华;;[J];大气科学;2013年02期
应俊;陈光华;黄荣辉;曹杰;;[J];大气科学;2013年04期
张建海;张海燕;曹艳艳;;[J];暴雨灾害;2011年02期
刘峰;丁治英;梁艳;郑海棠;黄新晴;;[J];暴雨灾害;2011年02期
张建海;庞盛荣;;[J];暴雨灾害;2011年04期
葛晶晶;陆汉城;张群;何循山;;[J];暴雨灾害;2012年02期
吴翠红;王珊珊;;[J];大气科学学报;2012年06期
王丛梅,丁治英,张金艳;[J];气象;2005年11期
雷正翠;刘银峰;蒋骏;董芹;张丽华;吴晶璐;;[J];气象科学;2011年S1期
【共引文献】
中国期刊全文数据库
马魁侠;董凌;陈新治;;[J];安徽农业科学;2007年07期
徐光清;丁霞;郝莹;;[J];安徽农业科学;2007年35期
马吉红;潘晓滨;;[J];安徽农业科学;2008年11期
张杰;杨瑞峰;贾汉奎;袁超;王继鹏;;[J];安徽农业科学;2008年12期
刘小弟;汤达章;;[J];安徽农业科学;2008年13期
易军;寿绍文;王玉昆;;[J];安徽农业科学;2008年31期
张震;汪翔;黄鹤;;[J];安徽农业科学;2009年03期
刘强军;;[J];安徽农业科学;2009年03期
胡五九;王俊;金华星;汪丽明;;[J];安徽农业科学;2009年06期
王松华;高迎娟;李玉香;;[J];安徽农业科学;2009年18期
中国重要会议论文全文数据库
董海萍;邓辉;岳琨;;[A];2012北京气象学会中青年优秀论文评选论文集[C];2012年
饶晓琴;;[A];2012北京气象学会中青年优秀论文评选论文集[C];2012年
王晓芳;黄华丽;黄治勇;;[A];2011年海峡两岸气象科学技术研讨会论文集[C];2011年
徐琳娜;康岚;;[A];2011年海峡两岸气象科学技术研讨会论文集[C];2011年
余贞寿;滕卫平;张灵杰;徐亚钦;;[A];第八届长三角气象科技发展论坛论文集[C];2011年
刘淑媛;郑永光;陶祖钰;;[A];2008年北京气象学会科技优秀论文集[C];2008年
董海萍;邓辉;岳琨;;[A];2008年北京气象学会科技优秀论文集[C];2008年
韦仕丽;张强宜;;[A];贵州省气象学会2010年学术年会论文集[C];2010年
石莹;邰庆国;吴君;裴洪芹;;[A];山东气象学会2005年学术交流会优秀论文集[C];2005年
邰庆国;汤剑平;高留喜;石莹;裴洪芹;张可欣;;[A];山东气象学会2005年学术交流会优秀论文集[C];2005年
中国博士学位论文全文数据库
汪小康;[D];南京大学;2011年
丛春华;[D];中国气象科学研究院;2011年
钟水新;[D];中国气象科学研究院;2011年
潘维玉;[D];南京大学;2011年
李春虎;[D];南京信息工程大学;2011年
姜勇强;[D];南京大学;2011年
邵海燕;[D];南京气象学院;2004年
梁俊明;[D];南京气象学院;2004年
覃丹宇;[D];南京气象学院;2004年
姚学祥;[D];南京气象学院;2004年
中国硕士学位论文全文数据库
袁有林;[D];兰州大学;2015年
【二级参考文献】
中国期刊全文数据库
;[J];Atmospheric and Oceanic Science L2010年01期
;[J];Atmospheric and Oceanic Science L2010年03期
;[J];Atmospheric and Oceanic Science L2011年05期
毛文书;彭骏;周强;张群;苏亚麒;;[J];成都信息工程学院学报;2008年06期
叶笃正,丑纪范,刘纪远,张增祥,王一谋,周自江,鞠洪波,黄签;[J];地理学报;2000年05期
张风,赵思雄;[J];Advances in Atmospheric S2004年05期
孙建华,张小玲,齐琳琳,赵思雄;[J];Advances in Atmospheric S2005年02期
Ghulam RASUL,Qamar-uz-Zaman CHAUDHRY,赵思雄,曾庆存,齐琳琳,张高英;[J];Advances in Atmospheric S2005年03期
Angkool WANGWONGCHAI,赵思雄,曾庆存;[J];Advances in Atmospheric S2005年03期
汤燕冰;甘晶晶;赵璐;高坤;;[J];Advances in Atmospheric S2006年05期
【相似文献】
中国期刊全文数据库
王晓君;马浩;;[J];地球科学进展;2011年11期
罗立辉;张耀南;周剑;潘小多;孙维君;;[J];冰川冻土;2013年03期
于灵雪;张树文;刘廷祥;卜坤;杨久春;;[J];地球信息科学学报;2013年04期
丁煌;陶树旺;肖子牛;程兴宏;;[J];高原气象;2013年02期
张南;龚建东;李泽椿;;[J];科技信息(科学教研);2007年17期
;[J];;年期
中国重要会议论文全文数据库
高丽;魏鸣;潘佳文;陈祺;;[A];第27届中国气象学会年会灾害天气研究与预报分会场论文集[C];2010年
刘春霞;齐义泉;梁建茵;;[A];中国气象学会2005年年会论文集[C];2005年
林惠娟;冀春晓;滕代高;;[A];第27届中国气象学会年会灾害天气研究与预报分会场论文集[C];2010年
何斌;;[A];中国气象学会2006年年会“灾害性天气系统的活动及其预报技术”分会场论文集[C];2006年
何斌;何锋;;[A];第五届长三角气象科技论坛论文集[C];2008年
李萍阳;;[A];推进气象科技创新加快气象事业发展——中国气象学会2004年年会论文集（下册）[C];2004年
何斌;;[A];第三届长三角气象科技论坛论文集[C];2006年
何斌;;[A];中国气象学会2006年年会“中尺度天气动力学、数值模拟和预测”分会场论文集[C];2006年
中国硕士学位论文全文数据库
刘洋;[D];中国海洋大学;2013年
&快捷付款方式
&订购知网充值卡
400-819-9993
《中国学术期刊（光盘版）》电子杂志社有限公司
同方知网数字出版技术股份有限公司
地址：北京清华大学 84-48信箱 知识超市公司
出版物经营许可证 新出发京批字第直0595号
订购热线：400-819-82499
服务热线：010--
在线咨询：
传真：010-
京公网安备75号君，已阅读到文档的结尾了呢~~
（2015大学论文）武汉区域分县可能最大降水预报研究,武汉年降水量,武汉降水量,武汉市降水量,武汉年平均降水量,武汉降水,武汉降水量查询,武汉市洪山区降水状况,武汉年均降水量,武汉工程降水
扫扫二维码，随身浏览文档
手机或平板扫扫即可继续访问
（2015大学论文）武汉区域分县可能最大降水预报研究
举报该文档为侵权文档。
举报该文档含有违规或不良信息。
反馈该文档无法正常浏览。
举报该文档为重复文档。
推荐理由：
将文档分享至：
分享完整地址
文档地址：
粘贴到BBS或博客
flash地址：
支持嵌入FLASH地址的网站使用
html代码：
&embed src='/DocinViewer--144.swf' width='100%' height='600' type=application/x-shockwave-flash ALLOWFULLSCREEN='true' ALLOWSCRIPTACCESS='always'&&/embed&
450px*300px480px*400px650px*490px
支持嵌入HTML代码的网站使用
您的内容已经提交成功
您所提交的内容需要审核后才能发布，请您等待！
3秒自动关闭窗口起报时次：
预报时效&&&&
:00 :00 03小时 :00 06小时 :00 09小时 :00 12小时 :00 15小时 :00 18小时 :00 21小时 :00 24小时 :00 27小时 :00 30小时 :00 33小时 :00 36小时 :00 39小时 :00 42小时 :00 45小时 :00 48小时 :00 51小时 :00 54小时 :00 57小时 :00 60小时 :00 63小时 :00 66小时 :00 69小时 :00 72小时当前位置： >
暖季试验简报总第十七期（2014年第9期）
发布时间：
10:47:47 & 作者：本站编辑 & 来源： 本站原创 & 浏览次数：
&&& 日上午，2014年度最后一期暖季试验暨第九期业务与科研结合试点工作暖季会商在中央气象台会商室西厅举行。来自中科院大气所、干部培训学院等单位的科学家、科研人员和中心业务台室预报员齐聚一堂，共同回顾并探讨8月16-18日发生在西北地区东部和华北等地的强对流天气过程的预报方法和形成机理。
&&&& 强天气预报中心预报员刘鑫华首先对实况进行了展示，并分析了非常规观测资料的特征（如图1、图2），之后给出了预报检验，指出16日、17日和18日08时未来12小时预报的预报均存在空报和漏报的现象。通过实况的分析指出，该过程是在500 hPa中国西部高压脊发展，脊前切断出低涡，低涡东移影响下产生的。接下来，刘鑫华从水汽、不稳定和触发条件三个方面进行了分析。
图1&&& &8月16日08时至17日08时（左图）和8月17日08时至18日08时（右图）强对流监测实况
图2& &8月16日17时红外云图（左图）和8月17日17时云顶亮温低于-52℃的卫星观测（右图）
&&&& 通过对比不同模式不同起报时刻对该次过程的预报，发现GRAPES-MESO对该次过程短时强降水的落区预报有一定的指示意义，降水量级预报偏小。而EC和T639都存在一定程度的位置偏差。EC的降水预报随着预报时效临近降水的位置有一个向西的调整，然而这种向西的调整十分微小，在实际预报中由于时间限制，预报员很难发现这种微小的调整。而NCEP的预报位置要比EC和T639要好，而量级仍然偏弱。在这种情况下，模式的降水预报结果在很大程度上会对预报员做出相对准确的短时强降水落区预报有一定的影响。
&&&& 针对8月16日的西北地区东部的强对流过程，以GFS为例来看，模式对于物理量的预报有很好的指示意义。从垂直风切变来看，西北地区的垂直风切变达到了18-20 m/s（图3）。在西北地区上空也报出了水汽舌（图4）。同时对流有效位能也超过了500 J/kg（图5）。大冰雹指数在该区域也超过了0.2（图6）。因此从诊断物理量方面来讲对于强对流的落区有很好的指示意义。然而，中高层水汽条件预报的差异，以及低层触发条件及水汽的预报偏差对于预报员不能准确判断强对流天气的落区有很大的影响。同时，值得注意的是，预报员原有的统计经验为，当大冰雹指数超过0.4-0.5以上就会考虑发生大冰雹的可能，而此次的预报来看大冰雹指数仅超过了0.2。这也是以后预报过程当中值得注意的问题。
&&&& 而对于西北地区东部和华北地区高低层温差超过30℃条件下，强对流天气落区预报的经验还十分欠缺，预报员还应该不断积累这方面的经验。
&&&& 总体而言，通过分析预报过程中存在如下问题和困惑：
&&&& 1 对于强对流短期预报时效如何有效根据模式结果进行订正
&&&& 2 目前预报时所做的分析都是基于环境背景场的分析，未来应该增加风暴尺度对流特征的物理量的应用
&&&& 3 如何使预报员能够在有效时间内快速锁定分类强对流落区（平台、流程、经验）
&&&&&&& 1）平台的优化和开发
&& &&&& 2）流程的优化
&&&&&&& 3）西北华北的预报经验
图3 GFS模式8月15日20时起报8月16日14时0-6 km垂直风切变
图4 GFS模式8月15日20时起报8月16日14时整层可降水量
图5 GFS模式8月15日20时起报8月16日14时CAPE
图6 GFS模式8月15日20时起报8月16日14时大冰雹指数
&&&&&& 国家气象中心章国材研究员也对此次过程进行了较为详细的分析，在详细分析了天气形势的演变之后，针对预报员提出的问题，从水汽条件、稳定度、垂直运动和地面辐合线四个方面详细分析了16日和17日强对流发生的原因，并给出了详细的物理量特征。认为，对于16日过程，银川850hPa由较湿转干，K指数小于30℃,不会出现强降水;但850-500hPa假相当位温差很大(12K)，CAPE=800-1000,应当出现对流性大风；宁夏中南部850、700hPa都维持湿区，湿层较厚，K指数很大(35-40℃)，有利于出现强降雨; 850-500hPa假相当位温差很大很大(12-16K)， 08时CAPE=1000,强烈不稳定，0、-20℃层高度合适，有利于出现冰雹；山西中北部850hPa维持低湿，仅700hPa出现湿区且增湿，K指数小，不会出现强降水； 850-500hPa假相当位温差维持较大数值(8)， 0、-20℃层高度合适，预报冰雹有难度；河北中南部850、700hPa都维持湿区，湿层较厚，K指数较大(维持35℃)，有利于出现强降雨， 850-500hPa假相当位温差激增(4 & 12K)，08时CAPE=1200,强烈不稳定，0、-20℃层高度合适，有利于出现冰雹。对于17日的过程，17日08时至14时河北内蒙交界区850hPa维持一个湿区(70-80%)，500 hPa 河北北部维持70-80％，湿层深厚，K指数从30激增到40，河北北部有利于产生强降雨; 850-500假相当位温08时已达8K，14时增至12K， 0℃、-20℃层高度,有利于产生冰雹。对于18日的强对流天气条件，从850 hPa 到500 hPa 湿度小，虽然存在影响系统，但难于产生强对流天气。
&&&&& 进而，章国材研究员强调一定要弄清楚实况，平台要给预报提供支撑。大气所孙建华研究员指出，EC是全球模式，强对流不是EC的重点，预报员要依赖中尺度模式，提炼中尺度模式结果。干部培训学院的俞小鼎教授从雷达气候学分型、500 hPa流型辨识、配料法等角度进行了分析和论述。许焕斌研究员从预报员与模式之间的关系角度谈了自己的看法。张小玲正研对于预报中存在的问题进行了分析和阐述。
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
主要撰写 刘鑫华
责任编辑 方& 翀 李萍阳
&&&&&&&&&&&&&
附件下载：