<div dir="ltr"><div class="gmail_default" style="font-size:small">Hi Craig,  yes, it's this kind of difference. </div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I did not find any option to specify an alternate earth radius in NCL when reading GRIB files. If someone knows, please kindly let me know.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Thanks,</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"><br clear="all"><div><div class="gmail_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><div class="gmail_quote">On Thu, Aug 2, 2018 at 12:32 PM, 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>
Are you talking about this sort of difference (using a Lambert CONUS 2.5 km grid)?<br>
<br>
Location: Cape Hatteras, NC<br>
lon = -75.528824<br>
lat =  35.250538<br>
<br>
alias mm5="proj +proj=lcc +lat_1=25 +lat_2=25 +lat_0=25 +lon_0=-95 +a=6370 +b=6370 +towgs84=0,0,0 +no_defs -f %.6f"<br>
<br>
mm5 << EOF<br>
-75.528824 35.250538<br>
EOF<br>
<br>
alias ndfd="proj +proj=lcc +lat_1=25 +lat_2=25 +lat_0=25 +lon_0=-95 +x_0=0 +y_0=0 +a=6371.229 +b=6371.229 +no_defs -f %.6f"<br>
<br>
ndfd << EOF<br>
-75.528824 35.250538<br>
EOF<br>
<br>
Results (km):<br>
x = 1791.188308   MM5/WRF spheroid, r = 6370000 m<br>
y = 1274.741690<br>
<br>
x = 1791.533892   NDFD grib octet 3.15 = 6 spheroid, r = 6371229 m<br>
y = 1274.987633<br>
<br>
Difference = 0.424165 km or 424.165 m<br>
<span class="HOEnZb"><font color="#888888"><br>
Craig<br>
</font></span><br>
P.S. I recently learned that you have to use cs2cs, instead of proj, when converting between different datums (not relevant in this case).<br>
<div class="HOEnZb"><div class="h5"><br>
> On Wed, 1 Aug 2018 13:57:55 -0600 Guoqing Ge - NOAA Affiliate <<a href="mailto:guoqing.ge@noaa.gov">guoqing.ge@noaa.gov</a>> wrote:<br>
> <br>
> We have a similar issue when using wgrib2 to read HRRR grib files. In this<br>
> situation, we have identified these inaccurate lat-lon values are due to<br>
> the use of different earth radius by WRF/geogrid.exe and wgrib2<br>
> respectively. Explicitly specify the geogrid-used earth radius (instead of<br>
> the default one) in the command line option, wgrib2 can reproduce the<br>
> original lat-lon values.  Based on similarity, we guess NCL should have the<br>
> same problem.<br>
> <br>
> Is there an option to specify earth radius in NCL when reading GRIB files?<br>
> <br>
> Thanks,<br>
> <br>
> Guoqing<br>
</div></div></blockquote></div><br></div>