Thermodynamics and quark susceptibilities: a Monte-Carlo approach to the PNJL model