Okay since I found that using the usual tensors according to xCoba is very messy, I thought of using CTensors.
So to be particular about what I am doing, I'll mention the whole problem and attach the semi-working code.
I have a AdS4 space, with global coordinates. I have a massive free scalar field coupled to the curved space. Now I can solve the equations of motion and find the solutions to the free scalar field, which I have already done.
Now I need to find the Energy-Momentum tensor for these scalar field solutions.
The scalar fields solutions turn out to be e^{-iwt}SphericalHarmonics[l,m,theta,phi](...)
where (...) is some Hypergeometric function.
Just to see if the code works fine, I am trying act act $\Nabla^2_{S^2}$ on SphericalHarmonicsY[l,m,theta,phi] which should give me the l(l+1)SphericalHarmonicsY[l,m,theta,phi] but I am not getting it. There are terms with basis with different colours, so what do I do ?
Having said that I am not even able to do a simple computation like, finding the covariant derivative on just one of the coordinate say r[]. That too is giving be some basis. What do I do ?