[ncl-talk] How to subset a variable for a particular area from unstructured gridded data

[Dennis Shea]

    latS  = ...    ; south boundary
    latN  = ...    ; north      "
    lonL  = ...   ; left boundary
    lonR  = ...  ; right        "

    f = addfile(...)
    LAT  = f->lat
    LON = f->lon
    ii  = ind(LAT.ge.latS .and. LAT.le.latN .and.  \
                LON.ge.lonL .and. LON.le.lonR)

    lat = LAT(ii)
    lon = LON(ii)
   prc = f->PRCC(:,ii)
   printVarSummary(lat)
*   printVarSummary(lon)*
   printVarSummary(prc)

 ======
*http://www.ncl.ucar.edu/Document/Functions/Built-in/ind.shtml*
<http://www.ncl.ucar.edu/Document/Functions/Built-in/ind.shtml>

On Thu, Jan 3, 2019 at 5:21 AM snehlata Tirkey <snehlata_tirkey at yahoo.com>
wrote:

> Dear NCL community,
>
> I have a global dataset with monthly temporal resolution on unstructured
> grid.
>
> Part of the ncdump is as follows-
>
> *dimensions:*
> * ncol = 48602 ;*
> * time = UNLIMITED ; // (1 currently)*
> * nbnd = 2 ;*
> * chars = 8 ;*
> * lev = 30 ;*
> * ilev = 31 ;*
> * psubcols = 4 ;*
> *variables:*
> * double lat(ncol) ;*
> * lat:long_name = "latitude" ;*
> * lat:units = "degrees_north" ;*
> * double lon(ncol) ;*
> * lon:long_name = "longitude" ;*
> * lon:units = "degrees_east" ;*
> * double area(ncol) ;*
> * area:long_name = "gll grid areas" ;*
> * double lev(lev) ;*
> * lev:long_name = "hybrid level at midpoints (1000*(A+B))" ;*
> * lev:units = "hPa" ;*
> * lev:positive = "down" ;*
> * lev:standard_name = "atmosphere_hybrid_sigma_pressure_coordinate" ;*
> * lev:formula_terms = "a: hyam b: hybm p0: P0 ps: PS" ;*
> *        float PRECC(time, ncol) ;*
> * PRECC:units = "m/s" ;*
> * PRECC:long_name = "Convective precipitation rate (liq + ice)" ;*
> * PRECC:cell_methods = "time: mean" ;*
>
> I am able to plot the above variable (PRECC) globally using sfXArray and
> sfYArray but how to subset the data for a smaller area for plotting or
> further processing?
>
> I could find information for 2D lat and lon type of grid in NCL but could
> not find information for such 1D lat/lon type of grid. Please help.
>
> Thanking in anticipation.
>
>
> Yours sincerely,
> Snehlata Tirkey
> e-mail:snehlata_tirkey at yahoo.com
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> 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/20190103/0eca58ec/attachment.html>