Conversion of Soren's matlab script for making the file final.log
which is suitable for input into graindex

Modules
		numpy.oldnumeric

Functions
		get_ds_string(g, ds_list) Attempt to emulate the graindex (hkl) syntax. Obviously you would not want a (10,0,0) peak! make_ds_list(cell, limit=2.0) Generates a list of d-spacings sqrt(...) sqrt(x)   Return the square root of x. write_graindex_gv(outfilename, gv, tth, eta, omega, intensity, unitcell) call with array of gvectors, tth, eta and omega vals   Using an indexing object to get the ring assignments

Data
		pi = 3.1415926535897931