forked from eyoungucla/chems
-
Notifications
You must be signed in to change notification settings - Fork 0
/
output_summary_atm_SiH4_xB_SAVE.txt
156 lines (155 loc) · 5.41 KB
/
output_summary_atm_SiH4_xB_SAVE.txt
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
1.65828e-01 # Red chi-square of fit
4.78000e+00 # Planet mass in Earth masses
3.50000e+03 # Temperature K
4.07285e+04 # Pressure in bar
1.18184e-02 # Mass fraction atmosphere
6.23029e-01 # Mass fraction silicate mantle
3.65152e-01 # Mass fraction metal
2.94695e-01 # xMgO melt
1.16618e-01 # xSiO2 melt
3.99445e-01 # xMgSiO3 melt
5.32126e-02 # xFeO melt
1.12754e-02 # xFeSiO3 melt
1.08318e-04 # xNa2O melt
5.61882e-03 # xNa2SiO3 melt
3.58143e-02 # xH2 melt
8.32130e-02 # xH2O melt
8.78208e-07 # xCO melt
4.97290e-09 # xCO2 melt
3.93450e-01 # xFe metal
8.72719e-02 # xSi metal
8.48176e-02 # xO metal
4.34462e-01 # xH metal
9.55206e-01 # xH2 gas
1.39962e-04 # xCO gas
2.74250e-07 # xCO2 gas
2.01485e-03 # xCH4 gas
2.19590e-10 # xO2 gas
1.50689e-02 # xH2O gas
1.54446e-04 # xFe gas
7.28092e-04 # xMg gas
4.95883e-03 # xSiO gas
5.23375e-04 # xNa gas
2.12035e-02 # xSiH4 gas
1.06615e+03 # Moles atm
2.74502e+03 # Moles silicate
4.08412e+03 # Moles metal
-2.17375e+00 # DIW condensed planet rel to IW
4.76008e+01 # wt percent SiO2 melt
4.20335e+01 # wt percent MgO melt
6.96089e+00 # wt percent FeO melt
5.33307e-01 # wt percent Na2O melt
1.08478e-01 # wt percent H2 melt
2.25231e+00 # wt percent H2O melt
3.69578e-05 # wt percent CO melt
3.28815e-07 # wt percent CO2 melt
8.38056e+01 # wt percent Fe metal
9.34846e+00 # wt percent Si metal
5.17575e+00 # wt percent O metal
1.67021e+00 # wt percent H metal
3.56061e+00 # Uncompressed density of metal g/cc
7.19678e+00 # Uncompressed density of pure Fe g/cc
2.82205e+01 # Experimental metal percent density deficit
3.25075e+00 # grams per mole atmosphere
5.92385e+01 # wt percent H2 atm
8.35098e+00 # wt percent H2O atm
9.16485e+00 # wt percent SiO atm
5.44375e-01 # wt percent Mg atm
2.65333e-01 # wt percent Fe atm
3.70138e-01 # wt percent Na atm
1.20598e-01 # wt percent CO atm
3.71287e-04 # wt percent CO2 atm
9.94178e-01 # wt percent CH4 atm
2.16153e-07 # wt percent O2 atm
-2.95992e+00 # DIW atmosphere rel to IW
2.24813e-03 # Model KEQ for rxn 1 : Na2SiO3 = Na2O + SiO2
2.34663e-02 # Model KEQ for rxn 2 : 1/2SiO2 + Fe = FeO + 1/2Si
8.60363e-02 # Model KEQ for rxn 3 : MgSiO3 = MgO + SiO2
1.42397e+02 # Model KEQ for rxn 4 : O + 1/2Si = 1/2SiO2
1.89737e-01 # Model KEQ for rxn 5 : 2Hmetal=H2,sil
5.50362e-01 # Model KEQ for rxn 6 : FeSiO3 = FeO + SiO2
5.72861e+04 # Model KEQ for rxn 7 : 2H2O,sil + Si=SiO2 + 2H2,sil
6.55209e-01 # Model KEQ for rxn 8 : CO,g + 1/2O2 = CO2,g
3.51563e+10 # Model KEQ for rxn 9 : CH4,g + 1/2O2 = 2H2,g + CO,g
5.27507e+00 # Model KEQ for rxn 10 : H2,g + 1/2O2 = H2O,g
3.53520e-01 # Model KEQ for rxn 11 : FeO = Fe,g + 1/2O2
3.00930e-01 # Model KEQ for rxn 12 : MgO = Mg,g + 1/2O2
5.17926e+00 # Model KEQ for rxn 13 : SiO2 = SiO,g + 1/2O2
1.25451e+04 # Model KEQ for rxn 14 : Na2O = 2Na,g + 1/2O2
9.20581e-07 # Model KEQ for rxn 15 : H2,g = H2,sil
1.40519e-06 # Model KEQ for rxn 16 : H2O,g = H2O,sil
1.54059e-07 # Model KEQ for rxn 17 : CO,g = CO,sil
4.45210e-07 # Model KEQ for rxn 18 : CO2,g = CO2,sil
8.44875e-12 # Model KEQ for rxn 19 : SiO + 2H2 = SiH4 + 1/2O2
INPUT FILE:
0.032 # MgO melt
0.001 # SiO2 melt
0.9211 #MgSiO3 melt
0.0005 #FeO melt, make small
0.0005 #FeSiO3 melt, 0.15 gives -1.55 for DIW at 0.9 xFe metal or 0.09 for -2
0.007 #Na2O melt
0.001 #Na2SiO3 melt
1.0e-9 #H2 melt
1.0e-9 #H2O melt
1.0e-9 #CO melt
1.0e-9 #CO2 melt
0.99 #Fe metal
1.0e-9 #Si metal
1.0e-9 #O metal
1.0e-9 #H metal
0.999 #H2 gas
0.001 #CO gas
1.0e-9 #CO2 gas
1.0e-9 #CH4 gas
1.0e-9 #O2 gas
1.0e-9 #H2O gas
1.0e-9 #Fe gas
1.0e-9 #Mg gas
1.0e-9 #SiO gas
1.0e-9 #Na gas
1.0e-13 #SiH4 gas
2300.0 #Moles atm, 250 is nominal
2000.0 #Moles silicate,
1800.0 #Moles metal
4.78 #Planet mass in Earth masses 0.5
3500.0 #surface T in K for 3000K at core-mantle boundary
300.0 #pressure penalty (weight)
1.5 # Weight factor for mass balance: wt_massbalance
0.5 # Weight factor for mole fraction sum constraints: wt_summing
1.0 #weight for atmosphere reactions: watm_m
1.0 #weight for solubility reactions: wsolub_m
5.0 #weight for intra-melt reactions: wmelt_m
5.0 # weight for evaporation reactions: wevap_m
1.0 #weight for f1: Na2SiO3 = Na2O + SiO2
1.0 #weight for f2: 1/2SiO2 + Fe = FeO + 1/2Si
1.0 #weight for f3: MgSiO3 = MgO + SiO2
1.0 #weight for f4: O + 1/2Si = 1/2SiO2
1.0 #weight for f5: 2Hmetal=H2,sil
1.0 #weight for f6: FeSiO3 = FeO + SiO2
2.0 #weight for f7: 2H2O,sil + Si=SiO2 + 2H2,sil
1.0 #weight for f8: CO,g + 1/2O2 = CO2,g
1.0 #weight for f9: CH4,g + 1/2O2 = 2H2,g + CO,g
5.0 #weight for f10: H2,g + 1/2O2 = H2O,g
1.0 #weight for f11: FeO = Fe,g + 1/2O2
1.0 #weight for f12: MgO = Mg,g + 1/2O2
1.0 #weight for f13: SiO2 = SiO,g + 1/2O2
1.0 #weight for f14: Na2O = 2Na,g + 1/2O2
0.5 #weight for f15: H2,g = H2,sil
1.0 #weight for f16: H2O,g = H2O,sil
2.0 #weight for f17: CO,g = CO,sil
1.0 #weight for f18: CO2,g = CO2,sil
1.0 #weight for f19: SiO + 3H2 = SiH4 + H2O
1.0 #weight for f20: Sum Si
1.0 #weight for f21: Sum Mg
5.0 #weight for f22: Sum O
1.0 #weight for f23: Sum Fe
5.0 #weight for f24: Sum H
1.0 #weight for f25: Sum Na
1.0 #weight for f26: Sum C
1.0 #weight for f27: melt mole fraction sum
1.0 #weight for f28: metal mole fraction sum
1.0 #weight for f29: atmosphere mole fraction sum
0 #integer seed, if 0 then seed is chosen randomly 36
2000 #of iterations at each temperature 20000
1.0e-06 #Offset in random positions of MCMC walkers, e.g. 1e-7 for xB or 1e-9 for xH2O
3500.0 # Core-mantle equilibration temperature