【已解决】求助！dim_avg_n_Wrap求wrfout垂直平均风场有区域空白

Mocca

如图，125E东侧无箭头，很奇怪。
用dim_avg_n_Wrap求再分析数据的同一时刻同一区域的垂直平均风场是有箭头的，另外分别画wrfout里面对应各个层次的风场也无缺测。相关代码如下：

        U = wrf_user_getvar(f1,"U",it)    ; 3D U at mass points               u = wrf_user_unstagger(U,U@stagger)
        V = wrf_user_getvar(f1,"V",it)    ; 3D V at mass points
        v = wrf_user_unstagger(V,V@stagger)

        pressure1 = (/700,600,550,500/)
        pressure2 = (/850,800,750,700/)
        u_850_ = wrf_user_intrp3d(u,p,"h", pressure2, 0., False)
        v_850_=  wrf_user_intrp3d(v,p,"h", pressure2, 0., False)
        u_500_ = wrf_user_intrp3d(u,p,"h", pressure1, 0., False)
        v_500_=  wrf_user_intrp3d(v,p,"h", pressure1, 0., False)

        u_850_!0                      = "lev"
        u_850_!1                      = "lat"
        u_850_!2                      = "lon"        
        u_850_&lev                    = pressure2
        u_850_@lon2d                    = lon
        u_850_@lat2d                    = lat
        u_850_@level_type            = "Isobaric surface (hPa)"
        u_850_@units                 = "m/s"
        u_850_@_FillValue            = 9999.
        copy_VarCoords(u_850_, v_850_)

        u_500_!0                      = "lev"
        u_500_!1                      = "lat"
        u_500_!2                      = "lon"        
        u_500_&lev                    = pressure1
        u_500_@lon2d                    = lon
        u_500_@lat2d                    = lat
        u_500_@level_type            = "Isobaric surface (hPa)"
        u_500_@units                 = "m/s"
        u_500_@_FillValue            = 9999.
        copy_VarCoords(u_500_, v_500_)


        u_500 = u_500_(lat|:,lon|:,lev|:)    ; reorder the dimensions
        v_500 = v_500_(lat|:,lon|:,lev|:)
        u_850 = u_850_(lat|:,lon|:,lev|:)
        v_850 = v_850_(lat|:,lon|:,lev|:)


        u_850_a = dim_avg_n_Wrap(u_850(:,:,:),2)
        v_850_a = dim_avg_n_Wrap(v_850(:,:,:),2)
        u_500_a = dim_avg_n_Wrap(u_500(:,:,:),2)
        v_500_a = dim_avg_n_Wrap(v_500(:,:,:),2)

啦啦海绵

楼主最后解决了吗？我是用的ERA资料 画出来的10m风数出来是负值，图片里不显示风...

Mocca

啦啦海绵 发表于 2021-8-13 16:58
楼主最后解决了吗？我是用的ERA资料 画出来的10m风数出来是负值，图片里不显示风...

啊..我后来发现是因为后面画图属性gsnAddCyclic=True了，应该是False。

啦啦海绵

Mocca 发表于 2021-8-14 11:08
啊..我后来发现是因为后面画图属性gsnAddCyclic=True了，应该是False。

我后来发现我读取数据的经纬度错了  自己太粗心了哈哈