def radialDistribution(this, bins, simulation = -1, states = "all", label = ''):
		print "Calculating radial distribution";
		simPath = this.simulations[simulation];
		length = this.stateNumbers[simulation];
		mdf.combine(simPath, length, );
		r, g = mda.calculateRadialDistribution(simPath, bins, length, states);
		mdplot.plotRadialDistribution(this.path, g, r, label = label);
		return r, g;
Пример #2
0
	def radialDistribution(this, bins, simulation = -1, states = "all", label = ''):
		print "Calculating radial distribution";
		simPath = this.simulations[simulation];
		length = this.stateNumbers[simulation];
		mdf.combine(simPath, length, );
		r, g = mda.calculateRadialDistribution(simPath, bins, length, states);
		
		
		datadirname = this.path + "/data/";
		if(not os.path.exists(datadirname)):
			os.makedirs(datadirname);
		#end
		
		np.save(datadirname + "/radial" + label + ".npy", r);
		np.save(datadirname + "/radialData" + label + ".npy", g);
		
		# mdplot.plotRadialDistribution(this.path, g, r, label = label);
		return r, g;