[GRASS-dev] Re: [GRASS GIS] #498: r.sun2 out of sync / broken svn history

[GRASS GIS]

#498: r.sun2 out of sync / broken svn history
---------------------+------------------------------------------------------
  Reporter:  hamish  |       Owner:  hamish          
      Type:  defect  |      Status:  assigned        
  Priority:  major   |   Milestone:  6.4.0           
 Component:  Raster  |     Version:  svn-develbranch6
Resolution:          |    Keywords:  r.sun           
  Platform:  Linux   |         Cpu:  x86-32          
---------------------+------------------------------------------------------
Comment (by hamish):

 I built and tried the old r.sun, there slopein= and aspin= are required
 and it has the same problems as rsun_aspslp.png's lower right pane.

 ----

 I ran the test case with step=0.0025 (9 sec) and step=0.01 (36 sec).
 They, together with a differences map between them are in the attached pic
 rsun36_9sec_diff.png.  Note the inner eastern wing in the 9sec image has a
 rather crisp edge- it would appear to be missing some near-sunset
 iterations??

 setup:
 {{{
 time r.sun -s elevin=gauss day=$DAY \
    beam_rad=rad_test.355.beam.9sec_step.try6 step=0.0025
 real    479m0.940s
 user    476m11.914s
 sys     1m12.181s


 time r.sun -s elevin=gauss day=$DAY \
    beam_rad=rad_test.355.beam.36sec_step.try7 step=0.01
 real    217m34.835s
 user    217m12.742s
 sys     0m23.561s


 r.colors -e col=bcyr map=...

 r.mapcalc diff36_9sec = \
    "rad_test.355.beam.36sec_step.try7 - rad_test.355.beam.9sec_step.try6"

 r.colors.stddev -z diff36_9sec
 }}}

 {{{
 plot_stuff()
 {
 r.colors $MAP -e color=bcyr
 d.rast "$MAP"
 d.vect gauss_200m_contours color=180:180:180
 echo "$1" | d.text color=black
 eval `r.univar -g $MAP`
 echo "sum: $sum" | d.text color=black at=1,5
 d.legend "$MAP" range=1300,1500
 }

 ###
 d.mon x2
 d.resize w=1024 h=768
 d.split.frame 4
 #
 d.frame uno
 MAP=rad_test.355.beam.9sec_step.try6
 plot_stuff "dt=0.0025 hr (9 sec), day 355"
 #
 d.frame dos
 MAP=rad_test.355.beam.36sec_step.try7
 plot_stuff "dt=0.01 hr (36 sec), day 355"
 #
 d.frame tres
 MAP=diff36_9sec
 d.rast $MAP
 d.legend $MAP
 echo "difference between dt=9,36 sec maps" | d.text color=black
 #
 d.frame cuatro
 MAP=rad_test.355.beam.3min_step.try6
 plot_stuff "dt=0.05 hr (3 minutes), day 355"
 ###
 d.frame full_screen
 #d.out.file out=rsun36_9sec_diff format=png
 }}}


 Hamish

-- 
Ticket URL: <https://trac.osgeo.org/grass/ticket/498#comment:11>
GRASS GIS <http://grass.osgeo.org>