Acceleration calculation

Calculating the stress divergence yields the acceleration (second time derivative of u). The divergence operation is just the conjugate of $\nabla^T$. For that I am using again subroutine del() with conj=1. This calculates the acceleration components from the stress tensor.