Update state with the second-order Guo force for TRT collision.
The first-order Hermite contribution is odd and therefore has to be weighted with the antisymmetric TRT relaxation rate. The second-order contribution is even and uses the symmetric relaxation rate.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| class(mus_source_op_type), | intent(in) | :: | fun | |||
| real(kind=rk), | intent(in) | :: | inState(:) | |||
| real(kind=rk), | intent(inout) | :: | outState(:) | |||
| integer, | intent(in) | :: | neigh(:) | |||
| real(kind=rk), | intent(in) | :: | auxField(:) | |||
| integer, | intent(in) | :: | nPdfSize | |||
| integer, | intent(in) | :: | iLevel | |||
| type(tem_varSys_type), | intent(in) | :: | varSys | |||
| type(tem_time_type), | intent(in) | :: | time | |||
| type(mus_convertFac_type), | intent(in) | :: | phyConvFac | |||
| type(mus_derVarPos_type), | intent(in) | :: | derVarPos(:) |