diff --git a/.cspell.json b/.cspell.json
index 4c1ef3d..83ffbff 100644
--- a/.cspell.json
+++ b/.cspell.json
@@ -36,6 +36,7 @@
     "pyproject.toml",
     "read_json.jl"
   ],
+  "ignoreWords": ["formfactor"],
   "language": "en-US",
   "languageSettings": [
     {
@@ -43,6 +44,5 @@
       "languageId": ["julia", "quarto", "text"]
     }
   ],
-  "version": "0.2",
-  "ignoreWords": ["formfactor", "pdfs"]
+  "version": "0.2"
 }
diff --git a/.cspell/julia.txt b/.cspell/julia.txt
index 4c59e8f..38a7c77 100644
--- a/.cspell/julia.txt
+++ b/.cspell/julia.txt
@@ -1,14 +1,17 @@
+colorbar
 compat
 eachslice
 guidefonthalign
 guidefontvalign
 joinpath
+latexstring
 lims
 minorticks
 mlims
 Mmap
 quadgk
 struct
+textrm
 xlab
 xlim
 ylab
diff --git a/.cspell/project.txt b/.cspell/project.txt
index 4f74cc9..fa1658e 100644
--- a/.cspell/project.txt
+++ b/.cspell/project.txt
@@ -1,8 +1,5 @@
-calv
-colorbar
-latexstring
 Mikhasenko
 Misha
 multline
-Semileptonic
-textrm
+pdfs
+semileptonic
diff --git a/docs/julia/lb2pkg.qmd b/docs/julia/lb2pkg.qmd
index f98dbd4..c9b6bde 100644
--- a/docs/julia/lb2pkg.qmd
+++ b/docs/julia/lb2pkg.qmd
@@ -50,7 +50,7 @@ end
 
 ## Deserialization of Objects to a Workspace
 
-Model components are deserialized from a JSON file into computational objects within a workspace for further manipulation. First, functions representing lineshapes and form factors are built. Following this, distributions are processed and added to the workspace.
+Model components are deserialized from a JSON file into computational objects within a workspace for further manipulation. First, functions representing lineshapes and form-factors are built. Following this, distributions are processed and added to the workspace.
 
 ```{julia}
 input = open(joinpath(@__DIR__, "..", "..", "models", "Lb2pKg.json")) do io
@@ -134,8 +134,8 @@ let
     _model = model
     f(σs) = Kibble(σs, ms^2) > 0 ? NaN : unpolarized_intensity(model, σs)
 
-    calv = [f(Invariants(ms; σ3, σ2)) for σ2 in σ2v, σ3 in σ3v]
-    heatmap(σ3v, σ2v, calv, colorbar=false, c=:speed)
+    values = [f(Invariants(ms; σ3, σ2)) for σ2 in σ2v, σ3 in σ3v]
+    heatmap(σ3v, σ2v, values, colorbar=false, c=:speed)
     plot!(border32(masses(model)),
         l=(3, :black),
         xlim=(σ3v[1], σ3v[end]), ylim=(σ2v[1], σ2v[end]),