MODULE array

  !USE mumps_bsplines
  use prec_const
  implicit none

  ! Arrays to store the rhs, for time integration
  COMPLEX(dp), DIMENSION(:,:,:,:,:), ALLOCATABLE :: moments_rhs_e ! (ip,ij,ikr,ikz,updatetlevel)
  COMPLEX(dp), DIMENSION(:,:,:,:,:), ALLOCATABLE :: moments_rhs_i ! (ip,ij,ikr,ikz,updatetlevel)

  ! To load collision matrix
  REAL(dp), DIMENSION(:,:), ALLOCATABLE :: iCe ! (ip,ij)
  REAL(dp), DIMENSION(:,:), ALLOCATABLE :: iCi ! (ip,ij)

  ! Intermediate steps in rhs of equations
  !COMPLEX(dp), DIMENSION(:,:,:), ALLOCATABLE:: moments_Apl, moments_Bpl, moments_Cpl, moments_Dpl,&
  !                                        moments_Epl, moments_Fpl, moments_Gpl, moments_Hpl  
END MODULE array