void ComputeTempAsphere::dof_compute()

{

  fix_dof = modify->adjust_dof_fix(igroup);

  adjust_dof_fix();

  // 6 dof for 3d, 3 dof for 2d

  // which dof are included also depends on mode

  MPI_Allreduce(&t,&scalar,1,MPI_DOUBLE,MPI_SUM,world);

  if (dynamic || tempbias == 2) dof_compute();

  if (tfactor == 0.0 && scalar != 0.0) 

  if (tfactor == 0.0 && atom->natoms != 0) 

    error->all(FLERR,"Temperature compute degrees of freedom < 0");

  scalar *= tfactor;

  return scalar;