mus_addForceToAuxField_MSL Subroutine

public subroutine mus_addForceToAuxField_MSL(fun, auxField, iLevel, time, varSys, phyConvFac, derVarPos)

This routine add body force to momentum in auxField for multispecies liquid model Refer to Appendix in PhD Thesis of K. Masilamani "Coupled Simulation Framework to Simulate Electrodialysis Process for Seawater Desalination"


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



integer, private :: dens_pos
integer, private :: mom_pos(3)
integer, private :: depField
real(kind=rk), private :: forceTerm(3)
real(kind=rk), private, dimension(varSys%nStateVars):: mass_dens
real(kind=rk), private, dimension(varSys%nStateVars):: massFrac
integer, private :: iElem
integer, private :: nElems
integer, private :: elemOff
integer, private :: nInputStates
integer, private :: iField
real(kind=rk), private :: forceField(fun%elemLvl(iLevel)%nElems*3)
type(mus_varSys_data_type), private, pointer:: fPtr