QuantumATK W-2024.09 version released on Sep 9, 2024
0 Members and 1 Guest are viewing this topic.
# Import an ElectrostaticDifferencePotential objectpotential = nlread('results.nc', ElectrostaticDifferencePotential)[0]# Calculate the meanv_z = numpy.apply_over_axes(numpy.mean,potential[:,:,:],[0,1]).flatten()v_z *= potential[:,:,:].unit()# Print out the resultdX, dY, dZ = potential.volumeElement()dz = dZ.norm()shape = potential.shape()for i in range(shape[2]): print dz*i, v_z