[ncl-talk] ncl script for NARR data

Dennis Shea shea at ucar.edu
Thu Jan 24 13:27:43 MST 2019


*http://www.ncl.ucar.edu/Document/Functions/Contributed/rcm2points_Wrap.shtml*
<http://www.ncl.ucar.edu/Document/Functions/Contributed/rcm2points_Wrap.shtml>

===

   lonDen = (/-104.867, -102.0/)
   latDen = (/ 39.75,     38.0/)

   f     = *addfile*
<http://www.ncl.ucar.edu/Document/Functions/Built-in/addfile.shtml>
("some_NARR_file.nc", "r")
   lat2d = f->xlat                         ; size = (nlat,nlon)
   lon2d = f->xlon                         ; size = (nlat,nlon)
   x     = f->X
   xDen  = *rcm2points_Wrap*
<http://www.ncl.ucar.edu/Document/Functions/Contributed/rcm2points_Wrap.shtml>
(lat2d, lon2d, x, lat, lon, 0)
   printVarSummary(xDen)

===


  dir_narr = "./"
  fil_narr = "narr-a_221_20110508_0000_000.grb"
  pth_narr = dir_narr + fil_narr
  f_narr   = addfile(pth_narr, "r")

;---Read variables
                                    ;          0               1
     2          3
  x_narr   = f_narr->TMP_221_ISBL   ; ( initial_time0_hours, lv_ISBL4,
gridx_221, gridy_221 )
  printVarSummary(x_narr)

  lat2d    = f_narr->gridlat_221
  lon2d    = f_narr->gridlon_221
  printMinMax(lat2d,0)
  printMinMax(lon2d,0)

;---Interpolate to user locations and levels for all times
  lat_pts = (/.../)
  lon_pts = (/.../)

  x_pts   = *rcm2points_Wrap*
<http://www.ncl.ucar.edu/Document/Functions/Contributed/rcm2points_Wrap.shtml>
(lat2d, lon2d, x_narr, lat, lon, 0)

  printVarSummary(x_pts)




On Thu, Jan 24, 2019 at 11:29 AM Mukul Tewari <mukultewari1 at gmail.com>
wrote:

> Hi,
>
> I would like to generate time series at a given lat/lon point from
> NARR. Does anyone has any suggestion, how to extract meteorological
> variables from NARR at a given point in order to generate time series
> over a year?
>
> thanks,
> Mukul
> _______________________________________________
> 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/20190124/1b75a2ac/attachment.html>


More information about the ncl-talk mailing list