Hi Elena,

you are reading from and writing to unallocated memory. You can see for yourself, either by using valgrind or the address sanitizer option for a newer gcc. With valgrind, you would do
make CXXFLAGS='-g -O0' test_boxrichards
valgrind ./test_boxrichards

You will get tons of warnings. They all come from your spatialparams file and improper use of your arrays
Scalar **xx;
Scalar *K_;
and so on. You are inconsistent about the fact that indexing in C arrays starts with 0 and ends with size-1. For example, you start to fill xx from index 1 instead of from index 0. You also don't initialize xx in the constructor (for example by xx(nullptr)), but immediately use it in the initialiseparam_ routine.

Please try to get rid of the valgrind warnings one after another. It would be safer and more convenient to use std::vector instead of pointers.

Kind regards
Bernd

On 02/19/2015 09:11 AM, Elena Shaporenko wrote:
Hi Bernd!
Sorry, I really forgot to attach input.file.
Need to calculate a little further up to 3500c
Best regards
Elena


2015-02-19 10:52 GMT+03:00 Bernd Flemisch <[email protected] <mailto:[email protected]>>:

    Hi Elena,

    can you please also provide your .input file? Are you using box or cc?

    I just tried box with your files and grid and an unmodified input
    file. It ran through until the prescribed end time of 3000s
    without problems.

    Kind regards
    Bernd


    On 02/17/2015 01:31 PM, Elena Shaporenko wrote:
    Hi Dumuxers!
    I have problems with Richards and 2p coupled.
    You will find an attached description of the problem and the
    necessary files.
    Tell me what can be done to solve this problem?
    Thanks a lot.
    Best regards
    Elena Shaporenko


    _______________________________________________
    Dumux mailing list
    [email protected]  <mailto:[email protected]>
    https://listserv.uni-stuttgart.de/mailman/listinfo/dumux


-- _______________________________________________________________

    Bernd Flemisch                         phone: +49 711 685 69162
    IWS, Universität Stuttgart             fax:   +49 711 685 60430
    Pfaffenwaldring 61            email:[email protected]  
<mailto:[email protected]>
    D-70569 Stuttgart            url:www.hydrosys.uni-stuttgart.de  
<http://www.hydrosys.uni-stuttgart.de>
    _______________________________________________________________


    _______________________________________________
    Dumux mailing list
    [email protected]
    <mailto:[email protected]>
    https://listserv.uni-stuttgart.de/mailman/listinfo/dumux




_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux


--
_______________________________________________________________

Bernd Flemisch                         phone: +49 711 685 69162
IWS, Universität Stuttgart             fax:   +49 711 685 60430
Pfaffenwaldring 61            email: [email protected]
D-70569 Stuttgart            url: www.hydrosys.uni-stuttgart.de
_______________________________________________________________

_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux

Reply via email to