wrfout的风场要怎么画成流线图呀？？

77777kkkkkk

https://www.ncl.ucar.edu/Applications/Scripts/wrf_gsn_8.ncl

edwardli

matplotlib框架下streamplot必须是均匀分布的严格递增阵列构成网格。Evenly spaced strictly increasing arrays to make a grid. If 2D, all rows of x must be equal and all columns of y must be equal; i.e., they must be as if generated by np.meshgrid(x_1d, y_1d).

使用wrf-python，把数值插值到等经纬度网格上即可。（思路肯定是这个）
参考：wrf.xy_to_ll（我猜的）
另，参考：https://wrf-python.readthedocs.i ... to-a-pressure-level

坎坷

我使用的插值函数，使用的是scipy.interpolate.griddata这个函数，经实测，边缘处有锯齿，其他处基本一致


from scipy.interpolate import griddata

class read_wrf_vars:
    def regrib_data(self, data, lon, lat):
        """
        插值获取网格化数据
        """
        # 截取边界
        min_lon, max_lon, min_lat, max_lat = lon[:, 0].max(), lon[:, -1].min(), lat[0, :].max(), lat[-1, :].min()

        dx = (max_lon - min_lon) / lat.shape[1]
        dy = (max_lat - min_lat) / lat.shape[0]

        resolution = 50

        # dx = np.fabs(lon[0,1] - lon[0,0])
        # dy = np.fabs(lat[0,1] - lat[0,0])

        order = 1
        while True:
            if dx*order >= 50:
                dx = dx*order // resolution * 50
                break
            else:
                order = order * 10
        dx = dx / order

        order = 1
        while True:
            if dy*order >= 50:
                dy = dy*order // resolution * 50
                break
            else:
                order = order * 10
        dy = dy / order

        min_lon = min_lon // dx * dx
        max_lon = max_lon // dx * dx
        min_lat = min_lat // dy * dy
        max_lat = max_lat // dy * dy

        grid_x, grid_y = np.arange(min_lon + dx, max_lon, dx), np.arange(min_lat + dy, max_lat, dy)

        grid_x2d, grid_y2d = np.meshgrid(grid_x, grid_y)

        grid_z0 = griddata(np.array([np.array(lat).flatten(), np.array(lon).flatten()]).T, np.array(data).flatten(), (grid_y2d.flatten(), grid_x2d.flatten()), method='linear')

        grid_z0 = grid_z0.reshape(grid_x2d.shape)

        return grid_z0, grid_x, grid_y