3. Examples for apatite saturation temperature (AST) calculation

Calculate apatite saturation temperature

The apatite saturation temperature is calculated using the empirical equation of Harrison & Watson (1984) https://doi.org/10.1016/0016-7037(84)90403-4

[1]:

## use the 3 lines below if you want to use the pyAp codes inthe subdirectory next to the current path
## comment them, if you want to use the pyAp installed elsewhere on your computer
import os, sys
if not os.path.exists('pyAp') and os.path.exists('../pyAp'):
    sys.path.insert(1, os.path.abspath('..'))

# import pyAp modules
import pyAp.ApSatTemp as calc_AST

Inputs are SiO $_{2}$ and P $_{2}$ O $_{5}$ concentrations in the melt.

Note that concetrations are in weight percent (between 0 and 100) instead of mass fractions (=concentration wt/100; between 0 and 1).

When mass fractions are used for calculation:

[2]:

print('>> AST (°C) =',calc_AST.AST(0.55,0.003))

>> Warning: Concentrations are expressed as mass fractions.
 Multiply them by 100 to convert into percentage!
>> AST (°C) = nan

When wt% concentrations are used for calculation:

[3]:

print('>> AST (°C) =', calc_AST.AST(55,0.3))

>> AST (°C) = 849.3629335400316

[ ]: