2Reference: Computer Physics Communications, 54, 295-305 (1989)