Commit b640c79d authored by Thomas Jahns's avatar Thomas Jahns 🤸
Browse files

Fix bug in handling of first element.

parent 7bf866a9
......@@ -319,7 +319,8 @@ CONTAINS
CALL assertion(np == SIZE(em), msg=context//': np /= SIZE(em)')
IF (pload(1) < 0.0_dp) THEN
DO ip=1,np
pload(1) = SUM(wload(1:em(1)))
DO ip=2,np
pload(ip)=SUM(wload(em(ip-1)+1:em(ip)))
ENDDO
ENDIF
......
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