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 @pdbfile@qual.psf
read coor card name @pdbfile@qual.crd

write psf card name @pdbfile@qual-water-donor.psf
write coor card name @pdbfile@qual-water-donor.crd

I want to know if anything more needs to be done to regenrate the psf?