[GRASS-user] Re: Trying python scripts.

[Nikos Alexandris]

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