p = [3, 2, 2];
else:
print "Unsupported number(" + str(mpi.cpu_count()) + ") of cpu cores, needs 8 or 12";
exit();
#end
c = 10;
b = usc.h2oDensityToCellSize(1);
print b;
c = np.array([c, c, c]);
s = b*c;
mdm = mdm();
mdm.compile();
mdm.setDir(path);
mdm.setProcessors(p);
mdm.setVolume(s);
mdm.setTemperature(300);
mdm.makeWater(c);
mdm.thermalize(therm, dt, 25*dt, 10, thermstep);
mdm.simulate(stab, dt, stabstep);
mdm.noseHoover(sim, dt, 25*dt, step);
mdm.makeVmdData();
# mdm.visualize();
mdm.radialDistribution(1024, label = "amorph", states = range(0,stableStates));
mdm.angularDistribution(1024, label = "amorph", states = range(0,stableStates));
