Hi, I would rather guess the images Saket gets are really full black,
which would mean that there are "NaN" values in the output of the
computation. This may result from some numerical instability. I suggest
to try first to halve the simulation time until the field gets finite.
Then one can determine where the instability starts.
Filip
On 5.6.2016 16:09, David Kleiven wrote:
Hi
by black images, do you mean greyscale images?
If so, you can change the colormap by adding the -Zc option. To use
the dkbluered map you can add /-Zc dkbluered/ as option to the h5topng
command.
More info on colormaps:
http://ab-initio.mit.edu/wiki/index.php/Color_tables_in_h5topng
This example works for me:
h5topng file.h5 -Zc dkbluered
Hope this helps.
David
On Sat, Jun 4, 2016 at 7:25 PM, saket wankhade
<saketwankh...@gmail.com <mailto:saketwankh...@gmail.com>> wrote:
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 <mailto:meep-discuss@ab-initio.mit.edu>
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
