专栏名称: 气象学家
【气象学家】公众号平台为您解读最新气象科研进展、分享气象实用编程技巧、追踪气象即时资讯。欢迎加入气象AI和Python交流群以及气象博士群!与5W+的专业人士一起交流互动!
今天看啥  ›  专栏  ›  气象学家

气象编程 | Python绘制全球涡度、散度和垂直速度

气象学家  · 公众号  ·  · 2024-07-13 16:11

文章预览

Python 绘制全球涡度、散度和垂直速度 作者:气象互助社-晓琛 邮箱:yuxichen1021@163.com 数据:nc数据 • 时间:2020年8月12日;水平分辨率:2.5° 2.5° • 层次(uv):1000hPa、850 hPa • 空间范围:30°S-60°N(37个格点),20°E-160°E (57个格点) 计算部分 01  计算部分 import numpy as np import xarray as xr import matplotlib.pyplot as plt import cartopy.crs as ccrs # 读取netCDF文件 data1 = xr.open_dataset( 'uwnd.20200812.nc' ) data2 = xr.open_dataset( 'vwnd.20200812.nc' ) data3 = xr.open_dataset( 'omega.20200812.1000hPa.nc' ) LAT = data1[ 'lat' ].values LON = data1[ 'lon' ].values ORGANnC1 = data1[ 'uwnd' ] ORGANnC2 = data2[ 'vwnd' ] ORGANnC3 = data3[ 'omega' ] # 提取 1000 hpa 层级的数据 u1000 = ORGANnC1.sel(time= "2020-08-12" ,level= '1e+03' ).values v1000 = ORGANnC2.sel(time= "2020-08-12" ,level= '1e+03' ).values # 提取 850 hpa 层级的数据 u850 = ORGANnC1.se ………………………………

原文地址:访问原文地址
快照地址: 访问文章快照
总结与预览地址:访问总结与预览