SageMath scripts used in the paper "MPDATA: Third-order accuracy for variable flows" (J. Comput. Phys. 359 2018).
- SageMath (known to work in the version 8.0)
- Tests also require pytest installed in the SageMath own version of Python.
You can install it using$ sage --python -m easy_install pytest
.
$ make
followed by $ make test
runs the tests.
The file mea.sage contains a routine to perform the modified equation analysis for a given flux function. The flux function of MPDATA is defined in schemes.sage. See test_upwind.sage and test_mpdata.sage for an example of combining the two.