| 87 | 
  | 
    HydroProps cr; | 
| 88 | 
  | 
    HydroProps cd; | 
| 89 | 
  | 
    calcHydroPropsAtCR(beads_, viscosity, temperature, cr); | 
| 90 | 
< | 
    calcHydroPropsAtCD(beads_, viscosity, temperature, cd); | 
| 90 | 
> | 
    //calcHydroPropsAtCD(beads_, viscosity, temperature, cd); | 
| 91 | 
  | 
    setCR(cr); | 
| 92 | 
  | 
    setCD(cd); | 
| 93 | 
  | 
     | 
| 156 | 
  | 
             | 
| 157 | 
  | 
            Xiott += Cij; | 
| 158 | 
  | 
            Xiotr += U[i] * Cij; | 
| 159 | 
< | 
            Xiorr += -U[i] * Cij * U[j] + (6 * viscosity * volume) * I;             | 
| 159 | 
> | 
            //Xiorr += -U[i] * Cij * U[j] + (6 * viscosity * volume) * I;     | 
| 160 | 
> | 
            Xiorr += -U[i] * Cij * U[j];  | 
| 161 | 
  | 
        } | 
| 162 | 
  | 
    } | 
| 163 | 
  | 
 | 
| 320 | 
  | 
             | 
| 321 | 
  | 
            Xitt += Cij; | 
| 322 | 
  | 
            Xitr += U[i] * Cij; | 
| 323 | 
< | 
            Xirr += -U[i] * Cij * U[j] + (6 * viscosity * volume) * I;             | 
| 323 | 
> | 
            //Xirr += -U[i] * Cij * U[j] + (6 * viscosity * volume) * I;             | 
| 324 | 
> | 
            Xirr += -U[i] * Cij * U[j]; | 
| 325 | 
  | 
        } | 
| 326 | 
  | 
    } | 
| 327 | 
  | 
 |