I am trying to get water water hydrogen bonds. I have added the donor list to TIP3 water model as described in the hbonds.inp. I get the same psf without any extra donors on regenerating the psf using this code:
read rtf card name top_all27.rtf
read para card name par_all27.prm
read psf card name @firstname.lastname@example.org
read coor card name @email@example.com
write psf card name @firstname.lastname@example.org
write coor card name @email@example.com
I want to know if anything more needs to be done to regenrate the psf?