Mercurial > repos > public > sbplib

diff +scheme/Euler1d.m @ 78:80948a4084f3

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Added the factor or 1/2 to the wall bc term as well.
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 25 Nov 2015 15:41:10 +0100
parents 0b07ff8a0a12
children 2102af217134
line wrap: on
line diff
--- a/+scheme/Euler1d.m	Wed Nov 25 15:19:55 2015 +0100
+++ b/+scheme/Euler1d.m	Wed Nov 25 15:41:10 2015 +0100
@@ -412,7 +412,7 @@
                 w_in = w_s(p_in);
                 w_ot = w_s(p_ot);
 
-                o = obj.Hi * tau * (w_in - R*w_ot);
+                o = 1/2*obj.Hi * tau * (w_in - R*w_ot);
             end
 
             closure = @closure_fun;