On Wed, 24 Oct 2007, shi lina wrote:
> I want to output the field at a fixed point from the beginning
> (time=0) to the end (time=end). I look up the tutorial and can't find the
> command. Can you give me some suggestions? Thank you very much!
There are several ways to do this.
For example, you can call "get-field-point" to get the field at a point,
and write a little step function to call this at each time and print it
out:
(define (my-print-point)
(print "point:, " (meep-time) ", "
(get-field-point Ez (vector3 1 2 3)) "\n"))
(run-until 100 my-print-point)
Modify as needed to print different field components, at a point different
from (1,2,3), etcetera.
You can also output to a single HDF5 file by using
(to-appended "my-point" (in-volume (volume (center 1 2 3)
(size 0 0 0))
(output-efield-z)))
and pass this to (run-until ...) or whatever. That will create a file
my-point.h5 consisting of a 1d dataset of the Ez field at (1,2,3) as a
function of time.
All of the above are described in detail by the manual, of course.
Regards,
Steven G. Johnson
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
