I will present the calculation of the generalised charm susceptibilities in 2+1 flavor QCD done on the lattice at several lattice spacings. I will show that, below the chiral crossover, these susceptibilities are well described by the hadron resonance gas (HRG) model if charmed hadrons not listed in tables of the Particle Data Group are included. However, the HRG description abruptly breaks down just above the chiral crossover. To understand this, we use a model for the charm pressure in which it is expressed as the sum of partial pressures from charmed baryons, charmed mesons, and charm quarks. I will present continuum estimates of these partial pressures and findings which show that, while the partial pressures of charmed mesons and baryons drop below their respective HRG predictions, the charm quark pressure becomes non-zero above the chiral crossover.