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;
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;