气象-探空资料绘图

'''
Author: your name
Date: 2021-08-10 12:00:04
LastEditTime: 2021-08-11 16:22:53
LastEditors: Please set LastEditors
Description: In User Settings Edit
FilePath: \Test\QHTJ\Python气象数据处理与绘图\TanKong.py
'''
import matplotlib.pyplot as plt
from mpl_toolkits.axes_grid1.inset_locator import inset_axes
import pandas as pd
import metpy.calc as mpcalc
from metpy.plots import Hodograph, SkewT
from metpy.units import units
import metpy
import numpy as np
from datetime import datetime
'''
探空站一般是为探测高空气象要素而建立的，通过探空气球来收集每天8点和20点的高空气象数据，遇到特殊天气(台风等)会进行加密观测。
可以获近地层、850、700、500、200百帕的温度、温度露点差、位势高度、风速风向等气象要素。
探空数据在天气预报有着重要的指示作用，可以分析出高空引导气流的位置、强度，及到达本地的时间和对当地天气的影响情况
'''
'''
T-lnP图中线条的物理意义
大气层结曲线：环境空气温度随高度的分布曲线。
大气状态曲线：气块作垂直运动时，温度随高度的变化曲线。

抬升凝结高度（LCL）：气块达到这个高度，水汽达到饱和，再向上会就凝结成水滴。
自由对流高度（LFC）：气块达到这个高度以后不需要任何外力就可以自由上升。
平衡高度（EL）：受力平衡点，垂直速度最大，加速度为0。
对流云顶高度（MPL）：气块能达到的最大高度，取决于EL层以上负不稳定能量刚好抵消CAPE时所达到的高度。

对流有效位能（CAPE）：正不稳定能量，是气块达到自由对流高度后可以获得继续上升的能量。
对流抑制能量（CIN）：负不稳定能量，是大气底层的气块要能达到自由对流高度，至少需从其它途径获得能量的下限
'''

def read_data(path):
    col_names = ['pressure', 'height', 'temperature', 'dewpoint', 'direction', 'speed']#气压、高度、气温、露点、风向、风速
    df = pd.read_fwf(path,skiprows=4, usecols=[0, 1, 2, 3, 6, 7], names=col_names)

    # df = pd.read_fwf('/home/kesci/work/