printf("total i: %f %f %f\n",ai->f.x,ai->f.y,ai->f.z);
if(aj==&(moldyn->atom[0]))
printf("total j: %f %f %f\n",aj->f.x,aj->f.y,aj->f.z);
printf("total i: %f %f %f\n",ai->f.x,ai->f.y,ai->f.z);
if(aj==&(moldyn->atom[0]))
printf("total j: %f %f %f\n",aj->f.x,aj->f.y,aj->f.z);