[ncl-talk] plot color fill problem
Chao Luo
chaoluo888 at gmail.com
Thu Dec 29 12:48:40 MST 2016
I am going to try it. Thanks!
On Thu, Dec 29, 2016 at 11:14 AM, David Brown <dbrown at ucar.edu> wrote:
> A small correction to Dennis' last message:
>
> The fill color of areas with missing data are controlled using the
> resource cnMissingValFillColor. By default this is set to "transparent"
> meaning that the background color will show up in these areas, but you can
> set it to any color you want.
> -dave
>
>
> On Wed, Dec 28, 2016 at 9:38 PM, Chao Luo <chaoluo888 at gmail.com> wrote:
>
>> Thanks, Dennis! This works.
>>
>> Chao
>>
>> On Wed, Dec 28, 2016 at 7:57 PM, Dennis Shea <shea at ucar.edu> wrote:
>>
>>> By default, 'gsn_csm_contour_map' uses 'gray' to fill land.
>>>
>>> See: http://www.ncl.ucar.edu/Applications/maponly.shtml
>>>
>>> To turn this off:
>>> res at mpFillOn = False
>>>
>>> ===
>>> The areas with no data (_FillValue) are set to the background color...
>>> white here
>>>
>>>
>>> On Wed, Dec 28, 2016 at 7:40 PM, Chao Luo <chaoluo888 at gmail.com> wrote:
>>>
>>>> Hi,
>>>>
>>>> I am using NCL script to plot map contour. The plot shows gray and
>>>> white color, which they are not defined in the script. I have tried some
>>>> map color fill options, it is still problem. Any suggestions and comments
>>>> are very appreciated.
>>>>
>>>> Thanks,
>>>>
>>>> Chao
>>>>
>>>> attachments are plot and script
>>>>
>>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>>>> ;load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
>>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>>>>
>>>> procedure metaUtil (z[*][*][*]:numeric, time[*], lat[*], lon[*])
>>>>
>>>> begin
>>>> z!0 = "time" ; name dimensions
>>>> z!1 = "lat"
>>>> z!2 = "lon"
>>>> z&time = time ; assign coordinate variables
>>>> z&lat = lat
>>>> z&lon = lon
>>>> end
>>>> begin
>>>> f = addfile ("data/IC_TO_CG.hdf", "r")
>>>> lat = f->Latitude
>>>> lon = f->Longitude
>>>> time = f->Day_of_year
>>>> ; print(f)
>>>> IC2CG = f->IC_TO_CG(:,:,:) ; IC to CG ratio
>>>> IC2CG at _FillValue = -1
>>>>
>>>> ;;;;average data to annualy data
>>>> ; AvgTime_IC_TO_CG = dim_avg_Wrap(IC2CG(Latitude|:, Longitude|:,
>>>> Day_of_year|:) )
>>>>
>>>> AvgTime_July = dim_avg_Wrap(IC2CG(Latitude|:, Longitude|:,
>>>> Day_of_year|181:211) )
>>>> AvgTime_Aug = dim_avg_Wrap(IC2CG(Latitude|:, Longitude|:,
>>>> Day_of_year|212:242) )
>>>>
>>>> printVarSummary(AvgTime_July)
>>>>
>>>> ;;;;; read mask data
>>>> f = addfile ("/home/cluo/landmask.nc","r")
>>>> lat_mask = f->lat
>>>> lon_mask = f->lon
>>>> land = f->land(:,:,:)
>>>> printVarSummary(land)
>>>> mask1 = land
>>>> printVarSummary(mask1)do i=0, 360-1
>>>> ii = 360+i
>>>> do j=0, 360-1
>>>> jj = 359-j
>>>> mask1(0,j,i) = land(0,jj,ii)
>>>> end do
>>>> end do
>>>>
>>>> do i=360, 720-1
>>>> ii = i-360
>>>> do j=0, 360-1
>>>> jj = 359-j
>>>> mask1(0,j,i) = land(0,jj,ii)
>>>> end do
>>>> end do
>>>>
>>>> do i = 0, 720-1
>>>> do j = 1, 360-1
>>>> if (mask1(0,j,i).ne.0) then
>>>> AvgTime_July(j,i) = 0
>>>> end if
>>>> end do
>>>> end do
>>>> ;*************create plots*******************
>>>> wks = gsn_open_wks ("ps", "narr") ; open workstation
>>>> gsn_define_colormap(wks,"BlAqGrYeOrRe") ; choose colormap
>>>> ; gsn_define_colormap (wks,"gui_default") ; choose color map
>>>> plot = new(1,graphic)
>>>>
>>>> res = True ; plot mods desired
>>>> for original grid
>>>>
>>>> res at gsnDraw = False ; Don't draw plot or advance the
>>>> res at gsnFrame = False ; frame after plot is created.
>>>>
>>>> res at cnFillOn = True ; color fill
>>>> res at cnLinesOn = False ; no contour lines
>>>> res at gsnSpreadColors = True ; use total colormap
>>>> res at gsnSpreadColorStart = 4
>>>> res at gsnSpreadColorEnd = -1
>>>> res at mpGridAndLimbOn = False
>>>> res at pmTickMarkDisplayMode = "Always" ; turn on tickmarks
>>>> res at tmXTOn = False
>>>> res at gsnAddCyclic = False ; regional data
>>>> res at mpOutlineOn = True ; turn on map outline
>>>> res at mpOutlineBoundarySets = "USStates" ; turn on state
>>>> boundaries
>>>> ; res at mpOutlineBoundarySets = "National"
>>>>
>>>>
>>>> ;;;; set boundary for subregion plot
>>>> res at mpLimitMode = "LatLon"
>>>> res at mpMinLonF = -125
>>>> res at mpMaxLonF = -70
>>>> res at mpMinLatF = 25
>>>> res at mpMaxLatF = 50
>>>> res at mpPerimOn = True
>>>> res at mpCenterLonF =-97
>>>> res at mpCenterLatF = 35
>>>> res at gsnCenterString = "Monthly mean IC/CG, August"
>>>> ; draw center subtitle
>>>> res at gsnCenterStringFontHeightF=25
>>>> res at gsnLeftString = "" ; draw left
>>>> subtitle
>>>> res at cnLevelSelectionMode = "ExplicitLevels" ; set explicit contour
>>>> levels
>>>> res at cnLevels = (/1,2.,3.,4.,5.,6.,8.,10./)
>>>> res at cnFillColors = (/10,20,30,40,50,65,75,85,96/)
>>>> res at gsnCenterString = "Climatology IC/CG, averaged in July"
>>>> plot(0) = gsn_csm_contour_map(wks,AvgTime_July(:,:),res)
>>>>
>>>> resP = True
>>>> resP at gsnPanelYWhiteSpacePercent = 5
>>>> resP at gsnPanelXWhiteSpacePercent = 5
>>>> resP at cnFillColors = (/5,10,20,30,40,50,65,75,85,96/)
>>>> resP at gsnMaximize = True
>>>> ; resP at gsnPanelLabelBar = True
>>>> resP at lbLabelBarOn = True
>>>> resP at gsnPaperOrientation = "portrait"
>>>> ; resP at lbOrientation = "vertical" ; vertical label bar
>>>> resP at lbOrientation = "horizontal"
>>>> resP at lbLabelAutoStride = True ; nice label bar labels
>>>> resP at lbTitleOn = True
>>>> resP at lbLabelFont = "helvetica"
>>>> ; resP at lbTitleString = "Sig lvl"
>>>> resP at lbTitlePosition = "Bottom"
>>>> resP at lbTitleFontHeightF = .022
>>>> resP at lbTitleDirection = "Across"
>>>> ; resP at txString = title
>>>> gsn_panel(wks,plot,(/1,1/),resP)
>>>> end
>>>>
>>>>
>>>> _______________________________________________
>>>> ncl-talk mailing list
>>>> ncl-talk at ucar.edu
>>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>>>
>>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20161229/2e5de7d9/attachment.html
More information about the ncl-talk
mailing list