From c50c487b547c18de063e7187133addb501fad773 Mon Sep 17 00:00:00 2001 From: ignis Date: Fri, 30 Aug 2019 09:28:06 +0900 Subject: [PATCH] RK coefficients moved into RK4 subroutine --- code/ysolve.f90 | 25 ++++++++++++------------- 1 file changed, 12 insertions(+), 13 deletions(-) diff --git a/code/ysolve.f90 b/code/ysolve.f90 index c2dba97..6c6743e 100644 --- a/code/ysolve.f90 +++ b/code/ysolve.f90 @@ -9,18 +9,6 @@ REAL, DIMENSION(:), ALLOCATABLE :: uxt,duxt REAL, DIMENSION(:), ALLOCATABLE :: dm - INTEGER :: istage - REAL, DIMENSION(5) :: a=(/ 970286171893./4311952581923., & - 6584761158862./12103376702013., & - 2251764453980./15575788980749., & - 26877169314380./34165994151039., & - 0. /), & - b=(/ 1153189308089./22510343858157., & - 1772645290293./4653164025191., & - -1672844663538./4480602732383., & - 2114624349019./3568978502595., & - 5198255086312./14908931495163. /) - CONTAINS !------------------------------------------------------------------------ SUBROUTINE parse @@ -331,7 +319,18 @@ END SUBROUTINE rhs END INTERFACE - REAL :: test_var + REAL, DIMENSION(5), PARAMETER :: a=(/ 970286171893./4311952581923., & + 6584761158862./12103376702013., & + 2251764453980./15575788980749., & + 26877169314380./34165994151039., & + 0. /), & + b=(/ 1153189308089./22510343858157., & + 1772645290293./4653164025191., & + -1672844663538./4480602732383., & + 2114624349019./3568978502595., & + 5198255086312./14908931495163. /) + + INTEGER :: istage istage=1; CALL substep(y1,y1,y2,yf) istage=2; CALL substep(y1,y2,y1,yf)