Field confinement

Dear All,

I have some questions regarding calculating the field confinement

1. I found this function (compute-energy-in-objects objects...) If I want to
calculate the energy in an object, for example an annulus
    (compute-energy-in-objects
                                  (make cylinder
                                  (center 0 0 0)
                                  (radius 1.0) (height infinity)
                                  (material 3.4))
                                  (center 0 0 0)
                                  (radius 0.5) (height infinity)
                                  (material air))
)
And I only want to know the energy in the region where the epsilon=3.4, can
I just add (compute-energy-in-dielectric 3.4 3.4) before that comment?
And also if I add (compute-dfield-energy) before computing the energy, does
it mean the calculated energy-in-object will be normalized by the total
energy in D?

2. When I use for example
(fix-efield-phase 1)
(output-efield-x 1)
the output file will be "e.k*.b1.x.te.h5"
How do I change the output file name?
for example change it to be "e.k*.b1_1.x.te.h5"?
Also when I use a loop

(do ((n 1 (+ n 1))) ((> n 10))
(fix-efield-phase n)
(output-efield-x n)
)

I only get one output h5 file, can anyone please tell me why?

Thanks a lot~
Mindy

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