python源码在此，请问怎么把nc文件数据在可视化时将lev从底层到高层排列？源码分享

王先生

# -*- coding: utf-8 -*-

"""
Created on Fri Mar  24 17:47:41 2019
@author: WYL
"""
from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
import numpy as np
import netCDF4 as nc
from netCDF4 import Dataset

obj = nc.Dataset("f:/merra2/QV/MERRA2_400.inst3_3d_asm_Np.20120601.SUB.nc4")
time_bnds = obj.variables['time_bnds'][:]
time = obj.variables['time'][0]
lev = obj.variables['lev'][:]
lat = obj.variables['lat'][:]
lon = obj.variables['lon'][220:280]
QV = obj.variables['QV'][0,:,100,220:280]
plt.figure(figsize=(16,8))
lons,levs = np.meshgrid(lon,lev)
x,y = lons,levs
lev=lon[::-1]
shade = plt.contourf(x,y,QV,20,cmap = 'jet')
plt.colorbar(shade)
print(help(plt.contourf))
plt.show()
如上传的图片，现在绘制的是个水汽（QV）垂直切面分布图。画图的时候发现nc4文件中lev项是从高空到地面的，求问写代码时候怎么反过来啊。。。这图是倒着的。源码分享给大家，希望会的帮看下应该怎么做更改，刚学了不久。
感谢@

Masterpiece

的指点，问题目前已经解决。代码放着给有需要的和我一样的新手。如果对代码有什么疑问或者代码错误，请及时指正，本人日常在线

Masterpiece

plt.gca().invert_yaxis()

王先生

自己给自己顶一下

王先生

Masterpiece 发表于 2019-3-25 21:26
plt.gca().invert_yaxis()

谢谢大神指点，成了。请问下这个语句只是将Y轴反转了一下是吗

Masterpiece

王先生 发表于 2019-3-25 21:34
谢谢大神指点，成了。请问下这个语句只是将Y轴反转了一下是吗

是的，正如你所见

LemonTree007

气压的坐标轴应该是取对数的，不是等距离分布的，要怎么弄。grads里是set zlog on

王先生

LemonTree007 发表于 2019-3-25 21:46
气压的坐标轴应该是取对数的，不是等距离分布的，要怎么弄。grads里是set zlog on

你问的这个我还真不知道，这里我的数据是位势高度，没考虑你说的这个问题，我现在去研究一下。感谢提问。有答案了我会及时回复你

王先生

LemonTree007 发表于 2019-3-25 21:46
气压的坐标轴应该是取对数的，不是等距离分布的，要怎么弄。grads里是set zlog on

https://blog.csdn.net/m0_37362454/article/details/81510617，在CSDN里找到了这个你看一下，然而我新手看不懂，觉得应该是他里面提到的操作方法。如果看会了回头指导我下

Masterpiece

LemonTree007 发表于 2019-3-25 21:46
气压的坐标轴应该是取对数的，不是等距离分布的，要怎么弄。grads里是set zlog on

import matplotlib.ticker as ticker
ax.set_yscale('log')
ax1 = plt.gca().yaxis
ax1.set_major_formatter(ticker.ScalarFormatter())
plt.minorticks_off()
plt.yticks([100,300,500,700,850,925,1000])

王先生

Masterpiece 发表于 2019-3-25 23:50
import matplotlib.ticker as ticker
ax.set_yscale('log')
ax1 = plt.gca().yaxis

ax和ax1区别是什么呢？不太懂，可以说的稍微详细点么谢谢