Re: ShadeLtContour not working!!

Thanks Adam,
now it works :)

-- AF

On Fri, Jun 28, 2013 at 3:19 AM, Adam Phillips <asphilli@ucar.edu> wrote:

> Hi Ahmad,
> My apologies, I forgot to add that you also need to tell NCL to draw the
> plot and advance the frame manually if you set gsnDraw and gsnFrame =
> False.
>
> So:
> res = True
>
> res@gsnDraw = False
> res@gsnFrame = False
> ....
>
> res@cnLevelSelectionMode = "ExplicitLevels"
> res@cnLevels = ispan(-10,10,1)
>
> plot = gsn_csm_contour(wks,T(1000,:,:),res)
> plot = ShadeLtContour(plot, 1.01, 3)
> draw(plot)
> frame(wks)
> end
>
> If that doesn't work let ncl-talk know.
> Adam
>
>
> On 06/27/2013 12:15 PM, Ahmad Farsyud wrote:
>
> Thank you Adam,
>
> But still not working !!
>
> 1. If I set gsnDraw and gsnFrame to 'False" i got a broken document.
> 2. If I just set gsnDraw to False, I got blank
> 3. I tried several ways to set the number which the contour will be shaded
> but still there is no effect.
>
> Please help me,
>
>
>
> 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/contributed.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
>
> begin
> f = addfile("ly.nc","r") ; Read in first time step and
> T = f->ly(:,:,:)
>
>
>
> T = where(T.lt.0, 0, T) ; make values less than 0 to be zero.
>
>
>
>
> ;*******************************************************************************
>
> wks = gsn_open_wks ("pdf", "annualzonalmean" ) ; open
> workstation
>
> res = True ; plot mods desired
> ;res@gsnDraw = False
> ;res@gsnFrame = False
>
> res@cnLineLabelsOn = True ; turn on line labels
> res@cnFillOn = False ; turn on color fill
> res@gsnLeftString = " "
> res@gsnRightString = " "
> res@gsnCenterString = " "
> res@gsnContourZeroLineThicknessF = 0
> ;res@cnLevelSpacingF = 1
> res@cnInfoLabelOn = False
>
> res@cnLevelSelectionMode = "ExplicitLevels"
> res@cnLevels = ispan(-10,10,1)
>
>
> plot = gsn_csm_contour(wks,T(1000,:,:),res)
> ;plot = ShadeLtContour(plot, 1, 3)
> ;plot = ShadeLtContour(plot, 1.51, 3)
> plot = ShadeLtContour(plot, 1.01, 3)
>
>
> end
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
> On Fri, Jun 28, 2013 at 1:22 AM, Adam Phillips <asphilli@ucar.edu> wrote:
>
>> Hi Ahmad,
>> Without knowing what NCL is picking as your contour levels, this is a bit
>> difficult to answer.
>> Note that in the ShadeLtContour documentation it states:
>> "Note that the shading begins at the first contour level less than *value
>> *, and not at the specified *value*. Check your plot to make sure that
>> the results are what you expected. "
>>
>> Again, not knowing what your contours are makes this a guess, but try
>> setting this:
>> plot = ShadeLtContour(plot, 1.51, 3)
>>
>> That will result in all areas less than the first contour less than 1.51
>> to be shaded.. If you are not happy with the results, try setting:
>> res@cnLevelSelectionMode = "ExplicitLevels"
>> res@cnLevels = ispan(-10,10,1)
>>
>> plot = gsn_csm_contour(wks,T(1000,:,:),res)
>> plot = ShadeLtContour(plot, 1.01, 3)
>>
>> Finally, and perhaps most importantly, you need to tell NCL to not draw
>> the plot upon the call to gsn_csm_contour and to not advance the frame, so
>> make sure you include the following lines no matter what:
>> res@gsnDraw = False
>> res@gsnFrame = False
>> In fact, you might want to try adding those two lines first to see if it
>> fixes the problem.
>>
>> If the above doesn't help, please let ncl-talk know.
>> Adam
>>
>>
>>
>> On 06/27/2013 08:48 AM, Ahmad Farsyud wrote:
>>
>> Dear NCL,
>>
>> I just wonder why the "ShadeLtContour" function in my code
>> is not working. I want to shade my plot if the contour less than 1.
>>
>> Here is my plot:
>>
>> 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/contributed.ncl"
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
>>
>> begin
>> f = addfile("ly.nc","r") ; Read in first time step and
>> T = f->ly(:,:,:)
>>
>> ;T = T(lev|:,time|:,lat|:)
>>
>> T = where(T.lt.0, 0, T)
>>
>>
>> ;printVarSummary(TEMP)
>>
>>
>>
>> ;*******************************************************************************
>>
>> wks = gsn_open_wks ("x11", "annualzonalmean" ) ; open
>> workstation
>>
>> res = True ; plot mods desired
>> res@cnLineLabelsOn = True ; turn on line labels
>> res@cnFillOn = False ; turn on color fill
>> res@cnMonoFillColor = True
>> res@gsnLeftString = " "
>> res@gsnRightString = " "
>> res@gsnCenterString = " "
>> res@gsnContourZeroLineThicknessF = 0
>> res@cnLevelSpacingF = 1
>> res@cnMonoFillPattern = True
>> res@cnInfoLabelOn = False
>>
>>
>> ;*******************************************************************************
>> ;changing aspect ratio
>>
>> ;res@vpXF = 0.13 ; change x-scale
>> ;res@vpWidthF = 0.75 ; change height and
>> width
>> ;res@vpHeightF = 0.45
>>
>>
>> plot = gsn_csm_contour(wks,T(1000,:,:),res)
>> plot = ShadeLtContour(plot, 1, 3)
>>
>> end
>>
>> thank you
>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>> --
>> ______________________________________________________________
>> Adam Phillips asphilli@ucar.edu
>> NCAR/Climate and Global Dynamics Division (303) 497-1726
>> P.O. Box 3000
>> Boulder, CO 80307-3000 http://www.cgd.ucar.edu/cas/asphilli
>>
>>
>
> --
> ______________________________________________________________
> Adam Phillips asphilli@ucar.edu
> NCAR/Climate and Global Dynamics Division (303) 497-1726
> P.O. Box 3000
> Boulder, CO 80307-3000 http://www.cgd.ucar.edu/cas/asphilli
>
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Jun 27 12:24:36 2013