<div dir="ltr"><div class="gmail_default" style="font-size:small">Stuart,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Were these plots generated with NCL? They look like Excel or similar to me.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Your VIC data has rectilinear lat / lon coordinate arrays which means you can easily slice the data by a lat / lon pair, like (45.4375,-96.4375) and do an XY plot.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">With the LIS data, however, the lat / lon arrays are 2D and they contain lots of missing values, so when you created the XY plot across 45.4375,-96.4375 using what looks like Excel, I'm not sure exactly what you were plotting. I see from your lat/lon data that index 0 of lat and 25 of lon is equal to 45.4375 and -96.4375. When I create plots at this time slice, however, they don't look the same to me (see timeseries.000005.png attached).</div><div class="gmail_default" style="font-size:small"><br class="gmail-Apple-interchange-newline"><span style="background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">Note: the VIC data goes from 1/1/1999 to 3/30/1999, but since there are only three LIS files, the data only goes from 1/1/1999 to 1/3/1999. You'll see this on the plots, where the orange curves only have three points.</span><br></div><div class="gmail_default" style="font-size:small"><span style="background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><br></span></div><div class="gmail_default" style="font-size:small">In order to understand your question, we have to better understand the original plots you created, and why you think they should be the same.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I modified both of your codes quite a bit and combined them into one, so you could better see the individual plots and the overlaid plots.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">See attached script and images. Hopefully this will help you figure things out.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">--Mary</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jul 12, 2018 at 9:07 AM, Smith, Stuart <span dir="ltr"><<a href="mailto:smit1770@purdue.edu" target="_blank">smit1770@purdue.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div lang="EN-US" link="#0563C1" vlink="#954F72"><div class="m_2959823101455758343WordSection1"><p class="MsoNormal">Hello,<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">I am using NCL to average simulated output variables over many grid cells, and evaluate the results between hydrologic models (LIS-VIC and VIC). The simulated outputs for both models should be the same. When I compare results of individual grid cells the results are the same, as shown below. However, when I average over the total number of grid cells (703 grid cells), the results are different which was unexpected (also shown below). I believe this is an error with my plotting script.<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">I wanted to ask for guidance on correcting the difference in calculation, as I can’t determine what is causing the differences. In the attached compressed file are examples of the LIS-VIC and VIC NetCDF output files along with the scripts used to calculate the average for each. I hope this file is not too large. I think the information within it is important for analyzing the differences. Thank you for your time. <u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><img width="481" height="289" style="width:5.0104in;height:3.0104in" id="m_2959823101455758343Picture_x0020_1" src="cid:image005.png@01D419D0.978C8E10"> <img width="476" height="286" style="width:4.9583in;height:2.9791in" id="m_2959823101455758343Picture_x0020_3" src="cid:image006.png@01D419D0.978C8E10"><u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Regards,<span class="HOEnZb"><font color="#888888"><u></u><u></u></font></span></p><span class="HOEnZb"><font color="#888888"><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">-Stuart <u></u><u></u></p><p class="MsoNormal"> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <u></u><u></u></p></font></span></div></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>