Homec4science
Diffusion COSOlver (improved_sugama) Browse

cosolver/src/array_mod.f90
cd1286e2d410improved_sugama

Last Change
Raw File

array_mod.f90
View Options

MODULE array
! contains array of pre-computed quantities
use prec_const
use FMZM
use basic
implicit none
! ! Electron collisional matrices
! DOUBLE ARRAYS
real(dp), dimension(:,:), allocatable :: CeiTdp,CeiFdp
real(dp), dimension(:,:), allocatable :: Ceedp
! global arrays
!! real(dp), dimension(:,:),allocatable :: CeiTdp_g,CeiFdp_g,Ceedp_g
!! TYPE(FM), dimension(:,:), allocatable :: CeiT,CeiF
!! TYPE(FM), dimension(:,:), allocatable :: Cee
! rows
TYPE(FM), dimension(:), allocatable :: rCeiT,rCeiF
TYPE(FM), dimension(:), allocatable :: rCee
! Ion collisional matrices
! DOUBLE ARRAYS
real(dp), dimension(:,:), allocatable :: CieTdp,CieFdp
real(dp), dimension(:,:), allocatable :: Ciidp
! ! FM ARRAYS
! TYPE(FM), dimension(:,:), allocatable :: CieT,CieF
! TYPE(FM), dimension(:,:), allocatable :: Cii
! rows
TYPE(FM), dimension(:), allocatable :: rCieT,rCieF
TYPE(FM), dimension(:), allocatable :: rCii
! Aux collisional rows for Field part of GK Sugama
TYPE(FM), dimension(:), allocatable :: roverlineVparaab
TYPE(FM), dimension(:), allocatable :: roverlineVperpab
TYPE(FM), dimension(:), allocatable :: roverlineWab
! Aux collisional rows for GK field part of Abel/Momentum conserving Lorentz
TYPE(FM), dimension(:), allocatable :: rUperp
TYPE(FM), dimension(:), allocatable :: rUpara
TYPE(FM), dimension(:), allocatable :: rQAbel
! Store non-null coefficients
! T4 coefficients
CHARACTER(len=NNLKPJlen),dimension(:), allocatable :: NNLKPJ
CHARACTER(len=T4len),dimension(:),allocatable :: T4STR
! T5 coefficients
CHARACTER(len=NNLKMPJlen),dimension(:), allocatable :: NNLKMPJ
CHARACTER(len=T5len),dimension(:),allocatable :: T5STR
! Tljpmf
CHARACTER(len=T5len),dimension(:),allocatable :: TljpmfSTR
CHARACTER(len=NNLJPMFlen),dimension(:),allocatable :: NNLJPMF
!Pipjm
!! CHARACTER(len=T5len),dimension(:),allocatable :: PipjmSTR
!! CHARACTER(len=NNPJMlen),dimension(:),allocatable :: NNPJM
CHARACTER(len=T5len),dimension(:,:,:),allocatable :: PipjmSTR
! Store speed functions
! Case C_{ab}
TYPE(FM),dimension(:),allocatable :: ekab_array
TYPE(FM), dimension(:), allocatable :: Erfkab_array
! Case C_{ba} (needed for Sugama)
TYPE(FM) ,dimension(:),allocatable :: ekba_array
TYPE(FM), dimension(:), allocatable :: Erfkba_array
! Case Improved Sugama GK
INTEGER, PARAMETER :: CurlI_IS_len = 100
CHARACTER(len=CurlI_IS_len), DIMENSION(:,:,:), ALLOCATABLE :: CurlIparaapjk_array_str
CHARACTER(len=CurlI_IS_len), DIMENSION(:,:,:), ALLOCATABLE :: CurlIperpapjk_array_str
CHARACTER(len=CurlI_IS_len), DIMENSION(:,:,:), ALLOCATABLE :: CurlIparabpjk_array_str
CHARACTER(len=CurlI_IS_len), DIMENSION(:,:,:), ALLOCATABLE :: CurlIperpbpjk_array_str
!
! Precomputed GK linear moment vector for GK Coulomb
CHARACTER(len=T5len), dimension(:,:,:,:,:), allocatable :: curlyNpjmn_STR
! Coefficient for spatial diffusion in GK test part (see. Abel, Sugama and Lorentz)
!
real(dp), dimension(:,:,:,:),allocatable :: curlyAshrg
!
END MODULE array
c4science · Help