温度平流计算结果离谱

liuna

请家园的各位大佬帮忙看看到底是哪里的问题呀，温度平流结果好奇怪
下面是代码和输出图
        dlon = (np.gradient(lon0) * np.pi / 180).reshape((1, -1))
        dlat = (np.gradient(lat0) * np.pi / 180).reshape((-1, 1))
        coslat = (np.array(np.cos((lat0) * np.pi / 180))).reshape((-1, 1))
        dx = a * coslat * dlon
        dy = a * dlat

        # 计算温度梯度
        t_np = np.array(t)
        dT_dx = np.gradient(np.array(t), axis=1) / dx
        dT_dy = np.gradient(np.array(t), axis=0) / dy

        temadvection0 = -(u0 * dT_dx + v0 * dT_dy)

        return lon0, lat0, u0, v0, temadvection0

贫道敬孔

可以直接使用mpcalc.advection

liuna

贫道敬孔 发表于 2025-7-6 09:29
可以直接使用mpcalc.advection

好的，感谢，但是为什么这样算出来结果不对嘞

tosaka

估计是dxdy计算有问题吧，再试试

liuna

tosaka 发表于 2025-7-6 23:00
估计是dxdy计算有问题吧，再试试

好的，谢谢啦

呆妹小霸王

没有图呀，把温度和流场画出图来自己判断下正负对不对。