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

GRASS GIS trac at osgeo.org
Wed Aug 12 06:30:33 EDT 2009

#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

 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

 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
 plot_stuff "dt=0.0025 hr (9 sec), day 355"
 d.frame dos
 plot_stuff "dt=0.01 hr (36 sec), day 355"
 d.frame tres
 d.rast $MAP
 d.legend $MAP
 echo "difference between dt=9,36 sec maps" | d.text color=black
 d.frame cuatro
 plot_stuff "dt=0.05 hr (3 minutes), day 355"
 d.frame full_screen
 #d.out.file out=rsun36_9sec_diff format=png


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

More information about the grass-dev mailing list