[GRASS-dev] computation bio_15 in r.bioclim

Paulo van Breugel p.vanbreugel at gmail.com
Tue Dec 15 05:56:48 PST 2015


Hi Markus,

With r.bioclim if the average rainfall across a year is 0, bio_15 (rainfall
seasonality) will be NA. This should, I think, be 0. Perhaps adding an if
condition in the mapcalc statement on line 505:

grass.mapcalc("$bio = if($precavg==0,0.0,round(100.0 * $precstddev /
$precavg))", bio = outpre + 'bio15', precstddev = precstddev, precavg =
precavg)

Cheers,

Paulo
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-dev/attachments/20151215/b15d19b6/attachment.html>


More information about the grass-dev mailing list