Thank you Steve
I could find their IDs by dump command.
I read the variable command deeper, I found
the variable formula (equal-style) x[atom ID] y[atom ID] z[atom ID],
so having the coordinates of atoms we can calculate their distance.
But newer problem is that my program will run in a loop.
So the corners atoms situation and IDs are changing in each loop.
I can define two groups containing only the two atoms I want their distance.
Is there any way to extract the IDs of the corner atoms which are the only
members of the defined group systematically?
Reza