模板堂网站源码,做网站推广和网络推广,做韩国外贸网站,简述网络营销的基本方法前言有些时候为了研究不同高度上气象要素之间的联系#xff0c;纬度-高度作为xy轴的图在一些SCI论文中比较常见。这是我研究的CMIP6数据中ua,va,wap,ta这几个气象要素在纬度-高度图上的不同。读取数据的地方就略去了#xff0c;需要注意的是需要进行平均#xff0c;温度转为…前言有些时候为了研究不同高度上气象要素之间的联系纬度-高度作为xy轴的图在一些SCI论文中比较常见。这是我研究的CMIP6数据中ua,va,wap,ta这几个气象要素在纬度-高度图上的不同。读取数据的地方就略去了需要注意的是需要进行平均温度转为摄氏度#W,V为垂直方向和纬向的速度WW.wap.mean((lon, time,),keep_attrsTrue,skipnaTrue)*100VV.va.mean((lon, time,),keep_attrsTrue,skipnaTrue)T_cT_all.ta-274.15TT_c.mean((lon, time,),keep_attrsTrue,skipnaTrue)# 这时数据格式为(plevel,lat)的格式数据的x轴为laty轴为plev。理解了数据的格式后转换一下思路就可以画图。之后就是画图因为垂直和纬向上风速的量级差的有点大所以我对这两个风速进行了均一化的处理。这样箭头就一样大了。操作plotax.contourf(V.lat, V.plev, V,levelsnp.arange(-0.06,0.060.01,0.01),cmapcmaps.GMT_panoply,extendboth# pivotmiddle,)cbfig.colorbar(plot,axax,orientationhorizontal, shrink0.8,aspect30,pad0.1,)cb.ax.tick_params(labelsize14)以上是V风的垂直填色图下面叠加风矢量图。 间隔取数值不然箭头太密集。xW.lat[::5] #间隔取数值不然箭头太密集yW.plevuV[:,::5]vW[:,::5]u_norm u / np.sqrt(u ** 2.0 v ** 2.0)v_norm v / np.sqrt(u ** 2.0 v ** 2.0)# 为什么要加符号是因为反转y轴的时候箭头方向不会反转所以要加符号ax.quiver(x, y, u_norm, -v_norm,#pivotmiddle,)ax.set_xticks([-90,-75,-50,-25,0,25,50,75,90])ax.set_xticklabels([r90$^\degree$S,r75$^\degree$S,r50$^\degree$S,r25$^\degree$S,r0$^\degree$,r25$^\degree$N, r50$^\degree$N,r75$^\degree$N, r90$^\degree$N])ax.set_ylim(100*100,1000*100)#绘制温度等值线,将负轮廓设置为实线而不是虚线plt.rcParams[contour.negative_linestyle] solidCSax.contour(T.lat,T.plev,T,30,colorsk,)ax.clabel(CS, inline1, fontsize10)ax.invert_yaxis()ax.set_yscale(symlog)#手动设置高度范围ax.set_yticks([100000,925*100,850*100,700*100,50000,20000,10000])ax.set_yticklabels([1000,925,850,700,500,200,100])好了图的效果在这里。e.g.最后需要注意的地方就是反转y轴以及改为对数坐标但这还远远不够。因为可能和我们要的效果不太一样还要手动修改一下ylabel这样和我们要的结果就一致了。
