The concept of the EDL began with Helmholtz in 1879, who described it as a simple capacitor formed by ions in the electrolyte aligning near a charged surface. In the early 1900s, Gouy and Chapman expanded this idea with the Gouy-Chapman (GC) model, introducing the concept of a diffuse layer where the ions remain mobile and do not adsorb directly onto the charged surface. Instead, the ion distribution within this layer is governed by a balance of electrostatic attraction and thermal motion, where the ion concentration decays exponentially over a distance known as the Debye length. In 1920, Stern combined the Helmholtz and GC models to create the Gouy-Chapman-Stern (GCS) model, which accounts for both an immobile Stern layer closest to the electrode surface and a mobile diffuse layer beyond it. However, these classical theories are generally ineffective in describing the electrode-electrolyte interface for concentrated 2 to 10 M aqueous solutions of e.g. LiCl, NaCl, and CsCl as has been experimentally demonstrated by Gaddam et al. and Smith et al. Most of the earlier models fall short in capturing the intricate interfacial behavior observed in ionic liquids (ILs) due to the interplay of van der Waals, Coulombic, dipolar, and solvophobic interactions. To this end, Kornyshev and co-workers in 2007 used the Poisson-Fermi equation to model concentrated ionic liquid interfaces employing a mean-field lattice-gas model to derive an analytical expression for the differential capacitance. Around the same time, have pores of different sizes, making it challenging to efficiently probe and differentiate the EDL structure in multiporous carbons using these theories. Modelling the EDL structure of water-in-salt electrolytes in the porous carbon using Goodwin-Kornyshev's mean-field theory including ion-pairing effects and parameters for pores of different sizes, would lead to a highly complex model with numerous elements.