Added HZ_calc.py - functional but ugly

Needs some tidying up