Ok. Problem solved by upgrading to meep-0.10. Regards, Fabian
Fabian Wilde wrote: > -----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
