Skip to content
Snippets Groups Projects
Commit d04b05b6 authored by Antoine Cyril David Hoffmann's avatar Antoine Cyril David Hoffmann :seedling:
Browse files

added explicit euler

parent cd2596f3
No related branches found
No related tags found
No related merge requests found
......@@ -70,6 +70,9 @@ CONTAINS
use parallel, ONLY: my_id
IMPLICIT NONE
SELECT CASE (numerical_scheme)
! Order 1 method
CASE ('EE')
CALL EE
! Order 2 methods
CASE ('RK2')
CALL RK2
......@@ -95,6 +98,22 @@ CONTAINS
IF (my_id .EQ. 0) WRITE(*,*) " Time integration with ", numerical_scheme
END SUBROUTINE set_numerical_scheme
!!! first order time schemes
SUBROUTINE EE
! Butcher coeff for Euler Explicit scheme
USE basic
USE prec_const
IMPLICIT NONE
INTEGER,PARAMETER :: nbstep = 1
CALL allocate_array(c_E,1,nbstep)
CALL allocate_array(b_E,1,nbstep,1,1)
CALL allocate_array(A_E,1,nbstep,1,nbstep)
ntimelevel = 1
c_E(1) = 0._xp
b_E(1,1) = 1._xp
A_E(1,1) = 0._xp
END SUBROUTINE EE
!!! second order time schemes
SUBROUTINE RK2
! Butcher coeff for clasical RK2 (Heun's)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment