- printf("[sic] creating atoms\n");
- //memset(&v,0,sizeof(t_3dvec));
- //r.y=0;
- //r.z=0;
- //r.x=0.23*sqrt(3.0)*LC_SI/2.0;
- //add_atom(&md,SI,M_SI,0,ATOM_ATTR_2BP,&r,&v);
- //r.x=-r.x;
- //add_atom(&md,SI,M_SI,0,ATOM_ATTR_2BP,&r,&v);
- create_lattice(&md,DIAMOND,LC_SI,SI,M_SI,ATOM_ATTR_2BP,0,10,10,10);
-
- /* set temperature */
- printf("[sic] setting temperature\n");
- set_temperature(&md,0.0);
+ //
+#ifdef ALBE
+ create_lattice(&md,DIAMOND,ALBE_LC_SI,SI,M_SI,
+ //create_lattice(&md,DIAMOND,ALBE_LC_C,C,M_C,
+#else
+ create_lattice(&md,DIAMOND,LC_SI,SI,M_SI,
+#endif
+ ATOM_ATTR_1BP|ATOM_ATTR_2BP|ATOM_ATTR_3BP|ATOM_ATTR_HB,
+ // ATOM_ATTR_2BP|ATOM_ATTR_HB,
+ 0,LCNTX,LCNTY,LCNTZ,NULL);
+ // 1,LCNTX,LCNTY,LCNTZ,NULL);
+
+ /* create zinkblende structure */
+ /*
+#ifdef ALBE
+ r.x=0.5*0.25*ALBE_LC_SIC; r.y=r.x; r.z=r.x;
+ create_lattice(&md,FCC,ALBE_LC_SIC,SI,M_SI,
+ ATOM_ATTR_1BP|ATOM_ATTR_2BP|ATOM_ATTR_3BP|ATOM_ATTR_HB,
+ 0,LCNTX,LCNTY,LCNTZ,&r);
+ r.x+=0.25*ALBE_LC_SIC; r.y=r.x; r.z=r.x;
+ create_lattice(&md,FCC,ALBE_LC_SIC,C,M_C,
+ ATOM_ATTR_1BP|ATOM_ATTR_2BP|ATOM_ATTR_3BP|ATOM_ATTR_HB,
+ 1,LCNTX,LCNTY,LCNTZ,&r);
+#else
+ r.x=0.5*0.25*TM_LC_SIC; r.y=r.x; r.z=r.x;
+ create_lattice(&md,FCC,TM_LC_SIC,SI,M_SI,
+ ATOM_ATTR_1BP|ATOM_ATTR_2BP|ATOM_ATTR_3BP|ATOM_ATTR_HB,
+ 0,LCNTX,LCNTY,LCNTZ,&r);
+ r.x+=0.25*TM_LC_SIC; r.y=r.x; r.z=r.x;
+ create_lattice(&md,FCC,TM_LC_SIC,C,M_C,
+ ATOM_ATTR_1BP|ATOM_ATTR_2BP|ATOM_ATTR_3BP|ATOM_ATTR_HB,
+ 1,LCNTX,LCNTY,LCNTZ,&r);
+#endif
+ */
+
+ /* check for right atom placing */
+ moldyn_bc_check(&md);
+
+ /* testing configuration */
+ //r.x=0.27*sqrt(3.0)*LC_SI/2.0; v.x=0;
+ //r.x=(TM_S_SI+TM_R_SI)/4.0; v.x=0;
+ //r.y=0; v.y=0;
+ //r.z=0; v.z=0;
+ //add_atom(&md,SI,M_SI,0,
+ // ATOM_ATTR_1BP|ATOM_ATTR_2BP|ATOM_ATTR_3BP|ATOM_ATTR_HB,
+ // ATOM_ATTR_2BP|ATOM_ATTR_HB,
+ // &r,&v);
+ //r.x=-r.x; v.x=-v.x;
+ //r.y=0; v.y=0;
+ //r.z=0; v.z=0;
+ //add_atom(&md,SI,M_SI,0,
+ // ATOM_ATTR_1BP|ATOM_ATTR_2BP|ATOM_ATTR_3BP|ATOM_ATTR_HB,
+ // ATOM_ATTR_2BP|ATOM_ATTR_HB,
+ // &r,&v);
+ //r.z=0.27*sqrt(3.0)*LC_SI/2.0; v.z=0;
+ //r.x=(TM_S_SI+TM_R_SI)/4.0; v.x=0;
+ //r.y=0; v.y=0;
+ //r.x=0; v.x=0;
+ //add_atom(&md,SI,M_SI,0,
+ // ATOM_ATTR_1BP|ATOM_ATTR_2BP|ATOM_ATTR_3BP|ATOM_ATTR_HB,
+ // ATOM_ATTR_2BP|ATOM_ATTR_HB,
+ // &r,&v);
+ //r.z=-r.z; v.z=-v.z;
+ //r.y=0; v.y=0;
+ //r.x=0; v.x=0;
+ //add_atom(&md,SI,M_SI,0,
+ // ATOM_ATTR_1BP|ATOM_ATTR_2BP|ATOM_ATTR_3BP|ATOM_ATTR_HB,
+ // ATOM_ATTR_2BP|ATOM_ATTR_HB,
+ // &r,&v);
+
+ /* set temperature & pressure */
+ set_temperature(&md,atof(argv[2])+273.0);
+ set_pressure(&md,BAR);
+
+ /* set amount of steps to skip before average calc */
+ set_avg_skip(&md,AVG_SKIP);
+
+ /* set p/t scaling */
+ //set_pt_scale(&md,0,0,T_SCALE_BERENDSEN,100.0);
+ //set_pt_scale(&md,P_SCALE_BERENDSEN,0.001,
+ // T_SCALE_BERENDSEN,100.0);
+ //set_pt_scale(&md,0,0,T_SCALE_DIRECT,1.0);
+ //set_pt_scale(&md,P_SCALE_BERENDSEN,0.001,0,0);