pub fn random_frac_xyz(zlow: f64, zhigh: f64) -> [f64; 3]
return random 3D fractional coords in xyz space x,y: (0, 1] z: (zlow, zhigh)