27 REAL(KIND=8),
DIMENSION(:),
ALLOCATABLE ::
buf_y1 28 REAL(KIND=8),
DIMENSION(:),
ALLOCATABLE ::
buf_ka 29 REAL(KIND=8),
DIMENSION(:),
ALLOCATABLE ::
buf_kb 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.
real(kind=8), dimension(:), allocatable buf_kb
Buffer B to hold tendencies.
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_ka
Buffer A to hold tendencies.
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.