[GRASS-dev] seemingly random errors in r.mapcalc

[Paulo van Breugel]

Hi,

In the last few revisions of GRASS GIS 7.0 (on Ubuntu 13.04) I am
experiencing this seemingly random errors with r.mapcalc. For example:

Running the following r.mapcalc computation gave me an ERROR:

GRASS 7.0.svn (AEA):~ > r.mapcalc --overwrite "glc_influence =
if(IIASA_crops==100,8.0, if((glc_change5 + glc_change6)==100, 4.0*IIASArev
+ 8.0*IIASA_crops/100.0, (glc_change1*IIASArev + 2.0*glc_change2*IIASArev +
3.0*glc_change3*IIASArev + 4.0*glc_change4*IIASArev + glc_div +
8.0*IIASA_crops) / 100.0))" --overwrite
ERROR: Error reading raster data
[Raster MASK present]

Rerunning the same, and it goes through fine:

GRASS 7.0.svn (AEA):~ > r.mapcalc --overwrite "glc_influence =
if(IIASA_crops==100,8.0, if((glc_change5 + glc_change6)==100, 4.0*IIASArev
+ 8.0*IIASA_crops/100.0, (glc_change1*IIASArev + 2.0*glc_change2*IIASArev +
3.0*glc_change3*IIASArev + 4.0*glc_change4*IIASArev + glc_div +
8.0*IIASA_crops) / 100.0))" --overwrite
 100%
[Raster MASK present]

This is not always the second time, sometimes it runs fine the first time,
but doesn't the second time, or it only runs fine the third time.

Initially I though it such errors are more likely to happen when the
equation includes double precision / floating rasters and integer rasters.
But in the example above all layers were converted to double precision
first.

I am really puzzled as to what is going on here. Any idea what I am doing
wrong here or what could be the cause?

Rgds

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