!       call reaget4(prompt,f1,f2,f3,f4)
!       Prompt for and get 4 reals; retry if
!       ill formed; return user originals if <cr>
        subroutine reaget4(prompt,f1,f2,f3,f4)
        implicit none
        real*4  f1,f2,f3,f4,vec(4)
        character*(*) prompt
        vec(1)=f1 ; vec(2)=f2 ; vec(3)=f3 ; vec(4)=f4
        call reavec(prompt,4,vec)
        f1=vec(1) ; f2=vec(2) ; f3=vec(3) ; f4=vec(4)
        return
        end