[GRASS-dev] Addons: r.sun2 segfault

[Markus Neteler]

Working on a large DEM, r.sun2 (addons) segfaults for me:

GRASS 6.4.svn (patUTM32):~ > g.region -p
projection: 1 (UTM)
zone:       32
datum:      wgs84
ellipsoid:  wgs84
north:      5157086
south:      5059541
west:       612488
east:       730098
nsres:      5
ewres:      5
rows:       19509
cols:       23522
cells:      458890698


GRASS 6.4.svn (patUTM32):~ > r.sun2 -s pat_dtm_5m horizon=horangle
horizonstep=15 \
     aspin=pat_dtm_5m.as slopein=pat_dtm_5m.sl day=180
insol_time=photoperiodo_d180
Mode 2: integrated daily irradiation
Segmentation fault

GRASS 6.4.svn (patUTM32):~ > gdb r.sun2
GNU gdb Red Hat Linux (6.5-37.el5_2.2rh)
...
This GDB was configured as "x86_64-redhat-linux-gnu"...Using host
libthread_db library "/lib64/libthread_db.so.1".

(gdb) r -s pat_dtm_5m horizon=horangle horizonstep=15
aspin=pat_dtm_5m.as slopein=pat_dtm_5m.sl day=180
insol_time=photoperiodo_d180
Starting program: /home/neteler/binaries/grass-6.4.svn/bin/r.sun2 -s
pat_dtm_5m horizon=horangle horizonstep=15 aspin=pat_dtm_5m.as
slopein=pat_dtm_5m.sl day=180 insol_time=photoperiodo_d180
[Thread debugging using libthread_db enabled]
Mode 2: integrated daily irradiation
[New Thread 47231340984336 (LWP 13269)]

Program received signal SIGSEGV, Segmentation fault.
[Switching to Thread 47231340984336 (LWP 13269)]
INPUT_part (offset=0, zmax=0x7fffc45840a8) at main.c:965
965                         horizonpointer[i] = (char)(rint(SCALING_FACTOR *
(gdb) where
#0  INPUT_part (offset=0, zmax=0x7fffc45840a8) at main.c:965
#1  0x0000000000405d6e in calculate (singleSlope=<value optimized
out>, singleAspect=<value optimized out>,
    singleAlbedo=<value optimized out>, singleLinke=<value optimized
out>, gridGeom=
      {xp = 2.3335382984350421e-310, yp = 2.3335374857606904e-310, xx0
= 6.9533063593499383e-310, yy0 = 2.3335374857180031e-310, xg0 =
2.0760920055667309e-317, yg0 = 4.4465908125712189e-323, stepx = 5,
stepy = 5, deltx = 117610, delty = 97545, stepxy = 5, sinlat =
2.3335382752661301e-310, coslat = 1.1670564030398969e-312}) at
main.c:1866
#2  0x00000000004071b1 in main (argc=9, argv=0x7fffc4584428) at main.c:745
(gdb)

As so often, I get lost in the gdb output :)
Any clue?

http://trac.osgeo.org/grass/browser/grass-addons/raster/r.sun_horizon/r.sun2

thanks,
Markus