[ncl-talk] Color filled contour disappears with contour_shade pattern

S Br sbr.climate at gmail.com
Fri May 3 00:15:59 MDT 2019 

Hi Adam,
It still doesn't work because I am using a regional domain. It shows the
following error.

warning:mpMinLonF is not a valid resource in correlation_NCL_contour at
this time
warning:mpMaxLonF is not a valid resource in correlation_NCL_contour at
this time
warning:mpMinLatF is not a valid resource in correlation_NCL_contour at
this time
warning:mpMaxLatF is not a valid resource in correlation_NCL_contour at
this time
warning:mpFillOn is not a valid resource in correlation_NCL_contour at this
time

Is there any other way, I can overlay the 2nd value with the patterns.

Thank you.
SB

On Thu, May 2, 2019 at 7:43 PM Adam Phillips <asphilli at ucar.edu> wrote:

> Hi SB,
> The plot that is going to be overlaid cannot be created by calling
> gsn_csn_contour_map. You should use gsn_csm_contour instead. So switch this:
>   plot2 =  gsn_csm_contour_map(wks,rain,res)
> with this:
>   plot2 =  gsn_csm_contour(wks,rain,res)
> Adam
>
>
> On Thu, May 2, 2019 at 11:59 AM S Br <sbr.climate at gmail.com> wrote:
>
>> Hi Rick,
>> You have correctly guessed my query. Now, I have checked that the new
>> variable 'rain' is correctly defined with the value gt/lt .39/-.39 but I am
>> still struggling to overlay this value with patterns with
>> another gsn_csm_contour_map. I am copying here the last part of my program.
>> Am I doing somewhere wrong? Else, do I need to send a second resources
>> (res2) for the 2nd plot?
>>
>>    plot1 =  gsn_csm_contour_map(wks,rain,res)
>>
>>   res at cnMonoFillPattern = True
>>   res at cnFillPattern = 12
>>   rain = where(rain.lt. .39 .and. rain.gt. -.39, rain at _FillValue, rain)
>>   plot2 =  gsn_csm_contour_map(wks,rain,res)
>>   overlay(plot1,plot2)
>>
>>   draw(plot1)
>>   frame(wks)
>>
>>
>>
>>
>>
>>
>>
>> On Thu, May 2, 2019 at 6:03 PM Rick Brownrigg <brownrig at ucar.edu> wrote:
>>
>>>
>>> So my (very limited) understanding from the code you sent, you drew a
>>> color-filled contour plot of the variable "rain", which was a time slice of
>>> the variable precip from the file. Then you intended to use gsn_csm_shade
>>> to hachure regions where rain was gt/lt .39/-.39 ?
>>>
>>> This is speculation:
>>>
>>> i) draw your color-filled plot as before
>>>
>>> ii) mask out the rain variable for only the values you want:
>>>
>>>    rain = where(rain.lt. .39 .and. rain.gt. -.39, rain at _FillValue, rain)
>>>
>>> iii) redraw the rain variable with gsn_csm_contour_map as before, but
>>> using cnFillPatterns instead of color
>>>
>>> iv) overlay the two plots as in the coneff_14 example (this requires you
>>> set res at gsnFrame = False and res at gsnDraw = False)
>>>
>>> HTH...
>>> Rick
>>>
>>>
>>>
>>> On Thu, May 2, 2019 at 10:30 AM S Br <sbr.climate at gmail.com> wrote:
>>>
>>>> Hi Rick,
>>>> Thanks for your reply but my region of interest cover a regional
>>>> lat/lon domain, so I can't use gsn_csm_contour(). Instead, I am using
>>>> gsn_csm_contour_map().
>>>> Actually, my main objective is to draw a correlation plot with contour
>>>> filled with color and then mark pattern to the region having significant
>>>> correlation value (see the map attached).
>>>> For example, I have calculated the threshold value of correlation to be
>>>> significant is 0.39. So, I want to draw patterns over the filled contour
>>>> for the region having the correlation value below -0.39 and above 0.39.
>>>> 'gsn_contour_shade' is working fine for this but unfortunately, I can't
>>>> use this with gsn_csm_contour_map()
>>>>
>>>> Is there any otherway to achieve this?
>>>>
>>>> Thanks.
>>>> SB
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> On Thu, May 2, 2019 at 4:45 PM Rick Brownrigg <brownrig at ucar.edu>
>>>> wrote:
>>>>
>>>>> The docs for gsn_contour_shade() state:
>>>>>
>>>>> You cannot mix and match color and shading with this routine. For an
>>>>> example of drawing shaded and filled contours, see example
>>>>> coneff_14.ncl <http://ncl.ucar.edu/Applications/coneff.shtml#ex14>.
>>>>>
>>>>> However, that example coneff_14 shows how to achieve the result via
>>>>> two calls to gsn_csm_contour():
>>>>>
>>>>> http://ncl.ucar.edu/Applications/coneff.shtml#ex14
>>>>>
>>>>> Hope that helps...
>>>>> Rick
>>>>>
>>>>> On Thu, May 2, 2019 at 5:39 AM S Br <sbr.climate at gmail.com> wrote:
>>>>>
>>>>>> Hi Users,
>>>>>> I am trying to plot a contour map filled with color and then apply
>>>>>> patterns to the values greater and smaller than the specified values using
>>>>>> the function gsn_contour_shade(plot,0.5,3.0,opt)
>>>>>>
>>>>>> However, when I use the gsn_contour_shade and overlay, the contour
>>>>>> filled with color disappears. I can only see the patterns without the
>>>>>> filled contour. Please have a look at the attached image.
>>>>>>  I am using NCL6.4.0 and my script is given below.
>>>>>>
>>>>>> Thank you, SB
>>>>>>
>>>>>> begin
>>>>>>
>>>>>>  in1 = addfile("precip.gpcp.mon.mean.nc","r")
>>>>>>   rain=in1->precip(0,:,:)
>>>>>>   printVarSummary(rain)
>>>>>> ;************************************
>>>>>> ;create plot
>>>>>> ;************************************
>>>>>>   wks = gsn_open_wks("pdf","overlay")             ; send graphics to
>>>>>> PNG file
>>>>>>   res                     = True
>>>>>>   res at gsnDraw           = False         ; Don't draw plots
>>>>>>   res at gsnFrame          = False         ; Don't advance frame
>>>>>>   res at cnFillOn          = True          ; Turn on contour fill
>>>>>>   res at cnLinesOn         = False         ; Turn off contour lines
>>>>>>   res at cnFillPalette       = "MPL_YlOrRd"   ; set color map
>>>>>>   res at tiMainString        = "Overlaying shaded contours on filled
>>>>>> contours"
>>>>>> ;  res at cnFillMode = "RasterFill"
>>>>>>   res at cnLevelSelectionMode =  "ManualLevels"
>>>>>>   res at cnMinLevelValF       = 0.0
>>>>>>   res at cnMaxLevelValF       = 5.0
>>>>>>   res at cnLevelSpacingF      = 0.5
>>>>>>
>>>>>> ;   res at gsnAddCyclic     = False
>>>>>>    res at mpMinLatF    =   -10.0               ; only plot 30S to 30N
>>>>>>    res at mpMaxLatF    =  54.0
>>>>>>    res at mpMinLonF    = 30.0              ; only plot 30S to 30N
>>>>>>    res at mpMaxLonF    = 100.0
>>>>>>
>>>>>>   plot = gsn_csm_contour_map(wks,rain,res)
>>>>>>
>>>>>>   opt = True
>>>>>>   opt at gsnShadeFillType = "pattern"      ; pattern fill
>>>>>>   opt at gsnShadeHigh = 2                  ; use pattern #2
>>>>>>   opt at gsnShadeLow = 17                  ; use pattern #17
>>>>>>   plot2 = gsn_contour_shade(plot,0.5,3.0,opt)
>>>>>>   overlay(plot,plot2)
>>>>>>
>>>>>>   draw(plot)
>>>>>>   frame(wks)
>>>>>> 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
>>
>
>
> --
> Adam Phillips
> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
> www.cgd.ucar.edu/staff/asphilli/   303-497-1726
>
> <http://www.cgd.ucar.edu/staff/asphilli>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190503/efeeec32/attachment.html>

More information about the ncl-talk mailing list