This routine add body force to velocity in auxField for weakly-compressible model.
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
class(mus_source_op_type), | intent(inout) | :: | fun | Description of method to update source |
||
real(kind=rk), | intent(inout) | :: | auxField(:) | output auxField array |
||
integer, | intent(in) | :: | iLevel | current level |
||
type(tem_time_type), | intent(in) | :: | time | current timing information |
||
type(tem_varSys_type), | intent(in) | :: | varSys | variable system definition |
||
type(mus_convertFac_type), | intent(in) | :: | phyConvFac | Physics conversion factor for current level |
||
type(mus_derVarPos_type), | intent(in) | :: | derVarPos(:) | position of derived quantities in varsys |
Type | Visibility | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|---|
integer, | private | :: | dens_pos | ||||
integer, | private | :: | vel_pos(3) | ||||
real(kind=rk), | private | :: | forceTerm(3) | ||||
integer, | private | :: | iElem | ||||
integer, | private | :: | nElems | ||||
integer, | private | :: | posInTotal | ||||
integer, | private | :: | elemOff | ||||
real(kind=rk), | private | :: | forceField(fun%elemLvl(iLevel)%nElems*3) | ||||
real(kind=rk), | private | :: | inv_rho |