allocate_body_eph Subroutine

    subroutine allocate_body_eph(me, nx, kx)
        class(body_eph),intent(inout) :: me
        integer,intent(in) :: nx
        integer,intent(in) :: kx
        me%nx = nx
        allocate(me%tx(      nx+kx,6))  ! columns are for each state element (rx,ry,rz,vx,vy,vz)
        allocate(me%bcoef(   nx,6))
        allocate(me%f(       nx,6))
    end subroutine allocate_body_eph