I
                        made a mistake with the output from
                        phenix.reduce, now it works well. In my first
                        attempt, the output from phenix.reduce was a
                        file whose only content was: "USER  MOD
                        reduce.3.24.130724 H: found=0, std=0, add=0,
                        rem=0, adj=0". It was solved when I renamed the
                        input pdb. I took erroneously the first file as
                        an input for phenix.refine and it produced an
                        error (obviously). Now I repeated the refinement
                        following exactly your instructions. I got the
                        same results (as expected). 
                     
                    Final:
                      r_work = 0.1755 r_free = 0.2044 bonds = 0.030
                      angles = 0.598
                    
                    
                    Final: r_work = 0.1663 r_free = 0.2017 bonds = 0.011
                    angles = 1.321