Matlab画SkewT-logP图（8月7日更新）

47779125

压高关系采用的是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等的值

期待版友们进一步互动

47779125

沙颖凯 发表于 2014-8-5 11:44
LCL和CAPE的程序我都有，等下我另外开个贴发出来。不过在这上面画CAPE有难度，因为不是每个case都有CAPE

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


模仿NCL的，如下

沙颖凯

改的很好啊，比我那个好看，我有空研究下
还有你把我名字打错了 - -。

freekiller

你们两都继续改，很完美了，鼓掌

Aires

希望你和沙颖凯 一起完善该程序，matlab版还没有这方便的内容，而且这个图用matlab画还是第一次，如果完善了话，希望你能发一个全面的教程，我们会适当加精的哦

tmmiw

支持一下这东东又学习了一次

哈斯巴人

47779125

沙颖凯 发表于 2014-8-3 17:57
改的很好啊，比我那个好看，我有空研究下
还有你把我名字打错了 - -。

不好意思 现在改过来了~

47779125

Aires 发表于 2014-8-4 08:23
希望你和沙颖凯 一起完善该程序，matlab版还没有这方便的内容，而且这个图用matlab画还是第一次，如果完善 ...

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

期待版友们进一步互动

沙颖凯

47779125 发表于 2014-8-5 11:33
不好意思 现在改过来了~

LCL和CAPE的程序我都有，等下我另外开个贴发出来。不过在这上面画CAPE有难度，因为不是每个case都有CAPE