<div dir="ltr"><div>You are correct, there are issues for the ERA_Interim data. Another user had issues with the JRA data set. Even though the original code was written by a NOAA person,  I see no reason why the function should only work with NCEP reanalysis data.<br><br></div><div>I cleaned up your code and used an updated &#39;epflux&#39; function. Nothing substantive was changed BUT I enabled much more printing. This enabled me to do a &#39;side-by-side&#39; with the NCEP results. Basically, all the variables look somewhat consistent *except* the <br><br>Vertical divergence : min=-252582   max=249489<br><br></div><div>which was *MUCH* larger than the NCEP values.<br><br></div><div>That is the current status.<br></div><div><br><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Dec 7, 2016 at 7:24 AM, Rick Brownrigg <span dir="ltr">&lt;<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>&gt;</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><div><div><div><div><div><div>Hi Josh,<br><br></div>I&#39;ve not had a chance to look at your script, but that &quot;Assignment type mismatch...&quot; error message always means one thing:  the variable on the lefthand side (lhs) has been defined elsewhere prior to the assignment, and it has a certain type, shape, and dimensions. The righthand side variable or expression must match in all regards -- type, shape, dimension. In your case, it doesn&#39;t.<br><br></div>If you simply want to reuse a variable name on the lhs, you can use the reassignment operator &quot;:=&quot;, which is an explicit indication by you the programmer that you don&#39;t care that the lhs gets stomped on; eg.:<br><br></div>     foo = (/ 1, 2, 3, 4 /)<br></div>     foo = &quot;A string, not an 1D integer array&quot;     ; assignment YIELDS AN ERROR, not same type, shape, or dimension<br>     foo := &quot;A string, not an 1D integer array&quot;    ; RE-assignment is OK.<br><br></div>If simple reuse is not what you want, try putting in some print statements to make sure you know what lhs and rhs are/should-be.<br><br></div>Hope that helps...<br></div>Rick<br></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On Tue, Dec 6, 2016 at 1:13 PM, Joshua Walston <span dir="ltr">&lt;<a href="mailto:joshwalston@gmail.com" target="_blank">joshwalston@gmail.com</a>&gt;</span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr"><div class="gmail_default" style="font-family:tahoma,sans-serif;color:rgb(0,0,0)"><div class="gmail_default" style="font-family:tahoma,sans-serif;color:rgb(0,0,0)">Hi NCL Users, <br><br></div><div class="gmail_default" style="font-family:tahoma,sans-serif;color:rgb(0,0,0)">I
 am attempting to use the Eliassen-Palm flux function via the 
ncl_contributed functions for the upcoming ncl_6.4 version. There seems 
to be an error in my calculations as the output figures seem odd. My 
script is running fine with no warnings, however the divergence and flux
 vectors seem incorrect. Can anyone see what I am doing wrong?<br><br></div><div class="gmail_default" style="font-family:tahoma,sans-serif;color:rgb(0,0,0)">Additionally,
 I have a 6 hourly data set that I would like to calculate the flux for 
(i.e. would prefer to not have a time average). Is it possible to loop 
through the data and get a return value of type &quot;list&quot; for each time, so
 I can plot the flux at each 6 hour time step. If so, how? When I 
attempt I receive a warning that states:<br><br>&quot;fatal:Assignment type mismatch, right hand side can&#39;t be coerced to type of left hand side&quot;  <br><br></div><div class="gmail_default" style="font-family:tahoma,sans-serif;color:rgb(0,0,0)">I have placed 3 figures (00z, 06z and the average over all times), my script, and the data (ERA40 Interim) here: <br><a href="http://goog_1241068196" target="_blank"><br></a></div><div class="gmail_default" style="font-family:tahoma,sans-serif;color:rgb(0,0,0)"><a href="ftp://pubfiles.dri.edu/pub/jwalston/" target="_blank">ftp://pubfiles.dri.edu/pub/jwa<wbr>lston/</a><br><br></div><div class="gmail_default" style="font-family:tahoma,sans-serif;color:rgb(0,0,0)">Thank you in advance for your help, <br></div>Josh<br><br><br><br><br><br><br><br></div></div>
<br></div></div>______________________________<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/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></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>