<div dir="ltr"><div dir="ltr"><div dir="ltr">Hi Rick,<div>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?</div><div><br></div><div><div>   plot1 =  gsn_csm_contour_map(wks,rain,res)</div><div><br></div><div>  res@cnMonoFillPattern = True<br></div><div>  res@cnFillPattern = 12</div><div>  rain = where(<a href="http://rain.lt/" target="_blank">rain.lt</a>. .39 .and. <a href="http://rain.gt/" target="_blank">rain.gt</a>. -.39, rain@_FillValue, rain)<br></div><div>  plot2 =  gsn_csm_contour_map(wks,rain,res)<br></div><div>  overlay(plot1,plot2)<br></div><div><br></div><div>  draw(plot1)</div><div>  frame(wks)</div></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, May 2, 2019 at 6:03 PM Rick Brownrigg <<a href="mailto:brownrig@ucar.edu">brownrig@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><br><div>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 ?</div><div><br></div><div>This is speculation:</div><div><br></div><div>i) draw your color-filled plot as before</div><div><br></div><div>ii) mask out the rain variable for only the values you want:</div><div><br></div><div>   rain = where(<a href="http://rain.lt" target="_blank">rain.lt</a>. .39 .and. <a href="http://rain.gt" target="_blank">rain.gt</a>. -.39, rain@_FillValue, rain)</div><div><br></div><div>iii) redraw the rain variable with gsn_csm_contour_map as before, but using cnFillPatterns instead of color</div><div><br></div><div>iv) overlay the two plots as in the coneff_14 example (this requires you set res@gsnFrame = False and res@gsnDraw = False)</div><div><br></div><div>HTH...</div><div>Rick<br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, May 2, 2019 at 10:30 AM S Br <<a href="mailto:sbr.climate@gmail.com" target="_blank">sbr.climate@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi Rick,<div>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(). </div><div>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). </div><div>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.</div><div>'gsn_contour_shade' is working fine for this but unfortunately, I can't use this with gsn_csm_contour_map()</div><div><br></div><div>Is there any otherway to achieve this?</div><div><br></div><div>Thanks.</div><div>SB</div><div><br></div><div><br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, May 2, 2019 at 4:45 PM Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div>The docs for gsn_contour_shade() state:</div><div><br></div><div>You cannot mix and match color and shading with this routine.
For an example of drawing shaded and filled contours, see
example <a href="http://ncl.ucar.edu/Applications/coneff.shtml#ex14" target="_blank">coneff_14.ncl</a>. <br></div><div><br></div><div>However, that example coneff_14 shows how to achieve the result via two calls to gsn_csm_contour():</div><div><br></div><div><a href="http://ncl.ucar.edu/Applications/coneff.shtml#ex14" target="_blank">http://ncl.ucar.edu/Applications/coneff.shtml#ex14</a></div><div><br></div><div>Hope that helps...</div><div>Rick<br></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, May 2, 2019 at 5:39 AM S Br <<a href="mailto:sbr.climate@gmail.com" target="_blank">sbr.climate@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr">Hi Users,<div>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)</div><div><br></div><div>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.</div><div><div> I am using NCL6.4.0 and my script is given below.</div></div><div><br></div><div>Thank you, SB</div><div><br></div><div><div>begin</div><div><br></div><div> in1 = addfile("<a href="http://precip.gpcp.mon.mean.nc" target="_blank">precip.gpcp.mon.mean.nc</a>","r")</div><div>  rain=in1->precip(0,:,:)</div><div>  printVarSummary(rain)</div><div>;************************************</div><div>;create plot</div><div>;************************************</div><div>  wks = gsn_open_wks("pdf","overlay")             ; send graphics to PNG file</div><div>  res                     = True<br></div><div>  res@gsnDraw           = False         ; Don't draw plots</div><div>  res@gsnFrame          = False         ; Don't advance frame</div><div>  res@cnFillOn          = True          ; Turn on contour fill</div><div>  res@cnLinesOn         = False         ; Turn off contour lines</div><div>  res@cnFillPalette       = "MPL_YlOrRd"   ; set color map</div><div>  res@tiMainString        = "Overlaying shaded contours on filled contours"</div><div>;  res@cnFillMode = "RasterFill"</div><div>  res@cnLevelSelectionMode =  "ManualLevels"</div><div>  res@cnMinLevelValF       = 0.0</div><div>  res@cnMaxLevelValF       = 5.0</div><div>  res@cnLevelSpacingF      = 0.5</div><div><br></div><div>;   res@gsnAddCyclic     = False</div><div>   res@mpMinLatF    =   -10.0               ; only plot 30S to 30N</div><div>   res@mpMaxLatF    =  54.0</div><div>   res@mpMinLonF    = 30.0              ; only plot 30S to 30N</div><div>   res@mpMaxLonF    = 100.0</div><div><br></div><div>  plot = gsn_csm_contour_map(wks,rain,res)</div><div><br></div><div>  opt = True</div><div>  opt@gsnShadeFillType = "pattern"      ; pattern fill</div><div>  opt@gsnShadeHigh = 2                  ; use pattern #2</div><div>  opt@gsnShadeLow = 17                  ; use pattern #17</div><div>  plot2 = gsn_contour_shade(plot,0.5,3.0,opt)</div><div>  overlay(plot,plot2)</div><div><br></div><div>  draw(plot)</div><div>  frame(wks)</div><div>end</div></div></div></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
</blockquote></div>
</blockquote></div>
</blockquote></div>
</blockquote></div>