s3d = Ga(
    'e_r e_theta e_phi', g=[
        1, r**2, r**2 * sin(th)**2], coords=X, norm=True)