<div dir="ltr">Dear ncl-talk,<div><br></div><div>I'm not sure of the best way forward on some code I am writing.</div><div><br></div><div>I have a fortran program that computes storm tracks of extratropical systems in an analysis dataset, as well as from GFS forecasts. The files output are a grads text file of the lat/lon positions and the storm numbers over the selected time periods.</div><div><br></div><div>My issue is that while I have the tracks for all the cyclones, I'm trying to decide the best way to match the GFS initial lat/lon positions to the analysis position so that I am computing the correct track and pressure error. That is, I want to find the lat/lon pair at the 00 hour that is closest to the analysis time for that run. My problem is that my latitude and longitude and pressure values in the analysis and forecast text files are one-dimensional. I have stored the data as a function of (storm number, time or forecast hour). Would it be best to convert the data, such as pressure, to be a function of (lat,lon) and use coordinate subscripting?</div><div><br></div><div>I tried to use closest_val_anyOrder but it is giving me different indices for lat and lon, which would not work for me. Attached is my code. You'll notice where I left off, the exit line on line 172.</div><div><br></div><div>Thanks,<br clear="all"><div><br></div>-- <br><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><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div>Andrew Kren<br>Assistant Scientist<br>CIMAS & NOAA/AOML<br>
325 Broadway, Boulder, CO 80305<br>
<a href="tel:%28303%29%20497-5418" value="+13034975847" target="_blank">(303) 497-5418</a><span><font color="#888888"><br>
</font></span></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
</div></div>