NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.
To learn about what's going on, see About the Archive Site.
NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.
To learn about what's going on, see About the Archive Site.
Hi Luis...
I'm not sure that "degrees" is a valid unit in this case, since it
really needs to be qualified (as in "degrees_north" or "deg_K").
Regardless, however, if your desire is to strip the units from this
computation, I suggest you use the "noUnit()" method....something like
this
def landSlideH(sr, si, sh, ts, tp, tpr, tpre):
ptot= noUnit(tp + tpr + tpre)
h = noUnit(sr*si*sh) * noUnit(ts + ptot)
return h
(here I assume that the units of "tp", "tpr" and "tpre" are
compatible...if not, then you'll need to use "noUnit()" on
each...likewise, "ts" and "ptot" are compatible...:)
tom
On Tue, Dec 7, 2010 at 6:59 PM, <luis.delcastillo@xxxxxxxxxxxx> wrote:
> For the Tp = no units and The other variables have been reclassified to
> values 0 - 5. Units degrees
> Correo electronico generado desde un BlackBerry de movistar.
>
> -----Original Message-----
> From: Tom Whittaker <whittaker@xxxxxxxx>
> Date: Tue, 07 Dec 2010 18:16:19
> To: Luis Del Castillo<luis.delcastillo@xxxxxxxxxxxx>
> Cc: <idvusers@xxxxxxxxxxxxxxxx>
> Subject: 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
>
--
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
idvusers archives: