Homec4science
Diffusion COSOlver (improved_sugama) Browse

cosolver/src/memory.f90
cd1286e2d410improved_sugama

Last Change
Raw File

memory.f90
View Options

subroutine memory
!
! initializes and allocate arrays
!
use array
use prec_const
use model
use basic
use grid
!
implicit none
IF(eicolls .and. is_e) THEN
!
! local arrays
CALL allocate_array(rCeiT,numbe(0,0),numbe(Pmaxe,Jmaxe))
CALL allocate_array(rCeiF,numbi(0,0),numbi(Pmaxi,Jmaxi))
CALL allocate_array(CeiTdp,lbare_s_l,lbare_e_l,numbe(0,0),numbe(Pmaxe,Jmaxe))
CALL allocate_array(CeiFdp,lbare_s_l,lbare_e_l,numbi(0,0),numbi(Pmaxi,Jmaxi))
!
!
ENDIF
IF(eecolls .and. is_self) THEN
! local array
CALL allocate_array(rCee,numbe(0,0),numbe(Pmaxe,Jmaxe))
CALL allocate_array(Ceedp,lbaree_s_l,lbaree_e_l,numbe(0,0),numbe(Pmaxe,Jmaxe))
ENDIF
IF(iecolls .and. is_i) THEN
! local array
CALL allocate_array(rCieT,numbi(0,0),numbi(Pmaxi,Jmaxi))
CALL allocate_array(rCieF,numbe(0,0),numbe(Pmaxe,Jmaxe))
CALL allocate_array(CieTdp,lbari_s_l,lbari_e_l,numbi(0,0),numbi(Pmaxi,Jmaxi))
CALL allocate_array(CieFdp,lbari_s_l,lbari_e_l,numbe(0,0),numbe(Pmaxe,Jmaxe))
!
!
ENDIF
IF(iicolls .and. is_self) THEN
CALL allocate_array(rCii,numbi(0,0),numbi(Pmaxi,Jmaxi))
CALL allocate_array(Ciidp,lbarii_s_l,lbarii_e_l,1,numbi(Pmaxi,Jmaxi))
! For GK Polarization Coulomb Operator
IF(IFPOL) THEN
CALL allocate_array(PipjmSTR,0,PMMAXX,0,JEmaxx,0,PMMAXX)
ENDIF
ENDIF
! Allocate error function arrays
CALL allocate_array(ekab_array,-5,ERFMAX)
CALL allocate_array(Erfkab_array,-5,ERFMAX)
IF(IFSUGAMA .or. IF0HSC) THEN
CALL allocate_array(ekba_array,-5,ERFMAX)
CALL allocate_array(Erfkba_array,-5,ERFMAX)
ENDIF
IF(IFIMPROVEDSUGAMA .and. (GKE .eq. 1 .or. GKI .eq. 1)) THEN
print *, 'Allocate arrays for CurlI(para/perp)(a/b)...'
CALL allocate_array(CurlIparaapjk_array_str,0,MAX(Pmaxe,Pmaxi),0,MAX(Jmaxe,Jmaxi),0,MAX(ImpSugamaJmax, ImpSugamaKmax))
CALL allocate_array(CurlIperpapjk_array_str,0,MAX(Pmaxe,Pmaxi),0,MAX(Jmaxe,Jmaxi),0,MAX(ImpSugamaJmax, ImpSugamaKmax))
CALL allocate_array(CurlIparabpjk_array_str,0,MAX(Pmaxe,Pmaxi),0,MAX(Jmaxe,Jmaxi),0,MAX(ImpSugamaJmax, ImpSugamaKmax))
CALL allocate_array(CurlIperpbpjk_array_str,0,MAX(Pmaxe,Pmaxi),0,MAX(Jmaxe,Jmaxi),0,MAX(ImpSugamaJmax, ImpSugamaKmax))
END IF
end subroutine memory
c4science · Help