fixed incorrect scaling of cos-acceleration viscosity
authorBerk Hess <hess@kth.se>
Wed, 8 May 2013 16:07:38 +0000 (18:07 +0200)
committerBerk Hess <hess@kth.se>
Wed, 8 May 2013 16:07:38 +0000 (18:07 +0200)
The cos-acceleration 1/viscosity output contained an anomalous factor
1/volume. This factor is now removed.
Fixes #1244

Change-Id: I9bf318b4e6557720683d50a0c2e887e306bacc3f

src/mdlib/mdebin.c

index 7fdfa7bba462a2e36fe343c8df787529c2cbdd16..0bd097945f2fc1de810db42ada767eb5b2f99c2b 100644 (file)
@@ -1091,7 +1091,7 @@ void upd_mdebin(t_mdebin       *md,
         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)