diff --git a/Source/Utils/TerrainMetrics.H b/Source/Utils/TerrainMetrics.H
index f52230a76..4f8e79ab3 100644
--- a/Source/Utils/TerrainMetrics.H
+++ b/Source/Utils/TerrainMetrics.H
@@ -344,17 +344,17 @@ amrex::Real
 Compute_Zrel_AtCellCenter (const int &i, const int &j, const int &k,
                            const amrex::Array4<const amrex::Real>& z_nd)
 {
-  const amrex::Real z_cc = 0.125*(  z_nd(i  ,j  ,k  ) + z_nd(i  ,j  ,k+1) +
-                                  + z_nd(i+1,j  ,k  ) + z_nd(i  ,j  ,k+1)
-                                  + z_nd(i  ,j+1,k  ) + z_nd(i  ,j+1,k+1)
-                                  + z_nd(i+1,j+1,k  ) + z_nd(i  ,j+1,k+1));
+    const amrex::Real z_cc = 0.125*(  z_nd(i  ,j  ,k  ) + z_nd(i  ,j  ,k+1) +
+                                      + z_nd(i+1,j  ,k  ) + z_nd(i  ,j  ,k+1)
+                                      + z_nd(i  ,j+1,k  ) + z_nd(i  ,j+1,k+1)
+                                      + z_nd(i+1,j+1,k  ) + z_nd(i  ,j+1,k+1));
 
-  // Note: we assume the z_nd array spans from the bottom to top of the domain
-  // i.e. no domain decomposition across processors in vertical direction
-  const amrex::Real z0_cc = 0.25*(  z_nd(i  ,j  ,0) + z_nd(i  ,j+1,0)
-                                  + z_nd(i+1,j  ,0) + z_nd(i+1,j+1,0));
+    // Note: we assume the z_nd array spans from the bottom to top of the domain
+    // i.e. no domain decomposition across processors in vertical direction
+    const amrex::Real z0_cc = 0.25*(  z_nd(i  ,j  ,0) + z_nd(i  ,j+1,0)
+                                      + z_nd(i+1,j  ,0) + z_nd(i+1,j+1,0));
 
-  return (z_cc - z0_cc);
+    return (z_cc - z0_cc);
 }
 
 /**