apendice_17_exercicio_nathalia

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

% Introdução à modelagem de bacias sedimentares
% Nathalia Barbosa Coelho - 2017101708
 
 % Calculation of the instantaneous syn-rift subsidence (Apêndice 21)
Ta = 1300; % temperatura na base da litosfera (ºC)
alpha = 0,00003; % coeficiente de expansão termal (1/K)
pm0 = 3330; % densidade do manto a 0 ºC (kg/m³)
pc0 = 2800; % densidade da crosta a 0 ºC (kg/m³)
ps = 2066; % densidade dos sedimentos (kg/m³)
yc = 32; % espessura inicial da crosta (km)
yl = 120; % espessura inicial do manto (km)
B = 1:4; % fator de estiramento - beta
% 1) Cálculo da subsidência com B variando de 1 a 4.
ys = yl*(1-(1/B))*((pm0-pc0)*(yc/yl)*(1-(alpha*Ta/2)*(yc/yl))-(pm0*alpha*Ta/2))/(pm0*(1-(alpha*Ta/2)-ps;
plot (ys,B,'r');
title 'Altura do bloco continental';
xlabel 'Densidade da crosta';
ylabel 'Altura';