Hi!

On Wed, Sep 15, 2010 at 10:15:20PM -0400, 月田翔 wrote:

> How should be done the output setting though it wants to know the
> time change of the electromagnetic field change by one point with
> the field as a numerical value when the electromagnetic field is
> analyzed with MEEP?

If I understand your question correctly, you want to plot E(\vec{r}, t) for a
fixed vector \vec{r} and all the times in your simulation.

One way to do it in the Scheme Meep is to define a function such as:

(define (my-output-trans)
    (print "transmitted field:\t " (meep-time) "\t "
    (get-field-point my-component meas-out-pt) "\n")
)

where my-component and meas-out-pt are the field component you are interested
in and the point in space (as a vector3) you want to "measure".

Then, when you run the simulation, you can run the simulation like:

(run-sources +
    (runtime)
    (at-beginning output-epsilon)
    (after-sources my-output-trans)
)

If you send Meep's output to a file, a 'grep transmitted' and a bit of
manipulation will give you a file with a column for t and a column for E(t).

There might be more efficient ways of doing it, though. The experts are
welcome! :-)

Bye & Good Luck!

Pablo B.


_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss

Reply via email to