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

[snehlata Tirkey]

Dear Dennis,
Thank you for the precise solution!
I am able to work with the data now.
Yours sincerely,Snehlata Tirkeye-mail: snehlata_tirkey at yahoo.com
 

    On Friday, 4 January 2019, 12:17:32 am GMT+5:30, Dennis Shea <shea at ucar.edu> wrote:  
 
     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

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 Tirkeye-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/20190104/b53fe0f5/attachment.html>