绘制土地利用率

onenuister

将wrf输出结果中的土地利用类型可视化，默认20种类型(见下图)
a2 = addfile("./wrfdata_output/geo_em.d03.nc","r")
LAT = wrf_user_getvar(a2, "lat", 0)
LON = wrf_user_getvar(a2, "lon", 0)
land1 = a2->LU_INDEX
land2 = a2->LANDMASK
land3 = a2->LANDUSEF


ds = dimsizes(LAT)
nrow = ds(0)-1
ncol = ds(1)-1


wks = gsn_open_wks("png", "./landuse3")
res = True
res@gsnDraw = False
res@gsnFrame = False
res@mpGridAndLimbOn = True
res@mpLandFillColor            = "white"
res@mpInlandWaterFillColor     = "white"
res@mpOceanFillColor           = "white"
; res@mpOutlineBoundarySets   = "National"
res@mpOutlineOn = True
res@mpDataSetName="Earth..4"
res@mpOutlineSpecifiers = (/"China: states","Taiwan"/)
res@mpDataBaseVersion = "MediumRes"
res@mpFillDrawOrder = "PostDraw"
res@mpFillOn = True
res@mpFillAreaSpecifiers = (/"water","land"/)
res@mpSpecifiedFillColors = (/"white","white"/)
res@mpSpecifiedFillColors = (/0,0/)
res@mpAreaMaskingOn         = True
res@mpMaskAreaSpecifiers    = (/"China: states","Taiwan"/)
res@cnLinesOn = False
res@cnLineLabelsOn = False
res@cnInfoLabelOn = True
res@cnFillOn = True
res@cnFillPalette = "default"
res@cnLevelSpacingF = 1
res@cnLevelSelectionMode = "ManualLevels"
res@pmTickMarkDisplayMode = "Always"
res@lbLabelBarOn = True
res@mpProjection = "LambertConformal"
res@mpLambertMeridianF = 108
res@mpLambertParallel1F = 30
res@mpLambertParallel2F = 60
res@mpLimitMode = "Corners"
res@mpLeftCornerLatF = LAT(0,0)
res@mpLeftCornerLonF = LON(0,0)
res@mpRightCornerLatF = LAT(nrow,ncol)
res@mpRightCornerLonF = LON(nrow,ncol)
res@tmXTOn = False
res@tmYROn = False
res@gsnAddCyclic = False
res@tfDoNDCOverlay = True
setvalues NhlGetWorkspaceObjectId()
"wsMaximumSize" : 100000000000000
end setvalues
plot = gsn_csm_contour_map(wks, land1(0,:,:), res)
draw(plot)
frame(wks)
end

onenuister

绘制得到的结果