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

Dennis Shea shea at ucar.edu
Wed Jan 7 08:35:23 MST 2015


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/20150107/a3de4618/attachment.html 


More information about the ncl-talk mailing list