#H in a box
#Tests free boundary conditions Poisson solver for PAW.

#Plane waves expansion:
pawecutdg 20 #much too low
ecut 10      #much too low


#kpts
kptopt 0
nkpt 1
kpt 0.0 0.0 0.0


# Self-consistent run to get the density
tolvrs     1.00d-4
nband 1
amu 1.00
nstep  10
ixc 1
icoulomb 1    #Isolated system

# H in a box
acell   3*7.0                    #This is too small
rprim      1.00   0.00   0.00
           0.00   1.00   0.00
           0.00   0.00   1.00
natom      1
xred      5.0000000000E-01  5.0000000000E-01  5.0000000000E-01
ntypat      1
typat      1
znucl     1.00

#Do a minimal test:
optforces  0 
optstress 0
prtwf  0
prtden 0
prteig 0

 pp_dirpath "$ABI_PSPDIR"
 pseudos "1h.paw"

#%%<BEGIN TEST_INFO>
#%% [setup]
#%% executable = abinit
#%% [files]
#%% files_to_test = 
#%%   t23.out, tolnlines = 15, tolabs = 4.300e-07, tolrel = 7.830e-05
#%% [paral_info]
#%% max_nprocs = 1
#%% [extra_info]
#%% authors = T. Rangel
#%% keywords =  PAW
#%% description = H in isolated boundary conditions (plane waves, PAW).
#%%   Tests free boundary conditions Poisson solver for PAW.
#%%<END TEST_INFO>