<div dir="ltr">Hi Liz,<div>In spite of the error message you were getting (which does need to be fixed), the signal/noise data is being plotted. There are several reasons it is not visible to you: 1) you are using the same color table as for the underlying data so it is very hard to distinguish; 2) the default dot size is very small -- you need to make the dots bigger; and 3) the AreaFill draw method requires continuous regions of non-missing data in order to draw -- it will not draw isolated points of good data in a sea of missing values. Also note that the second and third plots have no valid s/n data.</div><div><br></div><div>I have modified your script to set a single fill color (red -- so it is visible against both the predominately blue underlying plot and also against the black map outlines). I also made the dots bigger and set them closer together. And I turned off the constant fill message. Of course, these are all things you can tweak.</div><div><br></div><div>Further, my test script draws each panel member separately after the panel plot, so that you can see how more of the s/n data is visible if you make the plot bigger. Also I have a commented out a setting to change the FillMode to "cellfill". If you try that you will see the isolated points that have been eliminated when you use "areafill". Hope this helps give some perspective on the trade-offs involved in choosing how to plot your data.</div><div><br></div><div>I will create a JIRA ticket concerning the error message. I am attaching the modified version of your script.</div><div><br></div><div> -dave</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Sep 11, 2014 at 4:22 PM, Elizabeth Burakowski <span dir="ltr"><<a href="mailto:burakows@ucar.edu" target="_blank">burakows@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Thanks David. <div><br></div><div>Liz.</div><div><br></div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Sep 11, 2014 at 2:12 PM, David Brown <span dir="ltr"><<a href="mailto:dbrown@ucar.edu" target="_blank">dbrown@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Liz,<div>I am looking into this now. I will let you know what I find.</div><div> -dave</div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div>On Wed, Sep 10, 2014 at 4:09 PM, Elizabeth Burakowski <span dir="ltr"><<a href="mailto:burakows@ucar.edu" target="_blank">burakows@ucar.edu</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><div dir="ltr">Hi NCL-Talk, <div><br></div><div>I am trying to use cnConstFEnableFill to stiple regions in a constant fill color/pattern where signal/noise ratios are greater than two. However, I keep getting an error that says: </div><div><br></div><div>warning:ContourPlotInitialize: scalar field is constant; no contour lines will appear; use cnConstFEnableFill to enable fill</div><div><br></div><div>The plotting resources for my contour plot (res) and my signal:noise ratio plot (snres) using cnConstFEnableFill is below. The plot turns out perfectly fine, minus the desired stipling. </div><div><br></div><div>Is there a specific order I must put cnConstFEnableFill in my script that I am missing? The full script is pretty long, so I've snipped out the bulk of it and copied the plot resources information below. If the full script is desired, it is available here:</div><div><br></div><div>/glade/p/work/burakows/plot/NCL/LC_Diff/biascorr_diff_PD_DF_mult_mean_panel.ncl<br></div><div><br></div><div>Thanks, </div><div><br></div><div>Liz.</div><div><br></div><div><br></div><div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl"</div><div><br></div><div>;==========================================================</div><div>; The main code</div><div>;==========================================================</div><div><br></div><div>begin</div><div><br></div><div><snip> </div><div><br></div><div><br></div><div>;----------------------------------------------------------------------</div><div>; Plotting options section</div><div>;----------------------------------------------------------------------</div><div><br></div><div> pltType = "ps" ; plot type</div><div> pltDir = "./" ; plot directory</div><div> pltName = "Panel_biascorr_"+varw+"PD-DF_2011" <span style="white-space:pre-wrap"> </span>; plot name (ps file)</div><div> pltPath = pltDir+pltName ; plot path</div><div><br></div><div> wks = gsn_open_wks(pltType,pltPath) ; create workstation for ps file</div><div> </div><div> gsn_define_colormap(wks,"temp_diff_18lev")</div><div><br></div><div><span style="white-space:pre-wrap"> </span>res = True</div><div><br></div><div><span style="white-space:pre-wrap"> </span>res@gsnDraw<span style="white-space:pre-wrap"> </span>= False</div><div><span style="white-space:pre-wrap"> </span>res@gsnFrame<span style="white-space:pre-wrap"> </span>= False</div><div><br></div><div><span style="white-space:pre-wrap"> </span>res@cnFillOn = True ; color plot desired</div><div> <span style="white-space:pre-wrap"> </span>res@cnLinesOn = False ; turn off contour lines</div><div> <span style="white-space:pre-wrap"> </span>res@cnLineLabelsOn = False ; turn off contour labels</div><div> <span style="white-space:pre-wrap"> </span>res@cnInfoLabelOn<span style="white-space:pre-wrap"> </span> = False<span style="white-space:pre-wrap"> </span> ; turn off info label (top labels of indvid. plots)</div><div><span style="white-space:pre-wrap"> </span>res@cnFillMode = "RasterFill" ; turn raster on</div><div><span style="white-space:pre-wrap"> </span>res@cnLevelSelectionMode = "ManualLevels" ; Set contour levels manually</div><div><span style="white-space:pre-wrap"> </span>res@cnMinLevelValF = -4 ; minimum contour (degrees C, or mm)</div><div> res@cnMaxLevelValF = 4 ; maximum contour (degrees C, or mm)</div><div> res@cnLevelSpacingF = 0.5 ; contour interval</div><div><br></div><div> <span style="white-space:pre-wrap"> </span>res@lbLabelBarOn = False ; Will turn on in panel later</div><div><span style="white-space:pre-wrap"> </span>res@lbOrientation<span style="white-space:pre-wrap"> </span> = "Horizontal"<span style="white-space:pre-wrap"> </span> ; Horizontal label bar</div><div><br></div><div> <span style="white-space:pre-wrap"> </span>res@mpFillOn = False</div><div> <span style="white-space:pre-wrap"> </span>res@mpOutlineOn = True</div><div> <span style="white-space:pre-wrap"> </span>res@mpOutlineBoundarySets = "AllBoundaries"</div><div> <span style="white-space:pre-wrap"> </span>res@mpProjection = "CylindricalEquidistant"</div><div><br></div><div> <span style="white-space:pre-wrap"> </span>res@mpLimitMode = "LatLon" ; required</div><div> <span style="white-space:pre-wrap"> </span>res@mpMinLatF = minlat-bndadd</div><div> <span style="white-space:pre-wrap"> </span>res@mpMaxLatF = maxlat+bndadd</div><div> <span style="white-space:pre-wrap"> </span>res@mpMinLonF = minlon-bndadd</div><div> <span style="white-space:pre-wrap"> </span>res@mpMaxLonF = maxlon+bndadd</div><div> <span style="white-space:pre-wrap"> </span>res@mpCenterLonF = (minlon + maxlon)*0.5</div><div> <span style="white-space:pre-wrap"> </span>res@mpCenterLatF = (minlat + maxlat)*0.5</div><div><span style="white-space:pre-wrap"> </span></div><div><span style="white-space:pre-wrap"> </span>res@gsnLeftString = "" ; Turn off left subtitle</div><div> res@gsnRightString = "" ; Turn off right subtitle</div><div> <span style="white-space:pre-wrap"> </span>res@gsnMajorLatSpacing = 2</div><div> <span style="white-space:pre-wrap"> </span>;res@gsnMajorLonSpacing = 2</div><div> <span style="white-space:pre-wrap"> </span>res@gsnMinorLonSpacing = 2</div><div><br></div><div><span style="white-space:pre-wrap"> </span>res@gsnAddCyclic<span style="white-space:pre-wrap"> </span>= False<span style="white-space:pre-wrap"> </span>; regional grid (changes central meridian)/xwo</div><div><span style="white-space:pre-wrap"> </span></div><div><span style="white-space:pre-wrap"> </span></div><div> ;----Resources for stipling of significant signal/noise (>2)</div><div> snres <span style="white-space:pre-wrap"> </span>= True</div><div><span style="white-space:pre-wrap"> </span>snres@gsnDraw <span style="white-space:pre-wrap"> </span>= False</div><div><span style="white-space:pre-wrap"> </span>snres@gsnFrame<span style="white-space:pre-wrap"> </span>= False</div><div><span style="white-space:pre-wrap"> </span></div><div><span style="white-space:pre-wrap"> </span>snres@lbLabelBarOn<span style="white-space:pre-wrap"> </span>= False<span style="white-space:pre-wrap"> </span>; turn off label bar</div><div><span style="white-space:pre-wrap"> </span>snres@cnInfoLabelOn<span style="white-space:pre-wrap"> </span>= False<span style="white-space:pre-wrap"> </span>; turn off info label</div><div><span style="white-space:pre-wrap"> </span>snres@cnLineLabelsOn<span style="white-space:pre-wrap"> </span>= False<span style="white-space:pre-wrap"> </span>; turn off contour labels</div><div><span style="white-space:pre-wrap"> </span>snres@cnLinesOn<span style="white-space:pre-wrap"> </span>= False<span style="white-space:pre-wrap"> </span>; no contour lines</div><div><span style="white-space:pre-wrap"> </span>snres@cnConstFEnableFill= True<span style="white-space:pre-wrap"> </span>; allow constant values (1=sig) to fill</div><div><span style="white-space:pre-wrap"> </span>snres@cnFillOn<span style="white-space:pre-wrap"> </span>= True<span style="white-space:pre-wrap"> </span>; turn on fill</div><div><span style="white-space:pre-wrap"> </span>snres@cnMonoFillPattern = True<span style="white-space:pre-wrap"> </span>; single fill pattern</div><div><span style="white-space:pre-wrap"> </span>snres@cnFillColor<span style="white-space:pre-wrap"> </span>= "black"<span style="white-space:pre-wrap"> </span>; </div><div><span style="white-space:pre-wrap"> </span>snres@cnFillPattern<span style="white-space:pre-wrap"> </span>= 17<span style="white-space:pre-wrap"> </span>; stipling</div><div><span style="white-space:pre-wrap"> </span></div><div><span style="white-space:pre-wrap"><snip></span></div><div><span style="white-space:pre-wrap"> </span></div><div><span style="white-space:pre-wrap"> </span>plot_diff = gsn_csm_contour_map(wks,Wrfdiff_mean,res)</div><div><span style="white-space:pre-wrap"> </span>plot_sn = gsn_csm_contour(wks,Wrfdiff_sn2,snres)</div><div><span style="white-space:pre-wrap"> </span></div><div><span style="white-space:pre-wrap"> </span>overlay(plot_diff,plot_sn)</div><div><span style="white-space:pre-wrap"> </span></div></div><div><br></div><div><br></div></div>
<br></div></div>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br></div>
</div></div></blockquote></div><br></div>