Re: rectilinear grids to curvilinear grids

[1] I have no idea what you did w rgrid2rcm.
     Basically, you say it did not work. SO ... how is
     ncl-talk supposed to work with that information?

[2] Did you look at the ESMF Example script 29?
     Did you try understand what was being done?
     Did you try it?

On 5/22/14, 4:28 PM, Sha Feng wrote:
> Dennis,
>
> Sorry for typos. I did use "rgrid2rcm_Wrap". Also, I tried to use ESMF
> regrid functions, but I don't have weights for interpolation. Could you
> help me out, please?
>
> Thanks,
> Sha
>
>
> On Thu, May 22, 2014 at 3:23 PM, Dennis Shea <shea@ucar.edu> wrote:
>
>> You should use: rgrid2rcm or, better, rgrid2rcm_Wrap
>>
>> https://www.ncl.ucar.edu/Document/Functions/Built-in/rgrid2rcm.shtml
>>
>> https://www.ncl.ucar.edu/Document/Functions/Contributed/rgrid2rcm_Wrap.
>> shtml
>>
>> The 0.1deg x 0.1deg is a rectilinear grid ('rgrid(
>> The WRF is curvilinear.
>> ========
>>
>> It would be must faster to use the ESMF regrid functions.
>>
>> https://www.ncl.ucar.edu/Applications/ESMF.shtml
>>
>> Example 29 (Part B) shows how to interpolate the
>> rectilinear grid to a WRF grid.
>>
>>
>>
>>
>> On 5/22/14, 4:12 PM, Sha Feng wrote:
>>
>>> Hi Mary,
>>>
>>> I was trying to regrid data on 0.1deg x 0.1deg to a 12km x 12 km WRF
>>> domain
>>> using the function of "rcm2rgrid." The description of this function looks
>>> good match for this purpose. But after interpolation, the output doesn't
>>> look good to use (assuming nothing wrong in interpolation). After ncdump
>>> the variable, it shows a lot of
>>>
>>> ---- _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
>>> _,
>>> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
>>> _,
>>> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
>>> _,
>>> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
>>> _,
>>> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
>>> _,
>>> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
>>> _,
>>> _, _, _, _, _, _, _, _, _, _, _,
>>> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
>>> ---
>>>
>>> and ncview can't even run through. I attached my ncl codes as follows, and
>>> a plot for the source data that is carbon emission (I guess, rcm2rgrid
>>> might not be a good option for such scattering data). Could you take a
>>> look
>>> and help me for interpolation?
>>>
>>> ....
>>> f = addfile(dir+fname,"r")
>>> xin = f->E_Carbon
>>> latin = f->lat
>>> lonin = f->lon
>>> printMinMax(lonin,True)
>>> printMinMax(xin,True)
>>>
>>> a = addfile(dir+"wrfout_d"+dm+".nc","r")
>>> lat = wrf_user_getvar(a,"lat",0)
>>> lon = wrf_user_getvar(a,"lon",0)
>>> printMinMax(lon,True)
>>> data = rgrid2rcm_Wrap(latin,lonin,xin,lat,lon,0)
>>> printVarSummary(data)
>>> printMinMax(data,True)
>>> .....
>>>
>>> Many thanks!
>>> Sha
>>>
>>>
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu May 22 16:33:27 2014