[ncl-talk] Problem writing poisson_grid_fill results to a file

[Freddie Mpelasoka]

Many thanks Dennis. Indeed, "nan" in my input file was the  problem.

Kind regards, Freddie

On Thu, Jan 8, 2015 at 2:35 AM, Dennis Shea <shea at ucar.edu> wrote:

> The fundamental rule of data processing is *look at your data*
>
> Did you look at your data??? I don't think so.
>
> Also, you showed no attempt to write data out.
>
> [1]
> The file has 'nan' .... NCL does not understand 'nan'.
> You must change to _FillValue
> See: http://www.ncl.ucar.edu/Document/Functions/Built-in/isnan_ieee.shtml
>
> [2]
> After [1], use 'stat_dispersion' to examine the distribution of data.
> This is for *your* benefit ... not NCL's benefit
>
> [3]
> To write out a file see:
>     http://www.ncl.ucar.edu/Applications/list_io.shtml
>
> [4]
> I think your source data is messed up.
>
>
>
> ncl 0> f =
> addfile("corrected_extreme5Dayspr_bcc-csm1-1_rcp45_r1i1p1_2006-2035.nc","r")
>
> ncl 1> p = f->pr
> ncl 2> print(p)
> Variable: p
> Type: double
> Total Size: 655360 bytes
>             81920 values
> Number of Dimensions: 3
> Dimensions and sizes:   [T | 10] x [lat | 64] x [lon | 128]
> Coordinates:
>             T: [1..10]
>             lat: [-87.86379883923263..87.86379883923263]
>             lon: [   0..357.1875]
> Number Of Attributes: 0
> (0,0,0)  nan
> (0,0,1)  nan
> (0,0,2)  nan
> (0,0,3)  nan
> (0,0,4)  nan
> (0,0,5)  nan
> (0,0,6)  nan
> (0,0,7)  nan
> (0,0,8)  nan
> (0,0,9)  nan
> [snip]
> (0,0,43)        0.5284684300422668
> (0,0,44)        0.9607513546943665
> (0,0,45)        1.205325365066528
> (0,0,46)        1.090707898139954
> (0,0,47)        1.714292883872986
> (0,0,48)        1.940244793891907
> (0,0,49)        2.598035097122192
> (0,0,50)        2.348235607147217
> (0,0,51)        1.526988863945007
> (0,0,52)        0.7807736396789551
> (0,0,53)        0.1979905366897583
> (0,0,54)        0.1604289710521698
> [snip]
>
>
> On Wed, Jan 7, 2015 at 5:58 AM, Freddie Mpelasoka <
> freddie.mpelasoka at gmail.com> wrote:
>
>> 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.
>>
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150108/6aa8c90a/attachment.html