文章预览
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
………………………………