<p dir="ltr">Hi Ken,</p>
<p dir="ltr">> > ...<br>
> Before I supply a full script with external data, I wonder if I can make synthetic data to debug/recreate this more easily. I see that I can do "r.mapcalc 'x = 42'" to create a new vector raster with 1 point. Is there some easy method to create a 2D raster within grass? Or should I make a simple ASCII CSV file and read that in? If I can recreate the behavior with a small hand-made data set I think it would help me figure out what is wrong with the data.</p>
<p dir="ltr">r.mapcalc uses the region definition to create the new maps. So, if your region is one cell (you can check with g.region -p) you'll have a one cell raster as output. Therefore, if you want a 2D raster, just set a bigger region, something like g.region cols=2 rows=2 (check the manual for more options to set the region), and then the r.mapcalc operation you need.</p>
<p dir="ltr">Hope it helps,<br>
Vero</p>