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