<div dir="ltr">My solution is to use a smaller model domain with 'only' 16 million cells, with grid resolution of 10 x 10 m horizontally, and 1 m vertically.<div><br></div><div>is there anyway to remove the shade in the 3D view ?</div></div><br><div class="gmail_quote"><div dir="ltr">Le mar. 15 janv. 2019 à 08:31, Francois Chartier <<a href="mailto:fra.chartier@gmail.com">fra.chartier@gmail.com</a>> a écrit :<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="auto">i will try with a smaller model domain to see if the size of the domain the source of the problem. </div><br><div class="gmail_quote"><div dir="ltr">On Mon, Jan 14, 2019, 20:19 Francois Chartier <<a href="mailto:fra.chartier@gmail.com" target="_blank">fra.chartier@gmail.com</a> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr">it finally finished the process.  the files size that was created on the external disk was 50GB as you can see on the attached screenshot.</div><div>the interpolation took 49 min to complete.</div><div dir="ltr"><br><div><br></div><div><div>(Mon Jan 14 19:28:05 2019)                                                      </div><div>v.vol.rst input=Dec16CleanedJjan12@Toronto wcolumn=dbl_2 dmin=2 elevation=dec16jan14</div><div>280957 records selected from table</div><div>WARNING: Some points outside of region -- will ignore...</div><div>WARNING: Strip exists with insufficient data</div><div>Processing all selected output files will require 836000000 bytes of disk space for temp files</div><div>WARNING: There are points outside specified 2D/3D region--ignored 55094 points (total points: 280957)</div><div>WARNING: Points are more dense than specified 'DMIN'--ignored 156188 points (remain 124769)</div><div>WARNING: Unable to create 'cross_output' raster map without 'cross_input' raster map being specified</div><div>(Mon Jan 14 20:17:49 2019) Command finished (49 min 43 sec)     </div></div></div></div><div id="gmail-m_7934602069316905507m_-6629092632311730313DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2"><br>
<table style="border-top:1px solid rgb(211,212,222)">
        <tbody><tr>
        <td style="width:55px;padding-top:13px"><a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail" rel="noreferrer" target="_blank"><img src="https://ipmcdn.avast.com/images/icons/icon-envelope-tick-round-orange-animated-no-repeat-v1.gif" alt="" width="46" height="29" style="width: 46px; height: 29px;"></a></td>
                <td style="width:470px;padding-top:12px;color:rgb(65,66,78);font-size:13px;font-family:Arial,Helvetica,sans-serif;line-height:18px">Virus-free. <a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail" style="color:rgb(68,83,234)" rel="noreferrer" target="_blank">www.avast.com</a>
                </td>
        </tr>
</tbody></table><a href="#m_7934602069316905507_m_-6629092632311730313_DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2" width="1" height="1" rel="noreferrer"></a></div><br><div class="gmail_quote"><div dir="ltr">Le lun. 14 janv. 2019 à 20:13, Francois Chartier <<a href="mailto:fra.chartier@gmail.com" rel="noreferrer" target="_blank">fra.chartier@gmail.com</a>> a écrit :<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr">grass data is on the external drive, see screenshot attached.</div><div>freespace on external drive 1.77TB</div><div>what is interesting with the external drive is the temp file does not really explode in size like if files are on laptop:</div><div>it has been running for about 50 min now.</div><div><br></div><div dir="ltr">I changed the region, made the domain slightly smaller and used a lower resolution:100x100x5<br></div><div dir="ltr"><div>g.region</div><div><div>(Mon Jan 14 19:26:20 2019)                                                      </div><div>g.region n=4875000 s=4820000 e=680000 w=585000 t=350 b=50 nsres=100 ewres=100 tbres=5</div><div>(Mon Jan 14 19:26:21 2019) Command finished (0 sec) </div></div><div><br></div><div>for some reason it seems that i have 2 nsres and ewres</div><div><div>C:\Users\Francois Chartier>g.region -p3</div><div>projection: 1 (UTM)</div><div>zone:       17</div><div>datum:      nad83</div><div>ellipsoid:  grs80</div><div>north:      4875000</div><div>south:      4820000</div><div>west:       585000</div><div>east:       680000</div><div>top:        350.00000000</div><div>bottom:     50.00000000</div><div>nsres:      100</div><div>nsres3:     5</div><div>ewres:      100</div><div>ewres3:     5</div><div>tbres:      5</div><div>rows:       550</div><div>rows3:      11000</div><div>cols:       950</div><div>cols3:      19000</div><div>depths:     60</div><div>cells:      522500</div><div>cells3:     12540000000</div></div></div></div><div id="gmail-m_7934602069316905507m_-6629092632311730313gmail-m_288032140617036686m_-995905413045375258DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2"><br>
<table style="border-top:1px solid rgb(211,212,222)">
        <tbody><tr>
        <td style="width:55px;padding-top:13px"><a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail" rel="noreferrer" target="_blank"><img src="https://ipmcdn.avast.com/images/icons/icon-envelope-tick-round-orange-animated-no-repeat-v1.gif" alt="" width="46" height="29" style="width: 46px; height: 29px;"></a></td>
                <td style="width:470px;padding-top:12px;color:rgb(65,66,78);font-size:13px;font-family:Arial,Helvetica,sans-serif;line-height:18px">Virus-free. <a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail" style="color:rgb(68,83,234)" rel="noreferrer" target="_blank">www.avast.com</a>
                </td>
        </tr>
</tbody></table><a href="#m_7934602069316905507_m_-6629092632311730313_m_288032140617036686_m_-995905413045375258_DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2" width="1" height="1" rel="noreferrer"></a></div></div><br><div class="gmail_quote"><div dir="ltr">Le dim. 13 janv. 2019 à 08:52, Markus Neteler <<a href="mailto:neteler@osgeo.org" rel="noreferrer" target="_blank">neteler@osgeo.org</a>> a écrit :<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi,<br>
<br>
On Sat, Jan 12, 2019 at 11:56 PM Francois Chartier<br>
<<a href="mailto:fra.chartier@gmail.com" rel="noreferrer" target="_blank">fra.chartier@gmail.com</a>> wrote:<br>
><br>
> Hi,<br>
><br>
> I checked the temp file that is saved during the v.vol.rst interpolation process, and the file went up to over 50 GB and more until the space on my laptop (free space of 100GB) was exausted.<br>
<br>
I start to wonder if those tmp files created are in any way compressed<br>
(probably not?). This may be worth a bug ticket at<br>
<a href="https://trac.osgeo.org/grass/newticket" rel="noreferrer noreferrer" target="_blank">https://trac.osgeo.org/grass/newticket</a><br>
<br>
> So I copied all the grass data file onto an external hard drive with 1 Terabyte of space. started grass and browsed to the data that is<br>
> saved on the external hard drive.  when running vvolrst, the process also crashed due to not enough space....<br>
> the memory situation is much better during the process.<br>
> for some reason it does not seem to increase the size of the temp file on the hard drive...<br>
<br>
So the grassdata directory with the location and mapset in question is<br>
on the terabyte disk?<br>
<br>
> see error message below.<br>
> (Sat Jan 12 17:28:02 2019)<br>
> v.vol.rst --overwrite input=Dec16CleanedJjan12@Toronto wcolumn=dbl_2 dmin=1 elevation=dec16jan12<br>
> 280957 records selected from table<br>
> WARNING: Some points outside of region -- will ignore...<br>
> Processing all selected output files will require 2134724556 bytes of disk space for temp files<br>
> WARNING: There are points outside specified 2D/3D region--ignored 1 points (total points: 280957)<br>
> WARNING: Points are more dense than specified 'DMIN'--ignored 158199 points (remain 122758)<br>
> WARNING: Unable to create 'cross_output' raster map without 'cross_input' raster map being specified<br>
> ERROR: Not enough disk space - cannot write temp files<br>
<br>
I think we should try to improve the error message, telling on which<br>
device this occurs.<br>
<br>
The code in question is:<br>
vector/v.vol.rst/main.c, line 633<br>
<br>
            /* filling temp file with zeroes */<br>
            for (i = 0; i < n_rows; i++) {<br>
                if (!<br>
                    (fwrite<br>
                     (zero_array_cell, sizeof(FCELL), n_cols, Tmp_fd_cell))) {<br>
                    clean();<br>
                    G_fatal_error(_("Not enough disk space - cannot<br>
write temp files"));<br>
                }<br>
            }<br>
<br>
Does anyone know how to get a better error message?<br>
Would this help?<br>
<a href="https://linux.die.net/man/3/explain_fwrite" rel="noreferrer noreferrer" target="_blank">https://linux.die.net/man/3/explain_fwrite</a><br>
<br>
> (Sat Jan 12 17:54:14 2019) Command finished (26 min 12 sec)<br>
><br>
> below is the g.region -p3 printed, as you can see this is a large domain and high resolution.<br>
><br>
> projection: 1 (UTM)<br>
> zone:       17<br>
> datum:      nad83<br>
> ellipsoid:  grs80<br>
> north:      4924870<br>
> south:      4784003<br>
> west:       585000<br>
> east:       679717<br>
> top:        453.00000000<br>
> bottom:     2.16080000<br>
> nsres:      50.00603479<br>
> nsres3:     5.00007099<br>
> ewres:      50.00897571<br>
> ewres3:     5.00010558<br>
> tbres:      0.999643458980044<br>
<br>
Are the pixels/voxels non-square on purpose? This could be adjusted<br>
with g.region.<br>
<br>
> rows:       2817<br>
> rows3:      28173<br>
> cols:       1894<br>
> cols3:      18943<br>
> depths:     451<br>
> cells:      5335398<br>
> cells3:     240690193689<br>
<br>
That is still a lot (just guessing around):<br>
240,690,193,689 raster3d_cells * 64 bit (but it will be larger than<br>
32bit) = 1.540417e+13 bit<br>
<br>
So, in the first place we need to know where GRASS GIS really writes<br>
to in your system.<br>
<br>
Any devs having a suggestion here?<br>
<br>
Markus<br>
</blockquote></div>
<div id="gmail-m_7934602069316905507m_-6629092632311730313gmail-m_288032140617036686DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2"><br>
<table style="border-top:1px solid rgb(211,212,222)">
        <tbody><tr>
        <td style="width:55px;padding-top:13px"><a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail" rel="noreferrer" target="_blank"><img src="https://ipmcdn.avast.com/images/icons/icon-envelope-tick-round-orange-animated-no-repeat-v1.gif" alt="" width="46" height="29" style="width: 46px; height: 29px;"></a></td>
                <td style="width:470px;padding-top:12px;color:rgb(65,66,78);font-size:13px;font-family:Arial,Helvetica,sans-serif;line-height:18px">Virus-free. <a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail" style="color:rgb(68,83,234)" rel="noreferrer" target="_blank">www.avast.com</a>
                </td>
        </tr>
</tbody></table><a href="#m_7934602069316905507_m_-6629092632311730313_m_288032140617036686_DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2" width="1" height="1" rel="noreferrer"></a></div>
</blockquote></div>
</blockquote></div>
</blockquote></div>