[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