bug after r26901: rotational entropy missing from output
Gets Around
12:57:09 PM PDT - Sat, May 23rd 2015
In the two recent snapshot releases NWChem always reports rotational entropy as 0.0. The problem was a mistake made during refactoring src/vib/vib_wrtfreq.F to support single atoms. The inertprod value is now initialized to zero and never updated.

If you put
inertprod=inertia_A*inertia_B*inertia_C


back in before
if (dabs(inertprod).lt.1.0d-9) then


around line 238 of vib_wrtfreq.F, the code properly reports rotational entropy again.

Forum Vet
10:11:21 AM PDT - Tue, May 26th 2015
Committed to to the svn repository.
Thanks a lot
Quote:Mernst May 23rd 11:57 am
In the two recent snapshot releases NWChem always reports rotational entropy as 0.0. The problem was a mistake made during refactoring src/vib/vib_wrtfreq.F to support single atoms. The inertprod value is now initialized to zero and never updated.

If you put
inertprod=inertia_A*inertia_B*inertia_C


back in before
if (dabs(inertprod).lt.1.0d-9) then


around line 238 of vib_wrtfreq.F, the code properly reports rotational entropy again.


Forum >> NWChem's corner >> General Topics