diff --git a/sources/interpolations.F90 b/sources/interpolations.F90 index 3a4c406..3aef5f1 100644 --- a/sources/interpolations.F90 +++ b/sources/interpolations.F90 @@ -126,7 +126,7 @@ module interpolations ! flags for reconstruction corrections ! logical , save :: mlp = .false. - logical , save :: positivity = .false. + logical , save :: positivity = .true. logical , save :: clip = .false. ! interpolation coefficients @@ -208,7 +208,7 @@ module interpolations character(len=255) :: plimiter = "mm" character(len=255) :: climiter = "mm" character(len=255) :: mlp_limiting = "off" - character(len=255) :: positivity_fix = "off" + character(len=255) :: positivity_fix = "on" character(len=255) :: clip_extrema = "off" character(len= 16) :: stmp real(kind=8) :: cfl = 0.5d+00