[ferret_users] Ferret 5.81 Linux LET

[jimmyc@xxxxxxxxxxx]

Hi all-
I am doing some meteorological verification, computing equitable threat 
score. I want to use multiple thresh-holds. my jnl file is attached.
As of now, I specify the thresh-hold I want to compute the Equitable 
threat score. How can I repeat this calculation for a series of 
thresh-holds and plot the result as a line plot (value vs thresh-hold)?
Can you store the values one at a time in an array variable?
Comments or help would be appreciated.
Thanks

--
James Correia Jr.
Ph.D. Candidate in Agricultural Meteorology
Former President: ISU Graduate Meteorology Club
Dept. of Agronomy, Iowa State University
*********************************************
Email: jimmyc@iastate.edu
Phone: 515-294-9870
Web: http://bruce.agron.iastate.edu/jimmyc/
ISU GMC: http://www.stuorg.iastate.edu/gmc/
*********************************************
"There is a difference between knowing the path and walking the path"
-"Morpheus"

"Wisdom. Strength. Courage. Generosity. Each of us are born with one of
these. We must find the other three inside of us."
from "Into the West"
!obs file stage4
use $1
!model file from WRF
use $2

!define threshhold mask
let mask = IF apcp[i=1:199,j=1:199,k=1,d=1] GE $3 AND apcp[i=1:199,j=1:199,k=1,d=1] LT 500 then 1 else 0
!let mask1 = IF rainnc[i=1:200,j=1:200,l=1]+rainc[i=1:200,j=1:200,l=1] GE $3 then 1
let mask1 = IF e065[i=1:199,j=1:199,l=1,d=2] GE $3 then 1 else 0

!a is hits yes-yes
let a = IF mask EQ mask1 AND mask EQ 1 then 1
let a1 = a[i=1:199@sum,j=1:199@sum]

!b is forecast yes obs no
let b = IF mask1 EQ 1 AND mask EQ 0 then 1
let b1 = b[i=1:199@sum,j=1:199@sum]

!c is forecast no obs yes
let c = IF mask1 EQ 0 AND mask EQ 1 then 1
let c1 = c[i=1:199@sum,j=1:199@sum]

!d is no no
let d = IF mask EQ mask1 AND mask EQ 0 then 1
let d1 = d[i=1:199@sum,j=1:199@sum]

!ar is  a+b * a+c / a +b +c +d
let ar = ((a1+b1)*(a1+c1))/(a1+b1+c1+d1)

let ets = (a1-ar)/(a1+b1+c1-ar)
let bias = (a1+b1)/(a1+c1)

! bias and ETS are calculated
message/continue "ETS and bias are calculated"