*Jim<br><br><div>Trivially: did you do &quot;g.region rast=na_dem_sub_20&quot; before raster computations?<br>And r.mapcalc expression is without any quotes like r.mapcalc &quot;na_dem_sub_20 = na_dem_input - 20&quot; or r.mapcalc &#39;na_dem_sub_20 = na_dem_input - 20&#39;<br>


</div><br>------------------------------------<br>Alexander<br>