--- trunk/OOPSE/libmdtools/calc_dipole_dipole.F90 2003/08/27 16:16:01 727 +++ trunk/OOPSE/libmdtools/calc_dipole_dipole.F90 2003/08/27 16:25:11 730 @@ -176,15 +176,15 @@ contains t(3,atom2) = t(3,atom2) - ul2(1)*dudu2y + ul2(2)*dudu2x #endif + if (do_stress) then + #ifdef IS_MPI id1 = tagRow(atom1) id2 = tagColumn(atom2) #else id1 = atom1 id2 = atom2 -#endif - if (do_stress) then - +#endif if (molMembershipList(id1) .ne. molMembershipList(id2)) then ! because the d vector is the rj - ri vector, and