[GRASS-user] r.sun annual radiation

Salvatore Mellino salvatore.mellino at gmail.com
Sun Sep 4 15:18:26 EDT 2011


thanks for the script, it works. I want to be sure that the result is the global radiation map expressed in W/m2/yr (or it is W/pixel/yr?). I compare it with ones made with ArcGis and the values are differents. 
Another question…the albedo effect is excluded?



Il giorno 31/ago/2011, alle ore 17:43, Markus Neteler ha scritto:

> On Wed, Aug 31, 2011 at 9:59 AM, Markus Neteler <neteler a osgeo.org> wrote:
>> On Mon, Aug 29, 2011 at 2:08 PM, Salvatore Mellino
>> <salvatore.mellino a gmail.com> wrote:
>>> Dear Markus,
>>> I try to use rsun_crop.sh, but I received same error message, where I make a mistake?
>>> Computing radiation for day 075...
>>> Mode 2: integrated daily irradiation for a given day of the year
>>> /Applications/GRASS-6.4.app/Contents/MacOS/scripts/rsun_crop.sh: line 42: 27151 Segmentation fault: 11  r.sun elevin=SOLelev aspin=SOLaspect slopein=SOLslope lat="$lat" day="$i" beam_rad=b_rad.$DAY diff_rad=d_rad.$DAY refl_rad=r_rad.$DAY
>> wow, it triggers a  Segmentation fault.
>> I can replicate that.
>> I have filed it as ticket
>> http://trac.osgeo.org/grass/ticket/1431
> Hamish has sync'ed r.sun2 in 6.4svn with existing 6.5svn fixes in r47993.
> No more segfault, ticket has been closed.
> Lessons learned (from https://trac.osgeo.org/grass/ticket/498),
> - no need to define slope and aspect maps, done automatically by r.sun
> - no need to define latitude
> Find the updated script attached.
> Markus
> <rsun_crop_grass6.sh>

More information about the grass-user mailing list