-
Notifications
You must be signed in to change notification settings - Fork 68
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
Pratt isostasy #37
Comments
Hi Can I try tackle this issues? |
@michiboo of course, go right ahead 👍 |
Hi! I am attempting to solve this issue. To clarify, the output will be a grid of densities of ocean root and continent root? The input will be the assumed constant root thickness ? Does fatiando calculate an isostatic residual? |
Hi @andieie! I'm glad to have you onboard! I'll try to answer your questions.
Yes, but no necessarily a grid. The returned What @leouieda mentioned is that we should also enable the function to return an
Yes, the function will assume a constant depth for the bottom boundary of the crust, as Pratt does. This depth should be passed through the
Not specifically. I mean we don't plan to have a function that computes the isostatic residual anomaly entirely, but provide the tools to do so if anyone wants to. So, if I want to compute the gravitational field generated by the root of the crust (considering an Airy isostasy) or by the anomalous densities (Pratt isostasy), I would model it through prisms and then compute the gravitational field they generate (after #86 is merged, that will be very easy).
No problem! Feel free to ask anything! Hope my answers are helpful. If not, don't hesitate to write back. You can follow Heiskanen & Moritz (1967) if you need some background for the math and/or a nice figure to understand what the terms in the equation are. |
Description of the desired feature
PRs #17 and #36 added a function to calculate the Moho depth based on an Airy isostasy model. It would be great to have a Pratt isostasy calculation function as well. It would take the topography/bathymetry, respective densities, and a compensation depth and output the lateral density variation. The function should look something like this:
The body will be very similar to
isostasy_airy
to allow numpy arrays orxarray.DataArray
as input topographies. The docstring would benefit from having a figure to show what the terms in the equations are.A good example case would be calculate the density variation along the mid-atlantic ridge.
The text was updated successfully, but these errors were encountered: