How to plot H-Bond number vs extension (nm) for protein unfolding through SMD

GROMACS version: 2021
GROMACS modification: Yes/No
Here post your question
Dear all
I am facing problem in plotting graph for H-Bond versus Extension of protein after performing Steered MD simulation on my protein. Is there any command for the same or any script by which i can find this graph. Any help would be appreciated.