$\DeclareMathOperator\diag{diag}\DeclareMathOperator\SL{SL}$It is well-known that geodesic flow $g_t=\{\diag(e^t,e^{-t}) \}_{t>0}$ acts ergodically (actually mixing) on $\SL(2,\mathbb R)$ (Howe–Moore theorem) and thus for almost all lattice $\Lambda\in X:=\SL(2,\mathbb R)/{\SL(2,\mathbb Z)}$ and any $f\in C_c(X)$, we have
$$\lim_{T\to \infty} \frac{1}{T}\int_0^Tf(g_t.\Lambda)\,dt =\int_X f(x) \, d\mu(x),$$ where $\mu$ is the $\SL(2,\mathbb R)$ invariant probability Haar measure on $\SL(2,\mathbb R)/{\SL(2,\mathbb Z)}$.
Now what happens if we replace $(g_t)$ with $(h_t):=\{\diag(t^a,t^{-a})\}_{t>0}$ (assume $a>0$)? Do we still have for almost all $\Lambda \in X$ and compactly supported continuous function $f$ (in particular $f(h_t.\Lambda)$ is bounded on $[0,1]$ and thus the limit for the normalization by $T$ of this integral and the integral starting from $t=1$ are the same) that
$$\lim_{T\to \infty} \frac{1}{T}\int_0^Tf(h_t.\Lambda) \, dt=\int_X f(x) \, d\mu(x)? \tag{*}\label{star}$$
Or the equidistribution depends on the choice of $a>0$ or there is no equidistribution at all?
The difference between $g_t$ and $h_t$ is that $g_{s+t}=g_s g_t$ but $h_{st}=h_t h_s$ and thus not a one-parameter subgroup. Ideally, I wish someone had studied an analog of Weyl's criterion for the equidistribution of $\diag\{f(t),1/f(t)\}_{t>0} \cdot \Lambda$ on $X$.
By the way, I know one can perform a change of variable on $\frac{1}{T} \int_0^Tf(g_t.\Lambda) \, dt$ to reparametrize it to become $h_t$ action. But that is not what I meant. My question is about \eqref{star}, literally.
I hope those who claim that this is a poorly stated question could provide some solid evidence showing for example when $a=1$ the a.e. equidistribution is impossible.