Token GAMMA= of cast double (optional)

Linear scaling of radii to space filling volume. s(new) = a(s(old)+b), where ab=GAMMA(a-1) the average Wigner-Seitz radius.