On 06/03/2017 01:34 PM, Maximin Fotsing kuetche wrote:
I have in a .h5 file informations about a component Ez of electric
field. During simulations without (for incident field data) and with
structure I stored the field component in a line. I want to plot the
evolution of the field intensity for each point of the line.
Please can anyone guid? I also want to know how the fields are stored in
a .h5 file during the simulation. Is there any way to print them?
The Fourier-transformed fields computed when calling the "add-flux"
routine and output as an HDF5 file using "display-fluxes" is just a raw
dump of Meep's internal data. This data is not organized in a way which
is easily readable. We are planning to change this in a future release.
If you are interested in plotting a time evolution of the fields then
you just need to write your own run function, e.g.:
(define print-stuff
(lambda () (print "ez:, " (meep-time) ","
(get-field-point Ez (vector3 0 0 0)) "\n")))
(run-sources+ 100 (at-every 10 print-stuff))
This outputs the current simulation time and Ez at the center of the
computational cell at every 10 time units.
You can then use this information to compute the Fourier transform of
the fields in post processing. An example of this approach is shown in:
http://ab-initio.mit.edu/wiki/index.php/Meep_Tutorial/Multilevel-atomic_susceptibility
The Fourier-transformed fields are the steady-state fields which can be
used to obtain dependent quantities such as the field intensity.
_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
