diff --git a/src/interpolation.F90 b/src/interpolation.F90 index 2d7e9a2..42b0d61 100644 --- a/src/interpolation.F90 +++ b/src/interpolation.F90 @@ -69,7 +69,12 @@ module interpolation ds (i) = dvr(i) * dvl(i) if (ds(i) .gt. 0.0) then +#ifdef MINMOD + dv = sign(0.5, dvr(i)) * min(abs(dvr(i)), abs(dvl(i))) +#endif /* MINMOD */ +#ifdef LF dv = ds(i) / (dvr(i) + dvl(i)) +#endif /* LF */ vl(i) = vx(i) + dv vr(i) = vx(i) - dv