import loopfield
# compensation coils:
# r_in = 5mm, r_out = 30.25mm
# d_coil = 5.3mm
# w = 25mm + d_coil/2
# N_loop_1 = 9530
# N_loop_2 = 10623
# d_nomm = 100μm
# coil axis in z-direction
# sample midpoint in (0,0,0)
# untilted sample normal to (0,0,1)
field = loopfield.Field()
def save_3d_file(output_file, data, header):
fh = open(output_file, 'w')
fh.write(header + "\n")
shape = data.shape
for i in range(shape[0]):
block = data[i]
np.savetxt(fh, block, fmt="%.17g", delimiter="\t")
fh.write("\n")
fh.close()
def add_coil(pos,
N,
#!/usr/bin/env python3
#
# Maxwell coil plot example
#
import math
import loopfield as lf
import loopfield.plot as lfp
# field object
field = lf.Field(length_units=lf.cm, current_units=lf.A, field_units=lf.uT)
# Maxwell coil model with single current loops
R = 10
# center winding
c1 = lf.Loop([0, 0, 0], [1, 0, 0], R, 64)
# outer windings
c2 = lf.Loop([-R * math.sqrt(3. / 7.), 0, 0], [1, 0, 0],
R * math.sqrt(4. / 7.), 49)
c3 = lf.Loop([+R * math.sqrt(3. / 7.), 0, 0], [1, 0, 0],
R * math.sqrt(4. / 7.), 49)
# add windings to field
field.addLoop(c1)
field.addLoop(c2)
field.addLoop(c3)