Dear Meep users:
I cannot input TEM wave.
The following is consists of 2 infinite metalic planes,
and inputs Ey and Hx components.
It doesn't have Ex component, but it outputs that.
What's wrong?
Best regards,
Yuichi
------------------------------------------------
int simple_metallic(const char *mydirname) {
grid_volume v = vol2d(20.0, 11.0, 10);
structure s(v, eps, pml(1.0), identity()); //eps is -HUGE_VAL in
outside.
s.set_output_directory(mydirname);
fields f(&s);
// f.set_boundary(High, R, Metallic);
f.output_hdf5(Dielectric, v.surroundings());
continuous_src_time src(0.28, 0.0);
volume src_line(vec(5.0, 1.0), vec(5.0, 10.0));
volume src_boundary1(vec(0.0, 1.0), vec(20.0, 1.0));
volume src_boundary2(vec(0.0, 10.0), vec(20.0, 10.0));
f.add_volume_source(Ey, src, src_line, sin_w, -1.0); //sin wave
f.add_volume_source(Ex, src, src_boundary1, zero, 0.0);
f.add_volume_source(Ez, src, src_boundary1, zero, 0.0);
f.add_volume_source(Ex, src, src_boundary2, zero, 0.0);
f.add_volume_source(Ez, src, src_boundary2, zero, 0.0);
f.add_volume_source(Hz, src, src_line, sin_w, 0.0026546);
while (f.time() < para.total_time)
f.step();
}
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
