buffer1(i+imax*(1-1)) = buffer1(i+imax*(1-1)) + buffer2(1+(i-1)*kmax) * sphi_c(1,1+s_offset_c1)
