Running the NDSPMHD examples

NDSPMHD Documentation: home | installing | running the examples | basic tutorial | references

Running the NDSPMHD examples

See the main NDSPMHD page for details about the code.
See the install page for details about how to compile the code.

Preliminaries

After unpacking the tar file:

tar xvfz ndspmhd-examples.tar.gz

and entering the examples directory,

cd ndspmhd-examples/

you will see two subdirectories for hydrodynamical and MHD tests:

$ ls
HYDRO/	MHD/

For hydrodynamics (in the HYDRO directory), there are 8 tests:

test1-2Drandom-settling/	test3.1-2Dsodshock/		test6-2Dkhinstability/
test2-2Dpairing-instability/	test4-1Drarefaction/		test7-2Duniform-shear/
test3-1Dsodshock/		test5-1Dblastwave/		test8-3Dsedov-blast-wave/

Whilst for MHD there are 7:

test1-1.5Dbriowu-shock/		test3-1.5Dtoth-strong-mhdshock/	test5-2Dmhdrotor/		test7-2Dcurrent-loop-advection/
test2-1.75Dmhd-shock/		test4-2.5Dalfvenwave/		test6-2Dorszag-tang-vortex/

These are all standard test problems in the computational astrophysics literature and appropriate references can be found along side these tests in my papers on SPH and SPMHD algorithms. The reference list gives the most useful of these.

Running the tests

It is preferable to go through the examples in the order given. For each test there is a README-blah.txt file in the subdirectory describing how to run the code step-by-step, so things are fairly self-explanatory.

For example, entering the first test directory

cd test1-2Drandom-settling/

you will find the following files:

$ ls
Makefile		random.in		rpsph.in		splash.limits
README-settling.txt	randomav.in		splash.defaults

...so simply follow the instructions in the README-settling.txt file...

$ more README-settling.txt
This is a test showing how a random particle distribution will relax to an
ordered arrangement because of the Hamiltonian nature of SPH.

Type "make clean" and "make" to re-compile the 2DSPMHD binary.

Run the code using ./2DSPMHD random.in
Or with output to a file: ./2DSPMHD random.in >& random.out &

Plot the results using "nsplash random_0*.dat"
...

Note that you will need to have SPLASH installed to view the code output.

Next steps...

After running a few of the examples, you will have a rough idea of how the code works and what the important input options are. From there you can proceed to explore the tests further or continue towards setting up your own problem...

Daniel Price's SPH pages