【新提醒】Matlab画SkewT

本帖最后由 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等的值

期待版友们进一步互动