+
+ //md.atom[2].r.x=0.5*(a-1)*LC_SI;
+ //md.atom[2].r.y=0.5*(b-1)*LC_SI;
+ //md.atom[2].r.z=0;
+ //md.atom[2].element=C;
+ //md.atom[2].mass=M_C;
+
+ //md.atom[3].r.x=0.5*(a-1)*LC_SI;
+ //md.atom[3].r.y=0;
+ //md.atom[3].r.z=0;
+ //md.atom[3].element=SI;
+ //md.atom[3].mass=M_SI;