To characterize the BLG-methanol interfaces properly taking into account time-and space molecular fluctuations, we started by defining, at each MD time step, each graphene surface as a dense 2D grid formed by the atomic positions and by the midpoint between each pair of atoms belonging to the same hexagonal ring. With this instantaneous definition of the BLG layers, the density of methanol along the perpendicular direction is estimated considering the distance between the center of mass of each molecule and the closest graphene grid point, in this way accounting for graphene corrugation.