diff --git a/src/rk.f90 b/src/rk.f90 index 9182eafe..cd89d284 100644 --- a/src/rk.f90 +++ b/src/rk.f90 @@ -259,7 +259,7 @@ subroutine rk_scal(rkpar,n,dli,l,dzci,dzfi,grid_vol_ratio_f,alpha,dt,is_bound,u, dsdtrko => dsdtrko_t end if if(is_cmpt_wallflux) then - call cmpt_wallflux(n,is_bound,l,dli,dzci,dzfi,alpha,flux) + call cmpt_scalflux(n,is_bound,l,dli,dzci,dzfi,alpha,s,flux) f = (factor1*sum(flux(:)/l(:)) + factor2*sum(fluxo(:)/l(:))) fluxo(:) = flux(:) end if