real function zalignf(z)
      real z
      real p(6)
      data p/1606.3,-61.663,-993.49,15.555,212.31,13.193/
      zalignf=p(1)+p(2)*z+p(3)*z**2+p(4)*z**3+p(5)*z**4+p(6)*z**5
      end

      real function zalign(x)
      real x
      real par
      common /pawpar/par(2)
      real zalignf
      external zalignf
      zalign=par(1)*zalignf(x-par(2))
      end