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

add viscosity profile that is more consistent with original Steinberg… #5372

Merged

Conversation

jdannberg
Copy link
Contributor

…er and Calderwood formulation

This is based on a conversation with Bernhard Steinberger, who sent me this file, noting that it is more consistent with the profiles computed in his paper.

I have also updated the cookbook to show the difference between the profiles.

For all pull requests:

For new features/models or changes of existing features:

  • I have tested my new feature locally to ensure it is correct.
  • I have created a testcase for the new feature/benchmark in the tests/ directory.
  • I have added a changelog entry in the doc/modules/changes directory that will inform other users of my change.

@jdannberg jdannberg force-pushed the steinberger_continuous_viscosity branch 2 times, most recently from 67c89b9 to 778dbad Compare August 24, 2023 02:15
@jdannberg jdannberg force-pushed the steinberger_continuous_viscosity branch from 778dbad to 5be681b Compare August 24, 2023 13:37
Copy link
Member

@gassmoeller gassmoeller left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thank you for fixing this! The PR looks good, I just marked some small possible improvements. As far as I understand we prefer the new profile as a more faithful representation of the original paper. Should we make this more clear in the cookbook description so that users will preferentially use the new profile? Maybe make the new profile the default in the material model? (this would change user models, so maybe not?) Not sure what the best path is, let me know what you think.

optimize them based on observational data. During this optimization, the different
layers that make up the profiles (usually lithosphere, upper mantle, transition
zone, and lower mantle) are shifted left or right relative to each other.
The figures that show these optizimed profiles in the paper represent them as
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
The figures that show these optizimed profiles in the paper represent them as
The figures that show these optimized profiles in the paper represent them as

layers that make up the profiles (usually lithosphere, upper mantle, transition
zone, and lower mantle) are shifted left or right relative to each other.
The figures that show these optizimed profiles in the paper represent them as
piece-wise constant within 22 layers, and the profiles shown in Figures
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
piece-wise constant within 22 layers, and the profiles shown in Figures
piece-wise constant values within 22 layers, and the profiles shown in Figures

@@ -0,0 +1,5 @@
New: There is now a new viscosity profile data file that
is more consistent with how the profile is computed in
the original Steinberger and Calderwood paper.
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
the original Steinberger and Calderwood paper.
the original Steinberger and Calderwood (2006) paper.

@jdannberg
Copy link
Contributor Author

I addressed your comments.
For the default value, I checked this when making the PR, and the current default is set Radial viscosity file name = radial-visc.txt so it's still the profile with the boundary layers at the top and bottom. That's why I thought we probably do not want to change it.

But I added a sentence that we recommend the new profile.

@gassmoeller gassmoeller merged commit 82693fa into geodynamics:main Aug 29, 2023
7 checks passed
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.

None yet

2 participants