Hello,
We have an existing aspect model result that was run in a 3D chunk but was not output with spherical coordinates (it is not feasible to rerun the model for this purpose). We would like to use bounds in longitude, latitude and depth too crop the data that we are using in further analysis. We want to make sure we are using the same cartesian to spherical transformation used in Aspect. Can someone point us to where this information is (perhaps a file in deal.ii)? I tried searching myself but only found this for a 2D example (I think? - step75). Thank you for your help! - Magali
Hi Magali:
I think that this is the function you are looking for: aspect/source/utilities.cc at main · geodynamics/aspect · GitHub
Best
W.
Thank you!
Magali