1724       REAL(KIND=dp), DIMENSION(3)                        :: rcur  local
1733             rcur(:) = ri(:) - rj(:)
1734             CALL helium_pbc(helium, rcur)
1735             wnum(:) = wnum(:) + rcur(:)
1740          rcur(:) = ri(:) - rj(:)
1741          CALL helium_pbc(helium, rcur)
1742          wnum(:) = wnum(:) + rcur(:)
1835       REAL(KIND=dp), DIMENSION(3)                        :: rcur  local
1847             rcur(:) = ri(:) - rj(:)
1848             CALL helium_pbc(helium, rcur)
1849             wnum(:) = wnum(:) + rcur(:)
1861          rcur(:) = ri(:) - rj(:)
1862          CALL helium_pbc(helium, rcur)
1863          wnum(:) = wnum(:) + rcur(:)
1921       REAL(KIND=dp), DIMENSION(3)                        :: r1, r12, r2, rcur  local
1935             rcur(1) = r1(2)*r2(3) - r1(3)*r2(2)
1936             rcur(2) = r1(3)*r2(1) - r1(1)*r2(3)
1937             rcur(3) = r1(1)*r2(2) - r1(2)*r2(1)
1938             area(:) = area(:) + rcur(:)
1950          rcur(1) = r1(2)*r2(3) - r1(3)*r2(2)
1951          rcur(2) = r1(3)*r2(1) - r1(1)*r2(3)
1952          rcur(3) = r1(1)*r2(2) - r1(2)*r2(1)
1953          area(:) = area(:) + rcur(:)
2053       REAL(KIND=dp), DIMENSION(3)                        :: rcur, rsum  local
2065             rcur(1) = ri(2)*rj(3) - ri(3)*rj(2)
2066             rcur(2) = ri(3)*rj(1) - ri(1)*rj(3)
2067             rcur(3) = ri(1)*rj(2) - ri(2)*rj(1)
2068             rsum(:) = rsum(:) + rcur(:)
2080          rcur(1) = ri(2)*rj(3) - ri(3)*rj(2)
2081          rcur(2) = ri(3)*rj(1) - ri(1)*rj(3)
2082          rcur(3) = ri(1)*rj(2) - ri(2)*rj(1)
2083          rsum(:) = rsum(:) + rcur(:)
2109       REAL(KIND=dp), DIMENSION(3)                        :: r1, r12, r2, rcur  local
2124             rcur(1) = r1(2)*r2(3) - r1(3)*r2(2)
2125             rcur(2) = r1(3)*r2(1) - r1(1)*r2(3)
2126             rcur(3) = r1(1)*r2(2) - r1(2)*r2(1)
2127             area(:) = area(:) + rcur(:)
2143          rcur(1) = r1(2)*r2(3) - r1(3)*r2(2)
2144          rcur(2) = r1(3)*r2(1) - r1(1)*r2(3)
2145          rcur(3) = r1(1)*r2(2) - r1(2)*r2(1)
2146          area(:) = area(:) + rcur(:)
2203       REAL(KIND=dp), DIMENSION(3)                        :: com, r1, r12, r2, rcur  local
2219             rcur(1) = r1(2)*r2(2) + r1(3)*r2(3)
2220             rcur(2) = r1(3)*r2(3) + r1(1)*r2(1)
2221             rcur(3) = r1(1)*r2(1) + r1(2)*r2(2)
2222             moit(:) = moit(:) + rcur(:)
2234          rcur(1) = r1(2)*r2(2) + r1(3)*r2(3)
2235          rcur(2) = r1(3)*r2(3) + r1(1)*r2(1)
2236          rcur(3) = r1(1)*r2(1) + r1(2)*r2(2)
2237          moit(:) = moit(:) + rcur(:)
2341       REAL(KIND=dp), DIMENSION(3)                        :: com, rcur, ri, rj  local
2353             rcur(1) = ri(2)*rj(2) + ri(3)*rj(3)
2354             rcur(2) = ri(3)*rj(3) + ri(1)*rj(1)
2355             rcur(3) = ri(1)*rj(1) + ri(2)*rj(2)
2356             moit(:) = moit(:) + rcur(:)
2369          rcur(1) = ri(2)*rj(2) + ri(3)*rj(3)
2370          rcur(2) = ri(3)*rj(3) + ri(1)*rj(1)
2371          rcur(3) = ri(1)*rj(1) + ri(2)*rj(2)
2372          moit(:) = moit(:) + rcur(:)