强对流天气的中尺度分析_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
文档贡献者
评价文档：
&&¥2.00
&&¥2.00
&&¥2.00
&&¥2.00
&&¥2.00
&&¥1.00
&&¥3.00
喜欢此文档的还喜欢
强对流天气的中尺度分析
把文档贴到Blog、BBS或个人站等：
普通尺寸(450*500pix)
较大尺寸(630*500pix)
大小：31.88MB
登录百度文库，专享文档复制特权，财富值每天免费拿！
你可能喜欢散度、涡度与天气--《甘肃气象》1985年04期
散度、涡度与天气
【摘要】：正 县站的预报工作,自省局停止语言广播之后,123传真图便成了极其重要的工具,它包括散度图、涡度图和暴雨、冰雹围区预报图等等。如何正确理解和使用这些图,就成为一个十分重要的问题。下面就散度、涡度图的使用,作一简单介绍,供县站预报员参考。
【作者单位】：
【关键词】：
【正文快照】：
县站的预报工作,自省局停止语言广播之后,123传真图便成了极其重要的工具,它包括散度图、涡度图和暴雨、冰雹围区预报图等等。如何正确理解和使用这些图,就成为一个十分重要的问题。下面就散度、涡度图的使用,作一简单介绍,供县站预报员参考。散度形式相反的两组图。 水平散
欢迎：、、)
支持CAJ、PDF文件格式，仅支持PDF格式
【相似文献】
中国期刊全文数据库
张云惠;钱文新;;[J];新疆气象;2005年S1期
郑婧;朱星球;肖安;;[J];气象与减灾研究;2007年02期
张磊;李毓富;王贵朋;赵元富;;[J];科技情报开发与经济;2008年07期
黄庆国;覃彦英;罗思泽;;[J];气象研究与应用;2009年02期
邵振平;谷秀杰;张云平;郭献林;李定才;;[J];气象与环境科学;2009年S1期
梁明增;黄骏;许丽华;朱金菊;;[J];现代农业科技;2009年19期
李亚春;;[J];内蒙古气象;2009年05期
黄仁立;韦华红;;[J];气象研究与应用;2010年S2期
许新田;宁志谦;陶建玲;郭大梅;;[J];陕西气象;2007年04期
蒋元政;[J];气象;1984年01期
中国重要会议论文全文数据库
徐凤梅;李静;;[A];中国气象学会2008年年会极端天气气候事件与应急气象服务分会场论文集[C];2008年
孙莹;郑贤;;[A];中国气象学会2008年年会复杂地形影响下的天气与气候分会场论文集[C];2008年
李海燕;巴哈古丽;赵威;;[A];中国气象学会雷达气象学与气象雷达委员会第二届学术年会文集[C];2006年
曹俊元;查贲;倪琛玲;;[A];中国气象学会2006年年会“灾害性天气系统的活动及其预报技术”分会场论文集[C];2006年
王淑云;白春艳;周连科;;[A];中国气象学会2005年年会论文集[C];2005年
魏香;王华;;[A];中国气象学会2007年年会天气预报预警和影响评估技术分会场论文集[C];2007年
张青;叶飞;周月中;;[A];第七届长三角气象科技论坛论文集[C];2010年
沙文生;;[A];2004年滇桂黔三省区五州市气象协作区交流会特刊论文集[C];2004年
孙兴池;王文毅;王业宏;;[A];2006年海峡两岸气象科学技术研讨会论文集[C];2006年
孙兴池;王文毅;王业宏;;[A];中国气象学会2006年年会“灾害性天气系统的活动及其预报技术”分会场论文集[C];2006年
中国重要报纸全文数据库
章国材;[N];中国气象报;2007年
刘晓平 宫卫平;[N];中国气象报;2003年
纪晓玲　杨侃;[N];中国气象报;2010年
李红莉　林春泽;[N];中国气象报;2011年
陆铭;[N];中国气象报;2011年
王海;[N];中国气象报;2011年
中国博士学位论文全文数据库
王慧;[D];南京信息工程大学;2005年
慕建利;[D];南京信息工程大学;2009年
覃丹宇;[D];南京气象学院;2004年
姜勇强;[D];南京大学;2011年
中国硕士学位论文全文数据库
高维英;[D];兰州大学;2007年
刘学刚;[D];南京信息工程大学;2007年
王建鹏;[D];南京气象学院;2003年
关月;[D];中国科学院研究生院（大气物理研究所）;2007年
王勇;[D];南京信息工程大学;2007年
韦晋;[D];南京信息工程大学;2008年
武文辉;[D];南京信息工程大学;2007年
刘国忠;[D];南京信息工程大学;2007年
柳草;[D];中国气象科学研究院;2009年
丁做尉;[D];中国海洋大学;2010年
&快捷付款方式
&订购知网充值卡
400-819-9993
《中国学术期刊（光盘版）》电子杂志社有限公司
同方知网数字出版技术股份有限公司
地址：北京清华大学 84-48信箱 知识超市公司
出版物经营许可证 新出发京批字第直0595号
订购热线：400-819-82499
服务热线：010--
在线咨询：
传真：010-
京公网安备74号后使用快捷导航没有帐号？
只需一步，快速开始
扫一扫，访问微社区
只需一步, 快速开始
查看: 3761|回复: 73
matlab读取NCEP细网格资料及计算散度、涡度平流、温度平流、K指数等要素
在线时间 小时
登录后查看更多精彩内容~
才可以下载或查看，没有帐号？
本帖最后由 F117_ren_0 于
11:36 编辑
首先感谢气象家园各位元老，学习和借鉴了很多有用的东西，具体名字就不再点出来了，再次表示感谢！读取FNL数据首先安装read_grib工具箱和M_map工具箱，具体安装方法请百度！
1,读取FNL数据
grib_struct=read_grib('F:\FNL\fnl__00_c',-1); %%%读取NCEP再分析资料
& & read_grib('F:\FNL\fnl__00_c','inv'); %生成变量目录表文档
& & grib_struct(13).
& & Hight_500hPa=grib_struct(213).
%Hight_hPa=Hight_1000hPa_01;
for i=1:91;
& & for j=1:181;
& && &&&lon(i,j)=181-j;
& && &&&lat(i,j)=91-i;
& && &&&Hight(i,j)=Hight_500hPa((i-1)*360+j)/10.0;
这里是grib1格式，好像是07年以后是grib2格式，但读法一样，只是从排列发生变化，具体可以见%生成变量目录表文档%想读取的数据 Hight_500hPa=grib_struct(213).fltarray里，只要改动括号里的数字就ok了
2,计算散度
%%%%%%%%%%%%%%%%%%%%读取NFL
& & grib_struct=read_grib('F:\FNL\fnl__00_c',-1); %%%读取NCEP再分析资料
& & read_grib('F:\FNL\fnl__00_c','inv'); %生成变量目录表文档
& & %grib_struct(13).&&
& & U850=grib_struct(157).
& & V850=grib_struct(158).
for i=34:50;
& & for j=70:100;
& && &&&lon(i,j)=100-j;
& && &&&lat(i,j)=50-i;
& && &&&u850(i,j)=U850((i-1)*360+j);
& && &&&v850(i,j)=V850((i-1)*360+j);
u850=u850(34:end,70:end);
v850=v850(34:end,70:end);
& &R0=6371004;&&pi=3.1415926;
& &lat0=50;&&lon0=70;
& & for lati=50:-1:34
& && &latm(50-lati+1)= distance(lat0,lon0,lati,lon0)/180*pi*R0;
& &for loni=70:1:100
& && &lonm(loni-69)= distance(lat0,lon0,lat0,loni)/180*pi*R0;
& &[xx,yy]=meshgrid(lonm,latm);
& &div=divergence(xx,yy,u850,v850);&&
&&save D:\matt\div850 div xx yy
850hPa散度计算完成，这里的经纬度我选的是新疆的范围，如有需要可自行改动。
& & & & & & & & & & & & & & & &
3、计算涡度平流和温度平流
首先读取出UV，涡度，温度
lon=70:100;
lat=34:50;
PI=3.1415926;
R0=6371004;
[dtx850,dty850]=gradient(w850,1,1);
dx=gradient(lon,1)*PI/180;
dy=gradient(lat,1)*PI/180;
wp850=-((u850.*dtx850)./(dx.'*cos(lat*PI/180))'+v850.*(dty850'*diag(1./dy))')/R0;
以上w850为850hPa的涡度，wp850为计算出来的850hPa的涡度平流
[dtx850,dty850]=gradient(t850,1,1);
tp850=-((u850.*dtx850)./(dx.'*cos(lat.*PI/180))'+v850.*(dty850'*diag(1./dy))')/R0;
以上t850为850hPa的温度，tp850为计算出来的850hPa的温度平流。
& & & & & & & & & & & & & & & &
& & & & & & & & & & & & & & & &
4、计算K指数
%%%%%%%%%%%%%%%%%%K指数
% K=(T850-T500)+Td850-(T700-Td700)
%%%%%%%%%%%%%%%%%%
load D:\matt\Temp
load D:\matt\RH
%%%%t&0&&a=7.5; b=237.3;
%%%%t&=0 a=9.5 b=265.5
for i=1:17;
& & for j=1:31;
if t850(i,j)&0;
& & a=7.5;
& & b=237.3;
& & a=9.5;
& & b=265.5;
qq1(i,j)=a*t850(i,j);
qq2(i,j)=b+t850(i,j);
qq=qq1./qq2;
Es=E0*(10.^qq);%%%%%%%%饱和水汽压
E=Es.*rh850.*0.01;&&%%%该相对湿度下水汽压
T1=b*log10(E./Es);
T2=a-log10(E./Es);
Td850=T1./T2;
save D:\matt\Td850 Td850
load D:\matt\Temp
load D:\matt\RH
%%%%t&0&&a=7.5; b=237.3;
%%%%t&=0 a=9.5 b=265.5
for i=1:17;
& & for j=1:31;
if t700(i,j)&0;
& & a=7.5;
& & b=237.3;
& & a=9.5;
& & b=265.5;
qq1(i,j)=a*t700(i,j);
qq2(i,j)=b+t700(i,j);
qq=qq1./qq2;
Es=E0*(10.^qq);%%%%%%%%饱和水汽压
E=Es.*rh700.*0.01;&&%%%该相对湿度下水汽压
T1=b*log10(E./Es);
T2=a-log10(E./Es);
Td700=T1./T2;
save D:\matt\Td700 Td700
load D:\matt\Temp
load D:\matt\Td850
load D:\matt\Td700
K=(t850-t500)+Td850-(t700-Td700);
save D:\matt\Kzhishu K
首先要读出湿度和温度，我以上程序中load出来的就是存好的数据，后面就是计算了。
& & & & & & & & & & & & & & & &
load D:\matt\wdpll;
fnshp_P='D:\Program Files\toolbox\m_map\bou2_4m\bou2_4p.shp';% 省界
S = shaperead(fnshp_P, 'UseGeoCoords', true,...&&
&&'Selector',{@(NAME) strcmp(NAME,'新疆维吾尔自治区'), 'NAME'});
figure(1);
set (gca,'Position',[0.1,0.1,0.85,0.85], 'color','w')
lo=70:100;
[lo,la]=meshgrid(lo,la);
wp850=wp850.*10.^10;
wp850(isnan(wp850)==1)= 0;
[x,y]=meshgrid(70:0.1:100,34:0.1:50);
KK=griddata(lo(:),la(:),wp850(:),x,y,'v4');
m_proj('Equidistant Cylindrical','lon',[70,100],'lat',[34,50]);
[cs,h]=m_contour(x,y,KK,'LineWidth',2);
set(h,'ShowText','on','LevelList',[-500 -350 -200 -100 -50 0 50 100 200 350 500]);
clabel(cs,h,'fontsize',5);
m_grid('box','fancy');
%axis equal
m_plot(S.Lon,S.Lat,'-k','linewidth',2);
title('850涡度平流');
其他的也一样，就不一一表示了。
给自己设置一个个性的签名吧~
在线时间 小时
图片不让我上传，没图没真相啊，呵呵
给自己设置一个个性的签名吧~
在线时间 小时
print出来的图片不让贴出来？奇怪!
给自己设置一个个性的签名吧~
在线时间 小时
楼主好牛逼的样子，竟然知道grib数据的内部结构，真了不起
给自己设置一个个性的签名吧~
在线时间 小时
&成长值: 5415
为啥不能加图呢？
给自己设置一个个性的签名吧~
在线时间 小时
楼主好牛逼的样子，竟然知道grib数据的内部结构，真了不起
还是笨笨牛b，哈哈
给自己设置一个个性的签名吧~
在线时间 小时
&成长值: 3870
貌似楼主很牛逼的样子，狂支持，下载了！
给自己设置一个个性的签名吧~
在线时间 小时
这么详细啊 学学
给自己设置一个个性的签名吧~
在线时间 小时
为啥不能加图呢？
不知道，也许新疆这边不让上图吧，我已经让广州笨笨帮我上图了
给自己设置一个个性的签名吧~
在线时间 小时
高了半天终于插进来了，好费劲
给自己设置一个个性的签名吧~您所在位置： &
&nbsp&&nbsp&nbsp&&nbsp
数值天气预报产品在航空气象中的释用-气象出版社.pdf119页
本文档一共被下载：
次 ,您可免费全文在线阅读后下载本文档
文档加载中...广告还剩秒
需要金币：100 &&
数值天气预报产品在航空气象中的释用-气象出版社.pdf
你可能关注的文档：
··········
··········
世界气象组织技术手册
数值天气 预报产 品在航 空气象 中的释用
CA eM 航空气象先进技术应用工作小组编写
中国民航总局空管局气象处
本书较为全面地介绍了数值天气预报产品在航空气象中的释用,
内容包括全球模式、中
尺度模式、模式产品的直接使用、评估与检验,
未来展望等五个章节 。书后附有常用缩略语
本书可供航空气象预报人员参考使用。
气象出版社出版
北 京市 中关村 南大街 46
北京市 白河印刷厂印刷
开 本: 7 87× 1092 1/
16 印张: 7 .5 字数 : 185
2 0 01 年 1 1 月第 1 版 2 00 1 年 11 月第 1 次印刷
1― 110 0 定价 : 2 5.0 0 元
统一 书号 : 135 0 29 - 52 9 0
《数值天气预报产品在航空气象中的释用》一书是世界气象组织的技术文献,
本书没有太
多的理论研究和公式推导,
而从数值模式开始,
介绍了数值预报的产品在航空气象业务中的
正在加载中，请稍后...