Update state with source variable "force" with 1st order integration of force in lattice Boltzmann equation. Simuilar to derive routine but it updates the state whereas derive is used for tracking Refer to Appendix in PhD Thesis of K. Masilamani "Coupled Simulation Framework to Simulate Electrodialysis Process for Seawater Desalination"
This subroutine's interface must match the abstract interface definition proc_apply_source in derived/mus_source_type_module.f90 in order to be callable via applySrc function pointer.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| class(mus_source_op_type), | intent(in) | :: | fun | Description of method to apply source terms |
||
| real(kind=rk), | intent(in) | :: | inState(:) | input pdf vector |
||
| real(kind=rk), | intent(inout) | :: | outState(:) | output pdf vector |
||
| integer, | intent(in) | :: | neigh(:) | connectivity Array corresponding to state vector |
||
| real(kind=rk), | intent(in) | :: | auxField(:) | auxField array |
||
| integer, | intent(in) | :: | nPdfSize | number of elements in state Array |
||
| integer, | intent(in) | :: | iLevel | current level |
||
| type(tem_varSys_type), | intent(in) | :: | varSys | variable system |
||
| type(tem_time_type), | intent(in) | :: | time | Point in time at which to evaluate the variable. |
||
| 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 | |||
|---|---|---|---|---|---|---|---|
| type(mus_varSys_data_type), | private, | pointer | :: | fPtr | |||
| real(kind=rk), | private | :: | forceField(fun%elemLvl(iLevel)%nElems*3) | ||||
| integer, | private | :: | iElem | ||||
| integer, | private | :: | nElems | ||||
| integer, | private | :: | iDir | ||||
| integer, | private | :: | posInTotal | ||||
| integer, | private | :: | iField | ||||
| integer, | private | :: | depField | ||||
| integer, | private | :: | nScalars | ||||
| integer, | private | :: | |||||
| integer, | private | :: | nInputStates | ||||
| real(kind=rk), | private | :: | mass_dens(varSys%nStateVars) | ||||
| real(kind=rk), | private | :: | massFrac(varSys%nStateVars) | ||||
| real(kind=rk), | private | :: | forceTerm | ||||
| real(kind=rk), | private | :: | force_elem(3) | ||||
| real(kind=rk), | private, | dimension(3, varSys%nStateVars ) | :: | spcForce | |||
| integer, | private | :: | dens_pos | ||||
| integer, | private | :: | elemOff |