I have the following function (to calculate cell area) PI=3.1415926536 RES=0.0008333 R=6371007 r.mapcalc "test3 = (sin(y()+$RES/2) - sin(y()-$RES/2)) * ($RES * $PI/180) * $R^2"
This will give me a wrong results, apparently because R is an integer value. If I adapt the formula: r.mapcalc "test3 = (sin(y()+$RES/2) - sin(y()-$RES/2)) * ($RES * $PI/180) * float($R)^2" I get the correct result. From the help file I understand that if any of the arguments is a floating number the result should be a floating number. But that does not seem the case above so apparently I am missing something here. Can anybody explain the logic behind the above? Rgds Paulo
_______________________________________________ grass-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-dev
