ff1.make_xz_spectral()
ff2.make_xz_spectral()
#===================================================================#
#### ---- Stack velocity fields in the wall-normal direction ####
#===================================================================#
ff1.stack_ff_in_y()
ff2.stack_ff_in_y()
#===================================================================#
#### Calculate amplitude/phase difference ####
#===================================================================#
tolerance = 1e-9
print("==========================================")
print("Differences between " + args.File1[:-3] + " and " + args.File2[:-3])
print("==========================================")
message = "Amplitude difference between " + args.File1[:-3] + " and " + args.File2[:-3]
diff_A = Tests.difference(abs(ff1.velocityField), abs(ff2.velocityField), tolerance, message)
diff_P = mean(angle(ff1.velocityField,deg=True) - angle(ff2.velocityField,deg=True))
print("Amplitude diff \t\t%.2E" % diff_A)
print("")
print("Field 1 (deg) \t\t%.6f (avg)" % mean(angle(ff1.velocityField,deg=True)))
print("Field 2 (deg) \t\t%.6f (avg)" % mean(angle(ff2.velocityField,deg=True)))
print("Phase diff (deg) \t%.6f (avg)" % diff_P)
print("")
#===================================================================#
#### Euler form of Fourier transformed field ####
#===================================================================#
print("==========================================")
print("Reconstructed using Euler formula")
print("==========================================")
test_ff1 = abs(ff1.velocityField) * exp(1j * angle(ff1.velocityField))
test_ff2 = abs(ff2.velocityField) * exp(1j * angle(ff2.velocityField))
