爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 4800|回复: 0

[求助] 绘制时间-高度剖面图报错

[复制链接]
发表于 2022-12-10 16:19:18 | 显示全部楼层 |阅读模式

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

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

x
初学小白想绘制相对湿度的时间-高度剖面,这个错误不知道该如何处理,求助各位大佬~
代码:
import xarray as xr
import numpy as np
import datetime as dt
import cartopy.crs as ccrs
import cartopy.feature as cfeature
import cartopy.mpl.ticker as cticker
import matplotlib.pyplot as plt
import matplotlib.ticker as ticker
from decimal import Decimal
plt.rcParams['font.sans-serif']=['SimHei']
plt.rcParams['axes.unicode_minus']=False

ds=xr.open_dataset(r'F:\20200525.nc')
r=ds.r.loc['2020-05-24 20:00':'2020-05-26 20:00',:,27,98]
r = np.array(r).reshape((17,16))
time=ds.variables['time'][38:55]
level=ds.variables['level'][:]
lat=ds['latitude']
lon=ds['longitude']

fig = plt.figure(figsize=(9,6))
ax = fig.add_axes([0.1, 0.1, 0.8, 0.6])
ax.set_yscale('symlog')
ax.set_yticks([700,500,300,200,100])
ax.set_yticklabels(['700','500','300','200','100'])
ax.invert_yaxis()
ax.set_ylabel('Level (hPa)',fontsize=18)
ax.set_xlabel('time',fontsize=18)
c = ax.contourf(time,level,r.T,level=np.arange(0,100,10),extend = 'both',zorder=0, cmap=plt.cm.bwr)
position=fig.add_axes([0.65, 0.034, 0.3, 0.017])
fig.colorbar(c,cax=position,orientation='horizontal',format='%d',)
报错信息:OverflowError: int too big to convert

wind.txt

1.27 KB, 下载次数: 1, 下载积分: 金钱 -5

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

本版积分规则

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

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

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