求大佬，给地图规定区域添加颜色，其余地区空白

卑微江小澄

只想在青藏高原添加颜色，其它地区保持空白。求大佬
begin
all_files = systemfunc ("ls GLDAS_CLM_TP.*.nc")
fall      = addfiles (all_files, "r")                     
ST         = fall[:]->var85

u1    =ST(0:11,1:4,{24.500:41.500},{72.500:105.500})   ;从CLM中提取数据
u1_avg = dim_avg_n_Wrap(u1,0)                    ;1年的平均
u1_avg_1 = dim_avg_n_Wrap(u1_avg,0)              ;深度平均
printVarSummary(u1_avg_1)
u2    =ST(12:23,1:4,{24.500:41.500},{72.500:105.500})   ;从CLM中提取数据
u2_avg = dim_avg_n_Wrap(u2,0)                    ;1年的平均
u2_avg_1 = dim_avg_n_Wrap(u2_avg,0)      

u3    =ST(24:35,1:4,{24.500:41.500},{72.500:105.500})   ;从CLM中提取数据
u3_avg = dim_avg_n_Wrap(u3,0)                    ;1年的平均
u3_avg_1 = dim_avg_n_Wrap(u3_avg,0)

u4    =ST(36:47,1:4,{24.500:41.500},{72.500:105.500})   ;从CLM中提取数据
u4_avg = dim_avg_n_Wrap(u4,0)                    ;1年的平均
u4_avg_1 = dim_avg_n_Wrap(u4_avg,0)

u5    =ST(48:59,1:4,{24.500:41.500},{72.500:105.500})   ;从CLM中提取数据
u5_avg = dim_avg_n_Wrap(u5,0)                    ;1年的平均
u5_avg_1 = dim_avg_n_Wrap(u5_avg,0)

u6    =ST(60:71,1:4,{24.500:41.500},{72.500:105.500})   ;从CLM中提取数据
u6_avg = dim_avg_n_Wrap(u6,0)                    ;1年的平均
u6_avg_1 = dim_avg_n_Wrap(u6_avg,0)

u7    =ST(72:83,1:4,{24.500:41.500},{72.500:105.500})   ;从CLM中提取数据
u7_avg = dim_avg_n_Wrap(u7,0)                    ;1年的平均
u7_avg_1 = dim_avg_n_Wrap(u7_avg,0)

u8    =ST(84:95,1:4,{24.500:41.500},{72.500:105.500})   ;从CLM中提取数据
u8_avg = dim_avg_n_Wrap(u8,0)                    ;1年的平均
u8_avg_1 = dim_avg_n_Wrap(u8_avg,0)

u9    =ST(96:107,1:4,{24.500:41.500},{72.500:105.500})   ;从CLM中提取数据
u9_avg = dim_avg_n_Wrap(u9,0)                    ;1年的平均
u9_avg_1 = dim_avg_n_Wrap(u9_avg,0)

u10    =ST(108:120,1:4,{24.500:41.500},{72.500:105.500})   ;从CLM中提取数据
u10_avg = dim_avg_n_Wrap(u10,0)                    ;1年的平均
u10_avg_1 = dim_avg_n_Wrap(u10_avg,0)

STA=(/(u1_avg_1),(u2_avg_1),(u3_avg_1),(u4_avg_1),(u5_avg_1),(u6_avg_1),(u7_avg_1),(u8_avg_1),(u9_avg_1),(u10_avg_1)/)

STA@long_name="soil_temperature_avg"
STA!0="year"

STA!1="lat"

STA!2="lon"
copy_VarCoords(u1_avg_1(lat|:, lon|:), STA(0,:,:))

printVarSummary(STA)
rs  =   trend_manken(STA, False, 0)
copy_VarCoords(STA(lat|:, lon|:,year|:), rs(0, :,:))

printVarSummary(rs)
wks = gsn_open_wks("png", "987")
    gsn_define_colormap(wks, "BlueRed")
       
       
    res = True
    res@gsnFrame = False
    res@gsnDraw = False
    res@gsnAddCyclic = False   
        res@cnLevelSelectionMode = "ManualLevels"
        res@cnMissingValFillColor = "Gray"
    res@mpMaxLatF=40.125
    res@mpMinLatF=25.875
    res@mpMaxLonF=104.125
    res@mpMinLonF=74.875      
    res@mpDataSetName = "Earth..4"  
    res@mpDataBaseVersion = "MediumRes"
    res@mpOutlineBoundarySets = "National"         
    res@mpNationalLineThicknessF = 2
    res@mpOutlineOn = True        
    res@mpLandFillColor= "white"  
    res@mpInlandWaterFillColor="white"         
    res@cnFillOn = True
    res@cnLinesOn = False
    res@cnLineLabelsOn = False
    res@cnInfoLabelOn = False   
    res@cnLevelSelectionMode = "ExplicitLevels"
    res@lbLabelBarOn = True   
        res@mpDataBaseVersion     = "MediumRes"
    res@mpPerimOn =False
       
       
       
    res@gsnCenterString = "(a)"
    plot0 = gsn_csm_contour_map(wks, rs(0,:,:),res)
    res@gsnCenterString = "(b)"
    plot1 = gsn_csm_contour_map(wks, rs(1,:,:),res)


    resmp1                   =   True
    resmp1@gsLineColor       =   "black"
    resmp1@gsLineThicknessF  =   2.0


    resmp2                  =   True
    resmp2@gsLineColor      =   "black"
    resmp2@gsLineThicknessF =   3.0
    res@gsnCenterString = "Sensible heat net flux"



    plotmp1 = gsn_add_shapefile_polylines(wks, plot0, "shapefile/China_GuoJieShengJie_Polyline.shp", resmp1)
    plotmp2 = gsn_add_shapefile_polylines(wks, plot0, "shapefile/Qinghai-Tibet Plateau.shp", resmp2)

    plotmp11 = gsn_add_shapefile_polylines(wks, plot1, "shapefile/China_GuoJieShengJie_Polyline.shp", resmp1)
    plotmp22 = gsn_add_shapefile_polylines(wks, plot1, "shapefile/Qinghai-Tibet Plateau.shp", resmp2)


    res_panel = True
    panel_dims = (/1, 2/)
    res_panel@gsnMaximize = True
    gsn_panel(wks, (/plot0, plot1/), panel_dims, res_panel)
end

墒以光年

进此页面学习
http://www.ncl.ucar.edu/Applications/mask.shtml

卑微江小澄

墒以光年 发表于 2020-8-6 16:26
进此页面学习
http://www.ncl.ucar.edu/Applications/mask.shtml

我就是在这里学不明白才来提问的哭了ORZ

好运人

青藏高原那一块，还是西藏一个省？