diff --git a/src/forward.jl b/src/forward.jl
index e806051..77d15eb 100644
--- a/src/forward.jl
+++ b/src/forward.jl
@@ -394,7 +394,8 @@ function alifForward(   zit,
                         epsilonRec[i1,i2,i3,i4] = (alpha[i1,i2,i3,i4] * epsilonRec[i1,i2,i3,i4])
 
                         # compute epsilonRecA
-                        epsilonRecA[i1,i2,i3,i4] = (phi[i1,i2,i3,i4] * epsilonRec[i1,i2,i3,i4])
+                        epsilonRecA[i1,i2,i3,i4] = (rho[i1,i2,i3,i4] * 
+                                                        (phi[i1,i2,i3,i4] * epsilonRec[i1,i2,i3,i4]))
 
                         # compute avth
                         avth[i1,i2,i3,i4] = vth[i1,i2,i3,i4] + (beta[i1,i2,i3,i4] * a[i1,i2,i3,i4])