##数据来源:NCEP再分析数据
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap
import numpy as np
#import scipy.stats as ss
import netCDF4 as nc
import math
file1=r'C:\Users\59799\Desktop\uwnd.mon.mean.nc'
file2=r'C:\Users\59799\Desktop\vwnd.mon.mean.nc'
file3=r'C:\Users\59799\Desktop\air.mon.mean.nc'
uwnd=nc.Dataset(file1);vwnd=nc.Dataset(file2)
air=nc.Dataset(file3)
###通用量读取
lat=np.array(uwnd.variables['lat'])
lon=np.array(uwnd.variables['lon'])
level=np.array(uwnd.variables['level'])
###u,v,airtemp获取
u=np.array(uwnd.variables['uwnd'])
v=np.array(vwnd.variables['vwnd'])
airtemp=np.array(air.variables['air'])
##水平温度平流计算公式为 -(udt/dx+vdt/dy)
##对于离散数据,无法做微分会偏微分运算,采取计算方法中
#