--- a/src/SbpOperators/stencil.jl	Wed Dec 02 15:26:13 2020 +0100
+++ b/src/SbpOperators/stencil.jl	Wed Dec 02 17:10:18 2020 +0100
@@ -9,6 +9,18 @@
 end
 
 """
+    Stencil(weights::NTuple; center::Int)
+
+Create a stencil with the given weights with element `center` as the center of the stencil.
+"""
+function Stencil(weights::NTuple; center::Int)
+    N = length(weights)
+    range = (1, N) .- center
+
+    return Stencil(range, weights)
+end
+
+"""
     scale(s::Stencil, a)
 
 Scale the weights of the stencil `s` with `a` and return a new stencil.

--- a/test/testSbpOperators.jl	Wed Dec 02 15:26:13 2020 +0100
+++ b/test/testSbpOperators.jl	Wed Dec 02 17:10:18 2020 +0100
@@ -13,6 +13,10 @@
 
     @test eltype(s) == Float64
     @test SbpOperators.scale(s, 2) == SbpOperators.Stencil((-2,2), (2.,4.,4.,6.,8.))
+
+    @test SbpOperators.Stencil((1,2,3,4), center=1) == SbpOperators.Stencil((0, 3),(1,2,3,4))
+    @test SbpOperators.Stencil((1,2,3,4), center=2) == SbpOperators.Stencil((-1, 2),(1,2,3,4))
+    @test SbpOperators.Stencil((1,2,3,4), center=4) == SbpOperators.Stencil((-3, 0),(1,2,3,4))
 end
 
 # @testset "apply_quadrature" begin