On 05/11/2012 06:46 PM, schnetter(a)cct.lsu.edu wrote:
> @@ -478,7 +480,7 @@
> end do
> else !!$ Really implement C&W, page 197; which requires stencil 4.
> do i = 4, nx - 3
> - s=sign(1.d0, -dpress(i))
> + s=sign(one, -dpress(i))
> flatten = max(tilde_flatten(i), tilde_flatten(i+s))
> if (abs(1.d0 - flatten)> 0.d0) then
> trivial_rp(i-1) = .false.
Hi Erik,
you missed the line above: if (abs(1.d0 - flatten)> 0.d0) then
Cheers,
Bruno.