<div dir="ltr"><div dir="ltr"><div dir="ltr">A slight correction here, sorry.  It is CellFill that needs 4 adjacent grid points to draw each grid box.  RasterFill will draw a box for every point with non-missing data.  See the resource docs for cnFillMode.  RasterFill would be a good initial setting for Vanúcia to try, for debugging.</div><div dir="ltr"><br></div><div dir="ltr"><br></div><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Feb 17, 2020 at 7:56 PM Dennis Shea via ncl-talk <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>For a</div><div>      res@cnFillMode  = "RasterFill" <br></div><div>there must be 4 grid points that are adjacent to one another for a grid box to be drawn.<br></div><div>Likely, that is not the case for the <b>VERY SPARSELY POPULATED</b> arrays you have available.</div><div>In fact, NCL's  <a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/stat_dispersion.shtml" target="_blank"><b>stat_dispersion </b></a> function indicates that more than <b> 99.999+%</b> of the array values are missing [_FillValue]. Specifically, of the 946956000 values only 14 are not missing.</div><div><br></div><div>[time | 365] x [lat | 1410] x [lon | 1840]  ===> 946956000</div><div><br></div><div>=========================================</div><div>The script ou attached:</div><div>:<br></div><div>The following is <b>not</b> 'time'.<br></div><div>       time     = r(:,0,0)                      ; selecting just time</div><div>This assignment yields all the 'r' values over 365 times at the grid point with index values (0,0)</div><div>Try the following to see the values in 'time'<br></div><div>       print(time)                                ; values of array 'r' over all times at grid point (0,0)<br></div><div>       exit<br></div><div><br></div><div>'time' would be</div><div>       time = r<b>&</b>time                           ; extract all 'time' coordinates associated with variable 'r'</div><div><br></div><div>or, explicitly extract the variable 'time' from the file(s).</div><div>        time = dat[:]->time</div><div><br></div><div>The above are equivalent.</div><div>==============================</div><div>I have no idea why the following is being assigned:</div><div>         r@_FillValue = -999   <br></div><div>The _FillValue is already -999.0</div><div>===============================</div><div><br></div><div>Attached are two scripts :</div><div>[1] <b>tst_vanucia.ncl_v0  </b>uses stat_dispersion to illustrate its usage.</div><div>[2] <b>tst_vanucia.ncl_v1</b>: creates the attached plot</div><div>     <br></div><div>     It uses<a href="http://www.ncl.ucar.edu/Applications/polyg.shtml" target="_blank"><b> polymarkers </b></a>to illustrate the 'nfr' values and prints the desired values.<br></div><div>=====================<br></div><div>     YYYYMMDD   LAT    LON   VALUE<br></div><div>(0)        2017012112  -24.4  -57.01  1<br>(1)    2017012912  -19.58  -50.4  1<br>(2)    2017013012  -22.35  -57.97  1<br>(3)   2017020412  -17.45  -46.58  1<br>(4)   2017031112  -24.71  -57.28  1<br>(5)   2017100212  -20.21  -43.62  1<br>(6)   2017100612  -24.86  -57.87  1<br>(7)   2017101312  -22.87  -53.63  1<br>(8)   2017102712  -19.19  -57.16  1<br>(9)   2017111712  -22.51  -44.11  1<br>(10)  2017112412  -19.63  -54.81  3<br>(11)  2017121512  -18.2  -54.46  2<br>(12)   2017121612  -21.58  -42.03  1<br>(13)  2017122612  -17.79  -45.92  2</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Feb 17, 2020 at 7:27 AM Vanúcia Schumacher via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">




<div dir="ltr">
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
<span style="font-size:12pt;color:rgb(51,51,51)">Hi NCL's user</span></div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
<span style="font-size:12pt;color:rgb(51,51,51)">I would like to know how I can extract the date information (YYYYMMDD) and coords (LAT LON)  for each non-missing value from the variable obtained?</span><span><br>
</span>
<div><span style="font-size:12pt;color:rgb(51,51,51)">And how could I plot these points, because they are 14 valid pixels (cells-size) and when I plot only 1 pixel appears</span><br>
</div>
<div><br>
</div>
<span style="font-size:12pt;color:rgb(51,51,51)">Attached are my script and data which I am trying to do this.</span><br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
<span><br>
</span></div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
<span style="font-size:12pt;color:rgb(51,51,51)">I appreciate any help</span></div></div></blockquote></div></blockquote></div></div></div>