[GRASS-user] r.sun annual radiation

[Salvatore Mellino]

Hi,

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?

Thanks 

Salvatore 


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>