<div dir="ltr">Hi Tao:<div><br></div><div>I apologize, but I can't answer right now as I am traveling.</div><div><br></div><div>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.</div><div><br></div><div><div>;---Write wrf variables lat/lon to file for <a href="http://ncltalk.nc">ncltalk.nc</a> </div><div> system("rm -f <a href="http://file_for_ncltalk.nc">file_for_ncltalk.nc</a>")</div><div> fout = addfile("<a href="http://file_for_ncltalk.nc">file_for_ncltalk.nc</a>","c")</div><div> fout->xlat = xlat</div><div> fout->xlong = xlong</div></div><div><br></div><div>Barry</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jan 17, 2017 at 7:14 PM, Tao Lu <span dir="ltr"><<a href="mailto:hakufu.asano@gmail.com" target="_blank">hakufu.asano@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello Barry:<div><br></div><div>Thank you very much for your reply.</div><div>I am very appreciated that you shared your codes with us.</div><div><br></div><div>Some editions have been done in my situation.</div><div>However, I could not get contour after <span style="font-size:12.8px">grid = obj_anal_ic(lon_o,lat_o,obs_o,</span><span style="font-size:12.8px"><wbr>xlon_m,xlat_m,rscan,False).</span></div><div><span style="font-size:12.8px">I attached the edited code and run information, please help me check that.</span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">I also have some questions about the method for regridding:</span></div><div><span style="font-size:12.8px">1. There <span style="background-color:rgb(255,255,255)">is a item named "<font color="#000000">rightmost dimension or leftmost dimension</font>" in</span></span></div><div><span style="font-size:12.8px;background-color:rgb(255,255,255)"> <a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/obj_anal_ic.shtml" target="_blank">http://www.ncl.ucar.edu/<wbr>Document/Functions/Built-in/<wbr>obj_anal_ic.shtml</a></span></div><div><span style="font-size:12.8px;background-color:rgb(255,255,255)"> I googled it but have no ideal.</span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">2. Also according to the webpage in 1. </span><span style="font-size:12.8px">obs_o</span>@_FillValue will be ignored.</div><div> Does it mean that missing value won't be accounted for when interpolation?</div><div><br></div><div>3. "cband radar observations: min=0 max=220" before regridding</div><div> "var_zoom min=0 max=23.5354" after regridding</div><div> After regridding the radar density is much smaller than original maximum.</div><div> Is it normal? Maybe as you said rscan is important.</div><div><br></div><div>4. When you use</div><div> <span style="font-size:12.8px">grid = obj_anal_ic(lon_o,lat_o,obs_o,</span><span style="font-size:12.8px"><wbr>xlon_m,xlat_m,rscan,False)</span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px"> You calculate this</span></div><div><span style="font-size:12.8px"> ;---Get slice of lat/lon data at halfway points.</span></div><div><span style="font-size:12.8px"> nlat = dims2d(0)</span></div><div><span style="font-size:12.8px"> nlon = dims2d(1)</span></div><div><span style="font-size:12.8px"> xlat_m = xlat2d_m(:,nlon/2)</span></div><div><span style="font-size:12.8px"> xlon_m = xlon2d_m(nlat/2,:)</span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px"> Is it reasonable to calculate lat/lon at a quarter points?</span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">5. If use halfway points to regrid radar data, </span></div><div><span style="font-size:12.8px"><img src="cid:ii_159ad5c3151ea464" alt="Inline image 1" width="247" height="128" style="margin-right:0px"><br></span></div><div><span style="font-size:12.8px"> We will get a n(xlat_m) x n(xlon_m) regridded array. </span></div><div><span style="font-size:12.8px"> So, what will the array's lat/lon will be?</span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">6. This method is to interpolate radar data to wrf grid.</span></div><div><span style="font-size:12.8px"> Is the mean average method more suitable for radar data? Because radar have</span></div><div><span style="font-size:12.8px"> so many data points unlike rain gauge data?</span></div><div><span style="font-size:12.8px"><br></span></div><div><br></div><div>My ascii data converted from radar is too large, if you need it I will send it to you.</div><div>The data is something like below:</div><div><div># crain_composite_data:rain_<wbr>density:lat:lon:ndata=2406400:</div><div># 2015-09-09_00:00:00</div><div>-9999 47.32916666666667 137.01874999999998 </div><div>-9999 47.32916666666667 137.03125 </div><div>-9999 47.32916666666667 137.04375 </div><div>-9999 47.32916666666667 137.05625 </div><div>-9999 47.32916666666667 137.06875 </div><div>-9999 47.32916666666667 137.08124999999998 </div><div><br></div></div><div><br></div><div>Waiting for your reply.</div><div>Thank you very much.</div><div><br></div><div>Tao</div><div><span style="font-size:12.8px"> </span></div></div><div class="gmail_extra"><div><div class="h5"><br><div class="gmail_quote">On Tue, Jan 17, 2017 at 4:05 AM, Barry Lynn <span dir="ltr"><<a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Tao:<div><br></div><div>Here, I gridded the radar data to the WRF grid.</div><div><br></div><div>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).</div><div><br></div><div> grid = obj_anal_ic(lon_o,lat_o,obs_o,<wbr>xlon_m,xlat_m,rscan,False)<br></div><div><br></div><div>Please write if you have questions.</div><div><br></div><div>Barry</div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="m_3189478359095484208h5">On Mon, Jan 16, 2017 at 9:11 AM, Tao Lu <span dir="ltr"><<a href="mailto:hakufu.asano@gmail.com" target="_blank">hakufu.asano@gmail.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="m_3189478359095484208h5"><div dir="ltr">Dear ncl-talk<div><br></div><div>Hello.</div><div><br></div><div>I am new to regridding, if possible please give me some help and advice.</div><div><br></div><div><font size="4">What I am doing?: Compare two kind of grid data</font></div><div><b>I have two kind of data.</b></div><div>1. radar data: </div><div>resolution: 0.5 minute x 0.75 minute, as seen not regular grids</div><div>file format: ascii</div><div><img src="cid:ii_159a607a710bddc8" alt="Inline image 1" width="378" height="472" style="margin-right:0px"></div><div>2. output data from WRF:</div><div><div> dx = 30000,</div><div> dy = 30000,</div><div> map_proj = 'lambert',</div><div> ref_lat = 36.00,</div><div> ref_lon = 140.00,</div><div> truelat1 = 30.0,</div><div> truelat2 = 60.0,</div><div> stand_lon = 140.0,</div></div><div>file fomat: netcdf</div><div><font size="4"><img src="cid:ii_159a611b801e2cea" alt="Inline image 3" width="325" height="320" style="margin-right:0px"><br clear="all"></font></div><div><div><font size="4"><br></font></div><div><b>I have to regrid the two data to each other or to new grid.</b></div><div>1. Regrid wrf data grid to radar data grid?</div><div>2. Regrid radar data to wrf data?</div><div>3. Regrid both to another data?</div><div>Which is easiest?</div><div><br></div><div>I am really new to this, please give me some guide line, I will further operation on my work.</div><div><br></div><div>Thank you very much for your kind help.</div><div><br></div><div>Tao</div><div><br></div><div><br></div>-- <div class="m_3189478359095484208m_3288249013606581598m_-6628643434214747579gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div style="font-size:12.8px"><span style="font-size:small;font-family:arial">******************************</span><span style="font-size:small;font-family:arial"><wbr>************************</span><br></div><div><span style="font-size:small;font-family:arial">Tao Lu (M. Eng)</span></div><div><span style="font-family:arial,helvetica,sans-serif;font-size:12.8px">Laboratory of River Engineering and Hydrology,</span><br style="font-family:arial;font-size:small"><span style="font-size:small;font-family:arial">Dept. of Civil and Environmental Engineering,</span><br style="font-family:arial;font-size:small"><span style="font-family:arial,helvetica,sans-serif;font-size:12.8px">Graduate School of Science and Engineering,</span><br style="font-family:arial;font-size:small"><span style="font-size:small;font-family:arial">Chuo University</span><br style="font-family:arial;font-size:small"><span style="font-size:small;font-family:arial"> </span><br style="font-family:arial;font-size:small">1-13-27<span style="font-size:small;font-family:arial">,Kasuga,Bunkyo-ku,Tokyo</span><br style="font-family:arial;font-size:small"><span style="font-size:small;font-family:arial">112-8551, Japan</span><br style="font-family:arial;font-size:small"><span style="font-size:small;font-family:arial">TEL: </span>03-3817-1805<font size="2" face="arial, helvetica, sans-serif"><span lang="EN-US">; </span></font><span style="font-family:arial,helvetica,sans-serif;font-size:small">Phone: 070-2188-7509</span><br style="font-family:arial;font-size:small"><span style="font-size:12.8px">Email1: </span><a href="mailto:mail%3Amet.yamos@gmail.com" style="font-size:12.8px;color:rgb(17,85,204)" target="_blank">hakufu.asano@gmail.com</a><br style="font-size:12.8px"><span style="font-size:12.8px">Email2: </span><a href="mailto:mail%3Ayamoto@civil.chuo-u.ac.jp" style="font-size:12.8px;color:rgb(17,85,204)" target="_blank">lutao@civil.chuo-u.ac.<wbr>jp</a><br style="font-family:arial;font-size:small"><span style="font-size:small;font-family:arial">******************************</span><span style="font-size:small;font-family:arial"><wbr>***********************</span><span style="font-size:12.8px"><br></span></div></div></div></div></div></div></div></div>
</div></div>
<br></div></div>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="m_3189478359095484208m_3288249013606581598gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Barry H. Lynn, Ph.D<div><div>Senior Lecturer,</div><div><div><span style="color:rgb(136,136,136)">The Institute of the Earth Science, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">The Hebrew University of Jerusalem, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Givat Ram, Jerusalem 91904, Israel </span><br style="color:rgb(136,136,136)"></div><span style="color:rgb(136,136,136)">Tel: 972 547 231 170</span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Fax: (972)-25662581</span></div></div><div><span style="color:rgb(136,136,136)"><br></span></div><div>C.E.O, Weather It Is, LTD<br>Weather and Climate Focus<br><a href="http://weather-it-is.com" target="_blank">http://weather-it-is.com</a><br>Jerusalem, Israel<br>Local: 02 930 9525<br>Cell: 054 7 231 170<br>Int-IS: x972 2 930 9525<br>US <a href="tel:(914)%20432-3108" value="+19144323108" target="_blank">914 432 3108</a><br></div></div></div>
</div>
</blockquote></div><br><br clear="all"><div><br></div></div></div>-- <br><div class="m_3189478359095484208gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px">******************************<wbr>************************</span><br></div><div style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px">盧 涛 (ル タオ) </span><span style="font-family:arial,helvetica,sans-serif;font-size:small">修士課程2年</span><br style="font-size:12.8000001907349px"><p style="font-size:12.8000001907349px;margin:0px"><font size="2" face="arial, helvetica, sans-serif">〒<span lang="EN-US">112-8551 </span>東京都文京区春日<span lang="EN-US">1-13-27</span></font></p><p style="font-size:12.8000001907349px;margin:0px"><font size="2" face="arial, helvetica, sans-serif">中央大学理工学研究科都市環境学専攻</font></p><p style="font-size:12.8000001907349px;margin:0px"><font size="2" face="arial, helvetica, sans-serif">河川・水文<span lang="EN-US"><span lang="EN-US"><span lang="EN-US">研究</span></span><span lang="EN-US"><span lang="EN-US">室</span></span></span></font><span style="font-family:arial,helvetica,sans-serif;font-size:small">（山田正教授）</span></p></div><div style="font-size:12.8000001907349px"><p style="font-size:12.8000001907349px;margin:0px"><font size="2" face="arial, helvetica, sans-serif"><span lang="EN-US"><br></span></font></p><p style="font-size:12.8000001907349px;margin:0px"><font size="2" face="arial, helvetica, sans-serif"><span lang="EN-US">TEL: 03-3817-3406; </span></font><span style="font-family:arial,helvetica,sans-serif;font-size:small">Phone: 070-2188-7509</span></p>Email1: <a href="mailto:mail%3Amet.yamos@gmail.com" style="color:rgb(17,85,204);font-size:12.8000001907349px" target="_blank">hakufu.asano@gmail.com</a><br style="font-size:12.8000001907349px">Email2: <a href="mailto:mail%3Ayamoto@civil.chuo-u.ac.jp" style="color:rgb(17,85,204);font-size:12.8000001907349px" target="_blank">lutao@civil.chuo-u.ac.<wbr>jp</a><br style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px">******************************<wbr>*************************</span></div><span class=""><div style="font-size:12.8000001907349px"><span style="font-size:small;font-family:arial">******************************</span><span style="font-size:small;font-family:arial"><wbr>************************</span><br></div><div><span style="font-size:small;font-family:arial">Tao Lu (M. Eng)</span></div><div><span style="font-family:arial,helvetica,sans-serif;font-size:12.8px">Laboratory of River Engineering and Hydrology,</span><br style="font-family:arial;font-size:small"><span style="font-size:small;font-family:arial">Dept. of Civil and Environmental Engineering,</span><br style="font-family:arial;font-size:small"><span style="font-family:arial,helvetica,sans-serif;font-size:12.8px">Graduate School of Science and Engineering,</span><br style="font-family:arial;font-size:small"><span style="font-size:small;font-family:arial">Chuo University</span><br style="font-family:arial;font-size:small"><span style="font-size:small;font-family:arial"> </span><br style="font-family:arial;font-size:small">1-13-27<span style="font-size:small;font-family:arial">,Kasuga,Bunkyo-ku,Tokyo</span><br style="font-family:arial;font-size:small"><span style="font-size:small;font-family:arial">112-8551, Japan</span><br style="font-family:arial;font-size:small"><span style="font-size:small;font-family:arial">TEL: </span>03-3817-1805<font size="2" face="arial, helvetica, sans-serif"><span lang="EN-US">; </span></font><span style="font-family:arial,helvetica,sans-serif;font-size:small">Phone: 070-2188-7509</span><br style="font-family:arial;font-size:small"><span style="font-size:12.8px">Email1: </span><a href="mailto:mail%3Amet.yamos@gmail.com" style="font-size:12.8px;color:rgb(17,85,204)" target="_blank">hakufu.asano@gmail.com</a><br style="font-size:12.8px"><span style="font-size:12.8px">Email2: </span><a href="mailto:mail%3Ayamoto@civil.chuo-u.ac.jp" style="font-size:12.8px;color:rgb(17,85,204)" target="_blank">lutao@civil.chuo-u.ac.<wbr>jp</a><br style="font-family:arial;font-size:small"><span style="font-size:small;font-family:arial">******************************</span><span style="font-size:small;font-family:arial"><wbr>***********************</span><span style="font-size:12.8px"><br></span></div></span></div></div></div></div></div></div></div></div></div></div></div></div>
</div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Barry H. Lynn, Ph.D<div><div>Senior Lecturer,</div><div><div><span style="color:rgb(136,136,136)">The Institute of the Earth Science, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">The Hebrew University of Jerusalem, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Givat Ram, Jerusalem 91904, Israel </span><br style="color:rgb(136,136,136)"></div><span style="color:rgb(136,136,136)">Tel: 972 547 231 170</span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Fax: (972)-25662581</span></div></div><div><span style="color:rgb(136,136,136)"><br></span></div><div>C.E.O, Weather It Is, LTD<br>Weather and Climate Focus<br><a href="http://weather-it-is.com" target="_blank">http://weather-it-is.com</a><br>Jerusalem, Israel<br>Local: 02 930 9525<br>Cell: 054 7 231 170<br>Int-IS: x972 2 930 9525<br>US 914 432 3108<br></div></div></div>
</div>