Includes the mass loss driven by external photoevaporation into DustPy (Stammler and Birnstiel, 2022), following the implementation of Sellek et al.(2020) implementation, using the FRIED grid to calculate the mass loss rates (Haworth et al., 2018).
To setup the external photoevaporation module, add the following lines after initialization (the FRIED table is required to use this module).
from setup_externalPhotoevaporation import setup_externalPhotoevaporation_FRIED
setup_externalPhotoevaporation_FRIED(sim, fried_filename = "./friedgrid.dat", star_mass = 1., UV_flux = 1000.)
See the run_externalPhotoevaporation.py code for an example.
If you use this module, please cite Gárate et al. (in prep.)
This module is currently being implemented into the Dustpylib under the dustpylib.photoevaporation.external module.