se27xx / geotop Goto Github PK

Some variables in the code line:

1004 D = egy->Dlayer->co[l] + (egy->w_exc_ice->co[m]-egy->deltaw_exice-

seems not to be correctly allocated.

This needs to be adjusted!!

Please find below the debugging backtrace output.

E
.....
/home/ecor/activity/2016/rendena100_server/rendena100_aruba/rec/SnowTemperature0000L0010
/home/ecor/activity/2016/rendena100_server/rendena100_aruba/rec/SnowTemperature0000L0011
/home/ecor/activity/2016/rendena100_server/rendena100_aruba/rec/SnowTemperature0000L0012
3/6/1989 1:00 0.02% - Times: Elapsed (h:m:s) 0:02:58 Remaining (h:m) 230:38
3/6/1989 2:00 0.02% - Times: Elapsed (h:m:s) 0:03:00 Remaining (h:m) 229:50
3/6/1989 3:00 0.02% - Times: Elapsed (h:m:s) 0:03:04 Remaining (h:m) 230:18
3/6/1989 4:00 0.02% - Times: Elapsed (h:m:s) 0:03:07 Remaining (h:m) 228:19
3/6/1989 5:00 0.02% - Times: Elapsed (h:m:s) 0:03:15 Remaining (h:m) 233:36

Program received signal SIGSEGV, Segmentation fault.
SolvePointEnergyBalance (surfacemelting=, Tgd=Tgd@entry=-9999,
Tbottom=Tbottom@entry=-9999, EBd=, Convd=-9999,
surfacebalance=1, t=2.2083333333721384, Dt=Dt@entry=3600, i=i@entry=5667,
j=j@entry=4173, r=r@entry=50, c=c@entry=91, SL=SL@entry=0x23ce9f00,
SC=SC@entry=0x6eb9d0, V=V@entry=0xefe3620, egy=0x6ddab0, land=0x6dd2f0,
sl=0x6dd1f0, cnet=0x6dda20, top=0x6dd0f0, par=0x6dd3f0, ns=ns@entry=7,
ng=ng@entry=0, zmu=zmu@entry=10, zmT=zmT@entry=10, z0s=z0s@entry=0.0001,
d0s=d0s@entry=0, rz0s=rz0s@entry=0, z0v=0, d0v=0, rz0v=rz0v@entry=1,
hveg=0, v=v@entry=0.5, Ta=Ta@entry=-9.6966499867132825,
Qa=Qa@entry=0.0026893840615131374, P=P@entry=681.8266901097046, LR=7,
eps=eps@entry=0.97999999999999998, fc=fc@entry=0, LSAI=0, decaycoeff0=10,
Wcrn=Wcrn@entry=0x23fa63c8, Wcrnmax=0, Wcsn=0x23f202f8, Wcsnmax=0,
SWin=SWin@entry=63.807149110326094, LWin=LWin@entry=225.46345342289931,
SWv=0, LW=LW@entry=0x7fffffffd760, H=H@entry=0x7fffffffd7c0,
E=E@entry=0x7fffffffd7d0, LWv=LWv@entry=0x7fffffffd770,
Hv=Hv@entry=0x7fffffffd7e0, LEv=LEv@entry=0x7fffffffd7f0,
---Type to continue, or q to quit---
Etrans=Etrans@entry=0x7fffffffd800, Ts=Ts@entry=0x7fffffffd9c0,
Qs=Qs@entry=0x7fffffffd9d0, Eadd=Eadd@entry=0,
Hg0=Hg0@entry=0x7fffffffd810, Hg1=Hg1@entry=0x7fffffffd830,
Eg0=Eg0@entry=0x7fffffffd820, Eg1=Eg1@entry=0x7fffffffd840,
Qv=Qv@entry=0x7fffffffd9b0, Qg=Qg@entry=0x7fffffffd9e0,
Lob=Lob@entry=0x7fffffffd970, rh=rh@entry=0x7fffffffd920,
rv=rv@entry=0x7fffffffd930, rb=rb@entry=0x7fffffffd950,
rc=rc@entry=0x7fffffffd940, ruc=ruc@entry=0x7fffffffd960,
u_top=u_top@entry=0x7fffffffd9a0, decay=decay@entry=0x7fffffffd980,
Locc=Locc@entry=0x7fffffffd990, LWup_ab_v=LWup_ab_v@entry=0x7fffffffd7b0,
lpb=lpb@entry=0x7fffffffd9f0, dUsl=dUsl@entry=0x7fffffffd850)
at src/geotop/energy.balance.c:1004
1004 D = egy->Dlayer->co[l] + (egy->w_exc_ice->co[m]-egy->deltaw_exice->co[m])/par->excess_ice_density;
(gdb) backtrace
#0 SolvePointEnergyBalance (surfacemelting=,

Tgd=Tgd@entry=-9999, Tbottom=Tbottom@entry=-9999, EBd=<optimized out>, 
Convd=-9999, surfacebalance=1, t=2.2083333333721384, Dt=Dt@entry=3600, 
i=i@entry=5667, j=j@entry=4173, r=r@entry=50, c=c@entry=91, 
SL=SL@entry=0x23ce9f00, SC=SC@entry=0x6eb9d0, V=V@entry=0xefe3620, 
egy=0x6ddab0, land=0x6dd2f0, sl=0x6dd1f0, cnet=0x6dda20, top=0x6dd0f0, 
par=0x6dd3f0, ns=ns@entry=7, ng=ng@entry=0, zmu=zmu@entry=10, 
zmT=zmT@entry=10, z0s=z0s@entry=0.0001, d0s=d0s@entry=0, 
rz0s=rz0s@entry=0, z0v=0, d0v=0, rz0v=rz0v@entry=1, hveg=0, v=v@entry=0.5, 
Ta=Ta@entry=-9.6966499867132825, Qa=Qa@entry=0.0026893840615131374, 
P=P@entry=681.8266901097046, LR=7, eps=eps@entry=0.97999999999999998, 
fc=fc@entry=0, LSAI=0, decaycoeff0=10, Wcrn=Wcrn@entry=0x23fa63c8, 
Wcrnmax=0, Wcsn=0x23f202f8, Wcsnmax=0, SWin=SWin@entry=63.807149110326094, 
LWin=LWin@entry=225.46345342289931, SWv=0, LW=LW@entry=0x7fffffffd760, 
H=H@entry=0x7fffffffd7c0, E=E@entry=0x7fffffffd7d0, 
LWv=LWv@entry=0x7fffffffd770, Hv=Hv@entry=0x7fffffffd7e0, 
LEv=LEv@entry=0x7fffffffd7f0, Etrans=Etrans@entry=0x7fffffffd800, 
Ts=Ts@entry=0x7fffffffd9c0, Qs=Qs@entry=0x7fffffffd9d0, Eadd=Eadd@entry=0, 
Hg0=Hg0@entry=0x7fffffffd810, Hg1=Hg1@entry=0x7fffffffd830, 
Eg0=Eg0@entry=0x7fffffffd820, Eg1=Eg1@entry=0x7fffffffd840, 

---Type to continue, or q to quit---
Qv=Qv@entry=0x7fffffffd9b0, Qg=Qg@entry=0x7fffffffd9e0,
Lob=Lob@entry=0x7fffffffd970, rh=rh@entry=0x7fffffffd920,
rv=rv@entry=0x7fffffffd930, rb=rb@entry=0x7fffffffd950,
rc=rc@entry=0x7fffffffd940, ruc=ruc@entry=0x7fffffffd960,
u_top=u_top@entry=0x7fffffffd9a0, decay=decay@entry=0x7fffffffd980,
Locc=Locc@entry=0x7fffffffd990, LWup_ab_v=LWup_ab_v@entry=0x7fffffffd7b0,
lpb=lpb@entry=0x7fffffffd9f0, dUsl=dUsl@entry=0x7fffffffd850)
at src/geotop/energy.balance.c:1004
#1 0x0000000000411055 in PointEnergyBalance (i=i@entry=5667, r=50, c=91,

Dt=Dt@entry=3600, JDb=JDb@entry=726622.20833333337, 
JDe=JDe@entry=726622.25, L=L@entry=0x23ce9f00, C=C@entry=0x6eb9d0, 
S=S@entry=0x6ddfa0, G=G@entry=0x0, V=V@entry=0xefe3620, 
snowage=snowage@entry=0x6f0b50, A=A@entry=0x6dd050, E0=<optimized out>, 
Et=<optimized out>, Dtplot=Dtplot@entry=0, W=0, f=f@entry=0x0, 
SWupabove_v=SWupabove_v@entry=0x7fffffffdb10, 
Tgskin=Tgskin@entry=0x7fffffffdb20) at src/geotop/energy.balance.c:578

#2 0x0000000000414c6b in EnergyBalance (Dt=Dt@entry=3600,

JD0=JD0@entry=726622.20833333337, JDb=JDb@entry=726622.20833333337, 
JDe=JDe@entry=726622.25, L=L@entry=0x23ce9f00, C=C@entry=0x6eb9d0, 
S=S@entry=0x6ddfa0, G=G@entry=0x0, V=V@entry=0xefe3620, 

---Type to continue, or q to quit---
snowage=snowage@entry=0x6f0b50, A=A@entry=0x6dd050,
W=W@entry=0x7fffffffdc28) at src/geotop/energy.balance.c:147
#3 0x0000000000415923 in time_loop (A=A@entry=0x6dd050)

at src/geotop/geotop.c:309

#4 0x0000000000401aa5 in main (argc=2, argv=0x7fffffffddc8)

at src/geotop/geotop.c:167

(gdb)