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 <nete...@osgeo.org> wrote: >> On Mon, Aug 29, 2011 at 2:08 PM, Salvatore Mellino >> <salvatore.mell...@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> _______________________________________________ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user