<div dir="ltr">I dont know why I am getting this error. <div><br></div><div><div>Unknown command: u0.5=zinterp(u,height,0.5)</div><div>Unknown command: v0.5=zinterp(v,height,0.6)</div><div>Done.  Newly defined variable interp has 1 v-field.</div>
<div>Done.  Newly defined variable interp has 1.5 u-field.</div><div>DEFINE error:  Name too long; missing &#39;=&#39;</div><div>Done.  Newly defined variable interp has 1.5 v-field.</div><div>DEFINE error:  Name too long; missing &#39;=&#39;</div>
<div>Done.  Newly defined variable interp has 2 u-field.</div><div>Done.  Newly defined variable interp has 2 v-field.</div><div>Done.  Newly defined variable interp has 2.5 u-field.</div><div>DEFINE error:  Name too long; missing &#39;=&#39;</div>
<div>Done.  Newly defined variable interp has 2.5 v-field.</div><div>DEFINE error:  Name too long; missing &#39;=&#39;</div><div>Done.  Newly defined variable interp has 3 u-field.</div><div>Done.  Newly defined variable interp has 3 v-field.</div>
<div>Done.  Newly defined variable interp has 3.5 u-field.</div><div>DEFINE error:  Name too long; missing &#39;=&#39;</div><div>Done.  Newly defined variable interp has 3.5 v-field.</div><div>DEFINE error:  Name too long; missing &#39;=&#39;</div>
<div>Done.  Newly defined variable interp has 4 u-field.</div><div>Done.  Newly defined variable interp has 4 v-field.</div><div>Done.  Newly defined variable interp has 4.5 u-field.</div><div>DEFINE error:  Name too long; missing &#39;=&#39;</div>
<div>Done.  Newly defined variable interp has 4.5 v-field.</div><div>DEFINE error:  Name too long; missing &#39;=&#39;</div><div>Done.  Newly defined variable interp has 5 u-field.</div><div>Done.  Newly defined variable interp has 5 v-field.</div>
<div>Done.  Newly defined variable interp has 5.5 u-field.</div><div>DEFINE error:  Name too long; missing &#39;=&#39;</div><div>Done.  Newly defined variable interp has 5.5 v-field.</div><div>DEFINE error:  Name too long; missing &#39;=&#39;</div>
<div>Done.  Newly defined variable interp has 6 u-field.</div><div>Done.  Newly defined variable interp has 6 v-field.</div><div>Data Request Error:  File number out of range</div><div>  Variable = u0.5</div><div>  Error ocurred at column 6</div>
<div>DEFINE error:  Invalid expression.</div></div><div><br></div><div><br></div><div>Here&#39;s my code:</div><div><br></div><div><br></div><div><div>&quot;u0.5=zinterp(u,height,0.5)&quot;</div><div><br></div><div>&quot;v0.5=zinterp(v,height,0.6)&quot;</div>
<div><br></div><div>&quot;u1=zinterp(u,height,1)&quot;</div><div><br></div><div>&quot;define v1=&quot;zinterp(v,height,1)</div><div><br></div><div>&quot;define u1.5=&quot;zinterp(u,height,1.5)</div><div><br></div><div><br>
</div><div>&quot;define v1.5=&quot;zinterp(v,height,1.5)</div><div><br></div><div><br></div><div>&quot;define u2=&quot;zinterp(u,height,2)</div><div><br></div><div><br></div><div>&quot;define v2=&quot;zinterp(v,height,2)</div>
<div><br></div><div><br></div><div>&quot;define u2.5=&quot;zinterp(u,height,2.5)</div><div><br></div><div><br></div><div>&quot;define v2.5=&quot;zinterp(v,height,2.5)</div><div><br></div><div><br></div><div>&quot;define u3=&quot;zinterp(u,height,3)</div>
<div><br></div><div><br></div><div>&quot;define v3=&quot;zinterp(v,height,3)</div><div><br></div><div><br></div><div>&quot;define u3.5=&quot;zinterp(u,height,3.5)</div><div><br></div><div><br></div><div>&quot;define v3.5=&quot;zinterp(v,height,3.5)</div>
<div><br></div><div>&quot;define u4 =&quot;zinterp(u,height,4)</div><div><br></div><div><br></div><div>&quot;define v4=&quot;zinterp(v,height,4)</div><div><br></div><div><br></div><div>&quot;define u4.5=&quot;zinterp(u,height,4.5)</div>
<div><br></div><div><br></div><div>&quot;define v4.5=&quot;zinterp(v,height,4.5)</div><div><br></div><div><br></div><div>&quot;define u5=&quot;zinterp(u,height,5)</div><div><br></div><div><br></div><div>&quot;define v5=&quot;zinterp(v,height,5)</div>
<div><br></div><div><br></div><div>&quot;define u5.5=&quot;zinterp(u,height,5.5)</div><div><br></div><div><br></div><div>&quot;define v5.5=&quot;zinterp(v,height,5.5)</div><div><br></div><div><br></div><div>&quot;define u6=&quot;zinterp(u,height,6)</div>
<div><br></div><div>&quot;define v6=&quot;zinterp(v,height,6)</div><div><br></div><div>&#39;define umean=(u10+u0.5+u1+u1.5+u2+u2.5+u3.0+u3.5+u4+u4.5+u5+u5.5+u6)/13.0&#39;</div><div>&#39;define vmean=(v10+v0.5+v1+v1.5+v2+v2.5+v3.0+v3.5+v4+v4.5+v5+v5.5+v6)/13.0&#39;</div>
<div>&#39;define ushear=u5-u10&#39;</div><div>&#39;define vshear=v5-v10&#39;</div><div>&#39;define shear=sqrt(ushear*ushear+vshear*vshear)&#39;</div><div>&#39;define umotion=((umean+(7.5/(shear))*vshear))&#39;</div><div>&#39;define vmotion=((vmean-(7.5/(shear))*ushear))&#39;</div>
<div>&#39;define srh1=((u0.5-umotion)*(v10-vmotion)-(u10-umotion)*(v0.5-vmotion))&#39;</div><div>&#39;define srh2=((u1-umotion)*(v0.5-vmotion)-(u0.5-umotion)*(v1-vmotion))&#39;</div><div>&#39;define srh1km=srh1+srh2&#39;</div>
<div><br></div><div>&#39;set gxout shaded&#39;</div><div>&#39;set clevs 100 150 200 250 300 400 500 600 700 800&#39;</div><div>&#39;set ccols 0 14 4  11 5 3  7  12 8 2 6&#39;</div><div>&#39;d srh1km&#39;</div></div><div><br>
</div><div>I dont get what&#39;s wrong? </div><div><br></div><div>Thanks.</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Sat, Sep 7, 2013 at 4:50 PM, Nimrod Micael <span dir="ltr">&lt;<a href="mailto:nmicael@aggies.ncat.edu" target="_blank">nmicael@aggies.ncat.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">Just change it to height coordinates from that script? </div><div class="gmail_extra"><div><div class="h5">
<br><br><div class="gmail_quote">On Sat, Sep 7, 2013 at 4:49 PM, Nimrod Micael <span dir="ltr">&lt;<a href="mailto:nmicael@aggies.ncat.edu" target="_blank">nmicael@aggies.ncat.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">I see... so how would the shear from Surface to 1 km AGL be calculated? I cant find it online. <div><br>

</div><div>Thanks!</div></div><div><div><div class="gmail_extra"><br><br><div class="gmail_quote">On Sat, Sep 7, 2013 at 3:45 PM, Jeff Duda <span dir="ltr">&lt;<a href="mailto:jeffduda319@gmail.com" target="_blank">jeffduda319@gmail.com</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">Because in areas of elevated terrain, the top of the layer in which the shear is being computed is much closer to the surface than 1 km.  Said another way, they use the 900 mb level as a surrogate for the 1 km level.  If the actual surface pressure is close to 900 mb, however, the vertical distance between the surface and the surrogate 1 km level is going to be quite small.  Thus, the computed shear value will be expected to be lower than if you used the 1 km AGL instead.<span><font color="#888888"><div>



<br></div><div>Jeff</div></font></span></div><div><div><div class="gmail_extra"><br><br><div class="gmail_quote">On Sat, Sep 7, 2013 at 1:47 PM, Nimrod Micael <span dir="ltr">&lt;<a href="mailto:nmicael@aggies.ncat.edu" target="_blank">nmicael@aggies.ncat.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">@Jeff: <div><br></div><div>Ok because I find it weird that my 0-1 KM shear values at 16Z in GSO is about 1/2 of the value of the observation from a sounding taken at the time when using that shr1km(u,v) function. However, when I used Rob Hart&#39;s plotskew script to plot the sounding...the ERSH matched closer to the observation at that time. So I am wondering why such a huge error in that.  </div>




</div><div class="gmail_extra"><br><br><div class="gmail_quote">On Sat, Sep 7, 2013 at 2:08 PM, Jeff Duda <span dir="ltr">&lt;<a href="mailto:jeffduda319@gmail.com" target="_blank">jeffduda319@gmail.com</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">If you&#39;re calculating 0-1 km SRH over flat terrain near sea level, their computation is within reason.  However, if you&#39;re on higher terrain or sloped terrain (e.g., the Plains region of the US), you&#39;ll want to use a different method.<span><font color="#888888"><div>





<br>Jeff</div></font></span></div><div><div><div class="gmail_extra"><br><br><div class="gmail_quote">On Sat, Sep 7, 2013 at 11:28 AM, Nimrod Micael <span dir="ltr">&lt;<a href="mailto:nmicael@aggies.ncat.edu" target="_blank">nmicael@aggies.ncat.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">He stated it&#39;s 1 KM MSL.<div><br></div><div><br><div><br></div><div><span style="font-family:arial,sans-serif;font-size:13px">&quot;However, the question of AGL or MSL depends on your model data output -</span><br style="font-family:arial,sans-serif;font-size:13px">






<span style="font-family:arial,sans-serif;font-size:13px">if it compensates the pressure height. GFS IIRC does not do that, so it</span><br style="font-family:arial,sans-serif;font-size:13px"><span style="font-family:arial,sans-serif;font-size:13px">would be MSL.</span><br style="font-family:arial,sans-serif;font-size:13px">






<br style="font-family:arial,sans-serif;font-size:13px"><span style="font-family:arial,sans-serif;font-size:13px">To understand the calculation, you can check the original math from our</span><br style="font-family:arial,sans-serif;font-size:13px">






<span style="font-family:arial,sans-serif;font-size:13px">website:</span><br style="font-family:arial,sans-serif;font-size:13px"><a href="http://www.saakeskus.fi/grads-scripts/63-storm-relative-helicity-0-3km" style="font-family:arial,sans-serif;font-size:13px" target="_blank">http://www.saakeskus.fi/grads-scripts/63-storm-relative-helicity-0-3km</a>&quot;<br>






</div></div><div><br></div><div>So in order to calculate 0-1 KM shear it has to be 0-1 KM AGL?</div><div><br></div><div><br></div><div>Thanks.</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">






On Fri, Sep 6, 2013 at 1:30 PM, Jeff Duda <span dir="ltr">&lt;<a href="mailto:jeffduda319@gmail.com" target="_blank">jeffduda319@gmail.com</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">I don&#39;t use openGrads, so there may be some specifics I&#39;m not aware of.  uVAR and vVAR are just dummy names for the u- and v-components of the wind.  The inaccuracies are probably due to the exact specifics with how the 0-1 km shear is computed in the script.  One big question is how is the 1 km layer wind field being computed?  Is it just 1 km MSL or 1 km AGL?  Those are things you need to find out.  I recommend asking the authors whose email addresses are given on the linked page.<br>







<br>Jeff<br></div><div class="gmail_extra"><br><br><div class="gmail_quote"><div><div>On Fri, Sep 6, 2013 at 12:20 PM, Nimrod Micael <span dir="ltr">&lt;<a href="mailto:nmicael@aggies.ncat.edu" target="_blank">nmicael@aggies.ncat.edu</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><div dir="ltr"><h2 style="color:rgb(147,125,0);font-size:16px;font-family:&#39;bitstream vera sans&#39;,verdana,sans-serf">







<a name="140fa31109422b8a_140fa302f6f692dd_140f9f56f216ea7a_140f9c210c615c15_140f99c64d0502a2_140f945daa5b8bb2_140f4534ed44620b_140f44d21cee3a1f_srh1km_uvar_vvar_"><strong>srh1km</strong>(<em>uVAR,vVAR</em>)</a></h2><div>
What is uVAR and vVAR?  I dont have those variables. </div>



<div><br></div><div><a href="http://opengrads.org/doc/udxt/saakeskus/" target="_blank">http://opengrads.org/doc/udxt/saakeskus/</a><br></div><div><br></div><div>I tried replacing uVar and vVar with u and v wind but I am getting some inaccuracies with the 0-1 KM shear. Can someone help me out?</div>








<div><br></div><div>Thanks. </div><span><font color="#888888"><div><br></div>-- <br>Nimrod Micael
</font></span></div>
<br></div></div>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><span><font color="#888888"><br><br clear="all"><span><font color="#888888"><br>-- <br>Jeff Duda<br>Graduate research assistant<br>University of Oklahoma School of Meteorology<br>Center for Analysis and Prediction of Storms<br>







</font></span></font></span></div><span><font color="#888888">
<br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></font></span></blockquote></div><span><font color="#888888"><br><br clear="all"><div><br></div>-- <br>Nimrod Micael
</font></span></div>
<br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><span><font color="#888888"><div><br></div>-- <br>Jeff Duda<br>Graduate research assistant<br>University of Oklahoma School of Meteorology<br>Center for Analysis and Prediction of Storms<br>




</font></span></div><span><font color="#888888">
</font></span></div></div><span><font color="#888888"><br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></font></span></blockquote></div><span><font color="#888888"><br><br clear="all"><div><br></div>-- <br>Nimrod Micael
</font></span></div>
<br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br>Jeff Duda<br>Graduate research assistant<br>University of Oklahoma School of Meteorology<br>Center for Analysis and Prediction of Storms<br>
</div>
</div></div><br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br>Nimrod Micael
</div>
</div></div></blockquote></div><br><br clear="all"><div><br></div></div></div><span class="HOEnZb"><font color="#888888">-- <br>Nimrod Micael
</font></span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>Nimrod Micael
</div>