This python script runs a two-population dust evolution model according to Birnstiel, Klahr, Ercolano, A&A (2012). Parameters can be set by using the arguments, otherwise default parameters are used (see output). The parameters, their meaning and units can be seen by executing
This code is published on github.com/birnstiel.
For bug reports, questions, ... contact me via my website.
If you use this code in a publication, please cite at least Birnstiel, Klahr, Ercolano, A&A (2012), and possibly Birnstiel et al. (ApJL) 2015 if you use the size distribution reconstruction. I addition to that, it would be best practice to include the hash of the version you used to make sure results are reproducible, as the code can change.
Output is written in the folder
data/ by default (can be specified with option
The following files are created:
||Grain size grid||cm|
||drift-fragmentation limit on radial grid||cm|
||drift size limit on radial grid||cm|
||fragmentation limit on radial grid||cm|
||maximum particle size as function of radius and time||cm|
||lists several constants||see file contents|
||dust surface density as function of radius and time||g cm^-2|
||final dust surface density distribution (fct. of particle size and radius)||g cm^-2|
||gas surface density as function of radius and time||g cm^-2|
||times at which the snapshots were taken||s|
||small grain velocity as function of radius and time||cm s^-1|
||large grain velocity as function of radius and time||cm s^-1|
||gas velocity as function of radius and time||cm s^-1|
- [ ] proper integration of $da/dt$ instead of using exponential approximation.