<div dir="ltr">Many thanks Dennis. Indeed, "nan" in my input file was the problem.<div><br></div><div>Kind regards, Freddie</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jan 8, 2015 at 2:35 AM, Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu" target="_blank">shea@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>The fundamental rule of data processing is *look at your data*<br><br>Did you look at your data??? I don't think so.<br><br></div><div>Also, you showed no attempt to write data out.<br><br>[1]<br></div><div>The file has 'nan' .... NCL does not understand 'nan'.<br></div><div>You must change to _FillValue<br></div><div>See: <a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/isnan_ieee.shtml" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Built-in/isnan_ieee.shtml</a><br><br>[2]<br></div><div>After [1], use 'stat_dispersion' to examine the distribution of data.<br></div><div>This is for *your* benefit ... not NCL's benefit<br><br>[3] <br></div><div>To write out a file see:<br> <a href="http://www.ncl.ucar.edu/Applications/list_io.shtml" target="_blank">http://www.ncl.ucar.edu/Applications/list_io.shtml</a><br><br>[4] <br></div><div>I think your source data is messed up.<br></div><div><br><br></div><div><br>ncl 0> f = addfile("corrected_extreme5Dayspr_bcc-csm1-1_rcp45_r1i1p1_2006-2035.nc","r") <br>ncl 1> p = f->pr<br>ncl 2> print(p) <br>Variable: p<br>Type: double<br>Total Size: 655360 bytes<br> 81920 values<br>Number of Dimensions: 3<br>Dimensions and sizes: [T | 10] x [lat | 64] x [lon | 128]<br>Coordinates: <br> T: [1..10]<br> lat: [-87.86379883923263..87.86379883923263]<br> lon: [ 0..357.1875]<br>Number Of Attributes: 0<br>(0,0,0) nan<br>(0,0,1) nan<br>(0,0,2) nan<br>(0,0,3) nan<br>(0,0,4) nan<br>(0,0,5) nan<br>(0,0,6) nan<br>(0,0,7) nan<br>(0,0,8) nan<br>(0,0,9) nan<br></div>[snip]<br><div><div><div>(0,0,43) 0.5284684300422668<br>(0,0,44) 0.9607513546943665<br>(0,0,45) 1.205325365066528<br>(0,0,46) 1.090707898139954<br>(0,0,47) 1.714292883872986<br>(0,0,48) 1.940244793891907<br>(0,0,49) 2.598035097122192<br>(0,0,50) 2.348235607147217<br>(0,0,51) 1.526988863945007<br>(0,0,52) 0.7807736396789551<br>(0,0,53) 0.1979905366897583<br>(0,0,54) 0.1604289710521698<br></div><div>[snip]<br></div><div><br></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote"><span class="">On Wed, Jan 7, 2015 at 5:58 AM, Freddie Mpelasoka <span dir="ltr"><<a href="mailto:freddie.mpelasoka@gmail.com" target="_blank">freddie.mpelasoka@gmail.com</a>></span> wrote:<br></span><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class=""><div><div><div dir="ltr">I am trying to write the results of poisson_grid_fill function to a file instead of a plot by the attached simple script. I expected to generate a file with no missing values, the same as my input data. The examples only show the plots, but I need the numbers. Your help will be highly appreciated.</div>
</div></div><br></span>_______________________________________________<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>