name = "oned-test"
outdir = "result/%s" % name
prepare(outdir)
string = String(31)
string.source = source.HardSource(source.sin_oft, (f,), string.shape[0]/2)
hdf5 = h5py.File("%s/%s.hdf5" % (outdir, name),"w")
hdf5.attrs["name"] = name
hdf5.attrs["dt"] = dt
hdf5.attrs["dx"] = dx
hdf5.require_group("timeline")
for t in range(0,300):
string.update_dfield()
string.update_efield()
string.update_abc()
# string.update_source(t*dt)
string.efield[15] = source.sin_oft(t*dt, f)
string.update_bfield()
string.update_hfield()
# plot(string.efield, "/tmp/%.3d.png", t)
hdf5.require_group("timeline/"+str(t))
hdf5["timeline"][str(t)]["ex"] = string.efield
hdf5["timeline"][str(t)]["hy"] = string.hfield
print(t)
hdf5.close()
