I think we can split this function in three different functions:

1. drag, lift and pitch moment co-rotational computation
2. added mass force
3. tangent matrix