python气象绘图速成_气象数据可视化——利用Python绘制温度平流

进口作为pltfrom mpl_进口Basemapimport numpy一样np #进口ssimport netCDF4作为ncimport mathfile1 = r 'C:桌面\\ \\ 59799 \\ \\用户\u201Cfile2 = r 'C: \\ 59799 \\ \\用户桌面\\ ' file3 = r 'C:桌面\\ \\ 59799 \\ \\用户\u201D(file1) (file2) (file3) # # #普遍阅读量((\u201C纬度\u201D))((\u201C朗\u201D))((\u201C水平\u201D))# # # u, v, airtemp得到([' uwnd ']) ([' vwnd '])(['空气'])# #水平温度平流计算公式- (udt \/ dx + vdt \/ dy) # #对于离散数据,不做微分操作,使用# #中提到的差分法的计算方法。

但对于dy从赤道到北极是固定# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - # #首先计算dyR = 6371 * 1000 #地球半径6371公里()# #角弧度,经度和纬度区间的数据是°dy = R * #计算dxdxθ= []# # # R是小圆的半径90 n为0,因为我距离(len (lat[0:73])):打印(lat[我])(lat[我])R = R * (arg) #打印(R) (R * ()) (dx) # # dx 90 n的距离差为0每个纬度圈穿过°经度dx [0] = 0;

! ! ! ! ! !没有数据只有一个# # # # # # ! ! ! ! ! ! !没有数据只有一个# # # # # # ! ! !! !只有没有数据((73144))我的范围(len (lat)): j的范围(len(朗)):dTdt (i, j) = - (u (1、2、i, j) * dtdx (i, j) + v (1、2、i, j) * dtdy (i, j)) #图纸(X, Y)(经度[:],lat [1:72]) LON_ () LAT_() =映射技术(lon_0 = lon_0 lat_0 = lat_0 llcrnrlon =朗[0],llcrnrlat = lat [71], urcrnrlon =朗[143],urcrnrlat = lat [1]) () (X, Y, dTdt[1:72:],提出= coolwarm) (c)(1948年1月温度平流())