| Hi Alex, I believe you need to change your size in the y direction to infinity rather than no-size. (set! geometry (list (make block (center comb_width no-size (- 1 (/ comb_dept 2)) ) (size comb_width infinity comb_dept) (material (make dielectric (epsilon 1.8))) ) (make block (center 0.5 no-size (/ (- 1 comb_dept) 2) ) (size 1.0 infinity (- 1 comb_dept) ) (material (make dielectric (epsilon 1.8))) ) )) Attached is a .png of the permittivity with this change.  Best wishes, Tim Tim Saucer Physics 236/241/261 Lead GSI Sih Research Group PhD Candidate Office: 4409 Randall, Lab: 1277 Randall Department of Physics, University of Michigan Alex <to.sanja@...> writes: > > Hello. > I try to define geomentry of 2D structure with relative permittivity > 1.8. But then I look h5 output file by h5totxt I see one "1,1,1,1". > Ctl file: > (set! geometry-lattice (make lattice (size 1 no-size 1))); Specifies > the size of the unit cell > (define comb_dept (/ 9 10)) > (define comb_width (/ 5 10) ) > > (set! geometry (list > (make block (center comb_width no-size (- 1 (/ comb_dept 2)) ) (size > comb_width no-size comb_dept) (material (make dielectric (epsilon > 1.8))) ) > (make block (center 0.5 no-size (/ (- 1 comb_dept) 2) ) (size 1.0 > no-size (- 1 comb_dept) ) (material (make dielectric (epsilon 1.8))) ) > )) > (display geometry) > (set! resolution 32) > (set! k-point (vector3 1 0 0)) > (set! ensure-periodicity true) > (define-param fcen 0.5) > (define-param df 3.0) > (set! sources (list (make source (src (make gaussian-src (frequency > fcen) (fwidth df))) (component Hz) (center 0.1234 0)))) > > (set-param! k-point (vector3 0.4 0)) > (run-sources+ 300 (after-sources (harminv Hz (vector3 0.1234) fcen df))) > (define-param k-interp 50) > ;(run-k-points 300 (interpolate k-interp (list (vector3 0) (vector3 0.5)))) > (run-k-points 300 (append > (interpolate 15 (list (vector3 0) (vector3 0.1))) > (interpolate 5 (list (vector3 0.1) (vector3 0.4))) > (interpolate 15 (list (vector3 0.4) (vector3 0.5))) > ) ) > (run-until 200 (at-beginning output-epsilon) (at-end > output-efield-z)) > > OUTPUT: > [alex <at> nout meep]$ h5totxt 1d_walls-eps-000000.00.h5 > 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 > .... > > Alex Friman > |
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
