Added a nfar option to the d3_q2r input. When set to zero, it produces non-recentered force constants which are very similar to the ones from standard q2r, except that:
1. The cell index is R in crystal coordinates (in q2r it is R+(111))
2. The order of the Rs is different
3. Before every block of R there is the number of vectors in the block