destroy_wavefunctions Subroutine

subroutine destroy_wavefunctions(p)
type(wfns_t), pointer :: p

type(wfns_t), pointer :: q

do while (associated(p))
   call ps_clean_annotation(p%annotation)
   call destroy_pswf(p%wf)
   call delete(p%grid)
   q => p%next
   deallocate(p)
   p => q
enddo

end subroutine destroy_wavefunctions