diff --git a/src/motion/pint_pile.F b/src/motion/pint_pile.F index bf58419353..405dbbdbd0 100644 --- a/src/motion/pint_pile.F +++ b/src/motion/pint_pile.F @@ -89,7 +89,7 @@ CONTAINS pile_therm%g_fric(1) = 1.0_dp/pile_therm%tau END IF DO i = 2, p - pile_therm%g_fric(i) = 2.0_dp*pile_therm%lamb*normalmode_env%lambda(i) + pile_therm%g_fric(i) = 2.0_dp*pile_therm%lamb*SQRT(normalmode_env%lambda(i)) END DO DO i = 1, p ex = -dti2*pile_therm%g_fric(i) diff --git a/tests/Pimd/TEST_FILES b/tests/Pimd/TEST_FILES index b608c33cdd..b5276f6e8b 100644 --- a/tests/Pimd/TEST_FILES +++ b/tests/Pimd/TEST_FILES @@ -16,11 +16,11 @@ water-in-helium-lastforce.inp 9 7e-14 1.071 water-micro-helium.inp 9 1e-13 1.0713355696395585E-003 w512_pint_nose.inp 9 8e-14 9.5780711530210407 w512_pint_gle.inp 9 3e-13 -3.5646497002020765 -w512_pint_pile.inp 9 3e-12 0.64599197840093936 +w512_pint_pile.inp 9 3e-12 0.86626194213382846 w512_pint_piglet.inp 9 1e-11 74.229316490453456 centroid_velocity_init.inp 9 1.0E-14 0.32999751991202331 he32_density.inp 40 2e-14 -6.8167031351590079E-005 he32_only_worm.inp 40 1e-11 1.9115591620126446 he32_only_worm_restart.inp 40 1e-11 1.7411512073586795 -water-in-helium-worm.inp 9 1e-11 1.1453251192918812 +water-in-helium-worm.inp 9 1e-11 1.1453624372292848 #EOF