#   Ice VIII, GGA with low ecut and k point sampling
#   only one perturbation, file from DKlug.
#   Small ecut, and very small number of steps,
#   especially for RF

  ndtset 2

#Dataset 1 : SCF
  kptopt1 1
  tolwfr1 1.0d-22
   nstep1 10

#Dataset 2 : RF at q=0 0 0   (this is always needed for IFCs)

  getwfk2  1
  kptopt2  3
    nqpt2  1
     qpt2  0.0d0 0.0d0 0.0d0
 rfatpol2  11 11
   rfdir2  1 0 0
  rfphon2  1
  toldfe2  1.0d-10
   nstep2  2

# data common to all sets
   ngkpt  1 1 1
  nshiftk 1
  shiftk  0.0 0.0 0.0
   acell 3*8.0
 ecut   19.0
 ixc    11
 diemac 5.0
 natom  12
 ntypat  2
 nband  16
 occopt 1
 znucl 1 8
 typat   1 1 1 1 1 1 1 1 2 2 2 2
 rprim  -.5 .5 .72462 .5 -.5 .72462 .5 .5 -.72462

  xred     0.61537   0.19604   0.41933
           0.27671   0.19604   0.08067
           0.38463   0.80396   0.58067
           0.72329   0.80396   0.91933
           0.80396   0.38463   0.08067
           0.80396   0.72329   0.41933
           0.19604   0.61537   0.91933
           0.19604   0.27671   0.58067
           0.35639   0.10639   0.25000
           0.64361   0.89361   0.75000
           0.89361   0.64361   0.25000
           0.10639   0.35639   0.75000


 pp_dirpath "$ABI_PSPDIR"
 pseudos "PseudosGTH_pwteter/01h.pspgth, PseudosTM_pwteter/8o.pspnc"

#%%<BEGIN TEST_INFO>
#%% [setup]
#%% executable = abinit
#%% [files]
#%% files_to_test = 
#%%   t82.out, tolnlines = 1, tolabs = 1.100e-06, tolrel = 1.600e-04, fld_options = -medium
#%% [paral_info]
#%% max_nprocs = 1 
#%% [extra_info]
#%% authors = D. Klug
#%% keywords = NC, DFPT
#%% description = 
#%%   Ice VIII (H2O centered tetragonal), with only the Gamma point, and a
#%%   low ecut, and small number of steps. Test GS and RF (only one perturbation).
#%%   This is to check that local and non-local (so, varying number
#%%   of projectors) psps can be merged.
#%%   Same as test 81, except GGA.
#%%   Note that the underlying FFT grid does not respect the
#%%   symmetry of the lattice.
#%%   Contirbuted by 
#%%<END TEST_INFO>