いつも記事には大変お世話になっております。
　現在、DFTやTD-DFTで構造最適化した構造に対して、SAC-CI法を用いて
一点計算を行い、吸収・蛍光スペクトルを求めるという計算を行おうとしています。使用しているソフトウェアはGaussian16です。
質問したいことはSAC-CI法でのDirect法の使用に関してです。SAC-CI法では一時ファイルに情報を一旦保存しますが、I/Oにに非常に時間がかかるため、Direct法やIncore-SAC法を用いて必要なディスク容量を減らそうと考えています。
　しかし、インプットファイルでDirect法の使用やIncoreSACをキーワードで指定しても(下記インプット1,2)Direct法が適用されているようには見えず(かなりの容量の.rwfファイルが保存され、データの読み込みに時間がかかっているようです。）、Logファイルを見ますとSemi-Direct法が使用されているという記述があります(Semi-direct transformation下記アウトプット)
　どのようにすればDirect法やIncoreSAC法をSAC-CIでの計算に使用できるのか(インプットが間違っているのか、キーワードの指定法がまちがっているのか、等）ご教示いただけますと幸いです。また、もし可能であれば、ディスク容量を減らす手法が他にありましたら合わせて教えていただければ助かります。ご回答お待ちしています。

インプット1
%mem=8GB
%nprocshared=8
%nprocLinda=8
%chk=D:\testSACCI.chk
%rwf=D:\G16W\Scratch\testSACCI.rwf
%Int=D:\G16W\Scratch\testSACCI.int
%D2E=D:\G16W\Scratch\testSACCI.d2e
# 6-31+G(d,p) scrf=(cpcm,solvent=water) nosymm maxdisk=1000GB
#sac-ci=(full,singlet=nstate=8,direct,maxitdiag=999,maxitsac=9999,
maxitlin=999,maxr2op=10000000) scf=(direct,maxcyc=5000)

Title Card Required

0 1

インプット2
%nprocshared=4
%mem=32GB
# SAC-CI(InCoreSAC, MaxR2OP=4000000)/aug-cc-pvtz Nosymm
scrf=(cpcm,solvent=methanol)
scf=maxcyc=5000

Title Card Required

0 1

アウトプット
One-electron integrals computed using PRISM.
NBasis= 408 RedAO= T EigKep= 1.26D-06 NBF= 408
NBsUse= 404 1.00D-06 EigRej= 6.37D-07 NBFU= 404
ExpMin= 4.38D-02 ExpMax= 5.48D+03 ExpMxC= 8.25D+02 IAcc=3 IRadAn= 5 AccDes= 0.00D+00
Harris functional with IExCor= 205 and IRadAn= 5 diagonalized for initial guess.
HarFok: IExCor= 205 AccDes= 0.00D+00 IRadAn= 5 IDoV= 1 UseB2=F ITyADJ=14
ICtDFT= 3500011 ScaDFX= 1.000000 1.000000 1.000000 1.000000
FoFCou: FMM=F IPFlag= 0 FMFlag= 100000 FMFlg1= 0
NFxFlg= 0 DoJE=T BraDBF=F KetDBF=T FulRan=T
wScrn= 0.000000 ICntrl= 500 IOpCl= 0 I1Cent= 200000004 NGrid= 0
NMat0= 1 NMatS0= 1 NMatT0= 0 NMatD0= 1 NMtDS0= 0 NMtDT0= 0
Symmetry not used in FoFCou.
Requested convergence on RMS density matrix=1.00D-08 within5000 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on energy=1.00D-06.
No special actions if energy rises.
Inv3: Mode=1 IEnd= 18287883.
Iteration 1 A*A^-1 deviation from unit magnitude is 5.44D-15 for 375.
Iteration 1 A*A^-1 deviation from orthogonality is 2.63D-15 for 1819 1012.
Iteration 1 A^-1*A deviation from unit magnitude is 5.00D-15 for 375.
Iteration 1 A^-1*A deviation from orthogonality is 4.18D-15 for 1640 1081.
Error on total polarization charges = 0.01715
SCF Done: E(RHF) = -707.283534207 A.U. after 17 cycles
NFock= 17 Conv=0.57D-08 -V/T= 2.0015
Range of M.O.s used for correlation: 1 404
NBasis= 408 NAE= 61 NBE= 61 NFC= 0 NFV= 0
NROrb= 404 NOA= 61 NOB= 61 NVA= 343 NVB= 343

**** Warning!!: The largest alpha MO coefficient is 0.20194045D+03

Semi-Direct transformation.
ModeAB= 4 MOrb= 404 LenV= 1048284301
LASXX= 3346478955 LTotXX= 3346478955 LenRXX= 3346478955
LTotAB= 0 MaxLAS= 5470062030 LenRXY= 0
NonZer= 6692957910 LenScr= 13438411264 LnRSAI= 10913111408
LnScr1= 21911821312 LExtra= 556121766 Total= 50165944705
MaxDsk=134217728000 SrtSym= T ITran= 5
JobTyp=0 Pass 1: I= 1 to 404.