I would like to receive some input on how you would approach the problem I’m currently facing. I’ll start with a short description of what we are trying to do.
The goal is to create a real world representation (3D) using S40RTS for the mantle and the winterc model for the top 410km. From this we want to use the gravity post processor and we want to do a Stokes solve.
Current Idea to Approach this Problem
Obviously, for the gravity calculation to be as good as possible we need the density to be in line with the density given in the winterc dataset. My best guess is to read in the density from winterc as ascii data in the Initial composition model. S40RTS can however, (as far as I know) only be read in using the Initial temperature model. Reading the winterc temperature dataset directly is also an option. However, the density obtained from this approach is not in line with the density given in the winterc dataset, so it will result in a different calculated gravity field.
How best to deal with both datasets to obtain realistic values for both temperature and density?
Using the Initial composition model is only possible for the first 410km using the winterc ascii data, what to do with the rest of the mantle? Use the Initial temperature model?
How to get from temperature perturbations (perturbed from PREM?) to absolute temperatures for the S40RTS data?
Hope my problem statement is clear enough, otherwise please ask for clarification. Thanks in advance!