contourline color?

Hi,
Temperature line want to draw red color.But , it is black color.
* This is my code :
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
begin
; DATA OPEN  year = getenv("YY")  mon = getenv("MM")  day = getenv("DD")  hour = getenv("HH")
  date = getenv("DATE")
  DATADir = "../DATA/"  in_file = addfile(DATADir + "wrfout_d01_"+year+"-"+mon+"-"+day+"_"+hour+":00:00.nc", "r")
; WORKSTATION TYPE  type = "x11";  type = "ps"
; time dimension  times = wrf_user_list_times(in_file)  ntimes = dimsizes(times)
  do it = 0, 0;ntimes-1,6
; workstation open    if (it .lt. 10)      vtime = "0"+it    else      vtime = it    end if
    wks = gsn_open_wks(type,"D1_500_"+date+"+"+vtime+"_basic")    gsn_define_colormap(wks,"Rainbow")
; geopotential height [gpm]    PRESS = wrf_user_getvar(in_file,"pressure",it)    height = wrf_user_getvar(in_file,"z",it)    geo_h = wrf_user_intrp3d(height,PRESS,"h",500,0.,False)    geo_h_at_units = "gpm"    geo_h_at_description = "Geopotential Height"    wrf_smooth_2d(geo_h,100)
; temperature [C] at 700 hPa    tc = wrf_user_getvar(in_file,"tc",it) ; temperatuer [C]    tc_plane = wrf_user_intrp3d(tc,PRESS,"h",500.,0.,False) ; temperature [C] at 700 hPa    tc_plane_at_units = "C"    tc_plane_at_description = "Temperature"    wrf_smooth_2d(tc_plane,100)
; absolute vorticity    avo = wrf_user_getvar(in_file,"avo",it)    av_plane = wrf_user_intrp3d(avo,PRESS,"h",500.,0,False)
; vector    u_wind = wrf_user_getvar(in_file,"ua",it)    v_wind = wrf_user_getvar(in_file,"va",it)    u_plane = wrf_user_intrp3d(u_wind,PRESS,"h",500,0.,False)    v_plane = wrf_user_intrp3d(v_wind,PRESS,"h",500,0.,False)    wrf_smooth_2d(u_plane,100)    wrf_smooth_2d(v_plane,100)    u_plane             = u_plane*1.94386    v_plane             = v_plane*1.94386    u_plane_at_units       = "knots"    u_plane_at_description = "Vector";; Plot;    res = True    res_at_InitTime = True    res_at_Footer = False    res_at_TimeLabel = times(it)
    res_at_gsnMaximize = True
; geopotential height    res2 = res    res2_at_cnMonoLineColor = True    res2_at_cnLineColor = "blue"    res2_at_cnLineThicknessF = 2.2    res2_at_cnInfoLabelOn = False    res2_at_cnLineLabelsOn = True    res2_at_cnLineLabelBackgroundColor = "white"    res2_at_cnLevelSpacingF = 60.
    contour_geoh = wrf_contour(in_file,wks,geo_h,res2)
; temperature [C]    res3 = res    res2_at_cnMonoLineColor = True    res3_at_cnLineColors = "red"    res3_at_cnLineThicknessF = 2.2    res3_at_cnInfoLabelOn = False    res3_at_cnLineLabelBackgroundColor = "white"    res3_at_cnLevelSpacingF = 3.    res3_at_gsnContourPosLineDashPattern = 1
    contour_tc = wrf_contour(in_file,wks,tc_plane,res3)
; vector    res4                       = res    res4_at_vcLevelSelectionMode = "ManualLevels"    res4_at_vcGlyphStyle          = "WindBarb"    res4_at_vcRefMagnitudeF       = 10.    res4_at_vcRefLengthF          = 0.025    res4_at_vcRefAnnoOn           = True    res4_at_vcMinDistanceF        = 0.025    res4_at_vcWindBarbTickLengthF = 0.4    res4_at_vcWindBarbColor       = "black"
    vector                     = wrf_vector(in_file,wks,u_plane,v_plane,res4)
; vorticity    res5 = res    res5_at_cnFillOn = True    res5_at_cnLevelSpacingF = 5
    contour_vrt = wrf_contour(in_file,wks,av_plane,res5)
;    pltres                       = True    mpres                        = True    mpres_at_mpGeophysicalLineColor = "black"
    plot = wrf_map_overlays(in_file,wks,(/contour_vrt,vector,contour_geoh,contour_tc/),pltres,mpres)
  end do
end

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Sep 15 2009 - 18:58:32 MDT