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

Update continental extension cookbook #1837

Merged
merged 2 commits into from
Jul 6, 2017

Conversation

naliboff
Copy link
Contributor

@naliboff naliboff commented Jul 6, 2017

A few changes to the continental extension cookbook to fix the error in the radiogenic heating term (switched to compositional heating model) and a few additional minor parameter file adjustments.

@gassmoeller
Copy link
Member

/run-tests

@gassmoeller
Copy link
Member

Seems reasonable. Interesting that the results change quite a bit although the value should be the same. Do you have an explanation for that?

@naliboff
Copy link
Contributor Author

naliboff commented Jul 6, 2017

@gassmoeller -Are you looking at the differences between the high resolution images (strain rate at 10 Myr snapshot)? The shear band patterns and system overall starts evolving quite quickly once extension gets going. With the slightly higher geotherm, the system is simply at a different "state" and due to the non-linearity things change quickly. Honestly, I'm surprised it did not change more. Running the same setup with different codes (similar methods) but nearly identical initial conditions sometimes gives noticeably different results. I guess "the results change quite a bit" is somewhat up to interpretation. Working on trying to quantifying these types of changes at the moment :)

@gassmoeller
Copy link
Member

I guess my question was: Why are we at a slightly higher geotherm, if the radioactive heating before and now are the same (1.e5e-6 W/m^3)?
Is there a noticeable drop in heating over 10 Myr due to radioactive decay in the old version? Or is 1.5e-6 in the old version the initial value that is reduced by the assumed decay over 4.5 billion years?

@@ -7945,17 +7945,17 @@ \subsubsection{Continental extension}
\begin{figure}
\centering
\includegraphics[width=\textwidth]{cookbooks/continental_extension/continental_extension_cookbook_strainrate_5e6yr.png}
\caption{\it Strain rate ($s^{-1}$) after 5e6 years of extension. The black line marks the 500 $^{o}C$ isotherm.}
\caption{\it Strain rate ($s^{-1}$) after 5e6 years of extension. The black line marks the 550 $^{o}C$ isotherm.}
Copy link
Contributor

Choose a reason for hiding this comment

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

While you're there, can you change this to $550^\circ$C? That creates the proper "degree" symbol. Same for the other figures below.

Copy link
Contributor Author

Choose a reason for hiding this comment

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

👍

@naliboff
Copy link
Contributor Author

naliboff commented Jul 6, 2017

Ah, sorry for the confusion. The radioactive heating before was not correct as the the initial concentration should have been 1e6 instead of 1

@gassmoeller
Copy link
Member

Now it makes sense 😄. Ready to merge when Wolfgang's comment is addressed.

@gassmoeller gassmoeller merged commit d59e3b4 into geodynamics:master Jul 6, 2017
@naliboff naliboff deleted the update_extension_cookbook branch July 17, 2017 20:15
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

3 participants