SUBROUTINE ELEMENT_PROPERTIES(x, y, N,NCON, L, S, C, xbar, ybar)
IMPLICIT NONE

INTEGER                                 :: i
REAL                                    :: delX, delY
INTEGER                                 :: N

REAL (KIND=8)                   :: x(N), y(N), NCON(N,2)
REAL (KIND=8)                   :: L(N), S(N), C(N), xbar(N), ybar(N)

do i=1, N
        delX = x(NCON(i,2))-x(NCON(i,1))
        delY = y(NCON(i,2))-y(NCON(i,1))

        L(i) = (delX*delX + delY*delY)
        L(i) = DSQRT(L(i))

        c(i) = delX/ L(i)
        s(i) = delY/ L(i)

        xbar = (x(NCON(i,2)) + x(NCON(i,1))) * .5d0
        ybar = (y(NCON(i,2)) + y(NCON(i,1))) * .5d0
END DO
RETURN

END SUBROUTINE