45 |
|
|
46 |
|
!! @author Charles F. Vardeman II |
47 |
|
!! @author Matthew Meineke |
48 |
< |
!! @version $Id: doForces.F90,v 1.11 2005-03-08 21:05:46 gezelter Exp $, $Date: 2005-03-08 21:05:46 $, $Name: not supported by cvs2svn $, $Revision: 1.11 $ |
48 |
> |
!! @version $Id: doForces.F90,v 1.12 2005-03-21 20:51:06 chrisfen Exp $, $Date: 2005-03-21 20:51:06 $, $Name: not supported by cvs2svn $, $Revision: 1.12 $ |
49 |
|
|
50 |
|
|
51 |
|
module doForces |
847 |
|
if (PropertyMap(me_i)%is_Dipole) then |
848 |
|
|
849 |
|
mu_i = getDipoleMoment(me_i) |
850 |
< |
|
850 |
> |
|
851 |
|
!! The reaction field needs to include a self contribution |
852 |
|
!! to the field: |
853 |
|
call accumulate_self_rf(i, mu_i, eFrame) |