module quads
  implicit none
  real(kind=8) :: pc_pi, trd,frd,erd,it1,it2
  real(kind=8), allocatable, dimension(:,:) :: vr,vph,vth,rho
  integer :: lmax
  parameter(pc_pi = acos(-1.d0))
  parameter(trd = 2.d0/3.d0)
  parameter(frd = 4.d0/3.d0)
  parameter(erd = 8.d0/3.d0)
  parameter(it1 = 1.0d0/sqrt(4.0d0*pc_pi))
  parameter(it2 = 1.0d0/(4.0d0*pc_pi))
  lmax = 20


contains

  subroutine set_vars(vel1,vel2,vel3,den)
    real(kind=8), intent(in) :: vel1,vel2,vel3,den
    integer :: i,ast
    ALLOCATE (vr(0:lmax,-lmax:lmax), STAT = ast)
    ALLOCATE (vph(0:lmax,-lmax:lmax), STAT = ast)
    ALLOCATE (vth(0:lmax,-lmax:lmax), STAT = ast)
    ALLOCATE (rho(0:lmax,-lmax:lmax), STAT = ast)

    vr = vel1 ; vph = vel2; vth = vel3; rho = den;
  end subroutine set_vars