#Lattice constants
#*****************
    acell    3*1.0616530850E+01
    rprim    0.0 0.5 0.5
             0.5 0.0 0.5
             0.5 0.5 0.0

#Definition of the atoms
#***********************
    natom    2
   ntypat    2
    znucl   13  33
    typat    1   2
      ixc    1
     xred    0.0  0.0  0.0
             0.25 0.25 0.25

#k-point grid, SCF cycle and plane wave basis
#********************************************
    nstep    100
     ecut    3
  nshiftk    4
   shiftk    0.5 0.5 0.5
             0.5 0.0 0.0
             0.0 0.5 0.0
             0.0 0.0 0.5

nbdbuf 0

ndtset 4
jdtset 11
       21 22 23


# DATASET 11: GS calculations (WF in the iBZ)
# *******************************************
   prtden11   1
   kptopt11   1
   toldfe11   1.0d-12
    ngkpt11   2 2 2

# DATASETS 2? : Berry phase calculations
# **************************************
    getden2?    11
    getwfk2?    11
     nband2?    4
    tolwfr2?    1.0d-22
      iscf2?   -2
    kptopt2?    1
  berryopt2?   -1
 berrystep2?    5

     ngkpt21    4  2  2     rfdir21   1 0 0
     ngkpt22    2  4  2     rfdir22   0 1 0
     ngkpt23    2  2  4     rfdir23   0 0 1


  chksymbreak 0

 pp_dirpath "$ABI_PSPDIR"
 pseudos "PseudosTM_pwteter/13al.pspnc, PseudosTM_pwteter/33as.pspnc"

#%%<BEGIN TEST_INFO>
#%% [setup]
#%% executable = abinit
#%% [files]
#%% files_to_test = 
#%%   t20.out, tolnlines = 0, tolabs = 0.0, tolrel = 0.0
#%% [paral_info]
#%% max_nprocs = 1
#%% [extra_info]
#%% authors = Unknown
#%% keywords = NC, DFPT
#%% description = 
#%%   Test of berrystep
#%%   Al-As cristal, zinc-blende structure
#%%   Compute the polarization along strings
#%%   of k-points parallel to the primitive vectors of the reciprocal
#%%   lattice, and using multiple step berry phase calculation.
#%%   (Tests 21-30 are related to geometry optimization, molecular dynamics, etc)
#%%
#%%   Cannot be executed in parallel:
#%%   chkint_prt: ERROR -
#%%   Context : the value of the variable nproc is 2.
#%%   The value of the input variable berrystep is 5, while it must be
#%%   equal to 1
#%%<END TEST_INFO>