[ncl-talk] Interpolating CCSM grid to RAP/HRRR grid
Adam Phillips
asphilli at ucar.edu
Wed May 12 11:55:17 MDT 2021
Hi Malcolm,
One thing to check, is to make sure the input latitudes and longitudes are
in the same direction and within the same range as the regridded latitudes
and longitudes. Your input coordinates are:
lat: [20.26178..53.24607]
lon: [225..300]
Do your target grid longitudes fall within 225-300E? And do the latitudes
fall within 20-53N?
What is the output from :
printMinMax(lat2d,0)
printMinMax(lon2d,0)
Adam
On Wed, May 12, 2021 at 10:18 AM Wilson, Malcolm Tyree via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:
> Hi again,
>
> rgird2rcm ran fine without any errors, thank you! However, I'm still
> getting back nothing but filled values (every point is a fill value) in the
> returned, interpolated array. I even adjusted the data type of fi to double
> to see if that was causing issues, but the same issue appears. Below is the
> code that I'm using now, along with the data summary of both fi and the
> output (fi after interpolation) arrays. I'm wondering if this is supposed
> to happen, or is there something I'm missing? Thanks!
>
> ;Interpolation
>
> f = addfile ("~/small_hus_may_pert_2010-2050.nc","r")
> g = addfile ("~/hrrr.t00z.wrfprsf00.nc","r")
>
> lat = f->lat
> lon = f->lon
>
> lat2d = g->gridlat_0 ; (N,M)
> lon2d = g->gridlon_0
>
> ;LATD = todouble(lat2d) ; (N*M)
> ;LOND = (lon2d)
>
> fi = f->hus_pert
>
> hrrr_grid_data = rgrid2rcm_Wrap(lat,lon,fi,lat2d,lon2d, 0)
>
> printVarSummary(hrrr_grid_data)
>
> Variable: fi
> Type: float
> Total Size: 228384 bytes
> 57096 values
> Number of Dimensions: 4
> Dimensions and sizes: [time | 1] x [lev | 26] x [lat | 36] x [lon | 61]
> Coordinates:
> time: [ 0.. 0]
> lev: [0.9925562..0.00354464]
> lat: [20.26178..53.24607]
> lon: [225..300]
> Number Of Attributes: 1
> units : kg/kg
>
> Variable: hrrr_grid_data
> Type: float
> Total Size: 198134664 bytes
> 49533666 values
> Number of Dimensions: 4
> Dimensions and sizes: [time | 1] x [lev | 26] x [ygrid_0 | 1059] x
> [xgrid_0 | 1799]
> Coordinates:
> time: [ 0.. 0]
> lev: [0.9925562..0.00354464]
> Number Of Attributes: 3
> _FillValue : 9.96921e+36
> units : kg/kg
> ncl : rgrid2rcm used for interpolation
> ------------------------------
> *From:* 아나유즈바시츠 <ajuzbasic at pusan.ac.kr>
> *Sent:* Wednesday, May 12, 2021 8:06 AM
> *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,
>
> 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
> <https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.ncl.ucar.edu%2FDocument%2FFunctions%2FBuilt-in%2Frgrid2rcm.shtml&data=04%7C01%7Cmtw5270%40psu.edu%7C9f1ea306fc124903ca7308d9153e5a20%7C7cf48d453ddb4389a9c1c115526eb52e%7C0%7C1%7C637564179828314466%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000&sdata=dHzt0YSIMNPq61oYEX2kWf9Y%2FRKjSuhSqB11FZwZI1A%3D&reserved=0>.
> 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%7C9f1ea306fc124903ca7308d9153e5a20%7C7cf48d453ddb4389a9c1c115526eb52e%7C0%7C1%7C637564179828324453%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000&sdata=gqC2iltJNNr6mxVcQVabWp6xmzOj8i%2FbICeZ%2Fz9QRP4%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%7C9f1ea306fc124903ca7308d9153e5a20%7C7cf48d453ddb4389a9c1c115526eb52e%7C0%7C1%7C637564179828324453%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000&sdata=ogr4icTWYUSpDkU0ZLWAKBEu2kQvuJTYCsxC6yx90w4%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%7C9f1ea306fc124903ca7308d9153e5a20%7C7cf48d453ddb4389a9c1c115526eb52e%7C0%7C1%7C637564179828334443%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000&sdata=tS66RAR%2FoHbgmc0P4E%2FBjOsuUG9pX5jjbsza0WJu5Rk%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%7C9f1ea306fc124903ca7308d9153e5a20%7C7cf48d453ddb4389a9c1c115526eb52e%7C0%7C1%7C637564179828344441%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000&sdata=b2Q%2FEzHUPoBZAvWT0jgZNGC7MdZYbyT6bkUc%2F%2BkUnPE%3D&reserved=0>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
--
Adam Phillips
Associate Scientist, Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/ 303-497-1726
<http://www.cgd.ucar.edu/staff/asphilli>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20210512/0a5fbfe3/attachment.html>
More information about the ncl-talk
mailing list