Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Dust growth: mass-weighted interpolations of dust-gas quantities #377

Merged

Conversation

StephaneMichoulier
Copy link
Contributor

Dust growth: mass-weighted interpolations of dust-gas quantities are now default, option removed

Type of PR:
Modification to existing code

Description:
Option to interpolate dust gas quantities such as gas density or sound speed by mass or mass/density when using dust growth has been removed. Mass weighted interpolations are more accurate than mass/density interpolations.
It was not needed to keep the mass/density weighted interpolation in the code, and thus was removed.
Screenshot 2023-03-07 at 15 22 19
Left: Radial profiles of sound velocity interpolation for two different methods (red: mass weighted, blue: mass/density weighted). Right: Interpolation of the sound speed on each dust particle for two different methods. The profiles in the left panel are calculated from the distributions in the right panel by averaging the values ​​per bin of radius.

Testing:
wbymass option has been removed from force.f90, growth.f90 and test_dustgrowth.f90.
I ran the phantom test suite + test_dustgrowth. I tested the growingdisc setup with all the options.

Did you run the bots? no

@danieljprice danieljprice merged commit 86470eb into danieljprice:master Mar 7, 2023
@StephaneMichoulier StephaneMichoulier deleted the dgrowth_interpolation branch March 8, 2023 00:13
s-neilson pushed a commit to s-neilson/phantom that referenced this pull request Mar 18, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants