diff --git a/Source/radiation/filter.H b/Source/radiation/filter.H
index be6a023490..d50a2a35cd 100644
--- a/Source/radiation/filter.H
+++ b/Source/radiation/filter.H
@@ -540,129 +540,129 @@ namespace filter
         }
         else if constexpr (width == 3) {
             if (Er_stencil(-1) == -1.0_rt) {
-                lamfil = ff3b0(0) * lam(0) +
-                         ff3b0(1) * lam(1) +
-                         ff3b0(2) * lam(2) +
-                         ff3b0(3) * lam(3);
+                lamfil = ff3b0(0) * lam_stencil(0) +
+                         ff3b0(1) * lam_stencil(1) +
+                         ff3b0(2) * lam_stencil(2) +
+                         ff3b0(3) * lam_stencil(3);
             }
             else if (Er_stencil(-2) == -1.0_rt) {
-                lamfil = ff3b1(-1) * lam(-1) +
-                         ff3b1( 0) * lam( 0) +
-                         ff3b1( 1) * lam( 1) +
-                         ff3b1( 2) * lam( 2) +
-                         ff3b1( 3) * lam( 3);
+                lamfil = ff3b1(-1) * lam_stencil(-1) +
+                         ff3b1( 0) * lam_stencil( 0) +
+                         ff3b1( 1) * lam_stencil( 1) +
+                         ff3b1( 2) * lam_stencil( 2) +
+                         ff3b1( 3) * lam_stencil( 3);
             }
             else if (Er_stencil(-3) == -1.0_rt) {
-                lamfil = ff3b2(-2) * lam(-2) +
-                         ff3b2(-1) * lam(-1) +
-                         ff3b2( 0) * lam( 0) +
-                         ff3b2( 1) * lam( 1) +
-                         ff3b2( 2) * lam( 2) +
-                         ff3b2( 3) * lam( 3);
+                lamfil = ff3b2(-2) * lam_stencil(-2) +
+                         ff3b2(-1) * lam_stencil(-1) +
+                         ff3b2( 0) * lam_stencil( 0) +
+                         ff3b2( 1) * lam_stencil( 1) +
+                         ff3b2( 2) * lam_stencil( 2) +
+                         ff3b2( 3) * lam_stencil( 3);
             }
             else if (Er_stencil(1) == -1.0_rt) {
-                lamfil = ff3b0(3) * lam(-3) +
-                         ff3b0(2) * lam(-2) +
-                         ff3b0(1) * lam(-1) +
-                         ff3b0(0) * lam( 0);
+                lamfil = ff3b0(3) * lam_stencil(-3) +
+                         ff3b0(2) * lam_stencil(-2) +
+                         ff3b0(1) * lam_stencil(-1) +
+                         ff3b0(0) * lam_stencil( 0);
             }
             else if (Er_stencil(2) == -1.0_rt) {
-                lamfil = ff3b1( 3) * lam(-3) +
-                         ff3b1( 2) * lam(-2) +
-                         ff3b1( 1) * lam(-1) +
-                         ff3b1( 0) * lam( 0) +
-                         ff3b1(-1) * lam( 1);
+                lamfil = ff3b1( 3) * lam_stencil(-3) +
+                         ff3b1( 2) * lam_stencil(-2) +
+                         ff3b1( 1) * lam_stencil(-1) +
+                         ff3b1( 0) * lam_stencil( 0) +
+                         ff3b1(-1) * lam_stencil( 1);
             }
             else if (Er_stencil(3) == -1.0_rt) {
-                lamfil = ff3b2( 3) * lam(-3) +
-                         ff3b2( 2) * lam(-2) +
-                         ff3b2( 1) * lam(-1) +
-                         ff3b2( 0) * lam( 0) +
-                         ff3b2(-1) * lam( 1) +
-                         ff3b2(-2) * lam( 2);
+                lamfil = ff3b2( 3) * lam_stencil(-3) +
+                         ff3b2( 2) * lam_stencil(-2) +
+                         ff3b2( 1) * lam_stencil(-1) +
+                         ff3b2( 0) * lam_stencil( 0) +
+                         ff3b2(-1) * lam_stencil( 1) +
+                         ff3b2(-2) * lam_stencil( 2);
             }
             else {
-                lamfil = ff3(0,S) * lam(0) +
-                         ff3(1,S) * (lam(-1) + lam(1)) +
-                         ff3(2,S) * (lam(-2) + lam(2)) +
-                         ff3(3,S) * (lam(-3) + lam(3));
+                lamfil = ff3(0,S) * lam_stencil(0) +
+                         ff3(1,S) * (lam_stencil(-1) + lam_stencil(1)) +
+                         ff3(2,S) * (lam_stencil(-2) + lam_stencil(2)) +
+                         ff3(3,S) * (lam_stencil(-3) + lam_stencil(3));
             }
         }
         else if constexpr (width == 4) {
             if (Er_stencil(-1) == -1.0_rt) {
-                lamfil = ff4b0(0) * lam(0) +
-                         ff4b0(1) * lam(1) +
-                         ff4b0(2) * lam(2) +
-                         ff4b0(3) * lam(3) +
-                         ff4b0(4) * lam(4);
+                lamfil = ff4b0(0) * lam_stencil(0) +
+                         ff4b0(1) * lam_stencil(1) +
+                         ff4b0(2) * lam_stencil(2) +
+                         ff4b0(3) * lam_stencil(3) +
+                         ff4b0(4) * lam_stencil(4);
             }
             else if (Er_stencil(-2) == -1.0_rt) {
-                lamfil = ff4b1(-1) * lam(-1) +
-                         ff4b1( 0) * lam( 0) +
-                         ff4b1( 1) * lam( 1) +
-                         ff4b1( 2) * lam( 2) +
-                         ff4b1( 3) * lam( 3) +
-                         ff4b1( 4) * lam( 4);
+                lamfil = ff4b1(-1) * lam_stencil(-1) +
+                         ff4b1( 0) * lam_stencil( 0) +
+                         ff4b1( 1) * lam_stencil( 1) +
+                         ff4b1( 2) * lam_stencil( 2) +
+                         ff4b1( 3) * lam_stencil( 3) +
+                         ff4b1( 4) * lam_stencil( 4);
             }
             else if (Er_stencil(-3) == -1.0_rt) {
-                lamfil = ff4b2(-2) * lam(-2) +
-                         ff4b2(-1) * lam(-1) +
-                         ff4b2( 0) * lam( 0) +
-                         ff4b2( 1) * lam( 1) +
-                         ff4b2( 2) * lam( 2) +
-                         ff4b2( 3) * lam( 3) +
-                         ff4b2( 4) * lam( 4);
+                lamfil = ff4b2(-2) * lam_stencil(-2) +
+                         ff4b2(-1) * lam_stencil(-1) +
+                         ff4b2( 0) * lam_stencil( 0) +
+                         ff4b2( 1) * lam_stencil( 1) +
+                         ff4b2( 2) * lam_stencil( 2) +
+                         ff4b2( 3) * lam_stencil( 3) +
+                         ff4b2( 4) * lam_stencil( 4);
             }
             else if (Er_stencil(-4) == -1.0_rt) {
-                lamfil = ff4b3(-3) * lam(-3) +
-                         ff4b3(-2) * lam(-2) +
-                         ff4b3(-1) * lam(-1) +
-                         ff4b3( 0) * lam( 0) +
-                         ff4b3( 1) * lam( 1) +
-                         ff4b3( 2) * lam( 2) +
-                         ff4b3( 3) * lam( 3) +
-                         ff4b3( 4) * lam( 4);
+                lamfil = ff4b3(-3) * lam_stencil(-3) +
+                         ff4b3(-2) * lam_stencil(-2) +
+                         ff4b3(-1) * lam_stencil(-1) +
+                         ff4b3( 0) * lam_stencil( 0) +
+                         ff4b3( 1) * lam_stencil( 1) +
+                         ff4b3( 2) * lam_stencil( 2) +
+                         ff4b3( 3) * lam_stencil( 3) +
+                         ff4b3( 4) * lam_stencil( 4);
             }
             else if (Er_stencil(1) == -1.0_rt) {
-                lamfil = ff4b0(4) * lam(-4) +
-                         ff4b0(3) * lam(-3) +
-                         ff4b0(2) * lam(-2) +
-                         ff4b0(1) * lam(-1) +
-                         ff4b0(0) * lam( 0);
+                lamfil = ff4b0(4) * lam_stencil(-4) +
+                         ff4b0(3) * lam_stencil(-3) +
+                         ff4b0(2) * lam_stencil(-2) +
+                         ff4b0(1) * lam_stencil(-1) +
+                         ff4b0(0) * lam_stencil( 0);
             }
             else if (Er_stencil(2) == -1.0_rt) {
-                lamfil = ff4b1( 4) * lam(-4) +
-                         ff4b1( 3) * lam(-3) +
-                         ff4b1( 2) * lam(-2) +
-                         ff4b1( 1) * lam(-1) +
-                         ff4b1( 0) * lam( 0) +
-                         ff4b1(-1) * lam( 1);
+                lamfil = ff4b1( 4) * lam_stencil(-4) +
+                         ff4b1( 3) * lam_stencil(-3) +
+                         ff4b1( 2) * lam_stencil(-2) +
+                         ff4b1( 1) * lam_stencil(-1) +
+                         ff4b1( 0) * lam_stencil( 0) +
+                         ff4b1(-1) * lam_stencil( 1);
             }
             else if (Er_stencil(3) == -1.0_rt) {
-                lamfil = ff4b2( 4) * lam(-4) +
-                         ff4b2( 3) * lam(-3) +
-                         ff4b2( 2) * lam(-2) +
-                         ff4b2( 1) * lam(-1) +
-                         ff4b2( 0) * lam( 0) +
-                         ff4b2(-1) * lam( 1) +
-                         ff4b2(-2) * lam( 2);
+                lamfil = ff4b2( 4) * lam_stencil(-4) +
+                         ff4b2( 3) * lam_stencil(-3) +
+                         ff4b2( 2) * lam_stencil(-2) +
+                         ff4b2( 1) * lam_stencil(-1) +
+                         ff4b2( 0) * lam_stencil( 0) +
+                         ff4b2(-1) * lam_stencil( 1) +
+                         ff4b2(-2) * lam_stencil( 2);
             }
             else if (Er_stencil(4) == -1.0_rt) {
-                lamfil = ff4b3( 4) * lam(-4) +
-                         ff4b3( 3) * lam(-3) +
-                         ff4b3( 2) * lam(-2) +
-                         ff4b3( 1) * lam(-1) +
-                         ff4b3( 0) * lam( 0) +
-                         ff4b3(-1) * lam( 1) +
-                         ff4b3(-2) * lam( 2) +
-                         ff4b3(-3) * lam( 3);
+                lamfil = ff4b3( 4) * lam_stencil(-4) +
+                         ff4b3( 3) * lam_stencil(-3) +
+                         ff4b3( 2) * lam_stencil(-2) +
+                         ff4b3( 1) * lam_stencil(-1) +
+                         ff4b3( 0) * lam_stencil( 0) +
+                         ff4b3(-1) * lam_stencil( 1) +
+                         ff4b3(-2) * lam_stencil( 2) +
+                         ff4b3(-3) * lam_stencil( 3);
             }
             else {
-                lamfil = ff4(0,S) * lam(0) +
-                         ff4(1,S) * (lam(-1) + lam(1)) +
-                         ff4(2,S) * (lam(-2) + lam(2)) +
-                         ff4(3,S) * (lam(-3) + lam(3)) +
-                         ff4(4,S) * (lam(-4) + lam(4));
+                lamfil = ff4(0,S) * lam_stencil(0) +
+                         ff4(1,S) * (lam_stencil(-1) + lam_stencil(1)) +
+                         ff4(2,S) * (lam_stencil(-2) + lam_stencil(2)) +
+                         ff4(3,S) * (lam_stencil(-3) + lam_stencil(3)) +
+                         ff4(4,S) * (lam_stencil(-4) + lam_stencil(4));
             }
         }