There was a small output bug since Gromos96 bonds are modified to have
the reference value stored as (b0)^2 rather than just b0.
Fixes #1109.
Change-Id: Ie4516515f7cacb35d262c2059b5e461962340fec
b0 = 0;
switch (ftype) {
case F_BONDS:
b0 = 0;
switch (ftype) {
case F_BONDS:
b0 = idef->iparams[type].harmonic.rA;
b0 = idef->iparams[type].harmonic.rA;
+ break;
+ case F_G96BONDS:
+ b0 = sqrt(idef->iparams[type].harmonic.rA);
break;
case F_MORSE:
b0 = idef->iparams[type].morse.b0A;
break;
case F_MORSE:
b0 = idef->iparams[type].morse.b0A;