<div dir="ltr"><div>L.B.,<br>Set your dimension environment to a time varying span at the time you define the field. If you're building up the timeseries of fields as you go (which it sounds like you are), you'll have to create a new 2D field at each iteration that contains the temporal maximum of the field at all times leading up to that time. Perhaps the best way to do that is to simply plot the max between two fields - the current field and the field representing the max values leading up to the previous hour. Since GrADS doesn't have a function to do that, you'll have to use a clever combination of const and maskout to do that. You'll want to verify this for yourself, but your mask field should probably be the difference field between your current sig tor values and the previous max values.<br><br></div>Jeff<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, May 6, 2016 at 12:14 PM, 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 don't believe so no. How would I go about doing that? That makes more sense as to why its not working though. I just do an iteration and jump through timesteps and plot the values, so when I define new fields they are a fixed time.</div><div class="gmail_extra"><span class=""><br clear="all"><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>
<br></span><div><div class="h5"><div class="gmail_quote">On Fri, May 6, 2016 at 12:29 PM, Jeff Duda <span dir="ltr"><<a href="mailto:jeffduda319@gmail.com" target="_blank">jeffduda319@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"><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"><div><div>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></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">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></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" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><span><font color="#888888"><br><br clear="all"><br>-- <br><div>Jeff Duda<br>Graduate research assistant<br>University of Oklahoma School of Meteorology<br>Center for Analysis and Prediction of Storms<br></div>
</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" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br></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>