Homec4science
Diffusion COSOlver (new_coulomb) Browse

cosolver/src/dougherty_mod.f90
83b8f38f6c51new_coulomb

Last Change
Raw File

dougherty_mod.f90
View Options

module dougherty
! lineariezed dougherty collision operator (ony for like-species)
USE prec_const
USE basic
USE basic_mpi
USE model
USE FMZM
IMPLICIT NONE
PUBLIC
CONTAINS
subroutine doughertyTDK( l, k, rCabT)
! test component of DK dougherty
IMPLICIT NONE
!
INTEGER,INTENT(in) :: l,k
TYPE(FM), DIMENSION(:),INTENT(inout) :: rCabT
! local var
integer :: icol, irow
CALL FM_ENTER_USER_ROUTINE
irow = numba(l,k)
icol = irow
rCabT( icol) = rCabT( icol) - (2 *TO_FM( k ) + TO_FM(l))
CALL FM_EXIT_USER_ROUTINE
end subroutine doughertyTDK
subroutine doughertyFDK( l, k, rCabF)
! field component of DK dougherty
IMPLICIT NONE
!
INTEGER,INTENT(in) :: l,k
TYPE(FM), DIMENSION(:),INTENT(inout) :: rCabF
! loc. vars.
integer :: icol
CALL FM_ENTER_USER_ROUTINE
! momentum conserving term
IF( l .eq. 1 .and. k .eq. 0) THEN
icol = numba( 1, 0)
rCabF( icol) = rCabF(icol) + TO_FM(1)
ENDIF
! energy resotring term
IF( l .eq. 0 .and. k .eq. 1) THEN
icol = numba( 2, 0)
rCabF(icol) = rCabF(icol) - sqrt(TO_FM(2)) *2 /TO_FM(3)
icol = numba( 0, 1)
rCabF(icol) = rCabF(icol) + TO_FM(4) /TO_FM(3)
ENDIF
IF( l .eq. 2 .and. k .eq. 0) THEN
icol = numba( 2, 0)
rCabF(icol) = rCabF(icol) + 2 /TO_FM(3)
icol = numba( 0, 1)
rCabF(icol) = rCabF(icol) - sqrt( TO_FM(2))*2 /TO_FM(3)
ENDIF
CALL FM_EXIT_USER_ROUTINE
end subroutine doughertyFDK
subroutine doughertyselfDK( l, k, rCaa)
! self component of DK dougherty
IMPLICIT NONE
!
INTEGER,INTENT(in) :: l,k
TYPE(FM), DIMENSION(:),INTENT(inout) :: rCaa
CALL FM_ENTER_USER_ROUTINE
call doughertyTDK( l, k, rCaa)
call doughertyFDK( l, k, rCaa)
CALL FM_EXIT_USER_ROUTINE
end subroutine doughertyselfDK
end module dougherty
c4science · Help