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

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 183052|回复: 125

[源程序] Matlab画SkewT-logP图(8月7日更新)

  [复制链接]

新浪微博达人勋

发表于 2014-8-3 17:41:22 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 47779125 于 2014-8-7 10:59 编辑

压高关系采用的是62年美国标准大气(可在代码中选择热带大气)
legend:

粗线:
黑色实线:温度探空
蓝色实线:露点温度探空
红色虚线:状态曲线

细线:
浅红色:干绝热线(等位温线),单位:摄氏度
青绿色:湿绝热线(等theta_e线),单位:摄氏度
黄绿色:等饱和比湿线,单位:g/kg

灰白相间的线条:等温线,单位:摄氏度






以下是更新日志:


7月31日更新:

想编skewT-logP图的愿望很久,但因为忙于准备毕业论文而屡屡推迟。
最近看到版友沙颖凯发了SkewT-logP的代码,才想起心中有个心愿未了。
在仔细欣赏了沙颍凯的程序后,发现代码中还存在一些问题:
1.既然等温线是倾斜直线,必须要画
2.干绝热线、湿绝热线的数值?(原代码是根据最大值和最小值直接分成24等份,每张图都不一致,这不利于读者阅读)
3. 等theta_e的计算?
4. 气压与高度的对应(作者采用的是等温度递减率的大气压高公式?),200hpa才11公里是不是好像不太对?
5. 让等温线与过1000hPa还是1050hPa的平行线交点等于实际值?
6.x轴范围是不是过大,是否需要legend,图像的大小长宽比是不是应该由用户决定?
7.等等...

具体请见附件。其他文件同版友沙颖凯的代码
http://bbs.06climate.com/forum.php?mod=viewthread&tid=27066&extra=page%3D2

8月6日更新
1. 添加计算抬升凝结高度的代码
2. 添加状态曲线

8月7日更新
1.有网友提出,风向杆显示有问题。我仔细检查了一下,这个问题主要是对数坐标造成的。在现在的代码中,所有y轴的长度都乘以了log(y)的导数,即1/y。
2. 图片大小由用户设置
3. 调整了线条颜色粗线。


未解决的问题:
1.一些物理量如CAPE, PW等的值

期待版友们进一步互动


示意图

示意图

skewTlogP_plot.m

9.96 KB, 阅读权限: 10, 下载次数: 443, 下载积分: 金钱 -5

test_skewTlogP.m

622 Bytes, 阅读权限: 10, 下载次数: 283, 下载积分: 金钱 -5

windbarb_profile.m

6.5 KB, 阅读权限: 10, 下载次数: 306, 下载积分: 金钱 -5

soundings_RF6.mat

19.94 KB, 阅读权限: 10, 下载次数: 243, 下载积分: 金钱 -5

点评

感觉棒棒哒  发表于 2015-10-20 15:50

评分

参与人数 9金钱 +119 贡献 +22 收起 理由
杜超杰 + 10
ndycq + 5 很给力!
wlzhongouc + 20 + 2
陌小默的 + 10 + 2 很给力!
raul928 + 10 + 2 赞一个!
二爷名声在外 + 15 + 5
Aires + 14 + 4
mofangbao + 15 + 5
斥鷃 + 20 + 2 鼓励讨论!

查看全部评分

本帖被以下淘专辑推荐:

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

新浪微博达人勋

 楼主| 发表于 2014-8-5 11:48:19 | 显示全部楼层
沙颖凯 发表于 2014-8-5 11:44
LCL和CAPE的程序我都有,等下我另外开个贴发出来。不过在这上面画CAPE有难度,因为不是每个case都有CAPE

不需要画,只需要计算数值。在title下面显示就行。


模仿NCL的,如下

NCL中的T-logP

NCL中的T-logP
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2014-8-3 17:57:03 | 显示全部楼层
本帖最后由 沙颖凯 于 2014-8-3 17:58 编辑

改的很好啊,比我那个好看,我有空研究下
还有你把我名字打错了 - -。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-8-3 22:24:24 | 显示全部楼层
你们两都继续改,很完美了,鼓掌
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-8-4 08:23:03 | 显示全部楼层
希望你和沙颖凯 一起完善该程序,matlab版还没有这方便的内容,而且这个图用matlab画还是第一次,如果完善了话,希望你能发一个全面的教程,我们会适当加精的哦
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-8-4 09:00:16 | 显示全部楼层
支持一下这东东又学习了一次
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-8-5 11:09:15 | 显示全部楼层
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2014-8-5 11:33:32 | 显示全部楼层
沙颖凯 发表于 2014-8-3 17:57
改的很好啊,比我那个好看,我有空研究下
还有你把我名字打错了 - -。

不好意思 现在改过来了~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-8-5 11:34:20 | 显示全部楼层
Aires 发表于 2014-8-4 08:23
希望你和沙颖凯 一起完善该程序,matlab版还没有这方便的内容,而且这个图用matlab画还是第一次,如果完善 ...

是的,未解决的问题主要是以下俩:
1. 没有计算凝结高度、CAPE值等物理量
2. 没有画状态曲线

期待版友们进一步互动
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-8-5 11:44:16 | 显示全部楼层
47779125 发表于 2014-8-5 11:33
不好意思 现在改过来了~

LCL和CAPE的程序我都有,等下我另外开个贴发出来。不过在这上面画CAPE有难度,因为不是每个case都有CAPE
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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