155 |
A = 4.*d2 - d2p1 |
A = 4.*d2 - d2p1 |
156 |
B = 4.*d2p1 - d2 |
B = 4.*d2p1 - d2 |
157 |
C = d2 |
C = d2 |
158 |
D = d2p1; |
D = d2p1 |
159 |
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) |
160 |
A = 4.*d2m1 - d2 |
A = 4.*d2m1 - d2 |
161 |
B = 4.*d2 - d2m1 |
B = 4.*d2 - d2m1 |
162 |
C = d2m1 |
C = d2m1 |
163 |
D = d2; |
D = d2 |
164 |
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) |
165 |
!qMD = 0.5*( ( Qi + Qip ) - dp1h ) |
!qMD = 0.5*( ( Qi + Qip ) - dp1h ) |
166 |
qMD = 0.5*( ( 2.*Qi + DelIp ) - dp1h ) |
qMD = 0.5*( ( 2.*Qi + DelIp ) - dp1h ) |