<jupyter_start><jupyter_text>Preliminary Reference Earth Model (PREM) O PREM é um modelo unidimensional que representa as propriedades médias da Terra em função do raio planetário. Este modelo inclui a profundidade, densidade, velocidades sísmicas, atenuação e o grau de anisotropia nos limites das camadas que compõem o planeta.  Os dados serão baixados de uma biblioteca desenvolvida em Python chamada RockHound. Esta biblioteca possui uma série de conjunto de dados úteis para as mais diversas aplicações, todas em um só lugar. Ela faz parte de um projeto chamado [Fatiando a Terra](https://www.fatiando.org). **Referências**:* <NAME>., and <NAME>. 1981. *Preliminary reference Earth model*. Phys. Earth Plan. Int. 25:297-356.* <NAME>., and <NAME>. (2019). *Rockhound: Download geophysical models/datasets and load them in Python (Version 0.1.0)*. Zenodo. https://doi.org/10.5281/zenodo.3086002  Importando bibliotecas<jupyter_code>import rockhound as rh
import matplotlib.pyplot as plt<jupyter_output><empty_output><jupyter_text>Baixando os dados<jupyter_code>prem = rh.fetch_prem()
print (prem)<jupyter_output>radius   depth   density       Vpv       Vph      Vsv      Vsh  eta  \
0    6371.0     0.0   1.02000   1.45000   1.45000  0.00000  0.00000  1.0   
1    6370.0     1.0   1.02000   1.45000   1.45000  0.00000  0.00000  1.0   
2    6369.0     2.0   1.02000   1.45000   1.45000  0.00000  0.00000  1.0   
3    6368.0     3.0   1.02000   1.45000   1.45000  0.00000  0.00000  1.0   
4    6368.0     3.0   2.60000   5.80000   5.80000  3.20000  3.20000  1.0   
..      ...     ...       ...       ...       ...      ...      ...  ...   
194   400.0  5971.0  13.05366  11.23711  11.23711  3.65027  3.65027  1.0   
195   300.0  6071.0  13.06890  11.24809  11.24809  3.65794  3.65794  1.0   
196   200.0  6171.0  13.07979  11.25593  11.25593  3.66342  3.66342  1.0   
197   100.0  6271.0  13.08632  11.26063  11.26063  3.66670  3.66670  1.0   
198     0.0  6371.0  13.08850  11.26220  11.26220  3.66780  3.66780  1.0   

      Q_mu  Q_kappa  
0      0.0  57823.0  
1      0.0  57823.0  
2      0.0  57823.0  [...]<jupyter_text>Visualizando os dados de densidade e velocidades sísmicas<jupyter_code>fig, axes = plt.subplots(1, 2, figsize=(14, 10), sharey=True)
fig.suptitle("PREM: Preliminary Reference Earth Model")
ax = axes[0]
prem.plot("density", "depth", legend=False, ax=ax)
ax.invert_yaxis()
ax.set_xlabel("Densidade [g/cm³]")
ax.set_ylabel("Profundidade [km]")
ax.grid()
ax = axes[1]
for velocity in ["Vpv", "Vsv"]:
    prem.plot(velocity, "depth", legend=False, ax=ax, label=velocity)
ax.grid()
ax.legend()
ax.set_xlabel("Velocity [km/s]")
plt.show()<jupyter_output><empty_output>
{'fraction_non_alphanumeric': 0.11747211895910781, 'fraction_numerical': 0.18587360594795538, 'mean_word_length': 2.6762295081967213, 'pattern_counts': {'":': 0, '<': 16, '<?xml version=': 0, '>': 16, 'https://': 2, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 2, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
{'chain_length': '3', 'dir': 'jupyter-structured-clean-dedup', 'id': '327178', 'license': 'CC-BY-4.0', 'n_tokens_mistral': 1356, 'n_tokens_neox': 1089, 'n_words': 288, 'path': 'Content/codes/1. PREM.ipynb', 'repo_name': 'andrelreis/geofisica2', 'size': '2712'}