def BzAnalCircTCole(a, t, sigma):
wt, tbase, omega_int = setFrequency(t)
hz = HzanalCirc(sigma, omega_int/2/np.pi, 1., a, 'secondary')
# Treatment for inaccuracy in analytic solutions
ind = omega_int < 0.2
hz[ind] = 0.
hzTD, f0 = transFilt(hz, wt, tbase, omega_int, t)
return hzTD*mu_0
def dBzdtAnalCircTCole(a, t, sigma):
wt, tbase, omega_int = setFrequency(t)
hz = HzanalCirc(sigma, omega_int/2/np.pi, 1., a, 'secondary')
# Treatment for inaccuracy in analytic solutions
ind = omega_int < 0.2
hz[ind] = 0.
dhzdtTD = -transFiltImpulse(hz, wt, tbase, omega_int, t)
return dhzdtTD*mu_0