I have FP raster for which r.stats -c prints between others 687.742414-700.250492 90 if I run r.mapcalc "pok=if(dem_orig>687,if(dem_orig<701,dem_orig,null()),null())" output pok has all cells null. I guess that r.stats does not work. Radim