Take imaginary part of shape coefficients.
s1=imag(s)
This function replaces all shape coefficients with their imaginary parts.