<div dir="ltr"><div>A JIRA ticket has been opened:<br><br>Issue <a class="gmail-issue-created-key gmail-issue-link" href="https://jira.ucar.edu/browse/NCL-2764">NCL-2764 - setfileoption: specify earth radius in NCL when reading GRIB files</a> <br><br></div>Possibly:<br><b>setfileoption</b>("grb", "radius_of_earth",  ‎6.3781e6 )   ; if meters
<br>
<b>setfileoption</b>("grb", "radius_of_earth",  ‎6378.1     )      ; if km
</div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Aug 3, 2018 at 11:26 AM, Guoqing Ge - NOAA Affiliate <span dir="ltr"><<a href="mailto:guoqing.ge@noaa.gov" target="_blank">guoqing.ge@noaa.gov</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"><div class="gmail_default" style="font-size:small">HI Craig, thanks a lot for this work around.    </div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I still hope there is a direct option to easily set earth radius in NCL just as in WGRIB2.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Best,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Guoqing</div><div class="gmail_default" style="font-size:small"><br></div></div><div class="gmail_extra"><span class=""><br clear="all"><div><div class="m_-5553722840278939156gmail_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 dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">--</div><div>Guoqing Ge, PhD</div><div>CIRES Research Scientist<br></div><div>at NOAA/ESRL/GSD (Rm 2B502)<span><br></span></div><div dir="ltr"><span style="color:rgb(0,0,0)">325 Broadway, Boulder, CO 80305</span></div><div dir="ltr"><font color="#000000">(303)497-4163<br></font><div><span style="color:rgb(0,0,0)">&</span></div><div><span style="color:rgb(0,0,0)">Thursdays at</span></div><div><font color="#000000">Developmental</font><span style="color:rgb(0,0,0)"> Testbed Center (FL3-1039)</span></div><div>3450 Mitchell Ln, Boulder, CO 80301<span style="color:rgb(0,0,0)"><br></span></div><div>(303)497-8445</div><div><span style="color:rgb(0,0,0)"><br></span></div><div><span style="color:rgb(0,0,0)"><br></span></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
<br></span><span class=""><div class="gmail_quote">On Fri, Aug 3, 2018 at 11:22 AM, Craig Mattocks <span dir="ltr"><<a href="mailto:mattocks@mac.com" target="_blank">mattocks@mac.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Guoqing,<br>
<br>
I assume the NCL developers could make this setting available through their API.<br>
<br>
BTW, there are Fortran wrappers available for proj.4:<br>
<br>
<a href="https://github.com/mhagdorn/proj4-fortran" rel="noreferrer" target="_blank">https://github.com/mhagdorn/pr<wbr>oj4-fortran</a><br>
<br>
or you could use a Fortran system call, or an NCL systemfunc call:<br>
<br>
<a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/systemfunc.shtml" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/Docum<wbr>ent/Functions/Built-in/systemf<wbr>unc.shtml</a><br>
<br>
Best,<br>
Craig<br>
<br>
</blockquote></div><br></span></div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">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</a><br>
<br></blockquote></div><br></div>