
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
Q=Numeric.array( [ [astar*0.01*i+eps,0 ,0] for i in range(1,100)]    )

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)]
