@ -150,7 +150,7 @@ module io
outn = outn + 1
end do
if((outn-atom_num) /= atom_num) then
if((outn-node_num) /= atom_num) then
print *, "outn", (outn-node_num), " doesn't equal atom_num ", atom_num
end if