NCL作图问题

放逐流年

论坛里的各位大神，现在我想根据wrfout画下垫面土地利用数据，怎么根据wrfout文件中的LU_INDEX绘图呢？主要是怎么将不同的地类绘制成自己设置的颜色，用来区分不同的地类，而不是这种级别类的呢？

放逐流年

f   = addfile ("/myarray/DuHY/reachwrfout/usgsdiy20153/wrfout/wrfout_d03_2015-07-02_03:00:00", "r")
lu = wrf_user_getvar(f,"LU_INDEX",0)
wks = gsn_open_wks("x11" ,"landuse")
colors = (/"white","black","Red","Green","PaleGreen","lightseagreen","LightSkyBlue","Yellow","Blue"/)  
gsn_define_colormap(wks, colors)

res                  = True
res@tfDoNDCOverlay   = True
res@gsnDraw          = False
res@gsnFrame         = False
res@gsnRightString   = ""
res@gsnLeftString    = ""

fres                    = res
fres@cnFillOn           = True
fres@cnLinesOn          = False
fres@cnInfoLabelOn      = False
fres@cnFillMode="RasterFill"
fres@cnLevelSelectionMode = "ExplicitLevels"
;fres@cnExplicitLabelBarLabelsOn = True
;fres@ContourParameters = (/"city","crop","grass","tree","water","sand","bare"/)
;fres@ContourParameters = (/1,3,7,15,16,19,24/)
fres@cnLevels             = (/1,3,7,15,16,19,24/)
fres@cnLineThicknessF   = 2.0
fres@lbOrientation      = "Vertical"
fres@lbLabelFontHeightF = 0.01
fcontour                = gsn_csm_contour(wks,lu,fres)

mres = res
mres = wrf_map_resources(f,mres)    ; Add necessary resources for WRF map
map = gsn_csm_map(wks,mres)

overlay(map,fcontour)

draw(map)
frame(wks)
这是我用于绘图的代码

放逐流年

意思就是说怎么用户指定颜色呢？比如说我想让LU_INDEX=1的颜色为红色