NCL绘制风场问题请教

飞雪残剑

根据 @成信之子的程序绘制了micaps11类数据的风场图，但是图超出了我的plot边界了，请问如何屏蔽边界外的图呢，也就是不会值plot以外的风场。图如下，期待大佬们的解答，谢谢啦！代码如下：


begin
    wks = gsn_open_wks("png", get_script_prefix_name())
    ;-----------读取数据-----------------
    filePath = "uv.19061620.000"
    nrows = 36
    ncols = 61
    data = asciiread(filePath, (/2,nrows,ncols/), "float")
    ;-----------提取u风、v风-----------------
    u = data(0,:,:)  ;二维数组
    v = data(1,:,:)  ;二维数组
    ;-----------创建lon、lat二维数组-----------------
    lon = new((/nrows,ncols/),"float")
    lat = new((/nrows,ncols/),"float")
    do i = 0,nrows-1
      lat(i,:) = 80.0-i*1.714
    end do
    do j = 0,ncols-1
      lon(:,j) = 10.0+j*3.167
    end do
    shapeFile = "E:/gadm34_CHN_shp/ChinaAccurateMap/china-provinces.shp"
    res               = True
    res@gsnMaximize   = True
    res@gsnDraw       = False
    res@gsnFrame      = False
    res@mpOutlineOn   = True
    res@mpFillOn      = False  
    res@mpDataBaseVersion        = "MediumRes"
    res@pmTickMarkDisplayMode    = "Always"
    res@mpLimitMode              = "LatLon"
    res@mpMinLatF                = 15
    res@mpMaxLatF                = 55
    res@mpMinLonF                = 70
    res@mpMaxLonF                = 135
    map = gsn_csm_map(wks,res)
    pres             = True
    pres@gsLineColor = "blue"
    poly = gsn_add_shapefile_polylines(wks,map,shapeFile,pres)
    draw(map)
    ;wmsetp("col", 2)   
    wmsetp("wbs", .03)
    wmbarb(wks,lon,lat,-u*2.5,-v*2.5)
    frame(wks)
end