hacal function

Calibrate Hargreaves model