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

[Wilson, Malcolm Tyree]

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<mailto: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<mailto: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/0c32c486/attachment.html>