Re: contourline color?

Please try change this paragraph:

; 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)

To:

; temperature [C]
     res3 = res
     res3_at_cnMonoLineColor = True ;This line changed res2 -> res3
     res3_at_cnLineColor = "red" ; This line changed cnLineColors ->
cnLineColor
     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)

Wei Huang
huangwei_at_ucar.edu
VETS/CISL
National Center for Atmospheric Research
P.O. Box 3000 (1850 Table Mesa Dr.)
Boulder, CO 80307-3000 USA
(303) 497-8924

On Sep 15, 2009, at 6:58 PM, 감은아 wrote:

> 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

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Sep 16 2009 - 09:48:50 MDT