Dear Meep users,
I use meep to calculate the omega~k relation of a waveguide. And meep
0.10 can work out the omega~k very well. However, when I use
meep-0.20.3, it can not at this time. Even when I change "dielectric" to
"medium" in the ctl file, the output will stop at "harminv1" when meep
is still running.
I wonder whether you have encountered the same kind of question. The
following is my ctl file.
Zhanghua
***********************************************************************************
(define-param sz 5) ; size of cell in z direction
(define-param dpml 0.5) ; PML thickness
(set! geometry-lattice (make lattice (size no-size no-size sz)))
(set! eps-averaging? false)
(define diel1 (make dielectric (epsilon 1)
(polarizations
(make polarizability
(omega 1e-20) (gamma 0) (delta-epsilon 1e40)))))
(define diel2 (make dielectric (epsilon 2.09)))
(set! geometry
(list (make block (center 0 0 0) (size infinity infinity sz)
(material diel1))
(make block (center 0 0 0) (size infinity infinity 0.1)
(material diel2))
))
(set! pml-layers (list (make pml (direction Z) (thickness dpml))))
(set-param! resolution 20)
(define-param fcen 0.4)
(define-param df 2.0)
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ez) (center 0 0 0.015))))
(define-param k-interp 50)
(run-k-points 300 (interpolate k-interp (list (vector3 0 0 0) (vector3
5 0 0))))
**********************************************************************************
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
