Hello again,
Many thanks for all the suggestions about shell scripts [5801].
I've run into another problem with a floating point raster map in a Lat/Long location (1 degree resolution). The map plots fine provided I don't use a region which crosses the 180E meridian; if the region does cross 180E then I get an error - there's a transcript showing this below. A similar integer raster map (which the floating point map is calculated from) does not show a problem with any of these regions. The reference in the error message ("ERROR: cell_values_d: xdr_double failed for index 14.") to an index seems to vary depending on which column lies over 180E.
Any suggestions as to what might be going wrong?
David
Mapset <dorme> in Location <global>
GRASS 5.0.0 > g.region n=90 s=-90 w=-180 e=-150; d.eraseMapset <dorme> in Location <global>
GRASS 5.0.0 > d.rast TBWrich.corr
100%Mapset <dorme> in Location <global>
GRASS 5.0.0 > g.region n=90 s=-90 w=150 e=180; d.eraseMapset <dorme> in Location <global>
GRASS 5.0.0 > d.rast TBWrich.corr
100%Mapset <dorme> in Location <global>
GRASS 5.0.0 > g.region n=90 s=-90 w=165 e=-165; d.eraseMapset <dorme> in Location <global>
GRASS 5.0.0 > d.rast TBWrich.corr
ERROR: cell_values_d: xdr_double failed for index 14.Mapset <dorme> in Location <global>
GRASS 5.0.0 > g.region n=90 s=-90 w=151 e=-179; d.eraseMapset <dorme> in Location <global>
GRASS 5.0.0 > d.rast TBWrich.corr
ERROR: cell_values_d: xdr_double failed for index 0.Mapset <dorme> in Location <global>
GRASS 5.0.0 > g.region n=90 s=-90 w=179 e=-151; d.eraseMapset <dorme> in Location <global>
GRASS 5.0.0 > d.rast TBWrich.corr
ERROR: cell_values_d: xdr_double failed for index 28.