27 REAL(KIND=8),
DIMENSION(:),
ALLOCATABLE ::
buf_y1 28 REAL(KIND=8),
DIMENSION(:),
ALLOCATABLE ::
buf_f0 29 REAL(KIND=8),
DIMENSION(:),
ALLOCATABLE ::
buf_f1 39 IF (allocstat /= 0) stop
"*** Not enough memory ! ***" 49 REAL(KIND=8),
INTENT(IN) :: t
50 REAL(KIND=8),
DIMENSION(0:ndim),
INTENT(IN) :: y
51 REAL(KIND=8),
DIMENSION(0:ndim),
INTENT(OUT) :: res
60 SUBROUTINE step(y,t,dt,res)
61 REAL(KIND=8),
DIMENSION(0:ndim),
INTENT(IN) :: y
62 REAL(KIND=8),
INTENT(INOUT) :: t
63 REAL(KIND=8),
INTENT(IN) :: dt
64 REAL(KIND=8),
DIMENSION(0:ndim),
INTENT(OUT) :: res
integer ndim
Number of variables (dimension of the model)
subroutine, public step(y, t, dt, res)
Routine to perform an integration step (Heun algorithm). The incremented time is returned.
The equation tensor for the coupled ocean-atmosphere model with temperature which allows for an exten...
subroutine, public sparse_mul3(coolist_ijk, arr_j, arr_k, res)
Sparse multiplication of a tensor with two vectors: .
real(kind=8), dimension(:), allocatable buf_f0
Buffer to hold tendencies at the initial position.
Module with the integration routines.
subroutine tendencies(t, y, res)
Routine computing the tendencies of the model.
The model parameters module.
real(kind=8), dimension(:), allocatable buf_y1
Buffer to hold the intermediate position (Heun algorithm)
type(coolist), dimension(:), allocatable, public aotensor
- Tensor representation of the tendencies.
subroutine, public init_integrator
Routine to initialise the integration buffers.
real(kind=8), dimension(:), allocatable buf_f1
Buffer to hold tendencies at the intermediate position.