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

文章预览

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 ………………………………