add_ebin(md->ebin, md->ivcos, 1, &(ekind->cosacc.vcos), bSum);
/* 1/viscosity, unit 1/(kg m^-1 s^-1) */
tmp = 1/(ekind->cosacc.cos_accel/(ekind->cosacc.vcos*PICO)
- *dens*vol*sqr(box[ZZ][ZZ]*NANO/(2*M_PI)));
+ *dens*sqr(box[ZZ][ZZ]*NANO/(2*M_PI)));
add_ebin(md->ebin, md->ivisc, 1, &tmp, bSum);
}
if (md->nE > 1)