<div dir="ltr"><div><div><div>Dear Mary,<br><br></div>Thank you very much for pointing at the issue and the fix for it. Now, the plot could be corrected with no color on the ocean points. <br><br>But what about excluding the data values on the contours randomly on the land points? How can I disable this option? I have not set any resource for this, it is done somehow automatically. Please can you help me to avoid these values on the maps?<br><br></div>Best regards,<br></div>Ruksana <br><div><div><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Nov 21, 2016 at 12:37 PM, Mary Haley <span dir="ltr"><<a href="mailto:haley@ucar.edu" target="_blank">haley@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"><div class="gmail_default" style="font-size:small">Ruksana,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">It is very important that you look at your data before you assume that NCL is plotting it incorrectly. </div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">If you look at "ppt3" with printVarSummary and printMinMax, you will see that your missing values are not set correctly.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Here's the printVarSummary, which claims the missing value is -99.9:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><div class="gmail_default"><font face="monospace, monospace">Variable: ppt3</font></div><div class="gmail_default"><font face="monospace, monospace">Type: float</font></div><div class="gmail_default"><font face="monospace, monospace">Total Size: 3480 bytes</font></div><div class="gmail_default"><font face="monospace, monospace"> 870 values</font></div><div class="gmail_default"><font face="monospace, monospace">Number of Dimensions: 3</font></div><div class="gmail_default"><font face="monospace, monospace">Dimensions and sizes:<span class="m_-6263216940847785553gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>[time | 1] x [lat | 30] x [lon | 29]</font></div><div class="gmail_default"><font face="monospace, monospace">Coordinates: </font></div><div class="gmail_default"><font face="monospace, monospace"> time: [14400..14400]</font></div><div class="gmail_default"><font face="monospace, monospace"> lat: [15.25..29.75]</font></div><div class="gmail_default"><font face="monospace, monospace"> lon: [ 83.. 97]</font></div><div class="gmail_default"><font face="monospace, monospace">Number Of Attributes: 5</font></div><div class="gmail_default"><font face="monospace, monospace"> long_name :<span class="m_-6263216940847785553gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>daily precipitation analysis interpolated onto 0.5deg grids [mm/day]</font></div><div class="gmail_default"><font face="monospace, monospace"> units :<span class="m_-6263216940847785553gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>mm/day</font></div><div class="gmail_default"><font face="monospace, monospace"> _FillValue :<span class="m_-6263216940847785553gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>-99.9</font></div><div class="gmail_default"><font face="monospace, monospace"> missing_value :<span class="m_-6263216940847785553gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>-99.9</font></div><div class="gmail_default"><font face="monospace, monospace"> level_description :<span class="m_-6263216940847785553gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>Earth surface</font></div><div class="gmail_default"><br></div><div class="gmail_default">However, if you do a printMinMax on pp3, you will see a maximum value of 1e+20, which seems suspicious:</div><div class="gmail_default"><br></div><div class="gmail_default"><font face="monospace, monospace">daily precipitation analysis interpolated onto 0.5deg grids [mm/day] (mm/day) : min=-2497.37 max=1e+20</font></div><div class="gmail_default"><br></div><div class="gmail_default">My guess is that you need to set the missing value of ppt3 to 1e20:</div><div class="gmail_default"><br></div><div class="gmail_default"><font face="monospace, monospace">ppt3@_FillValue = 1e20</font></div><div class="gmail_default"><br></div><div class="gmail_default">--Mary</div><div class="gmail_default"><br></div><div class="gmail_default"> </div></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On Sat, Nov 19, 2016 at 2:35 PM, Ruksana Abedin <span dir="ltr"><<a href="mailto:ruksana.abedin@gmail.com" target="_blank">ruksana.abedin@gmail.com</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 class="h5"><div dir="ltr"><div><div><div>Hi, <br><br></div>I have plotted a panel of two bias maps. The plot, input files and the script is attached. Problem is although my data is only on land points, on the first bias map, the ocean is colored as orange, which should not be the case. And the maps show few values of the contours on the land points irregularly, but I do not want any values to be shown. Can you please help me to do these two corections? <br><br>Thank you for your time and kind efforts.<br><br></div>Best regards,<br></div>Ruksana <br><div><div><div><br><br><br>begin<br><br>; read in data<br><br>f3=addfile(".........../WAH_AP<wbr>HRO_1998_2007_ANN_sum_<a href="http://precip_bias.nc" target="_blank">precip_<wbr>bias.nc</a>", "r");<br>f4=addfile(".........../WAH_GP<wbr>CC_1998_2007_ANN_sum_precip_<a href="http://bias.nc" target="_blank">bi<wbr>as.nc</a>", "r");<br><br>ppt3 = f3->precip<br>time = f3->time<br>lat = f3->lat<br>lon = f3->lon<br><br>ppt4 = f4->precip<br>time = f4->time<br>lat = f4->lat<br>lon = f4->lon<br><br>wks = gsn_open_wks("eps","ppt1") <br>plot = new(2,graphic) <br><br>printVarSummary(ppt3)<br>printVarSummary(ppt4)<br><br> gsn_define_colormap(wks,"perc<wbr>ent_11lev") ; choose colormap for bias map<br> ;gsn_reverse_colormap(wks)<br> res = True<br> res@gsnDraw = False<br> res@gsnFrame = False<br> res@cnInfoLabelOn = False <br> res@cnFillOn = True ; turn on color<br> res@cnLinesOn = False ; turn off contour lines<br> res@cnLevelSpacingF = 0.10 ; contour interval<br> res@cnFillDrawOrder = "PreDraw" ; draw contours first<br> res@lbLabelStride = 2 ; stride on label bar<br> res@gsnSpreadColors = True ; use full colormap<br> res@gsnSpreadColorEnd = -3 ; -3 don't use land color<br> res@gsnAddCyclic = False ; regional data <br> res@gsnStringFontHeightF = 0.01 <br> res@tiMainString = ""<br> res@gsnRightString = ""<br> ;res@mpProjection = "LambertConformal"<br> ;res@gsnMaskLambertConformal = True<br> res@gsnLeftString = ""<br> res@mpMaxLatF = 30;max(lat) ; zoom in on region<br> res@mpMinLatF = 15;min(lat)<br> res@mpMinLonF = 83;min(lon)<br> res@mpMaxLonF = 97; max(lon)<br> res@mpCenterLonF = 10 ; def is zero<br> res@mpLandFillColor = "Transparent"<br> res@cnLevelSelectionMode= "ManualLevels"<br> <br> res@cnMinLevelValF = -2150.00 ;for bias maps <br> res@cnMaxLevelValF = 4450;100<br> res@cnLevelSpacingF = 150 ; 150<br> <br> res@lbLabelBarOn = False<br> res@lbLabelAutoStride = True<br> ppt3@_FillValue = 0.0<br> ppt3@_FillValue = -999<br><br>plot(0) = gsn_csm_contour_map(wks,ppt3(0<wbr>,:,:),res)<br>plot(1) = gsn_csm_contour_map(wks,ppt4(0<wbr>,:,:),res)<br><br>;*****************************<wbr>******************************<br>; create panel<br>;*****************************<wbr>******************************<br><br>resP = True <br>resP@gsnFrame = False <br>resP@gsnMaximize = True<br>resP@gsnPanelLabelBar = True<br>resP@lbLabelFontHeightF = 0.01<br>gsn_panel(wks,plot,(/2,2/),res<wbr>P) ; now draw as one plot <br><br>end<br><br></div></div></div></div>
<br></div></div>______________________________<wbr>_________________<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/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br></div>