# | Line 289 | Line 289 | FVECT sd, vd, vu; | |
---|---|---|
289 | d = DOT(sd,vu)/sqrt(1.0-d*d); | |
290 | if (d >= 1.0) | |
291 | tau = 0.0; | |
292 | + | else if (d <= -1.0) |
293 | + | tau = 180.0; |
294 | else | |
295 | tau = acos(d) * (180./PI); | |
296 | return( exp( sigma*( (35.2 - tau*.31889 - 1.22*exp(-.22222*tau))*1e-3 |
– | Removed lines |
+ | Added lines |
< | Changed lines |
> | Changed lines |