[GRASS5] Re: r.sun questions

[Jaro Hofierka]

Hamish wrote:

> 
> * I moved the output of r.sun_out.txt into the "history" of each output file.
> This was an issue when running mode 2 for each day of the year and creating 
> statistics on the results with r.series. Each day would overwrite the last.
> It can be viewed with r.info or parsed using the map's file in the 
> $MAPSET/hist/ directory. The appended command history at the bottom is from a
> library function G_command_history() I wrote for src/libes/gis/history.c that
> isn't in CVS yet.
> 

It's a good idea, thanks.

> 
> Questions:
> - Declination:  Meaningless in mode 2?

no, it is still necessary in mode 2.

> - Latitude min-max:  I'm not sure if the PROJ updates changed this. Are the
> recent speedups caused by it only transforming one position (center of the 
> map*) or does it calculate for each cell, and what gets reported above is 
> simply the last instance? And why is la_max=0 (ie unset after init)?
> ????

la_max=0 is wrong, it should be -90.


>    * if so, 'g.region -l' shows "-45.3119" near the SE corner, not the middle.
>      Maybe this a common southern hemisphere bug: (-45 + 15/60) != -45.25  ?
> - sunrise & sunset times are infact slightly different for bigger maps; ok.
> 
> 
> I'm working on:
> - currently seg faults with incidout=
> - currently lets you mix mode 1 and mode 2 only options without error

Yes, this fix would be helpful as well.
Many thanks.

Jaro