Global kinetic ballooning mode simulations in BOUT++