missing documentation code gmx_potential
Summary
I would like to calculate the electrostatic potential profil for a membrane after realising a selection on x and y, conserving all z. I could not directly calculate the potential with gmx_potential because the selection I need gives me incompatible trajectory files (different number of atoms). So I would like to calculate it from a 3D maps of charge density obtain with Volmap - VMD.
I try to tune my script with a homogeneous membrane of DLPC, comparing the results I obtain with those obtained with gmx_potential. I followed the instruction on the manual to apply Poisson's equation, I looked at the MD units but when I realised the integrations there is a difference with the gmx_potential profil.
I try to find informations about the code or how gmx_potential works in more details but it was unsuccessful.
Could you please give me more information about how gmx_potential works?