[GRASS-dev] computation bio_15 in r.bioclim

GRASS GIS GRASS GIS developer
1

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

2

Hi Paolo,

On Tue, Dec 15, 2015 at 2:56 PM, Paulo van Breugel
<p.vanbreugel@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