When I look into .h5, I think I get the point. The reason is the difference between the epsilon of Ag and air is too small.
----- Original Message ---- From: Andrew Andrew <[EMAIL PROTECTED]> To: meep-discuss@ab-initio.mit.edu Sent: Thursday, April 19, 2007 3:14:44 PM Subject: [Meep-discuss] why png is black? Dear Steven, My codes is based on an example which is used to solve the bandstructure of holey waveguide http://ab-initio.mit.edu/wiki/index.php/Meep_Tutorial/Band_diagram%2C_resonant_modes%2C_and_transmission_in_a_holey_waveguide What I do is changing Si to Ag. The codes is as follows, ------------------------------------------------------------------------------------------------- ; Some parameters to describe the geometry: (define-param eps 13) ; dielectric constant of waveguide (define-param w 1.2) ; width of waveguide (define-param r 0.36) ; radius of holes ; The cell dimensions (define-param sy 12) ; size of cell in y direction (perpendicular to wvg.) (define-param dpml 1) ; PML thickness (y direction only!) (set! eps-averaging? false) (set! geometry-lattice (make lattice (size 1 sy no-size))) (define-param Ag (make dielectric (epsilon 1) (polarizations (make polarizability (omega 1e-20) (gamma 0) (delta-epsilon 1e40)) ))) (set! geometry (list (make block (center 0 0) (size infinity w infinity) (material Ag)) (make cylinder (center 0 0) (radius r) (height infinity) (material air)))) (set! pml-layers (list (make pml (direction Y) (thickness dpml)))) (set-param! resolution 40) (define-param fcen 0.25) ; pulse center frequency (define-param df 1.5) ; pulse freq. width: large df = short impulse (set! sources (list (make source (src (make gaussian-src (frequency fcen) (fwidth df))) (component Hz) (center 0.1234 0)))) (set! symmetries (list (make mirror-sym (direction Y) (phase -1)))) (define-param kx false) ; if true, do run at specified kx and get fields (define-param k-interp 19) ; # k-points to interpolate, otherwise (if kx (begin (set! k-point (vector3 kx)) (run-sources+ 300 (at-beginning output-epsilon) (after-sources (harminv Hz (vector3 0.1234 0) fcen df))) (run-until (/ 1 fcen) (at-every (/ 1 fcen 20) output-hfield-z))) (run-k-points 300 (interpolate k-interp (list (vector3 0) (vector3 0.5))))) ------------------------------------------------------------------------------------------------------------ Why holey-wvg-bands-eps-000000.00.png is black? If I use PEC, I can get the structure image. Thank you very much! mario Ahhh...imagining that irresistible "new car" smell? Check out new cars at Yahoo! Autos. _______________________________________________ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss __________________________________________________ Do You Yahoo!? Tired of spam? Yahoo! Mail has the best spam protection around http://mail.yahoo.com
_______________________________________________ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss