<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">If you look at the output from the "printMinMax" commands, you will see that you have very large negative values for the minimum:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><div class="gmail_default" style="font-size:small">daily mean temperature analysis [degC] (degC) : min=-1.07374e+09 max=1.78964</div><div class="gmail_default" style="font-size:small">daily mean temperature analysis [degC] (degC) : min=-1.07374e+09 max=4.93295</div><div class="gmail_default" style="font-size:small">daily mean temperature analysis [degC] (degC) : min=-1.07374e+09 max=2.17995</div><div class="gmail_default" style="font-size:small">daily mean temperature analysis [degC] (degC) : min=-1.07374e+09 max=1.05083</div><div style="font-size:small"><br></div><div style="font-size:small">This indicates to me that your missing values are not correct on the file. You can easily fix this by setting this minimum value to missing for all four values:</div><div style="font-size:small"><br></div><div><div>tmin = min(TmeanM0)</div><div>TmeanM0 = where(TmeanM0.eq.tmin,TmeanM0@_FillValue,TmeanM0)</div><div>TmeanM1 = where(TmeanM1.eq.tmin,TmeanM1@_FillValue,TmeanM1)</div><div>TmeanM2 = where(TmeanM2.eq.tmin,TmeanM2@_FillValue,TmeanM2)</div><div>TmeanM3 = where(TmeanM3.eq.tmin,TmeanM3@_FillValue,TmeanM3)</div><div><br></div><div><br></div><div>--Mary</div><div><br></div><div><br></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jan 25, 2017 at 3:45 AM, Ruksana Abedin <span dir="ltr"><<a href="mailto:ruksana.abedin@gmail.com" target="_blank">ruksana.abedin@gmail.com</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><div><div><div>Dear Mary and NCL users,<br><br></div>I have got similar problem while making bias maps for temperature. The script, variable summary, min-max, the panel plot temp_mean_diff(3).jpg. and the input files are given. This time missing value has been set properly but still I am getting colour on the Ocean points. <br><br>These biases are calculated by subtracting observation from the climate model. However, if I do the opposite, bias = observation-model; I get no problem with the same script and the panel plot is temp_mean_diff(4).jpg. </div>Can you please help me to get rid of this problem. I need to have the bias maps using the first method. <br>Thank you very much for your kind support, this really keeps me going on with using ncl. <br><br></div>Best regards,<br></div>Ruksana <br><div><div><div><div><div><div><br>begin<br><br>; read in data for bias maps <br><br>fM0=addfile("/.../WAH_APHRO_<wbr>TMEAN_1988_2007_DJF_bias.nc", "r");<br>fM1=addfile("/.../WAH_APHRO_<wbr>TMEAN_1988_2007_MAM_bias.nc", "r");<br>fM2=addfile("/.../WAH_APHRO_<wbr>TMEAN_1988_2007_JJAS_bias.nc", "r");<br>fM3=addfile("/.../WAH_APHRO_<wbr>TMEAN_1988_2007_ON_bias.nc", "r");<br>;*****************************<wbr>******************************<wbr>******************************<wbr>******************************<wbr>********<br>; read variable<br><br>TmeanM0 = fM0 ->tave(:,:,:)<br>timeM0 = fM0 ->time<br>latM0 = fM0 ->lat<br>lonM0 = fM0 ->lon<br><br>TmeanM1 = fM1 ->tave(:,:,:)<br>timeM1 = fM1 ->time<br>latM1 = fM1 ->lat<br>lonM1 = fM1 ->lon<br><br>TmeanM2 = fM2 ->tave(:,:,:)<br>timeM2 = fM2 ->time<br>latM2 = fM2 ->lat<br>lonM2 = fM2 ->lon<br><br>TmeanM3 = fM3 ->tave(:,:,:)<br>timeM3 = fM3 ->time<br>latM3 = fM3 ->lat<br>lonM3 = fM3 ->lon<br><br>printVarSummary(TmeanM0)<br>printVarSummary(TmeanM1)<br>printVarSummary(TmeanM2)<br>printVarSummary(TmeanM3)<br><br>printMinMax(TmeanM0, 0)<br>printMinMax(TmeanM1, 0)<br>printMinMax(TmeanM2, 0)<br>printMinMax(TmeanM3, 0)<br><br>;*****************************<wbr>******************************<br>; create plots<br>;*****************************<wbr>******************************<br> <br> wks = gsn_open_wks("eps","temp_mean_<wbr>diff") <br> plot = new(4,graphic) <br> gsn_define_colormap(wks,"temp_<wbr>diff_18lev") ; choose colormap<div><div class="h5"><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></div></div> res@cnLineLabelsOn = False<br> res@lbLabelAutoStride = True<br> res@lbLabelBarOn = False<br> TmeanM1@_FillValue = 0.0<br> TmeanM1@_FillValue = -99.9<br> res@cnMinLevelValF = -5.<br> res@cnMaxLevelValF = 5.; 1<br> res@cnLevelSpacingF = 1; 1<br> res@lbLabelAutoStride = True<br> <br>plot(0) = gsn_csm_contour_map(wks,<wbr>TmeanM0(0,:,:),res)<br>plot(1) = gsn_csm_contour_map(wks,<wbr>TmeanM1(0,:,:),res)<br>plot(2) = gsn_csm_contour_map(wks,<wbr>TmeanM2(0,:,:),res)<br>plot(3) = gsn_csm_contour_map(wks,<wbr>TmeanM3(0,:,:),res)<span class=""><br><br>;*****************************<wbr>******************************<br>; create panel<br>;*****************************<wbr>******************************<br>resP = True <br>resP@gsnFrame = False <br>resP@gsnMaximize = True<br></span>resP@gsnPanelLabelBar = True ; add common labelbar <br>;resP@lbLabelFontHeightF = 0.005<br>;resP@gsnPanelMainString =""<br>;resP@txString =""<br>resP@gsnPanelFigureStrings =(/"DJF","MAM","JJAS","ON"/) <wbr> ; add strings to panel <br>gsn_panel(wks,plot,(/4,4/),<wbr>resP) ; now draw as one plot <br><br>end<br><br>;*****************************<wbr>******************************<wbr>******************************<wbr>******************************<wbr>********<br><u>Variable information<br></u><br>Variable: TmeanM0<span class=""><br>Type: float<br>Total Size: 3480 bytes<br> 870 values<br>Number of Dimensions: 3<br>Dimensions and sizes: [time | 1] x [lat | 30] x [lon | 29]<br>Coordinates:<br></span> time: [19980101..19980101]<span class=""><br> lat: [15.25..29.75]<br> lon: [ 83.. 97]<br>Number Of Attributes: 5<br></span> long_name : daily mean temperature analysis [degC]<br> units : degC<span class=""><br> _FillValue : -99.9<br> missing_value : -99.9<br> level_description : Earth surface<br><br></span>Variable: TmeanM1<span class=""><br>Type: float<br>Total Size: 3480 bytes<br> 870 values<br>Number of Dimensions: 3<br>Dimensions and sizes: [time | 1] x [lat | 30] x [lon | 29]<br>Coordinates:<br></span> time: [19971016..19971016]<span class=""><br> lat: [15.25..29.75]<br> lon: [ 83.. 97]<br>Number Of Attributes: 5<br></span> long_name : daily mean temperature analysis [degC]<br> units : degC<span class=""><br> _FillValue : -99.9<br> missing_value : -99.9<br> level_description : Earth surface<br><br></span>Variable: TmeanM2<span class=""><br>Type: float<br>Total Size: 3480 bytes<br> 870 values<br>Number of Dimensions: 3<br>Dimensions and sizes: [time | 1] x [lat | 30] x [lon | 29]<br>Coordinates:<br></span> time: [19980201..19980201]<span class=""><br> lat: [15.25..29.75]<br> lon: [ 83.. 97]<br>Number Of Attributes: 5<br></span> long_name : daily mean temperature analysis [degC]<br> units : degC<span class=""><br> _FillValue : -99.9<br> missing_value : -99.9<br> level_description : Earth surface<br><br></span>Variable: TmeanM3<span class=""><br>Type: float<br>Total Size: 3480 bytes<br> 870 values<br>Number of Dimensions: 3<br>Dimensions and sizes: [time | 1] x [lat | 30] x [lon | 29]<br>Coordinates:<br></span> time: [19980501..19980501]<span class=""><br> lat: [15.25..29.75]<br> lon: [ 83.. 97]<br>Number Of Attributes: 5<br></span> long_name : daily mean temperature analysis [degC]<br> units : degC<span class=""><br> _FillValue : -99.9<br> missing_value : -99.9<br> level_description : Earth surface<br></span>(0) daily mean temperature analysis [degC]: min=-1.07374e+09 max=1.78964<br>(0) daily mean temperature analysis [degC]: min=-1.07374e+09 max=4.93295<br>(0) daily mean temperature analysis [degC]: min=-1.07374e+09 max=2.17995<br>(0) daily mean temperature analysis [degC]: min=-1.07374e+09 max=1.05083<br><br><br></div></div></div></div></div></div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Nov 29, 2016 at 5: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">Sorry for the delay. I was out of the office last week.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I assume you are talking about the contour line labels. Try setting:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">res@cnLineLabelsOn = False</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I suggest that you visit our examples page at:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="http://test.www.ncl.ucar.edu/Applications/" target="_blank">http://test.www.ncl.ucar.edu/A<wbr>pplications/</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">You will see lots of topics on this page, like "contour labels" and "contour effects". If you click on any one of the topics, you will see a page with a list of examples, each with a sample NCL script and graphic. This is one of the best ways to learn about the various graphical resources in NCL and how to set them.</div><div class="gmail_default"><br></div><div class="gmail_default">Good luck,</div><div class="gmail_default"><br></div><div class="gmail_default">--Mary</div><div class="gmail_default"><br></div></div><div class="m_-6007865439731416054HOEnZb"><div class="m_-6007865439731416054h5"><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Nov 21, 2016 at 6:55 AM, Ruksana Abedin <span dir="ltr"><<a href="mailto:ruksana.abedin@gmail.com" target="_blank">ruksana.abedin@gmail.com</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><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="m_-6007865439731416054m_-6748738048797173202HOEnZb"><div class="m_-6007865439731416054m_-6748738048797173202h5"><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_-6007865439731416054m_-6748738048797173202m_-8135060990140217472m_-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_-6007865439731416054m_-6748738048797173202m_-8135060990140217472m_-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_-6007865439731416054m_-6748738048797173202m_-8135060990140217472m_-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_-6007865439731416054m_-6748738048797173202m_-8135060990140217472m_-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_-6007865439731416054m_-6748738048797173202m_-8135060990140217472m_-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_-6007865439731416054m_-6748738048797173202m_-8135060990140217472m_-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="m_-6007865439731416054m_-6748738048797173202m_-8135060990140217472h5">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="m_-6007865439731416054m_-6748738048797173202m_-8135060990140217472h5"><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_b<wbr>ias.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>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div>