Dear Gromacs users I want to calculate the angle between P-N vector in the head group of DPPC molecules with the z-axis (bilayer normal) using gangle command.I know that for the z axis I should use (-g2 z) but do not know how to define a vector for P-N. Your help would be greatly appreciated. Regards Delara