爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 6040|回复: 0

[求助] 探空图绘制LCL语句出错

[复制链接]

新浪微博达人勋

发表于 2022-8-20 12:41:33 | 显示全部楼层 |阅读模式

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

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

x
import matplotlib.pyplot as plt
from mpl_toolkits.axes_grid1.inset_locator import inset_axes
import numpy as np
import pandas as pd

import metpy.calc as mpcalc
from metpy.plots import Hodograph, SkewT
from metpy.units import units

df=pd.read_csv(r'F:/短时科/zhuzhu/57993_2022061300.csv',na_values={'9999'},keep_default_na=False)
#print(df)

#df = df.dropna(subset=('temperature', 'dewpoint', 'direction', 'speed', 'u_wind', 'v_wind'), how='all').reset_index(drop=True)


p = df['pressure'].values[1:150] * units.hPa
T = df['temperature'].values[1:150] * units.degC
Td = df['dewpoint'].values[1:150] * units.degC
wind_speed = df['speed'].values[1:150]*2.5 * units.knots
wind_dir = df['direction'].values[1:150] * units.degrees
u, v = mpcalc.wind_components(wind_speed, wind_dir)
#print(Td)
interval = np.logspace(2, 3, num=25) * units('hPa')
idx = mpcalc.resample_nn_1d(p, interval)

fig = plt.figure(figsize=(9, 9))
skew = SkewT(fig, rotation=30)
skew.plot(p, T, 'r', linewidth=2)
skew.plot(p, Td, 'g', linewidth=2)
skew.plot_barbs(p[idx], u[idx], v[idx])
skew.plot_dry_adiabats()
skew.plot_moist_adiabats()
skew.plot_mixing_lines()
skew.ax.axvline(0, color='c', linestyle='--', linewidth=2)


skew.ax.set_ylim(1000, 100)
skew.ax.set_xlim(-40, 60)

lcl_pressure, lcl_temperature = mpcalc.lcl(p[0], T[0], Td[0])
学习网上大家的语句绘制探空图,但是计算lcl的语句则报错 unsupported operand type(s) for -: 'str' and 'float',请问有大神能帮解答吗。

密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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