[ncl-talk] About regridding from
Dennis Shea
shea at ucar.edu
Tue Jan 17 11:15:12 MST 2017
[1]
A general rule for regridding is to go from high-resolution to
low-resolution. *If* the variable is 'smoothly varying' virtually any
interpolation method (eg: bilinear) will do. If the data is fractal (eg:
hourly/3hrly/daily precipitation), then other methods may be more
appropriate.
Interpolating from low-resolution to high-resolution, does *not* provide
any additional information. However, depending upon your objective and use,
this may be what you need to do.
[2]
Look at ESMF Regridding: http://www.ncl.ucar.edu/Applications/ESMF.shtml
Look at Example 21
[a] Take the radar data and make one-dimensional
lat_radar/lon_radar/value_radar arrays. These are your input ('Source').
Create a weight file (say)
Opt at WgtFileName= "WGT.Radar_to_WRF.nc" ; default is "weights_file.nc"
[b] The target grid is defined by the WRF: XLAT, XLONG variables. There
are numerous examples on the ESMF page to illustrate how to do this.
[c] Run the script
[d] The weight file can be used subsequently for other regriddings.
=============
NOTE:
It is your responsibility to look at and study the ESMF URL. Read it
carefully before proceeding. Responding to ncl-talk questions, especially,
ESMF can be quite time consuming and we have our own jobs to do.
Good Luck
On Tue, Jan 17, 2017 at 10:41 AM, Barry Lynn <barry.h.lynn at gmail.com> wrote:
> Hi Tao:
>
> I apologize, but I can't answer right now as I am traveling.
>
> However, in order to provide ncl-talkers with enough information/data, you
> can write out the data grid (lat/lon) you are trying to plot on from WRF,
> and upload one time step of it (or two). You should also upload the radar
> data.
>
> ;---Write wrf variables lat/lon to file for ncltalk.nc
>
>
> system("rm -f file_for_ncltalk.nc")
> fout = addfile("file_for_ncltalk.nc","c")
> fout->xlat = xlat
> fout->xlong = xlong
>
> Barry
>
> On Tue, Jan 17, 2017 at 7:14 PM, Tao Lu <hakufu.asano at gmail.com> wrote:
>
>> Hello Barry:
>>
>> Thank you very much for your reply.
>> I am very appreciated that you shared your codes with us.
>>
>> Some editions have been done in my situation.
>> However, I could not get contour after grid =
>> obj_anal_ic(lon_o,lat_o,obs_o,xlon_m,xlat_m,rscan,False).
>> I attached the edited code and run information, please help me check that.
>>
>> I also have some questions about the method for regridding:
>> 1. There is a item named "rightmost dimension or leftmost dimension" in
>> http://www.ncl.ucar.edu/Document/Functions/Built-in/obj_anal_ic.shtml
>> I googled it but have no ideal.
>>
>> 2. Also according to the webpage in 1. obs_o at _FillValue will be ignored.
>> Does it mean that missing value won't be accounted for when
>> interpolation?
>>
>> 3. "cband radar observations: min=0 max=220" before regridding
>> "var_zoom min=0 max=23.5354" after regridding
>> After regridding the radar density is much smaller than original
>> maximum.
>> Is it normal? Maybe as you said rscan is important.
>>
>> 4. When you use
>> grid = obj_anal_ic(lon_o,lat_o,obs_o,xlon_m,xlat_m,rscan,False)
>>
>> You calculate this
>> ;---Get slice of lat/lon data at halfway points.
>> nlat = dims2d(0)
>> nlon = dims2d(1)
>> xlat_m = xlat2d_m(:,nlon/2)
>> xlon_m = xlon2d_m(nlat/2,:)
>>
>> Is it reasonable to calculate lat/lon at a quarter points?
>>
>> 5. If use halfway points to regrid radar data,
>> [image: Inline image 1]
>> We will get a n(xlat_m) x n(xlon_m) regridded array.
>> So, what will the array's lat/lon will be?
>>
>> 6. This method is to interpolate radar data to wrf grid.
>> Is the mean average method more suitable for radar data? Because
>> radar have
>> so many data points unlike rain gauge data?
>>
>>
>> My ascii data converted from radar is too large, if you need it I will
>> send it to you.
>> The data is something like below:
>> # crain_composite_data:rain_density:lat:lon:ndata=2406400:
>> # 2015-09-09_00:00:00
>> -9999 47.32916666666667 137.01874999999998
>> -9999 47.32916666666667 137.03125
>> -9999 47.32916666666667 137.04375
>> -9999 47.32916666666667 137.05625
>> -9999 47.32916666666667 137.06875
>> -9999 47.32916666666667 137.08124999999998
>>
>>
>> Waiting for your reply.
>> Thank you very much.
>>
>> Tao
>>
>>
>> On Tue, Jan 17, 2017 at 4:05 AM, Barry Lynn <barry.h.lynn at gmail.com>
>> wrote:
>>
>>> Hi Tao:
>>>
>>> Here, I gridded the radar data to the WRF grid.
>>>
>>> You will need to adapt what I did to your situation, but the most
>>> important step is to decide on the degree of smoothing (which is set by
>>> rscan).
>>>
>>> grid = obj_anal_ic(lon_o,lat_o,obs_o,xlon_m,xlat_m,rscan,False)
>>>
>>> Please write if you have questions.
>>>
>>> Barry
>>>
>>> On Mon, Jan 16, 2017 at 9:11 AM, Tao Lu <hakufu.asano at gmail.com> wrote:
>>>
>>>> Dear ncl-talk
>>>>
>>>> Hello.
>>>>
>>>> I am new to regridding, if possible please give me some help and advice.
>>>>
>>>> What I am doing?: Compare two kind of grid data
>>>> *I have two kind of data.*
>>>> 1. radar data:
>>>> resolution: 0.5 minute x 0.75 minute, as seen not regular
>>>> grids
>>>> file format: ascii
>>>> [image: Inline image 1]
>>>> 2. output data from WRF:
>>>> dx = 30000,
>>>> dy = 30000,
>>>> map_proj = 'lambert',
>>>> ref_lat = 36.00,
>>>> ref_lon = 140.00,
>>>> truelat1 = 30.0,
>>>> truelat2 = 60.0,
>>>> stand_lon = 140.0,
>>>> file fomat: netcdf
>>>> [image: Inline image 3]
>>>>
>>>> *I have to regrid the two data to each other or to new grid.*
>>>> 1. Regrid wrf data grid to radar data grid?
>>>> 2. Regrid radar data to wrf data?
>>>> 3. Regrid both to another data?
>>>> Which is easiest?
>>>>
>>>> I am really new to this, please give me some guide line, I will further
>>>> operation on my work.
>>>>
>>>> Thank you very much for your kind help.
>>>>
>>>> Tao
>>>>
>>>>
>>>> --
>>>> ******************************************************
>>>> Tao Lu (M. Eng)
>>>> Laboratory of River Engineering and Hydrology,
>>>> Dept. of Civil and Environmental Engineering,
>>>> Graduate School of Science and Engineering,
>>>> Chuo University
>>>>
>>>> 1-13-27,Kasuga,Bunkyo-ku,Tokyo
>>>> 112-8551, Japan
>>>> TEL: 03-3817-1805; Phone: 070-2188-7509
>>>> Email1: hakufu.asano at gmail.com <mail%3Amet.yamos at gmail.com>
>>>> Email2: lutao at civil.chuo-u.ac.jp <mail%3Ayamoto at civil.chuo-u.ac.jp>
>>>> *****************************************************
>>>>
>>>> _______________________________________________
>>>> ncl-talk mailing list
>>>> ncl-talk at ucar.edu
>>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>>>
>>>
>>>
>>> --
>>> Barry H. Lynn, Ph.D
>>> Senior Lecturer,
>>> The Institute of the Earth Science,
>>> The Hebrew University of Jerusalem,
>>> Givat Ram, Jerusalem 91904, Israel
>>> Tel: 972 547 231 170
>>> Fax: (972)-25662581
>>>
>>> C.E.O, Weather It Is, LTD
>>> Weather and Climate Focus
>>> http://weather-it-is.com
>>> Jerusalem, Israel
>>> Local: 02 930 9525
>>> Cell: 054 7 231 170
>>> Int-IS: x972 2 930 9525
>>> US 914 432 3108 <(914)%20432-3108>
>>>
>>
>>
>>
>> --
>> ******************************************************
>> 盧 涛 (ル タオ) 修士課程2年
>>
>> 〒112-8551 東京都文京区春日1-13-27
>>
>> 中央大学理工学研究科都市環境学専攻
>>
>> 河川・水文研究室(山田正教授)
>>
>>
>> TEL: 03-3817-3406; Phone: 070-2188-7509
>> Email1: hakufu.asano at gmail.com <mail%3Amet.yamos at gmail.com>
>> Email2: lutao at civil.chuo-u.ac.jp <mail%3Ayamoto at civil.chuo-u.ac.jp>
>> *******************************************************
>> ******************************************************
>> Tao Lu (M. Eng)
>> Laboratory of River Engineering and Hydrology,
>> Dept. of Civil and Environmental Engineering,
>> Graduate School of Science and Engineering,
>> Chuo University
>>
>> 1-13-27,Kasuga,Bunkyo-ku,Tokyo
>> 112-8551, Japan
>> TEL: 03-3817-1805; Phone: 070-2188-7509
>> Email1: hakufu.asano at gmail.com <mail%3Amet.yamos at gmail.com>
>> Email2: lutao at civil.chuo-u.ac.jp <mail%3Ayamoto at civil.chuo-u.ac.jp>
>> *****************************************************
>>
>
>
>
> --
> Barry H. Lynn, Ph.D
> Senior Lecturer,
> The Institute of the Earth Science,
> The Hebrew University of Jerusalem,
> Givat Ram, Jerusalem 91904, Israel
> Tel: 972 547 231 170
> Fax: (972)-25662581
>
> C.E.O, Weather It Is, LTD
> Weather and Climate Focus
> http://weather-it-is.com
> Jerusalem, Israel
> Local: 02 930 9525
> Cell: 054 7 231 170
> Int-IS: x972 2 930 9525
> US 914 432 3108 <(914)%20432-3108>
>
> _______________________________________________
> 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/20170117/e59119ab/attachment-0001.html
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 44666 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170117/e59119ab/attachment-0003.png
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 4669 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170117/e59119ab/attachment-0004.png
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 392777 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170117/e59119ab/attachment-0005.png
More information about the ncl-talk
mailing list