Cylindrical droplets analysis

Python script to analyze the density profile of cylindrical droplets. It analyzes 2D density files '.dat', generated by GROMACS tool 'gmx densmap', in which density has been averaged over one spatial direction and time, to extract relevant physical quantities. First it divides the droplet density plot in three regions: an external one, an internal one, which correspond to the droplet bulk, and an interfacial one. Then it divides in interfacial region in three new parts: the water-vapor interface, the water-substrate interface and the contact line region. Analyzing them separately it calculates the water and solute (if present) bulk densities and, using their profiles in the interfacial regions, also the droplets boundaries defined ad the equimolar Gibbs dividing surface. Using them, it calculates the solute adsorption at the two interfaces and at the three-phase contact line, and the contact angle with the substrate