--- forces-edip.1.3.f	Fri Apr 19 10:42:49 2002
+++ forces-edip.1.4.f	Fri Apr 19 10:43:03 2002
@@ -2,7 +2,7 @@
 C   forces-edip.f
 C   -------------
 
-C   Version 1.3f
+C   Version 1.4f
 
 C   Force and Energy Calculation with the
 C   Environment-Dependent Interatomic Potential
@@ -358,10 +358,9 @@
                   
 C   COORDINATION AND NEIGHBOR FUNCTION par_c<r<par_b
         
-                  if(r .lt. par_b) then
-                    if(r .lt. par_c) then
+                  if(r .lt. par_c) then
                     Z = Z + 1.0d0
-                   else
+                  else if(r .lt. par_b) then
                     nz = nz + 1
                     numz(nz) = j
                     xinv = bmc/(r-par_c)
@@ -376,10 +375,8 @@
                     sz_dy(nz) = dy
                     sz_dz(nz) = dz
                      sz_r(nz) = r
-                   end if 
-C  r < par_C
                   end if 
-C  r < par_b
+C  r < par_c < par_b
                   end if 
 C  fixZ .eq. 0
                 end if