153 |
A = 4.*d2 - d2p1 |
A = 4.*d2 - d2p1 |
154 |
B = 4.*d2p1 - d2 |
B = 4.*d2p1 - d2 |
155 |
C = d2 |
C = d2 |
156 |
D = d2p1; |
D = d2p1 |
157 |
dp1h = max(min(A,B,C,D),0. _d 0)+min(max(A,B,C,D),0. _d 0) |
dp1h = max(min(A,B,C,D),0. _d 0)+min(max(A,B,C,D),0. _d 0) |
158 |
A = 4.*d2m1 - d2 |
A = 4.*d2m1 - d2 |
159 |
B = 4.*d2 - d2m1 |
B = 4.*d2 - d2m1 |
160 |
C = d2m1 |
C = d2m1 |
161 |
D = d2; |
D = d2 |
162 |
dm1h = max(min(A,B,C,D),0. _d 0)+min(max(A,B,C,D),0. _d 0) |
dm1h = max(min(A,B,C,D),0. _d 0)+min(max(A,B,C,D),0. _d 0) |
163 |
!qMD = 0.5*( ( Qi + Qip ) - dp1h ) |
!qMD = 0.5*( ( Qi + Qip ) - dp1h ) |
164 |
qMD = 0.5*( ( 2.*Qi + DelIp ) - dp1h ) |
qMD = 0.5*( ( 2.*Qi + DelIp ) - dp1h ) |