[gradsusr] Max Function with user defined fields?

[L.B.]

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.

*-------*
William (L.B.) LaForce IV
Meteorologist
www.tornadoinsanity.com

On Fri, May 6, 2016 at 12:29 PM, Jeff Duda <jeffduda319 at gmail.com> wrote:

> L.B.,
> 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.
>
> Jeff
>
> On Fri, May 6, 2016 at 11:12 AM, L.B. <bcbass2989 at gmail.com> wrote:
>
>> 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:
>> 'define stp = cterm*ltermA*hterm*stermA'
>>
>> 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.
>>
>> 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.
>>
>> 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?
>>
>> Here is a snippet of the code that I am trying to work with:
>>
>> ...
>>
>>     if ('hgtprs(z=botzlev)' <= hgt0c)
>>     'wt=0'
>>     endif
>>     if ('hgtdiffZero' > 0 & 'hgtdiffTop' < 0)
>>     'wt=hgtdiffZero/(HGTtop0C-HGT0C)'
>>     endif
>>     if ('hgtdiffTop' >= 0)
>>     'wt=1'
>>     endif
>>     nextZlev = botzlev+1
>>     'define levDiff = hgtprs(z='nextZlev')-hgtprs(z='botzlev')'
>>     'hailsum=hailsum+(wt*E*levDiff)'
>>
>>     botzlev = botzlev + inc
>>     endwhile
>>
>>   'set gxout shaded'
>>   'SHI=hailsum*0.1'
>>   'MEHS=0.1*(2.54*pow(SHI,0.5))'
>> endif
>>
>> 'define hs = max(smth9(MEHS/2.54),t=1,t=24)'
>>
>> *-------*
>> William (L.B.) LaForce IV
>> Meteorologist
>> www.tornadoinsanity.com
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>>
>
>
> --
> Jeff Duda
> Graduate research assistant
> University of Oklahoma School of Meteorology
> Center for Analysis and Prediction of Storms
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20160506/00e1f048/attachment.html