diff --git a/sources/evolution.F90 b/sources/evolution.F90 index 4de16dd..d5a8e47 100644 --- a/sources/evolution.F90 +++ b/sources/evolution.F90 @@ -231,7 +231,6 @@ module evolution evolve => evolve_ssprk34 order = 3 registers = 2 - cfl = 2.0d+00 * cfl name_int = "3ʳᵈ order 4-step SSPRK(4,3)" case ("ssprk(5,3)", "SSPRK(5,3)", "rk3.5", "RK3.5") @@ -239,7 +238,6 @@ module evolution evolve => evolve_ssprk35 order = 3 registers = 3 - cfl = 2.65062919143939d+00 * cfl name_int = "3ʳᵈ order 5-step SSPRK(5,3)" case ("ssprk(5,4)", "SSPRK(5,4)", "rk4.5", "RK4.5") @@ -247,7 +245,6 @@ module evolution evolve => evolve_ssprk45 order = 4 registers = 3 - cfl = 1.50818d+00 * cfl name_int = "4th order 5-step SSPRK(5,4)" case ("ssprk(10,4)", "SSPRK(10,4)", "rk4.10", "RK4.10") @@ -255,7 +252,6 @@ module evolution evolve => evolve_ssprk410 order = 4 registers = 3 - cfl = 6.0d+00 * cfl name_int = "4th order 10-step SSPRK(10,4)" case ("ssprk(m,2)", "SSPRK(m,2)", "ssprk2(1)m", "SSPRK2(1)M") @@ -265,7 +261,6 @@ module evolution order = 2 registers = 3 stages = max(2, min(9, stages)) - cfl = (stages - 1) * cfl write(name_int, "('2ⁿᵈ-order ',i0,'-step embedded " // & "SSPRK(',i0,',2) method')") stages, stages @@ -298,7 +293,6 @@ module evolution order = 3 registers = 3 stages = 4 - cfl = 2.0d+00 * cfl name_int = "3ʳᵈ-order 4-step embedded SSP3(2)4 method" betas(:) = [ 5.5d-01, -2.7d-01, 5.0d-02 ] @@ -321,7 +315,6 @@ module evolution end do n = n - 1 stages = max(4, n**2) - cfl = (n - 1) * n * cfl write(name_int, "('3ʳᵈ-order order ',i0,'-step embedded SSPRK3(2)n² method')") stages betas(:) = [ 5.5d-01, -2.7d-01, 5.0d-02 ]