DO j=1,jmax
  DO i=1,imax
    DO k=1,kmax
      buffer1(i+imax*(j-1)) = &
      buffer1(i+imax*(j-1)) + &
      buffer2(k+(i-1)*kmax) * &
      sphi_c(k,j+s_offset_c1)
    ENDDO
  ENDDO
ENDDO
