<div dir="ltr"><div>Hi all,</div><div><br></div><div>I have written an involved script (which I run automatically), in part, calling r.basin to get the stream channel slope (preferably at the basin outlet). For all basin outlet locations, r.basin defines the watershed correctly, but for about 10% of the locations the *_parameters.csv file is not created and I get errors:</div><div><br></div><div>===========================<br>Hypsometric | quantiles<br>===========================<br>360 | 0.025<br>358 | 0.05<br>356 | 0.1<br>351 | 0.25<br>341 | 0.5<br>330 | 0.75<br>321 | 0.9<br>316 | 0.975<br><br><br>Done!<br>../src/gtk/dcclient.cpp(250): assert "Assert failure" failed in wxFreePoolGC(): Wrong GC<br>../src/gtk/dcclient.cpp(250): assert "Assert failure" failed in wxFreePoolGC(): Wrong GC<br>../src/gtk/dcclient.cpp(250): assert "Assert failure" failed in wxFreePoolGC(): Wrong GC<br>../src/gtk/dcclient.cpp(250): assert "Assert failure" failed in wxFreePoolGC(): Wrong GC<br>------------------------------<br>------------------------------<br> 100%<br>Tot. cells 1886.0<br>Tot. area 188600.0<br>Max distance 771.033539<br>Traceback (most recent call last):<br>  File "/home/teaiii/.grass7/addons/scripts/r.width.funct", line 132, in <module><br>    sys.exit(main())<br>  File "/home/teaiii/.grass7/addons/scripts/r.width.funct", line 87, in main<br>    prc[4,0] , prc[4,1] = findint(kl,0.5) , 0.5<br>  File "/home/teaiii/.grass7/addons/scripts/r.width.funct", line 125, in findint<br>    z1, z2, f1, f2 = kl[int(Xf[0])][0], kl[int(Xf[0]-1)][0], kl[int(Xf[0])][1], kl[int(Xf[0]-1)][1]<br>TypeError: only length-1 arrays can be converted to Python scalars<br><br>------------------------------<br><br>An ERROR occurred running r.basin<br>Please check for error messages above or try with another pairs of outlet<br>coordinates<br>Removing region definition <original><br>grep: mjk473_dem_parameters.csv: No such file or directory<br><br>Channel average slope (percent): <br>grep: mjk473_dem_parameters.csv: No such file or directory<br><br>Drainage Area (sq km): <br><br></div><div><br></div><div>I know the basin boundary is generated correctly because I generate a map automatically.</div><div><br></div><div>So, the above is a problem; I have sent a lot of time trying to understand where I'm going wrong. I just re-built GRASS 7.4.0 and that had no effect in correcting the problem.</div><div><br></div><div>r.basin generates waaay more than I need. Is there a simpler way the get the average main channel slope? I've looked through the available modules and can not find another way to do this. I just need a single value...</div><div><br></div><div>Best regards,</div><div>Tom<br></div><div><br></div><div><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div><br></div></div></div></div></div></div>
</div></div>