Dear meep users,
I am trying to compute electric field distribution at
wavelength 550nm.My problem is that I get black images after applying
h5topng command.Does any one have idea how to deal with it ?
Here is my .ctl file;;;; a=1um ,lambda=550nm
************************************************************************************
(define-param resl 200);
(define-param dpml 0.5) ; thickness of PML layers
(define-param sx 0.07) ;
(define-param sz 0.07)
(define-param rad 0.03)
(set-param! k-point (vector3 0 0 0))
;(define sx0 (+ sx (* 2 dpml)))
(define sy 1);
(define sy0 (+ sy (* 2 dpml))) ;
(set! ensure-periodicity true)
(define background (make dielectric (epsilon 1.000 )))
(set! eps-averaging? false)
(set! force-complex-fields? false)
(set! Courant 0.5)
(set! output-single-precision? true)
(define gold (make dielectric (epsilon 1)
(polarizations
(make polarizability
(omega 1e-20) (gamma 0.042747) (sigma 4.0314e+41))
(make polarizability
(omega 0.33472) (gamma 0.19438) (sigma 11.363))
(make polarizability
(omega 0.66944) (gamma 0.27826) (sigma 1.1836))
(make polarizability
(omega 2.3947) (gamma 0.7017) (sigma 0.65677))
(make polarizability
(omega 3.4714) (gamma 2.0115) (sigma 2.6455))
(make polarizability
(omega 10.743) (gamma 1.7857) (sigma 2.0148))
)))
(set! geometry-lattice (make lattice (size sx sy0 sz)))
(set! geometry (list
(make block (center 0 0 0) (size sx sy sz) (material
background))
(make sphere (center 0 0 0) (radius rad) (material gold))
))
(set! pml-layers (list
(make pml (thickness dpml) (direction Y))
))
(set! sources (list
(make source
(src (make continuous-src (frequency 1.8181)))
(component Ex)
(center 0 0.3 0) (size sx 0 sz)
)
)
)
(set! resolution resl)
(define monitor-xy
(volume
(center 0 0 0)(size sx sy 0)
)
)
(define monitor-xz
(volume
(center 0 0 0)(size sx 0 sz)
)
)
(run-until 200
(in-volume monitor-xy (at-beginning output-epsilon))
(in-volume monitor-xy (to-appended "ex" (at-every 0.6 output-efield-x)))
)
*******************************************************************************
Any suggestion would be greatly appreciated.
Thank you!
_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
