Берете координаты в небесной (ECI) экваториальной системе координат , поворачиваете вокруг оси Z на угол, равный звездному времени, и получаете координаты в земной системе, в этой самой ECEF.
Чтобы точно, надо небесные координаты брать не средние на J2000.0 , а истинные на дату. Но раз вам не надо точно, то можете брать J2000.0, ошибка будет порядка 1/3 градуса.
Опять же, для расчета звездного времени можно использовать UTC вместо UT1, разница не более 0.6 сек времени, то есть 9" дуги.