Example_EXOPLANET_WASP43b_1.m#

 1% -------------------------------------------------------------------------
 2% EXAMPLE: EXOPLANET WASP-43b - METALLICITY 1
 3%
 4% Compute equilibrium vertical composition with a metallicity 1 of WASP-43b
 5%   
 6% URL RESULTS TEA:
 7% https://github.com/dzesmin/RRC-BlecicEtal-2015a-ApJS-TEA/tree/master/Fig6/WASP43b-solar
 8%
 9% @author: Alberto Cuadra Lara
10%          PhD Candidate - Group Fluid Mechanics
11%          Universidad Carlos III de Madrid
12%                  
13% Last update Oct 12 2022
14% -------------------------------------------------------------------------
15
16LS = {'C2H2_acetylene', 'C2H4', 'C', 'CH4', 'CO2', 'CO', 'H2', 'H2O', 'H2S', 'H', 'HCN', 'He', 'HS_M', 'N2', 'N', 'NH3', 'O', 'S'};
17%% INITIALIZE
18self = App(LS);
19%% INITIAL CONDITIONS
20metallicity = 1;
21Fuel = {'H', 'He', 'C', 'N', 'O', 'S'};
22Ni_abundances = abundances2moles(Fuel, 'abundances.txt', metallicity)';
23T = linspace(100, 4000, 300);
24p = logspace(-5, 2, 300);
25
26self.PD.S_Fuel = Fuel;
27self.PD.N_Fuel = Ni_abundances;
28self = set_prop(self, 'TR', 300, 'pR', 1);
29self = set_prop(self, 'TP', T, 'pP', p);
30%% SOLVE PROBLEM
31self = solve_problem(self, 'TP');
32%% POSTPROCESSING CONFIGURATION
33self.Misc.config.label_type = 'long';
34%% DISPLAY RESULTS (PLOTS)
35plot_molar_fractions(self, self.PS.strP, 'Xi', 'p', 'ydir', 'reverse', 'xscale', 'log');