Dear users,
I wanted to create a Python script that can compute an adiabatic profile similar to what ASPECT does in the Steinberger model, and then use it in ASPECT as an ASCII file.
I want to do this to have more control over my temperature profile (Python is easier than muparser), but I didn’t understand how to do that.
The manual says that the adiabatic temperature is calculated based on Perplex’s file for the Steinberger material model. The same applies to the Burnman compressible model, but it uses Burnman’s file, which is created from Perplex’s file after choosing a reference temperature. This might not be suitable for every occasion because you don’t choose an adiabatic surface temperature (except when you create the file, am I wrong)?
From what I understood, when ASPECT uses the Steinberger model, it computes the adiabatic temperature based on the surface adiabatic temperature using Perplex’s file. I tried to create an adiabatic profile using Python but failed. Unlike Burnman’s ASCII file (where density, temperature, gravity, etc. are well plotted), in the original Perplex file, there is every case of temperature vs. pressure, etc. So, I didn’t understand which values I should keep fixed to have an adiabatic profile that aligns with density, etc. I tried to keep the entropy fixed, but the results are not as good as the Burnman ASCII file:
Left: Perpplex (density is incorrect, so the result must be wrong, with a fixed entropy of 2000). Right Burnman (the density is correct)
Thank you in advance for the help