how to get the averaged value from the mutlilevel MultiFabs?

[ztdepztdep]

I want to get the total averaged value from the mutlilevel MultiFabs? Which function does this job?

[WeiqunZhang]

There is a function called volumeWeightedSum in AMReX_MultiFabUtil.H that probably does what you want. It returns a volume weighted sum of all "leaf" cells.

[ztdepztdep]

I got this error:"error: invalid initialization of reference of type ‘const amrex::Vector<const amrex::MultiFab*>&’ from expression of type ‘amrex::Vector<amrex::MultiFab, std::allocator >’ Real average=volumeWeightedSum(Rho,1,geom,nsubsteps);" I declar the first parameter as amrex::Vector Rho;