-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1
Hello Steven and everybody!
I would like to let meep calculate the modes of disk/ringresonator in
cylindrical dimensions. 2D calculation runs fine, but when i try to add
a z-dimension, I always get the following output:
[EMAIL PROTECTED] ~/fdtd/disk_3d/mode_01 $ meep m=15 disk_3d.ctl
command-line param: m=15
- -----------
Initializing structure...
Working in Cylindrical dimensions.
block, center = (0,0,0)
size (1,1e+20,0.25)
axes (1,0,0), (0,1,0), (0,0,1)
dielectric constant epsilon = 12.25
meep: chunks[0] intersects with chunks[2]
What am I doing wrong?
The ctl-file looks like this:
(set! dimensions CYLINDRICAL)
(set-param! m 15)
(define-param n 3.5) ; index of waveguide
(define-param r 1.0) ; radius of disk
(define-param h 0.25) ; height of disk
(define-param pad 10) ; padding between waveguide and edge of PML
(define-param padz 5) ; padding between waveguide and edge of PML
(define-param dpml 2) ; thickness of PML
(define latticesizex (+ r pad dpml)) ; cell size
(define latticesizez (* 2 (+ (/ h 2) padz dpml))); cell size z-direction
(set! geometry-lattice (make lattice (size latticesizex no-size
latticesizez)))
;(set! geometry-lattice (make lattice (size latticesizex no-size no-size)))
(set! geometry (list
(make block (center 0 0 0) (size r infinity h) (material
(make dielectric (index n))))))
;(set! geometry (list
; (make block (center (/ r 2) 0 0) (size r infinity
infinity) (material (make dielectric (index n))))))
(set! pml-layers (list (make pml (thickness dpml))))
(set-param! resolution 10)
(define-param fcen 2094395102 ) ; pulse center frequency
(define-param df 209439510 ) ; pulse width (in frequency)
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ez) (center (- r 0.1) 0 0))))
(run-sources+ 300
(at-beginning output-epsilon)
(after-sources (harminv Ez (vector3 (- r 0.1)) fcen df)))
;(run-until (/ 1 fcen) (at-every (/ 1 fcen 20) output-efield-z))
(run-until (/ 1 fcen)
(in-volume (volume (center 0) (size (* 2 latticesizex)))
(at-beginning output-epsilon)
(to-appended "ez"
(at-every (/ 1 fcen 20) output-efield-z))))
Best regards,
Fabian
- --
_______________________________________
Fabian Wilde
Universitaet Hamburg
Institut fuer Angewandte Physik
Jungiusstrasse 11
20355 Hamburg
T.: +49 (40) 42838 - 2932
Fax: +49 (40) 42838 - 2931
email: [EMAIL PROTECTED]
_______________________________________
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.5 (GNU/Linux)
Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org
iD8DBQFE7HKhAxP20ZH7w5IRApUEAJ0VyhNZROPPFc0Wo1KcBOg0YMPG7ACeJzpY
CQBn11cEv5lliofqHGvw4DU=
=k9P/
-----END PGP SIGNATURE-----
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
