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