+ /* 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_mean_skip(&md,1000);
+
+ /* 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);