Thanks very much. I got some reasonable results.

Unexpectedly, I had to set nres.
Subst.doc led me to believe
it would automatically be set after the psf was read.