C -*- Mode: Fortran; -*-
c-----------------------------------------------------------------------
      module alphas
      contains
      subroutine SetAlphas(utilde,alpha,durl,evl,evr,ilo,ihi)
      use mesh
      use mesh_common
      implicit none

      double precision:: utilde(ILO:IHI,nvar)
      double precision:: alpha(ILO:IHI,nvar)
      double precision:: durl(ILO:IHI,nvar)
      double precision:: evl(ILO:IHI,nvar,nvar)
      double precision:: evr(ILO:IHI,nvar,nvar)
      integer:: i,j,k,l,m,ilo,ihi
c     
      do i=ilo,ihi,1
         alpha(i,:)=matmul(evl(i,:,:),durl(i,:))
      enddo
c
      return
      end subroutine SetAlphas
      end module alphas