subroutine rmean(x,m,n)
c
c remove offset from data array
c a.tento 4/02 v.1.0
c
       implicit none
       integer m, n
       real x
       dimension x(m,3)
c
        integer i, j
        double precision zero, aved,  bufd
        parameter ( zero = 0.0d+0 )
        dimension aved(3)
c
        aved(1) = zero
        aved(2) = zero
        aved(3) = zero
        do 1 i = 1, n
        do 2 j = 1, 3
        aved(j) = aved(j) + dble(x(i,j))
2       continue
1       continue
c
        do 3 j = 1, 3        
        aved(j) = aved(j) / dble(n)
3       continue
c
        do 4 i = 1, n
        do 5 j = 1, 3
        bufd = dble(x(i,j)) - aved(j)
        x(i,j) = sngl(bufd)
5       continue
4       continue
c
c
        return
c
        end