hello meep user

I need some help with h5topng
I do FDTD simulation by Meep with photonic crystal structure structure. This is 
my ctl file 

; some parameters for geometry
(define notomi (make dielectric (index 12.96)))
(define-param Nx 30) ; separuh nilai genap dari kolom silinder
(define-param Ny 50) ; separuh nilai genap dari baris silinder 
(define-param r 0.35) ; jari-jari silinder

; some parameters for computational cell
(define-param tepiy 10); y-length between PML and slab
(define-param tepix 5); x-length between PML and slab
(define-param dpml 1); PML thickness
(define-param sy (+ (* 2 tepiy) (* 2 Ny (/ (sqrt 3) 2)) dpml)) ; cell size in y 
direction

(define-param sx (+ (* 2 tepix) Nx dpml dpml)) ; cell size in x direction

(set! geometry-lattice (make lattice (size sx sy no-size)))

;struktur definition
(set! default-material air) ; everything is air
(set! geometry 
   (append 
                 (geometric-objects-duplicates (vector3 0 (sqrt 3) 0) 0 Ny
                        (geometric-objects-duplicates (vector3 1 0 0) 0 Nx
                                 (list
                                  (make cylinder (center (+ (+ (* (/ sx 2) -1) 
dpml) tepix) (+ (+ (* (/ sy 2) -1) dpml) tepiy) 0) (radius r) (height infinity) 
  
                               (material notomi)))))

             (geometric-objects-duplicates (vector3 0 (sqrt 3) 0) 0 (- Ny 1)
                         (geometric-objects-duplicates (vector3 1 0 0) 0 (- Nx 
1)
                                 (list
                                       (make cylinder (center (+ (+ (+ (* (/ sx 
2) -1) dpml) 0.5) tepix) (+ (+ (+ (* (/ sy 2) -1) dpml) tepiy) (/ (sqrt 3) 2)) 
0) (radius r) (height infinity)   
                               (material notomi)))))))

;PML and resolution
(set! pml-layers (list (make pml (thickness dpml))))
(set! resolution 20)

;source
(define-param fcen 0.50)
(define-param df 0.2) 
(define-param w 4) 
(set! sources (list
               (make source
                 (src (make gaussian-src (frequency fcen) (fwidth df)))
                 (component Ey)
                 (center 0 0)
                 (size 0 w))))

(run-until 200
           (at-beginning output-epsilon)
           (at-end output-efield-z)))  


Then I try to see the picture by command " h5topng -t 0:30 -R -Zc dkbluered -a 
yarg -A cobamodif1-eps-000000.00.h5 cobamodif1-ez.h5"
 and I only got pictures of totally black rectangular. Is there something wrong 
with my ctl file or my command. I urgently need your suggestion. Thank you for 
the help





      
____________________________________________________________________________________
Be a better friend, newshound, and 
know-it-all with Yahoo! Mobile.  Try it now.  
http://mobile.yahoo.com/;_ylt=Ahu06i62sR8HDtDypao8Wcj9tAcJ
_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss

Reply via email to