Running the dust settling test

This test is described in Price & Laibe 2015

Set up the problem as usual with SETUP=dustsettle:

mkdir runs/dustsettle
cd runs/dustsettle
~/phantom/scripts/writemake.sh dustsettle > Makefile
make setup
make
make moddump

To run the test you first need to relax the disc atmosphere:

./phantomsetup relax

I recommend you do NOT add dust at this stage (i.e. set dtg=0). Just run the relaxation with gas only:

./phantom relax.in

PL15 recommended you relax the disc for ~15 orbits. By default each dump file is output every 0.1 orbits, so 15 orbits corresponds to dump file number 150.

Then use moddump to add dust:

./phantommoddump relax_00150 disc_00000 0.

and run the dust settling test:

./phantom disc.in

If you keep the same dtmax as in relax.in then every 10 orbits corresponds to disc_00100, disc_00200 etc. which are the times shown in Figures 7-9 of the PL15 paper.