
hcp_a=2.51
hcp_c=1.624*hcp_a


astar=4*math.pi/hcp_a/math.sqrt(3.0)
cstar=2*math.pi/hcp_c
eps=0.0011

(v1,v2,v3)=cella.Brillvectors

print " I VETTORI DELLO SPAZIO RECIPROCO SONO "
print cella.Brillvectors



Q=Numeric.array( [ (v1*(0.005*i+eps) + 0.0*v2 + 0.0*v3)  for i in range(1,200)]    )

debyerange=[-1,-1,-1]
debye=astar/4
SMcellrange=[2,2,2]
cellrange= [-1,-1,-1]
Kcellrange=[-1,-1,-2]
sigmacellrange=[-1,-1,-1]
sigmacharge=3
IS_IT_SCREENED=0
USE_DFM=1
selectors=[(1.0,1.0,0.0), (1.0,-1.0,0.0), (1.0,0.0,1.0), (-1.0,1.0,0.0)]
