[ncl-talk] Interpolating CCSM grid to RAP/HRRR grid

[아나유즈바시츠]

Hi Malcom,

Apologies, I mixed it up. The function you need is rgrid2rcm, which is, as
you said, opposite of rcm2rgrid. This one:
http://www.ncl.ucar.edu/Document/Functions/Built-in/rgrid2rcm.shtml. This
should do what you want, without giving you all the missing values.



On Wed, May 12, 2021 at 9:00 PM Wilson, Malcolm Tyree <mtw5270 at psu.edu>
wrote:

> Hi,
>
> So, I'm trying to re-grid from a rectilinear grid (CCSM) to a curvilinear
> grid (HRRR). The CCSM had one-dimensional coordinates arrays while the HRRR
> has 2-dimensional coordinate arrays. So, the opposite of what rcm2rgrid is
> doing, which is why I tried using linint2_points_Wrap instead. Thanks!
>
> -Malcolm W.
> ------------------------------
> *From:* 아나유즈바시츠 <ajuzbasic at pusan.ac.kr>
> *Sent:* Tuesday, May 11, 2021 10:42 PM
> *To:* Wilson, Malcolm Tyree <mtw5270 at psu.edu>
> *Cc:* ncl-talk at ucar.edu <ncl-talk at ucar.edu>
> *Subject:* Re: [ncl-talk] Interpolating CCSM grid to RAP/HRRR grid
>
> Hi Malcom,
>
> It seems to me that you are trying to regrid from curvilinear grid (with
> 2D lat and lon) to rectilinear grid. You should rcm2rgrid function for
> that, not linint, as the linint one is meant to interpolate from one
> rectilinear grid to another.
>  Try to use this
> https://www.ncl.ucar.edu/Document/Functions/Contributed/rcm2rgrid_Wrap.shtml
> <https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.ncl.ucar.edu%2FDocument%2FFunctions%2FContributed%2Frcm2rgrid_Wrap.shtml&data=04%7C01%7Cmtw5270%40psu.edu%7C35391c78894144d8f84008d914ef9fc5%7C7cf48d453ddb4389a9c1c115526eb52e%7C0%7C1%7C637563841698774881%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=c8VOtCN6C6k3Bh7AadVKPva%2B8tO127MQhGIblceCsLE%3D&reserved=0>
>
> On Wed, May 12, 2021 at 10:51 AM Wilson, Malcolm Tyree via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
> Hello Everyone,
>
> I've been trying to interpolate humidity and temperature data from CCSM to
> the RAP/HRRR grid using NCL (specifically using the function linint2_points
> and linint2_points_Wrap). I created a script that should interpolate the
> data from CCSM to the HRRR grid, but the script returned an array that
> contains nothing but filled values (9.96921e+36) and no data of the
> variable field I'm interpolating. I don't understand this happened, and I
> used both the linint2_points and linint2_points_Wrap functions. Is this
> supposed to happen, or and I'm missing something? Below is the code I used.
> Any advice would be appreciated. Thanks!
>
> -Malcolm W.
>
>   f = addfile ("~/small_hus_may_pert_2010-2050.nc
> <https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fsmall_hus_may_pert_2010-2050.nc%2F&data=04%7C01%7Cmtw5270%40psu.edu%7C35391c78894144d8f84008d914ef9fc5%7C7cf48d453ddb4389a9c1c115526eb52e%7C0%7C1%7C637563841698774881%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=nZF82pZlPtdXFSwL8fKEx7y7M1JHUXMpDykFqoPacj0%3D&reserved=0>","r")
>
>   g = addfile ("~/hrrr.t00z.wrfprsf00.nc
> <https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fhrrr.t00z.wrfprsf00.nc%2F&data=04%7C01%7Cmtw5270%40psu.edu%7C35391c78894144d8f84008d914ef9fc5%7C7cf48d453ddb4389a9c1c115526eb52e%7C0%7C1%7C637563841698784876%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=DqpNWw6mKEiPm8AD6MphsTN81pjixrN9500YT8rtlO4%3D&reserved=0>
> ","r")
>
>   lat = f->lat      ;(length=36)
>   lon = f->lon    ;(length=61)
>
>   lat2d = g->gridlat_0              ; (shape: (1059,1799))
>   lon2d = g->gridlon_0            ; (shape: (1059,1799))
>
>   LAT1D = ndtooned(lat2d)
>   LON1D = ndtooned(lon2d)
>
>   fi = f->hus_pert ;(humidity data I'm interpolating)
>
>   fo = linint2_points_Wrap(lon,lat,fi, False, LON1D,LAT1D, 0)
>
>   hrrr_data = onedtond( ndtooned(fo), (/1,26,1059,1799/) ) ;(conatins
> nothing but filled values)
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
> <https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fmailman.ucar.edu%2Fmailman%2Flistinfo%2Fncl-talk&data=04%7C01%7Cmtw5270%40psu.edu%7C35391c78894144d8f84008d914ef9fc5%7C7cf48d453ddb4389a9c1c115526eb52e%7C0%7C1%7C637563841698784876%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=5%2FTIRa09iW7Si%2BVoUNYAP2%2FMpTXzUTkl79HqyOnifxc%3D&reserved=0>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20210512/2e5f66bd/attachment.html>