In [11]:

    

'%3.3e' % (1.268e7 - 2e6)


    

    
Out[11]:





'1.068e+07'

In [13]:

    

def mass(number, molar_mass):
    return (number / (6.022e23)) * molar_mass


    

In [17]:

    

55.845 + 32.066


    

    
Out[17]:





87.911

In [38]:

    

atoms = [1e6, 1.02e6, 3.93e5, 4.45e5]
molar_masses = [60.0843, 40.3044, 71.8444, 87.9111]

tot = 0
for a, m in zip(atoms, molar_masses):
    tot += mass(a, m)

rock_mass = (tot * (10**6) * 1.1) / (1.66 * 10 **(-24))


    

In [31]:

    

'%3.3e' % (1.068E7 - 7.08E6)


    

    
Out[31]:





'3.600e+06'

In [32]:

    

'%3.3e' % (.5 * 1.95E6)


    

    
Out[32]:





'9.750e+05'

In [36]:

    

15.9994 + 1.0794 + 1.0794


    

    
Out[36]:





18.1582

In [37]:

    

# N2, CO, H2O
atoms = [9.750E5, 7.08E6, 3.6E6]
molar_masses = [28.01348, 28.0101, 18.1582]

for a, m in zip(atoms, molar_masses):
    tot = mass(a, m)
    print (tot * (10**6)) / (1.66 * 10 **(-24))


    

    




2.73226512826e+13
1.98380544429e+14
6.53922765122e+13

In [42]:

    

print ((6.53922E13 + 1.98380544429e+14 + 2.73226512826e+13) / rock_mass)


    

    




1.5695043936

In [43]:

    




    

    




1.854696278e+14

In [45]:

    

14.00674 + 3*(1.01794)


    

    
Out[45]:





17.060560000000002

In [46]:

    

# H2O, CH4, NH3 
atoms = [1.068E7, 7.08E6, 1.95E6]
molar_masses = [18.1582, 16.0825, 17.06]

for a, m in zip(atoms, molar_masses):
    tot = mass(a, m)
    print (tot * (10**6)) / (1.66 * 10 **(-24))


    

    




1.93997086986e+14
1.13903738501e+14
3.32785809462e+13

In [50]:

    

print ((1.93997086986e+14 + 1.13903738501e+14 + 3.32785809462e+13) / rock_mass)


    

    




1.83954327444

In [50]:

    




    

In [50]:

    




    

In [52]:

    

'%3.3e' % (1.086E7 + (.1 * 1.41E7))


    

    
Out[52]:





'1.227e+07'

In [ ]: