Homec4science
Diffusion COSOlver (master) Browse

cosolver/src/auxval_mod.f90
c583b6e1861amaster

Last Change
Raw File

auxval_mod.f90
View Options

MODULE auxval
!
! evaluates intermediate quantities used in the collision model operators
!
USE PREC_CONST
USE BASIC
USE MODEL
USE ARRAY
USE speed_functions
use abel
use lorentz
use sugama
use improved_sugama
use coulomb_pol
USE gkcoulomb_coeffs
use coeff
!
IMPLICIT NONE
!======================================================================!
CONTAINS
SUBROUTINE evaluate_auxval_test
!
implicit none
!
call set_kernel_function
call evaluate_auxval_erf
!
IF(IFIMPROVEDSUGAMA .and. (GKE .eq. 1 .or. GKI .eq. 1)) THEN
!
call evaluate_CurlIpjk
!
ENDIF
!
END SUBROUTINE evaluate_auxval_test
!
!----------------------------------------------------------------------
!
SUBROUTINE evaluate_auxval_erf
! precompute error functions
implicit none
call evaluate_ekab
call evaluate_Erfkab
!
!
IF(IFSUGAMA .or. IF0HSC) THEN
!
call invert_model
call evaluate_ekba
call evaluate_Erfkba
!
call invert_model
!
ENDIF
!
END SUBROUTINE evaluate_auxval_erf
!
!----------------------------------------------------------------------
!
SUBROUTINE evaluate_auxval_field
! pre-compute useful quantities/coefficients
!
IMPLICIT NONE
!
call set_kernel_function
call evaluate_auxval_erf
!
IF(IFSUGAMA .and. (GKE .eq. 1 .or. GKI .eq. 1) .and. .not. IFIMPROVEDSUGAMA ) THEN
!
CALL evaluate_overlineVparaab
CALL evaluate_overlineVperpab
CALL evaluate_overlineWab
!
ENDIF
IF(IFIMPROVEDSUGAMA .and. (GKE .eq. 1 .or. GKI .eq. 1)) THEN
call evaluate_CurlIpjk
IF(ADD_IS_TO_OS) THEN
CALL evaluate_overlineVparaab
CALL evaluate_overlineVperpab
CALL evaluate_overlineWab
ENDIF
ENDIF
IF(IFABEL .and. ADDBACK) THEN
IF(IFGK .and. .not. DKBACK) THEN
CALL evaluate_rUperpAbelGK
CALL evaluate_rUparaAbelGK
CALL evaluate_rQAbelGK
ELSE
CALL evaluate_rUparaAbelDK
CALL evaluate_rQAbelDK
ENDIF
ENDIF
IF(ISELF .eq. 8 .or. ESELF .eq. 8) THEN
IF(addback) THEN
IF(IFGK .and. .not. DKBACK) THEN
CALL evaluate_rUparaLorentzGK
CALL evaluate_rUperpLorentzGK
ELSE
CALL evaluate_rUparaLorentzDK
ENDIF
ENDIF
ENDIF
!
END SUBROUTINE evaluate_auxval_field
!
!---------------------------------------------------
!
SUBROUTINE evaluate_ekab
! computes ekab for k=0,k=Jmax(e,i)
IMPLICIT NONE
!
INTEGER :: ii
!
do ii=-5,ERFMAX
ekab_array(ii) = ekab(ii)
end do
!
END SUBROUTINE evaluate_ekab
!
!---------------------------------------------------
!
SUBROUTINE evaluate_Erfkab
! computes Erfkab for k=0,k=Jmax(e,i)
integer :: ii
!
do ii=-5,ERFMAX
Erfkab_array(ii) = Erfkab(ii)
end do
!
END SUBROUTINE evaluate_Erfkab
!
!---------------------------------------------------
!
SUBROUTINE evaluate_ekba
! computes ekba for k=0,k=Jmax(e,i)
IMPLICIT NONE
!
INTEGER :: ii
!
do ii=-5,ERFMAX
ekba_array(ii) = ekab(ii)
end do
!
END SUBROUTINE evaluate_ekba
!
!---------------------------------------------------
!
SUBROUTINE evaluate_Erfkba
! computes Erfkba for k=0,k=Jmax(e,i)
integer :: ii
!
do ii=-5,ERFMAX
Erfkba_array(ii) = Erfkab(ii)
end do
!
END SUBROUTINE evaluate_Erfkba
!
!---------------------------------------------------
!
SUBROUTINE invert_ab
! inverts the species labels a->b and b->a and related quantities/routines
!
IMPLICIT NONE
!
call invert_model
call invert_basic
call evaluate_auxval_field
!
END SUBROUTINE invert_ab
!
!---------------------------------------------------
!
END MODULE auxval
c4science · Help