<div dir="ltr"><div>The problem seems to be fixed in the grass7.5-svn I last checked out from SVN and compiled about 3 weeks ago.<br><br></div>Hum...<br></div><div class="gmail_extra"><br><div class="gmail_quote">On 16 March 2018 at 11:22, Pierre Roudier <span dir="ltr"><<a href="mailto:pierre.roudier@gmail.com" target="_blank">pierre.roudier@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>OK, something additional I noticed: I obviously wrapped the calculations for each day of the year in a script.<br><br></div><div>Somehow, it seems I can re-process manually a day that was hanging during the script.<br><br></div>Below is the script I used:<br><div><br>for DAY in $(seq 1 365);<br>do<br>    DOY=$(printf "%03d\n" $((DAY)))<br>    <br>    echo "Processing day ${DOY} ..."<br>    <br>    eval `g.findfile element=cell file="glob_rad_${DOY}"`<br>    <br>    if [ ! "$file" ]<br>    then<br>    <br>        r.sun elev=dem_1k aspect=aspect_1k slope=slope_1k day=$DAY glob_rad=glob_rad_${DOY}<br>    <br>        sleep 15<br>    <br>    fi<br>done<br><br><br></div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On 16 March 2018 at 11:04, Pierre Roudier <span dir="ltr"><<a href="mailto:pierre.roudier@gmail.com" target="_blank">pierre.roudier@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div>Hi,<br><br></div>Just bumped into the weirdest bug in r.sun on GRASS 7.4 as installed from the Ubuntu repository.<br><br></div>I am using r.sun in Mode 2 to compute the daily average global radiation. The GRASS location is using EPSG:3031: this is data over the entire Antarctica.<br><br></div>For some specific days, that seem to be in winter, when most of the continent doesn't recieve any solar radiation, the command would hang forever at 9%:<br><br>For example I had to kill this command:<br><br>[Raster MASK present]<br>GRASS 7.4.0 (modis_lst_2015):/media/pierre<wbr>/1TB_ext4/code > time r.sun elev=dem_1k aspect=aspect_1k slope=slope_1k day=159 glob_rad=glob_rad_159 --o --verbose<br>Number of threads <1><br>Mode 2: integrated daily irradiation for a given day of the year<br>Using Linke constant: 3.000000<br>Using albedo constant: 0.200000<br>Using slope map <slope_1k><br>Using aspect map <aspect_1k><br>^C 9%<br><br>But for the day before and the day after the command was successful in less than a minute:<br><br>GRASS 7.4.0 (modis_lst_2015):/media/pierre<wbr>/1TB_ext4/code > time r.sun elev=dem_1k aspect=aspect_1k slope=slope_1k day=160 glob_rad=glob_rad_160 --o --verbose<br>Number of threads <1><br>Mode 2: integrated daily irradiation for a given day of the year<br>Using Linke constant: 3.000000<br>Using albedo constant: 0.200000<br>Using slope map <slope_1k><br>Using aspect map <aspect_1k><br> 100%<br><br>real    0m30.019s<br>user    0m29.602s<br>sys    0m0.372s<br><br></div><div>It also happened on day 151... I'm a bit stuck, so any pointers would be appreciated.<br><br></div>Cheers,<br><br></div>Pierre<br></div>
</blockquote></div><br></div>
</div></div></blockquote></div><br></div>