[GRASS-dev] computation bio_15 in r.bioclim

Markus Metz markus.metz.giswork at gmail.com
Wed Dec 16 01:25:51 PST 2015


Hi Paolo,

On Tue, Dec 15, 2015 at 2:56 PM, Paulo van Breugel
<p.vanbreugel at gmail.com> wrote:
> 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)
>
Fixed in r67163.

Markus M

> Cheers,
>
> Paulo
>


More information about the grass-dev mailing list