<div dir="auto"><div>Hi,<br><div class="gmail_extra"><br><div class="gmail_quote">On Sep 29, 2017 3:16 AM,  <<a href="mailto:balotti.and@gmail.com">balotti.and@gmail.com</a>> wrote:<br type="attribution"><blockquote class="quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
  
    
  
  <div text="#000000" bgcolor="#FFFFFF">
    <p>Hi Anna,<br>
      Here the test on NC dataset, I reduced a little bit the resolution
      to speedup computation:<br>
    </p>
    <ol>
      <li>g.region raster=elevation res=30 -p</li>
      <li>r.resamp.interp input=elevation@PERMANENT output=dem30<br>
      </li>
      <li>r.slope.aspect elevation=dem30@PERMANENT slope=slo aspect=asp</li>
      <li>r.sun.daily -t elevation=dem30@PERMANENT aspect=asp@PERMANENT
        slope=slo@PERMANENT start_day=91 end_day=304 glob_rad=sol_glob
        glob_rad_basename=sol nprocs=4</li>
      <li>t.rast.series input=sol@PERMANENT method=sum output=sol_sum<br>
      </li>
    </ol>
    <p>And the output querying on a random pixel is:</p>
    <blockquote>
      <p>est, nord: 636522.196262, 219573.598131<br>
        sol_glob@PERMANENT: value: 28270.75<br>
        sol_sum@PERMANENT: value: 1632588.35400391</p>
    </blockquote>
    <p>The only thing that could create problem maybe is that during the
      execution of r.sun.daily I received many warnings like this
      (translated: can't rename file NULL/CELL)<br>
    </p>
    <blockquote>
      <p><font size="-1">WARNING: Non riesco a rinominare il file null
          'C:\...\GRASS_data\nc_basic_<wbr>spm_grass7\PERMANENT\.tmp/<wbr>unknown\4860.1'
          in
'C:\...\GRASS_data\nc_basic_<wbr>spm_grass7\PERMANENT\cell_<wbr>misc\sol_glob\null':
          File exists</font><br>
        <font size="-1">WARNING: Non riesco a rinominare il file cella
          'C:\...\GRASS_data\nc_basic_<wbr>spm_grass7\PERMANENT\.tmp/<wbr>unknown\4860.0'
          in
          'C:\...\GRASS_data\nc_basic_<wbr>spm_grass7\PERMANENT\fcell\<wbr>sol_glob':
          File exists</font></p>
    </blockquote>
    <p>Since sum by r.sun.daily is less than sum by t.rast.series, maybe
      the summed file is not properly summed and represent just the sum
      of the last days of the period?<br>
      P.S. I am on Windows with GRASS 7.2.1 <br></p></div></blockquote></div></div></div><div dir="auto"><br></div><div dir="auto">That is helpful, thanks, I think I know what might be wrong looking at the code. It apparently works on Linux, but not on Windows. I don't have time today, but will fix it as soon as possible. Thank you!</div><div dir="auto"><br></div><div dir="auto">Anna</div><div dir="auto"><div class="gmail_extra"><div class="gmail_quote"><blockquote class="quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div text="#000000" bgcolor="#FFFFFF"><p>
    </p>
    <pre class="m_-625531586379638853moz-signature" cols="72">---
Andrea Balotti</pre><div class="elided-text">
    <div class="m_-625531586379638853moz-cite-prefix">Il 28/09/2017 19:00, Anna Petrášová ha
      scritto:<br>
    </div>
    <blockquote type="cite">
      <pre>Hi,


On Thu, Sep 28, 2017 at 5:30 AM,  <a class="m_-625531586379638853moz-txt-link-rfc2396E" href="mailto:balotti.and@gmail.com" target="_blank"><balotti.and@gmail.com></a> wrote:
</pre>
      <blockquote type="cite">
        <pre>Hi,
I am using  r.sun.daily to perform global radiation maps over a range of
dates but I am a little bit confused about the outputs of this module.
I created two outputs: the accumulated map of the whole time period and a
temporal dataset with all the daily maps.
I have tried to aggregate the daily maps in the temporal dataset with
t.rast.series (average and sum) supposing that one of them was the same with
the accumulated output of r.sun.daily but I was wrong.
Here the query on a pixel: the first line is the accumulated map from
r.sun.daily, the last two line are the aggregated maps from t.rast.series.
sol_glob_veg@PERMANENT: value: 24291.95
sol_veg_sum@PERMANENT: value: 1495683.87133789
sol_veg_avg@PERMANENT: value: 6989.<a href="tel:(769)%20688-6865" value="+17696886865" target="_blank">17696886865</a>

My question now is which kind of accumulation perform r.sun.daily?
TIA
</pre>
      </blockquote>
      <pre>it just sums the maps from each day, so t.rast.series with sum should
give you the same results. I tested it in NC sample datatset and
indeed it gave me identical results. Could you test it in the sample
dataset and perhaps send the commands you are using?

Thank you,

Anna

</pre>
      <blockquote type="cite">
        <pre>--
Andrea Balotti


______________________________<wbr>_________________
grass-user mailing list
<a class="m_-625531586379638853moz-txt-link-abbreviated" href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a>
<a class="m_-625531586379638853moz-txt-link-freetext" href="https://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">https://lists.osgeo.org/<wbr>mailman/listinfo/grass-user</a>
</pre>
      </blockquote>
    </blockquote>
    <br>
  </div></div>

</blockquote></div><br></div></div></div>