等高线图一直只是0上的一个点出现，但是当我手动计算等高线时，它们是变化的并且是实质性的。不确定我的代码或数学有问题吗？ ：
Re = 6371000
rho = np.linspace(-20*Re, 20*Re, 2000)
z1 = np.linspace(-10*Re, 10*Re, 2000)
z = z1[:, np.newaxis]
Be = 3120*10**(-9)
F = (rho**2)*Be*(Re/np.sqrt(rho**2 + z**2))**3
plt.contour(rho, z1, F)