Fork of NITSOL project by H.F.Walker: Taken from here Changes: common blocks eliminated, finite-difference precision can be specified etc.