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

Dennis Shea shea at ucar.edu
Thu Oct 8 20:23:55 MDT 2020


   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/20201008/d3c5048d/attachment.html>


More information about the ncl-talk mailing list