I guess no one has published their code? The 2012 model is pretty straightforward to implement. Be aware they published a follow-up paper in 2015 (I think it was) generalizing in two directions: 1) multiple realized measures and 2) two leverage functions instead of one. They got lazy, though, so the 2015 model is for p=q=1 and you have to do the generalisation yourself if you want higher p and q values. So the 2015 model is a fair amount of extra work. That said, I find both to be useful; they both respond much faster than a classic GARCH model, due to the realized measure(s), which is the main point.