Re: [Meep-discuss] Missing ex values

[Zheng Li]

王健 wrote: Zheng Li，您好！       Thank you for your quick reply!        The attachfile is my ctl file and some results i have got.       Au_no.ctl : ctl file run without any structure       Au_no.out : output file(you can see that the total time step is 100.01 meep time unit, which corresponding to a 20002 outputs of ex )       Au_ex.hdf : Ex values is (16979)       .jpg      : fourior transformation of ex which shows a peak at 680nm is not 700nm as i have specified in the variable of fcen=700nm in                   my ctl file. 　　   ======== 2008-06-10 11:11:25 您在来信中写道： ========   王健 wrote: Hi, meepgroup, I have used meep for metal hole calculation. the resolution is set  to 100. the whole time for caculation is 100.01 meep time. And i want to get a point Ex value  at every 0.005 meep time. So the total Ex values will be 20002, but when i used matlab to investigate the hdf5 file,  ex is noly a 1* 16979 vectors. it is seems strange that (20002-16979) values of ex have been missing. In order to in vestigate this,  i have define a "my-time" value to output meep-time at every time step. It shows that the meep-time is right at every time step. So I think that at some time there is no output of ex. how does this happen? can anyone help me with this? Best wishes!                   Best regards [EMAIL PROTECTED] Jian Wang National Lab for Infrared Physics, Shanghai Institute of Technical Physics, Chinese Academy of Sciences, 500 Yutian Road, Shanghai, China 200083 _______________________________________________ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss Hi, I cannot give further helps on your topic without your script file. By the way, you may check your points-missing problem via a function ----------------------- (get-field-point c pt) ----------------------- Zheng Li 2008-6-10 = = = = = = = = = = = = = = = = = = = = = = 　　　　　　　　致 礼！   　　　　　　　　　　　　　　王健 　　　　　　　　　　　　　　 [EMAIL PROTECTED] 　　　　　　　　　　　　　　　2008-06-10   Hi， It is clear that the "time=" terms actually appear 16979 times in your "Au_no.out".  So It seems there are some float-round-errors occurring in the function *at-every* when you use the minimal time step.  In your case,  do not use the *at-every* function. Zheng Li 2008-6-10

_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss