Two major changes: 1. New formula (6) in apply_AK.py 2. Handle profiles defined on pressure boundaries