diff --git a/ARBTrajecExample.py b/ARBTrajecExample.py deleted file mode 100644 index fbd1a16..0000000 --- a/ARBTrajecExample.py +++ /dev/null @@ -1,29 +0,0 @@ -from __future__ import division -import numpy as np -from ARBTools.ARBInterp import tricubic -#from ARBTools.ARBTrajec -from A5_03032019 import trajectories - -############################### Skookum ############################### - -if __name__ == '__main__': - mAr = 39.948 # Mass Ar - N = 100 # number of atoms - T = 0.01 # temperature in K - - moment = 3 # magnetic moment in Bohr magnetons - mJ x g-factor (neglecting nuclear spin, use mF if I!=0) - # testing with Ar in 3P2 state: mJ = 2, low-field-seeking - - tmax = 0.01 # end time of simulation - timestep = 1e-6 # acceptable timestep depends on shape of potential! - - fieldname = "ExampleVectorField" - field = np.genfromtxt(fieldname+'.csv', delimiter=',') # field to be interpolated - - ArTest = trajectories(field) # creates class to iterate particle trajectories - ArSample = ArTest.atoms(T,N, mAr) # creates random sample of Ar atoms - np.savetxt('Rand_In.csv', ArSample, delimiter=',') # Save initial sample to file - - ArTest.Iterate(ArSample, tmax, timestep, mAr, moment) # tracks motion of particles through field for some time - - np.savetxt('Rand_Out.csv', ArSample, delimiter=',') # Save final sample to file