Dear Shea<div>Thank you for sharing scripts, it's really helpful. One more question, if I want to calculate the meission amount x(i,j) into emission intensity, should I calculate area manually, as this scipts：</div><div><div>rad = tofloat(4.0*atan(1.0)/180.0)</div><div> re = 6.3712200E+03</div><div> rr = tofloat(re*rad)</div></div><div> area(i,j)=rr*(lat(i+1)-lat(i))*rr*(lon(j+1)-lon(j))*cos(rad*lat(i))</div><div>Or any command can return the area of certain grid directly?</div><div>Thank you</div><div><br></div><div>Ziyuan</div><div><br><br><br><blockquote name="replyContent" class="ReferenceQuote" style="padding-left:5px;margin-left:5px;border-left:#b6b6b6 2px solid;margin-right:0">-----原始邮件-----<br>
<b>发件人:</b><span id="rc_from">"Dennis Shea" <shea@ucar.edu></span><br>
<b>发送时间:</b><span id="rc_senttime">2017-03-05 00:09:00 (星期日)</span><br>
<b>收件人:</b> "Ziyuan Wang" <1301111586@pku.edu.cn><br>
<b>抄送:</b> "ncl-talk@ucar.edu" <ncl-talk@ucar.edu><br>
<b>主题:</b> Re: [ncl-talk] Question about regrid 1) different commands get different results; 2)input emission amount or emission intensity<br><br><div dir="ltr"><div><div><div><div><div><div>I have nor run your script.<br><br>I am not sure what you mean "results are different." The three methods are different. <br><br>(a) linint2 takes the 4 nearest source grid locations and does a bilinear interpolation. <br></div>(b) does some form of 'global' area conservation<br></div>(c) does a local area average<br><br></div>Hence, I am not surprised that they return different results.<br><br>===<br></div>Attached is a script that uses ESMF regridding via bilinear, patch and conserve methods.<br><br></div> If you run the script you will also get a nc fil with the interpolated grids<br><br></div>Good Luck<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Mar 3, 2017 at 2:18 PM, Ziyuan Wang <span dir="ltr"><<a href="mailto:1301111586@pku.edu.cn" target="_blank">1301111586@pku.edu.cn</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div link="#0563C1" vlink="#954F72" lang="ZH-CN"><div class="m_-1248244829172655769WordSection1"><p class="MsoNormal"><span lang="EN-US">Hi ncl-talk members<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">Now I have an emission amount inventory (lat 33:73, lon -12:45), with resolution 0.1*0.1, the emission variable E_HFC_1234 is the emission amount of each grid with unit Mg, as 'HFC.nc' show. I want to interpolate to lower resolution 1*1.<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">1. I tried 3 different commands: "</span><span style="font-size:10.0pt;font-family:Courier" lang="EN-US">linint2", "</span><strong><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:#852d85" lang="EN-US"><a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/area_conserve_remap_Wrap.shtml" target="_blank"><span style="color:#852d85">area_conserve_remap_Wrap</span></a>", and "</span></strong><span lang="EN-US">area_hi2lores_Wrap" But three return results are different (try in regrid.ncl)<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">2) The input "fi" of each command should be the emission amout? or the emission intensity of each grid? Do I need to calculate E_HFC_1234 into emission intensity, by dividing the area of each grid? If so, which command can I use to calculate area?<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">Thanks<br><br>Best<span class="HOEnZb"><font color="#888888"><u></u><u></u></font></span></span></p><span class="HOEnZb"><font color="#888888"><p class="MsoNormal"><span lang="EN-US">Ziyuan<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US"><u></u> <u></u></span></p></font></span></div></div><br>______________________________<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/<wbr>mailman/listinfo/ncl-talk</wbr></a><br>
<br></wbr></blockquote></div><br></div>
</blockquote><br><span><br><p class="MsoNormal" style="font-family: arial, helvetica, clean, sans-serif; font-size: 13px; margin: 0cm 0cm 0pt; line-height: 19px;"><font style="font-size: 12px; font-family: verdana;"><span lang="EN-US">王梓元(WANG Ziyuan)</span></font></p><p class="MsoNormal" style="font-family: arial, helvetica, clean, sans-serif; font-size: 13px; margin: 0cm 0cm 0pt; line-height: 19px;"><span class="yui-tag-span yui-tag" tag="span" style="font-size: 12px; font-family: verdana;">中国北京大学环境科学与工程学院<span lang="EN-US"><br> College of Environmental Science and Engineering <br> Peking University, China<br></span>● 中国北京海淀区颐和园路<span lang="EN-US">5</span>号 邮编：<span lang="EN-US">100871<br> No.5, Yiheyuan Road, Haidian Dist. <br> Beijing 100871, P.R. China <br></span>● 手机<span lang="EN-US"> Cell</span>：<span lang="EN-US"> (+86)132-6999-2003</span></span></p><p class="MsoNormal" style="font-family: arial, helvetica, clean, sans-serif; font-size: 13px; margin: 0cm 0cm 0pt; line-height: 19px;"><font style="font-size: 12px; font-family: verdana;">● 电子邮件<span lang="EN-US"> Email</span>：<a href="mailto:wangziyuan91@163.com" target="_blank">wangziyuan91@163.com</a>, <a href="mailto:ziyuanwang@oku.edu.cn" target="_blank">ziyuanwang@pku.edu.cn</a></font></p></span></div>