#   C in diamond structure; Very rough  2x2x2  special point grid; low ecut.
#   Only 3 q-points

 elph2_imagden 0.1 eV
 ngkpt 3*1
 ngqpt 3*1

 ndtset 4

# Spin polarized system
 nsppol 2
 nspden 2
 occopt 2
 occ 7*1 3*0 5*1 5*0 # 10 bands in total
 spinat 0.0 0.0 1.0 0.0 0.0 1.0

#This should not change with convergence studies
 nshiftk 1
 shiftk 0.0 0.0 0.0

# Ground state density
 tolvrs1 1.0d-16

# NSCF at k
 getden2 1
 iscf2 -2
 nband2 10 10
 nbdbuf2 2
 tolwfr2 1.0d-8 # Far too low for realistic values !

# For a q != 0 calculation, one should also do NSCF at k+q
 getden3 1
 getwfk3 2
 iscf3 -2
 nqpt3 1
 nband3 10 10
 nbdbuf3 2
 tolwfr3 1.0d-8 # Far too low for realistic values !
 qpt3 0.0 0.0 0.0 # change for the qpt !

# RESPFN calculation
 getwfk4 2
 getwfq4 3
 qpt4 0.0 0.0 0.0
 nband4 10 10
 nbdbuf4 2
 ieig2rf4 5 # Produce GKK files
 bdeigrf4 10 # 8+2 buff
 smdelta4 1
 kptopt4 3 # Need full k-point set for finite q response
 rfphon4 1 # Do phonon response
 rfatpol4 1 2 # Treat displacements of all atoms
 rfdir4 1 1 1 # Do all directions (symmetry will be used)
 tolwfr4 1.0d-8

# Cell dependant parameters
 acell 3*10
 rprim 1 0 0 0 1 0 0 0 1
 nsym 1
 natom 2
 typat 1 1
 xcart -1.1746436874 0 0
        1.1746436874 0 0
 nband 10 10
 ntypat 1
 znucl 8
 ecut 10 # Far too low !
 enunit 2
 nstep 50
 istwfk *1

# test_chain = t83.in, t84.in  to be added when python script will be enabled


 pp_dirpath "$ABI_PSPDIR"
 pseudos "PseudosTM_pwteter/8o.pspnc"

#%%<BEGIN TEST_INFO>
#%% [setup]
#%% executable = abinit
#%% [shell] 
#%% post_commands = 
#%% [files]
#%% files_to_test = 
#%%   t83.out, tolnlines = 7, tolabs = 1.1e-3, tolrel = 2.0e-4, fld_options = -medium
#%% [paral_info]
#%% max_nprocs = 2
#%% [extra_info]
#%% authors = Y. Gillet 
#%% keywords = NC, DFPT, EPH_OLD
#%% description = 
#%%  O2
#%%  Oxygen spin-temperature-dependent of the band energies
#%%<END TEST_INFO>