Re: [idvusers] IDV formulas

Luis -- What are the units of each of your fields?

tom

On Tue, Dec 7, 2010 at 2:57 PM, Luis Del Castillo
<luis.delcastillo@xxxxxxxxxxxx> wrote:
> Tom after you send me your email I did this:
>
> (newUnit(Sr, "ries", "R")*newUnit(SI, "ries", "R")*newUnit(Sh, "ries", 
> "R"))*(newUnit(Ts, "ries", "R") + (newUnit(P, "ries", "R") + newUnit(P1, 
> "ries", "R") + newUnit(P2, "ries", "R")))
>
> This is the only way that I saw a correct evaluation of the formula and I see 
> it when I plot each layer values as Value plots.
> If I just run (Sr*SI*Sh)*(Ts + Tp) I get a bad value something like 
> (4*5*2)*(3 + 50)  = 4 and since I use precipitation accumulated 72 hr I must 
> get a high value.
>
> So could be true about the incompatibility between Ts and Tp so am wondering 
> what could be the best way to reassign unit here because the output like 250 
> or 450 is going to be in hight risk, low risk.
>
> Bay the way I am trying to write a procedure like this one
>
> def     landSlideH(sr, si, sh, ts, tp, tpr, tpre):
>        ptot= tp + tpr + tpre
>        h = (sr*si*sh)*(ts + ptot)
>        return h
>
> Can I use newUnit inside it?
>
> Thanks
> Luis
>
> -----Original Message-----
> From: Tom Whittaker [mailto:whittaker@xxxxxxxx]
> Sent: martes, 07 de diciembre de 2010 03:23 p.m.
> To: Luis Del Castillo
> Cc: idvusers@xxxxxxxxxxxxxxxx
> Subject: Re: [idvusers] IDV formulas
>
> Hi Luis...
>
> It would seem to me that if there is a unit compatibility issue with
> this formla:
>
> equation is  H = (Sr*SI*Sh)*(Ts + Tp) .
>
> it would have to be between Ts and Tp.  Is that true?  Or is there
> some other issue here?
>
> Best wishes,
>
> tom
>
> --
> Tom Whittaker
> University of Wisconsin-Madison
> Space Science & Engineering Center (SSEC)
> Cooperative Institute for Meteorological Satellite Studies (CIMSS)
> 1225 W. Dayton Street
> Madison, WI  53706  USA
> ph: +1 608 262 2759
>



-- 
Tom Whittaker
University of Wisconsin-Madison
Space Science & Engineering Center (SSEC)
Cooperative Institute for Meteorological Satellite Studies (CIMSS)
1225 W. Dayton Street
Madison, WI  53706  USA
ph: +1 608 262 2759



  • 2010 messages navigation, sorted by:
    1. Thread
    2. Subject
    3. Author
    4. Date
    5. ↑ Table Of Contents
  • Search the idvusers archives: