inconsistance with radiation force on each dipole along different axis

[GoogleCodeExporter]

Hallo,
I am trying to calculate the force on each dipole of a sphere due to gaussian 
beam and I have used he following command:
 adda -m 1.2 0 -lambda 0.3865 -beam lminus 1.2 0 0 0  -eq_rad 1 -Cpr_mat -store_force

and the out put stored on the file name VisFrp-Y , consist six column files(if 
I am right) the first three represent the coordinate of each dipole and the 
last three column for the corresponding force on each dipole. For such a sphere 
located at the beam waist and with the default beam propagation direction along 
z-axis (0,0,1), the net forces along the x, and y direction should be zero and 
the one along the beam axis should be different from zero but when I sum up the 
net forces on each dipole the force along z and y are close to zero where as 
the the fx is not. 
I would like to know if any had the same problem?
I apologize for long mail.

Besira 

Original issue reported on code.google.com by besirat...@gmail.com on 28 Jun 2011 at 12:14

[GoogleCodeExporter]

Besira, please check again your processing of VisFrp-Y.dat.

I have run the same simulation and then summed corresponding columns in 
VisFrp-Y.dat. I have obtained:
Sum{Fx,Fy,Fz}={1.90169E-16,-1.09397E-16,0.01142} (as it should be), contrary to 
your statement. Moreover, it agrees (qualitatively) with 
Cpr = (7.721251857e-017,1.873792634e-016,0.2870429707) (from file CrossSec-Y).

Original comment by yurkin on 29 Jun 2011 at 5:12

[GoogleCodeExporter]

Thank you Yurkin,

It is right, it works fine. I made mistake while summation.
Thanks. 

Besira

Original comment by besirat...@gmail.com on 29 Jun 2011 at 9:25

[GoogleCodeExporter]

Original comment by yurkin on 29 Jun 2011 at 10:15

• Changed state: Invalid