The output file is printing black images. I converted h5 to txt to check, the field values are NaN. Help me please.

Here is my code:

(set! geometry-lattice (make lattice (size 10 10 no-size)))
(set! pml-layers (list (make pml (thickness 1))))
(set-param! resolution 30)

(set! geometry
   (list
(make block (material (make dielectric (epsilon-diag 1.31 1.52 1.52) (epsilon-offdiag 1 0 0)))
        (center 0 0 0)
        (size 2.63 2.63 0)
        (e1 1 1 0) (e2 -1 1 0) (e3 0 0 1))
(make block (material (make dielectric (epsilon-diag 11.37 11.45 0.14) (epsilon-offdiag 11.37 0 0)))
        (center 0 0 0)
        (size 2 2 0)
        (e1 0.35 0.93 0)(e2 0.35 -0.93 0)(e3 0 0 1))
    (make block (material (make dielectric (epsilon 1)))
        (center 0 0 0)
        (size 1 1 0)
        (e1 1 1 0) (e2 1 -1 0) (e3 0 0 1))
))

(define-param fcen 0.005)

(set! sources (list (make source
                      (src (make continuous-src (frequency fcen) ))
                      (component Hz) (size 0 8 0)
            (center -3.9 0 0))))

(run-until 50
 (at-beginning output-epsilon)
    (to-appended "ez" (at-every .1 output-hfield-z))
)


_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss

Reply via email to