r.x=(-0.5+0.25+0.125)*ALBE_LC_SI;
r.y=(-0.5+0.25+0.125)*ALBE_LC_SI;
r.z=(-0.5+0.25)*ALBE_LC_SI;
md->atom[4372].r.x=(-0.5+0.125+0.125)*ALBE_LC_SI;
md->atom[4372].r.y=(-0.5+0.125+0.125)*ALBE_LC_SI;
r.x=(-0.5+0.25+0.125)*ALBE_LC_SI;
r.y=(-0.5+0.25+0.125)*ALBE_LC_SI;
r.z=(-0.5+0.25)*ALBE_LC_SI;
md->atom[4372].r.x=(-0.5+0.125+0.125)*ALBE_LC_SI;
md->atom[4372].r.y=(-0.5+0.125+0.125)*ALBE_LC_SI;