Fix a bug in the function 'compute_vz' of the c_functions module

When converting the velocity vector from spherical to cartesian
coordinates, a factor "sin(phi)" was erroneously used instead of
the correct "sin(theta)". This lead to wrong pulsationally-driven
broadening profiles computed with the function 'macroscopic_broadening_profile'.