Marcello Gorini: > ...let's go step by step: > ls -l myscript.py > -rwxr-xr-x 1 marcello marcello 1227 2010-05-28 08:13 myscript.py
looks fine > g.gisenv > GISDBASE=/home/marcello/grassdata > LOCATION_NAME=Ocean_floor > MAPSET=level4 > MONITOR=x0 > GRASS_GUI=tcltk just out of curiosity (and not that, I think, it makes any differene): are you running (grass and) the script via the gui's command line or in text-mode? > g.region -p > projection: 99 (Mercator) > zone: 0 > datum: wgs84 > ellipsoid: wgs84 > north: 3947900 > south: 3633400 > west: -7366700 > east: -7022600 > nsres: 3700 > ewres: 3700 > rows: 85 > cols: 93 > cells: 7905 looks ok (below also) > r.info map1 ... > | Range of data: min = 0.000000 max = 10000.000000 ... > r.info map2 ... > | Range of data: min = 1000.000000 max = 1000.000000 ... > Result of "r.info outmap": > ERROR: Raster map <outmap> not found Hmmm... I don't know what could be wrong here. I guess running directly in- grass what the python script is expected to do ( r.mapcalc "outmap = map1 + map2" ) is ok, isn't it? Could you also please copy-paste the exact (python) command you use within the grass-session to get the desired "outmap"? Nikos _______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
