Sure you can get the rhombohedral lattice parameters by peak matching in Fullprof program. But Fullprof will give you only hexagonal lattice parameters which you have to trasform to the rhombohedral cell. But the transformation is easy.
In case you know the phase and only need a refinement PowderCell is still a nice tool since theer you can also change easily into the different settings or subgroups.