[ncl-talk] Time series at each grid point from a netcdf file

[Zilore Mumba]

Thanks! Dennis, I will try the code as soon as I have access to my computer.
For Rick, I want the output in text (txt) format. If I can have each grid
point on a line, then I can run fortran on the data outside of the ncl
script. I am assuming this is what the code from Dennis does. I will try it

On Fri, Oct 9, 2020 at 4:24 AM Dennis Shea <shea at ucar.edu> wrote:

>    prc(time,lat,lon)
> ---------------------
>    dimprc = dimsizes(prc)
>    ntim = dimprc(0)
>    nlat  = dimprc(1)
>    mlon = dimprc(2)
>
> ; loop over each grid point
>
>    do nl=0,nlat-1
>       do ml=0,mlon-1
>            prcts = prc(:,nl,ml)      ; prc(tim)   <=== time series
>            :
>           [use NCL or invoke fortran via a shared object]
>            :
>       end do    ; ml
>    end do       ; nl
>
> On Thu, Oct 8, 2020 at 1:58 PM Zilore Mumba via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
>> I have a netcdf file containing daily precipitation data over 30 years on
>> a grid of 96 points in x and 99 points in y. I would like to work with the
>> data on a grid-point by grid-point basis. I think this would be best
>> achieved by having a time series of the data at each grid point. I can then
>> use, eg. fortran to make some computation, such as the start/end of the wet
>> session.
>> I wonder if someone would have a script that can generate such time
>> series. My search in the ncl documentation and archives achieves nothing.
>> I will appreciate any help or advice on this.
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at mailman.ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20201009/9abd655f/attachment.html>