The codes for generating the self consistent solutions can be found in the directoryCodes for generating self-consistent solution.
They are written in the Julia programming language.
The codes for obtaining the plots are in the directory Codes for plotting. You will need Mathematica and IJulia (a Julia kernel for Jupyter notebooks) to run them.
The data files for average curvature obtained by running these codes are stored in OpenDocument Spreadsheet (.ods) format in the directory data.