Commit 64006d88 authored by Sven Harig's avatar Sven Harig
Browse files

small fix in energy calculation

parent 7d8a66de
......@@ -1255,7 +1255,8 @@ subroutine compute_energy
use PARAMETERS
!
implicit none
integer :: el, elnodes(3), depH, a, b
integer :: el, elnodes(3)
real(kind=wp) :: depH, a, b
real(kind=wp), parameter :: one_third=1._wp/3._wp, plop0=1026.0_wp
!$OMP MASTER
......@@ -1281,12 +1282,12 @@ subroutine compute_energy
enddo
!$OMP END DO
!$OMP MASTER
Epot = 0.5_wp*g*plop0* Epot
Ekin = 0.5_wp* plop0* Ekin
!$OMP END MASTER
end subroutine compute_energy
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment