double precision sign, den, anum

         anum = sngl(b-discri)

         if (abs((t-ta)*anum).lt.(tupper-tlower)*abs(den)) then

            t = t + anum*(ta-t)/den