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
add notebook for grain size plunge #5376
add notebook for grain size plunge #5376
Conversation
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Very nice, thanks for cleaning up the notebook! I have a few small comments, but afterwards this is good to go.
{ | ||
"cell_type": "markdown", | ||
"source": [ | ||
"This is the same field boundary grain size that is plotted in Figure 4 of Mulyokiva & Bercovici (1.6e-5 m). Note that we still had to convert from roughness to grain size after calculating $r_f$ from their equation (8), so their naming scheme is a bit misleading (one would expect it to be field boundary roughness, not grain size)." |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
"This is the same field boundary grain size that is plotted in Figure 4 of Mulyokiva & Bercovici (1.6e-5 m). Note that we still had to convert from roughness to grain size after calculating $r_f$ from their equation (8), so their naming scheme is a bit misleading (one would expect it to be field boundary roughness, not grain size)." | |
"This is the same field boundary grain size that is plotted in Figure 4 of Mulyukova & Bercovici (1.6e-5 m). Note that we still had to convert from roughness to grain size after calculating $r_f$ from their equation (8), so their naming scheme is a bit misleading (one would expect it to be field boundary roughness, not grain size)." |
"output_type": "stream", | ||
"name": "stdout", | ||
"text": [ | ||
"residual: 4164445.1831699545\n", |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Please remove the outputs from the jupyter notebook.
{ | ||
"cell_type": "markdown", | ||
"source": [ | ||
"Again, our viscosity (2.89e18 Pa s) is similar but slightly smaller than what is plotted in the figure (approximately 3.4e18 Pa s). This is expected, since we are in the diffusion creep regime, where a lower grain sizes reduces the viscosity." |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
"Again, our viscosity (2.89e18 Pa s) is similar but slightly smaller than what is plotted in the figure (approximately 3.4e18 Pa s). This is expected, since we are in the diffusion creep regime, where a lower grain sizes reduces the viscosity." | |
"Again, our viscosity (2.89e18 Pa s) is similar but slightly smaller than what is plotted in the figure (approximately 3.4e18 Pa s). This is expected, since we are in the diffusion creep regime, where a lower grain size reduces the viscosity." |
{ | ||
"cell_type": "markdown", | ||
"source": [ | ||
"We can see that the results given in the paper only approximately fuflfil the equation.\n", |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
"We can see that the results given in the paper only approximately fuflfil the equation.\n", | |
"We can see that the results given in the paper only approximately fulfil the equation.\n", |
Thank you for finding all my typos! I've fixed those, but I have left the output in the notebook, since that allows it to view the correct test results in the browser just by opening the file, without actually having to run the output. I thought this would be useful for easily comparing the notebook and the ASPECT results. |
/rebuild |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Looks good to go.
This is a follow-up PR to #5363 that adds documentation. Specifically, it adds a jupyter notebook that computes the analytical solution that the benchmark reproduces.
For all pull requests: