<div dir="ltr"><div>L.B.,<br>Is your SigTor parameter defined as a time varying field? If not, then that's why the max() function isn't working the way you want.<br><br></div>Jeff<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, May 6, 2016 at 11:12 AM, L.B. <span dir="ltr"><<a href="mailto:bcbass2989@gmail.com" target="_blank">bcbass2989@gmail.com</a>></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'm not sure if I have my dimension settings wrong or if the max function is not capable of calculating the max of user defined fields. Take for instance the significant tornado parameter:<div>'define stp = cterm*ltermA*hterm*stermA'</div><div><br></div><div>To get the max function to work correctly, I cannot just do max(stp,t=1,t=24), I have to expand out each of the terms in that stp algorithm down to the actual fields available in the model data such as CAPEsfc for instance.</div><div><br></div><div>I have a new script that I created that does some if statements and such that simply cannot be broken down into one single line or their individual fields from the model data and thus, I cannot get the max function to work. It only will calculate one timestep.</div><div><br></div><div>I have tried setting t to vary, such as set t 1 24, etc but nothing works. Is there a way around this, or am I missing something?</div><div><br></div><div>Here is a snippet of the code that I am trying to work with:</div><div><br></div><div><div>...</div><div><br></div><div> if ('hgtprs(z=botzlev)' <= hgt0c)</div><div> 'wt=0'</div><div> endif</div><div> if ('hgtdiffZero' > 0 & 'hgtdiffTop' < 0)</div><div> 'wt=hgtdiffZero/(HGTtop0C-HGT0C)'</div><div> endif</div><div> if ('hgtdiffTop' >= 0)</div><div> 'wt=1'</div><div> endif</div><div> nextZlev = botzlev+1</div><div> 'define levDiff = hgtprs(z='nextZlev')-hgtprs(z='botzlev')'</div><div> 'hailsum=hailsum+(wt*E*levDiff)'</div><div><br></div><div> botzlev = botzlev + inc</div><div> endwhile</div><div><br></div><div> 'set gxout shaded'</div><div> 'SHI=hailsum*0.1'</div><div> 'MEHS=0.1*(2.54*pow(SHI,0.5))'</div><div>endif</div><div><br></div><div>'define hs = max(smth9(MEHS/2.54),t=1,t=24)'</div><div><br></div><div><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><b>-------</b></div>William (L.B.) LaForce IV<br>Meteorologist<div><a href="http://www.tornadoinsanity.com" target="_blank">www.tornadoinsanity.com</a></div></div></div></div></div></div></div></div>
</div></div>
<br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature">Jeff Duda<br>Graduate research assistant<br>University of Oklahoma School of Meteorology<br>Center for Analysis and Prediction of Storms<br></div>
</div>