Started implementing zhao2012 in python