diff -Naur lammps-14Sep09/lib/meam/meam_dens_init.F lammps-16Sep09/lib/meam/meam_dens_init.F
--- lammps-14Sep09/lib/meam/meam_dens_init.F	2009-08-18 11:02:09.000000000 -0600
+++ lammps-16Sep09/lib/meam/meam_dens_init.F	2009-09-16 08:58:03.000000000 -0600
@@ -158,7 +158,7 @@
               cikj = (2.d0*(xik+xjk) + a - 2.d0)/a
               Cmax = Cmax_meam(elti,eltj,eltk)
               Cmin = Cmin_meam(elti,eltj,eltk)
-              if (cikj.ge.Cmax.or.cikj.lt.0.d0) then
+              if (cikj.ge.Cmax) then
                 goto 10
 c     Note that we never have 0<cikj<Cmin here, else sij=0 (rejected above)
               else
@@ -368,7 +368,7 @@
         cikj = (2.d0*(xik+xjk) + a - 2.d0)/a
         Cmax = Cmax_meam(elti,eltj,eltk)
         Cmin = Cmin_meam(elti,eltj,eltk)
-        if (cikj.ge.Cmax.or.cikj.lt.0.d0) then
+        if (cikj.ge.Cmax) then
           goto 10
         else if (cikj.le.Cmin) then
           sij = 0.d0