diff --git a/docs/models/PDE/wave.md b/docs/models/PDE/wave.md
index 80cd8692e..f01cc2df8 100644
--- a/docs/models/PDE/wave.md
+++ b/docs/models/PDE/wave.md
@@ -5,8 +5,8 @@ wave
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/wave/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/wave/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/wave/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/wave/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/PDE/wave2d.md b/docs/models/PDE/wave2d.md
index 410ba69b6..2b38ae225 100644
--- a/docs/models/PDE/wave2d.md
+++ b/docs/models/PDE/wave2d.md
@@ -5,8 +5,8 @@ wave2d
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/wave2d/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/wave2d/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/wave2d/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/wave2d/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/article/d3q19_heat_adj_art.md b/docs/models/article/d3q19_heat_adj_art.md
index 346bfc69b..cb3f59617 100644
--- a/docs/models/article/d3q19_heat_adj_art.md
+++ b/docs/models/article/d3q19_heat_adj_art.md
@@ -5,8 +5,8 @@ d3q19_heat_adj_art
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q19_heat_adj_art/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q19_heat_adj_art/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q19_heat_adj_art/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q19_heat_adj_art/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/electrokinetic/d2q9_npe_guo.md b/docs/models/electrokinetic/d2q9_npe_guo.md
index 082b8fe96..038df2c0d 100644
--- a/docs/models/electrokinetic/d2q9_npe_guo.md
+++ b/docs/models/electrokinetic/d2q9_npe_guo.md
@@ -5,8 +5,8 @@ d2q9_npe_guo
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_npe_guo/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_npe_guo/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_npe_guo/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_npe_guo/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/electrokinetic/d2q9_poison_boltzmann.md b/docs/models/electrokinetic/d2q9_poison_boltzmann.md
index 52fc4eb2c..7025ce42e 100644
--- a/docs/models/electrokinetic/d2q9_poison_boltzmann.md
+++ b/docs/models/electrokinetic/d2q9_poison_boltzmann.md
@@ -5,8 +5,8 @@ d2q9_poison_boltzmann
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_poison_boltzmann/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_poison_boltzmann/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_poison_boltzmann/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_poison_boltzmann/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/experimental/d2q9_diff.md b/docs/models/experimental/d2q9_diff.md
index 9ac32fbc8..9e9eb29fa 100644
--- a/docs/models/experimental/d2q9_diff.md
+++ b/docs/models/experimental/d2q9_diff.md
@@ -5,8 +5,8 @@ d2q9_diff
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_diff/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_diff/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_diff/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_diff/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/experimental/d2q9_hb.md b/docs/models/experimental/d2q9_hb.md
index 8735039b8..3ce29a739 100644
--- a/docs/models/experimental/d2q9_hb.md
+++ b/docs/models/experimental/d2q9_hb.md
@@ -5,8 +5,8 @@ d2q9_hb
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_hb/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_hb/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_hb/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_hb/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/experimental/d2q9_inc.md b/docs/models/experimental/d2q9_inc.md
index 23d59c452..264358cae 100644
--- a/docs/models/experimental/d2q9_inc.md
+++ b/docs/models/experimental/d2q9_inc.md
@@ -5,8 +5,8 @@ d2q9_inc
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_inc/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_inc/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_inc/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_inc/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/flow/d2q9/d2q9.md b/docs/models/flow/d2q9/d2q9.md
index eeba83622..23c1bdc1d 100644
--- a/docs/models/flow/d2q9/d2q9.md
+++ b/docs/models/flow/d2q9/d2q9.md
@@ -5,8 +5,8 @@ The `d2q9` model is a basic 2D Lattice Boltzmann Method model for flow simulatio
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/flow/d2q9/d2q9_BC.md b/docs/models/flow/d2q9/d2q9_BC.md
index 19bc23d2f..c162fb727 100644
--- a/docs/models/flow/d2q9/d2q9_BC.md
+++ b/docs/models/flow/d2q9/d2q9_BC.md
@@ -5,8 +5,8 @@ The `d2q9` model is a basic 2D Lattice Boltzmann Method model for flow simulatio
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_BC/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_BC/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_BC/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_BC/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/flow/d2q9/d2q9_BC_autosym.md b/docs/models/flow/d2q9/d2q9_BC_autosym.md
index 18b01da08..60f6dac2b 100644
--- a/docs/models/flow/d2q9/d2q9_BC_autosym.md
+++ b/docs/models/flow/d2q9/d2q9_BC_autosym.md
@@ -5,8 +5,8 @@ The `d2q9` model is a basic 2D Lattice Boltzmann Method model for flow simulatio
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_BC_autosym/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_BC_autosym/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_BC_autosym/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_BC_autosym/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/flow/d2q9/d2q9_autosym.md b/docs/models/flow/d2q9/d2q9_autosym.md
index fe2a05c51..dd273d075 100644
--- a/docs/models/flow/d2q9/d2q9_autosym.md
+++ b/docs/models/flow/d2q9/d2q9_autosym.md
@@ -5,8 +5,8 @@ The `d2q9` model is a basic 2D Lattice Boltzmann Method model for flow simulatio
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_autosym/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_autosym/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_autosym/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_autosym/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/flow/d2q9_SRT.md b/docs/models/flow/d2q9_SRT.md
index 096d3e465..90a06a9a5 100644
--- a/docs/models/flow/d2q9_SRT.md
+++ b/docs/models/flow/d2q9_SRT.md
@@ -5,8 +5,8 @@ d2q9_SRT
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_SRT/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_SRT/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_SRT/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_SRT/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/flow/d2q9_cumulant.md b/docs/models/flow/d2q9_cumulant.md
index 59760d1ca..efc488246 100644
--- a/docs/models/flow/d2q9_cumulant.md
+++ b/docs/models/flow/d2q9_cumulant.md
@@ -5,8 +5,8 @@ d2q9_cumulant
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_cumulant/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_cumulant/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_cumulant/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_cumulant/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/flow/d2q9_les.md b/docs/models/flow/d2q9_les.md
index fd3f53e3e..112a8c5b2 100644
--- a/docs/models/flow/d2q9_les.md
+++ b/docs/models/flow/d2q9_les.md
@@ -5,8 +5,8 @@ d2q9_les
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_les/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_les/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_les/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_les/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/flow/d2q9_thin_film.md b/docs/models/flow/d2q9_thin_film.md
index ead457667..ace469247 100644
--- a/docs/models/flow/d2q9_thin_film.md
+++ b/docs/models/flow/d2q9_thin_film.md
@@ -5,8 +5,8 @@ d2q9_thin_film
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_thin_film/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_thin_film/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_thin_film/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_thin_film/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/flow/d3q19.md b/docs/models/flow/d3q19.md
index e00cfc332..f322423a1 100644
--- a/docs/models/flow/d3q19.md
+++ b/docs/models/flow/d3q19.md
@@ -5,8 +5,8 @@ d3q19
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q19/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q19/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q19/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q19/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/flow/d3q19_les.md b/docs/models/flow/d3q19_les.md
index 596bdeb1f..8a6791747 100644
--- a/docs/models/flow/d3q19_les.md
+++ b/docs/models/flow/d3q19_les.md
@@ -5,8 +5,8 @@ d3q19_les
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q19_les/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q19_les/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q19_les/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q19_les/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/flow/d3q27.md b/docs/models/flow/d3q27.md
index 60d7f41f5..02b297325 100644
--- a/docs/models/flow/d3q27.md
+++ b/docs/models/flow/d3q27.md
@@ -5,8 +5,8 @@ d3q27
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/flow/d3q27_cumulant/d3q27_cumulant.md b/docs/models/flow/d3q27_cumulant/d3q27_cumulant.md
index 3b9fa21fb..5eb9cc0ee 100644
--- a/docs/models/flow/d3q27_cumulant/d3q27_cumulant.md
+++ b/docs/models/flow/d3q27_cumulant/d3q27_cumulant.md
@@ -5,8 +5,8 @@ d3q27_cumulant
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_cumulant/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_cumulant/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_cumulant/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_cumulant/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/flow/d3q27_cumulant/d3q27_cumulant_AVG.md b/docs/models/flow/d3q27_cumulant/d3q27_cumulant_AVG.md
index 11b54bf88..9cc752d05 100644
--- a/docs/models/flow/d3q27_cumulant/d3q27_cumulant_AVG.md
+++ b/docs/models/flow/d3q27_cumulant/d3q27_cumulant_AVG.md
@@ -5,8 +5,8 @@ d3q27_cumulant_AVG
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_cumulant_AVG/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_cumulant_AVG/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_cumulant_AVG/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_cumulant_AVG/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/flow/d3q27_cumulant/d3q27_cumulant_AVG_IB.md b/docs/models/flow/d3q27_cumulant/d3q27_cumulant_AVG_IB.md
index 76f1ce4c0..e2fc36398 100644
--- a/docs/models/flow/d3q27_cumulant/d3q27_cumulant_AVG_IB.md
+++ b/docs/models/flow/d3q27_cumulant/d3q27_cumulant_AVG_IB.md
@@ -5,8 +5,8 @@ d3q27_cumulant_AVG_IB
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_cumulant_AVG_IB/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_cumulant_AVG_IB/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_cumulant_AVG_IB/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_cumulant_AVG_IB/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/flow/d3q27_cumulant/d3q27_cumulant_AVG_IB_SMAG.md b/docs/models/flow/d3q27_cumulant/d3q27_cumulant_AVG_IB_SMAG.md
index 83f2396ac..6a3f8d94f 100644
--- a/docs/models/flow/d3q27_cumulant/d3q27_cumulant_AVG_IB_SMAG.md
+++ b/docs/models/flow/d3q27_cumulant/d3q27_cumulant_AVG_IB_SMAG.md
@@ -5,8 +5,8 @@ d3q27_cumulant_AVG_IB_SMAG
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_cumulant_AVG_IB_SMAG/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_cumulant_AVG_IB_SMAG/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_cumulant_AVG_IB_SMAG/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_cumulant_AVG_IB_SMAG/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/flow/d3q27_cumulant/d3q27_cumulant_AVG_SMAG.md b/docs/models/flow/d3q27_cumulant/d3q27_cumulant_AVG_SMAG.md
index 78787a9a5..1b16aab1f 100644
--- a/docs/models/flow/d3q27_cumulant/d3q27_cumulant_AVG_SMAG.md
+++ b/docs/models/flow/d3q27_cumulant/d3q27_cumulant_AVG_SMAG.md
@@ -5,8 +5,8 @@ d3q27_cumulant_AVG_SMAG
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_cumulant_AVG_SMAG/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_cumulant_AVG_SMAG/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_cumulant_AVG_SMAG/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_cumulant_AVG_SMAG/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/flow/d3q27_cumulant/d3q27_cumulant_IB.md b/docs/models/flow/d3q27_cumulant/d3q27_cumulant_IB.md
index eaab7cb49..511155c9b 100644
--- a/docs/models/flow/d3q27_cumulant/d3q27_cumulant_IB.md
+++ b/docs/models/flow/d3q27_cumulant/d3q27_cumulant_IB.md
@@ -5,8 +5,8 @@ d3q27_cumulant_IB
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_cumulant_IB/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_cumulant_IB/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_cumulant_IB/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_cumulant_IB/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/flow/d3q27_cumulant/d3q27_cumulant_IB_SMAG.md b/docs/models/flow/d3q27_cumulant/d3q27_cumulant_IB_SMAG.md
index 9c465847a..47f4f74ea 100644
--- a/docs/models/flow/d3q27_cumulant/d3q27_cumulant_IB_SMAG.md
+++ b/docs/models/flow/d3q27_cumulant/d3q27_cumulant_IB_SMAG.md
@@ -5,8 +5,8 @@ d3q27_cumulant_IB_SMAG
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_cumulant_IB_SMAG/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_cumulant_IB_SMAG/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_cumulant_IB_SMAG/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_cumulant_IB_SMAG/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/flow/d3q27_cumulant/d3q27_cumulant_SMAG.md b/docs/models/flow/d3q27_cumulant/d3q27_cumulant_SMAG.md
index f7831797e..202c91068 100644
--- a/docs/models/flow/d3q27_cumulant/d3q27_cumulant_SMAG.md
+++ b/docs/models/flow/d3q27_cumulant/d3q27_cumulant_SMAG.md
@@ -5,8 +5,8 @@ d3q27_cumulant_SMAG
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_cumulant_SMAG/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_cumulant_SMAG/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_cumulant_SMAG/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_cumulant_SMAG/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/flow/experimental/d3q27_BGK.md b/docs/models/flow/experimental/d3q27_BGK.md
index 3a34fb174..11c547d3e 100644
--- a/docs/models/flow/experimental/d3q27_BGK.md
+++ b/docs/models/flow/experimental/d3q27_BGK.md
@@ -5,8 +5,8 @@ d3q27_BGK
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_BGK/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_BGK/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_BGK/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_BGK/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/flow/experimental/d3q27_BGK_galcor.md b/docs/models/flow/experimental/d3q27_BGK_galcor.md
index cf0873d55..7b27a4061 100644
--- a/docs/models/flow/experimental/d3q27_BGK_galcor.md
+++ b/docs/models/flow/experimental/d3q27_BGK_galcor.md
@@ -5,8 +5,8 @@ d3q27_BGK_galcor
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_BGK_galcor/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_BGK_galcor/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_BGK_galcor/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_BGK_galcor/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/flow/qibb/d3q27_cumulant_qibb_small.md b/docs/models/flow/qibb/d3q27_cumulant_qibb_small.md
index 895d351ec..c569ff64a 100644
--- a/docs/models/flow/qibb/d3q27_cumulant_qibb_small.md
+++ b/docs/models/flow/qibb/d3q27_cumulant_qibb_small.md
@@ -5,8 +5,8 @@ d3q27_cumulant_qibb_small
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_cumulant_qibb_small/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_cumulant_qibb_small/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_cumulant_qibb_small/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_cumulant_qibb_small/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/heat/d3q19_heat.md b/docs/models/heat/d3q19_heat.md
index 58f097581..e927851db 100644
--- a/docs/models/heat/d3q19_heat.md
+++ b/docs/models/heat/d3q19_heat.md
@@ -5,8 +5,8 @@ d3q19_heat
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q19_heat/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q19_heat/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q19_heat/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q19_heat/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/heat/experimental/d2q9_heat.md b/docs/models/heat/experimental/d2q9_heat.md
index affcd940a..be3370c7b 100644
--- a/docs/models/heat/experimental/d2q9_heat.md
+++ b/docs/models/heat/experimental/d2q9_heat.md
@@ -5,8 +5,8 @@ d2q9_heat
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_heat/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_heat/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_heat/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_heat/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/moving/d2q9_plate.md b/docs/models/moving/d2q9_plate.md
index bf7f4e550..874007a1b 100644
--- a/docs/models/moving/d2q9_plate.md
+++ b/docs/models/moving/d2q9_plate.md
@@ -5,8 +5,8 @@ d2q9_plate
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_plate/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_plate/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_plate/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_plate/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/multiphase/d2q9_csf.md b/docs/models/multiphase/d2q9_csf.md
index 9f9cad14d..194142f96 100644
--- a/docs/models/multiphase/d2q9_csf.md
+++ b/docs/models/multiphase/d2q9_csf.md
@@ -5,8 +5,8 @@ d2q9_csf
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_csf/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_csf/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_csf/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_csf/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/multiphase/d2q9_kuper.md b/docs/models/multiphase/d2q9_kuper.md
index f942be1e8..0ea220cc5 100644
--- a/docs/models/multiphase/d2q9_kuper.md
+++ b/docs/models/multiphase/d2q9_kuper.md
@@ -5,8 +5,8 @@ d2q9_kuper
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_kuper/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_kuper/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_kuper/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_kuper/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/multiphase/d2q9_pf.md b/docs/models/multiphase/d2q9_pf.md
index 5c4602dda..054c831af 100644
--- a/docs/models/multiphase/d2q9_pf.md
+++ b/docs/models/multiphase/d2q9_pf.md
@@ -5,8 +5,8 @@ d2q9_pf
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_pf/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_pf/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_pf/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_pf/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/multiphase/d2q9_pf_pressureEvolution.md b/docs/models/multiphase/d2q9_pf_pressureEvolution.md
index 06d939910..35a46eea8 100644
--- a/docs/models/multiphase/d2q9_pf_pressureEvolution.md
+++ b/docs/models/multiphase/d2q9_pf_pressureEvolution.md
@@ -5,8 +5,8 @@ d2q9_pf_pressureEvolution
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_pf_pressureEvolution/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_pf_pressureEvolution/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_pf_pressureEvolution/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_pf_pressureEvolution/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/multiphase/d2q9_pf_velocity/d2q9_pf_velocity.md b/docs/models/multiphase/d2q9_pf_velocity/d2q9_pf_velocity.md
index 973ec6901..1a1994fd8 100644
--- a/docs/models/multiphase/d2q9_pf_velocity/d2q9_pf_velocity.md
+++ b/docs/models/multiphase/d2q9_pf_velocity/d2q9_pf_velocity.md
@@ -5,8 +5,8 @@ The 'd2q9_pf_velocity' model is a multiphase 2D(planar) lattice Boltzmann model
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_pf_velocity/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_pf_velocity/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_pf_velocity/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_pf_velocity/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/multiphase/d2q9_pf_velocity/d2q9_pf_velocity_GF.md b/docs/models/multiphase/d2q9_pf_velocity/d2q9_pf_velocity_GF.md
index 4f8749425..455439ebc 100644
--- a/docs/models/multiphase/d2q9_pf_velocity/d2q9_pf_velocity_GF.md
+++ b/docs/models/multiphase/d2q9_pf_velocity/d2q9_pf_velocity_GF.md
@@ -5,8 +5,8 @@ The 'd2q9_pf_velocity' model is a multiphase 2D(planar) lattice Boltzmann model
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_pf_velocity_GF/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_pf_velocity_GF/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_pf_velocity_GF/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_pf_velocity_GF/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/multiphase/d2q9_pf_velocity/d2q9_pf_velocity_GF_autosym.md b/docs/models/multiphase/d2q9_pf_velocity/d2q9_pf_velocity_GF_autosym.md
index 9918157db..fcf58ed17 100644
--- a/docs/models/multiphase/d2q9_pf_velocity/d2q9_pf_velocity_GF_autosym.md
+++ b/docs/models/multiphase/d2q9_pf_velocity/d2q9_pf_velocity_GF_autosym.md
@@ -5,8 +5,8 @@ The 'd2q9_pf_velocity' model is a multiphase 2D(planar) lattice Boltzmann model
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_pf_velocity_GF_autosym/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_pf_velocity_GF_autosym/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_pf_velocity_GF_autosym/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_pf_velocity_GF_autosym/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/multiphase/d2q9_pf_velocity/d2q9_pf_velocity_Outflow.md b/docs/models/multiphase/d2q9_pf_velocity/d2q9_pf_velocity_Outflow.md
index bfc63c44e..f322f4031 100644
--- a/docs/models/multiphase/d2q9_pf_velocity/d2q9_pf_velocity_Outflow.md
+++ b/docs/models/multiphase/d2q9_pf_velocity/d2q9_pf_velocity_Outflow.md
@@ -5,8 +5,8 @@ The 'd2q9_pf_velocity' model is a multiphase 2D(planar) lattice Boltzmann model
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_pf_velocity_Outflow/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_pf_velocity_Outflow/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_pf_velocity_Outflow/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_pf_velocity_Outflow/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/multiphase/d2q9_pf_velocity/d2q9_pf_velocity_Outflow_autosym.md b/docs/models/multiphase/d2q9_pf_velocity/d2q9_pf_velocity_Outflow_autosym.md
index a8ee001c3..bdfe3a3af 100644
--- a/docs/models/multiphase/d2q9_pf_velocity/d2q9_pf_velocity_Outflow_autosym.md
+++ b/docs/models/multiphase/d2q9_pf_velocity/d2q9_pf_velocity_Outflow_autosym.md
@@ -5,8 +5,8 @@ The 'd2q9_pf_velocity' model is a multiphase 2D(planar) lattice Boltzmann model
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_pf_velocity_Outflow_autosym/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_pf_velocity_Outflow_autosym/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_pf_velocity_Outflow_autosym/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_pf_velocity_Outflow_autosym/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/multiphase/d2q9_pf_velocity/d2q9_pf_velocity_autosym.md b/docs/models/multiphase/d2q9_pf_velocity/d2q9_pf_velocity_autosym.md
index cbd1676bd..2724e1122 100644
--- a/docs/models/multiphase/d2q9_pf_velocity/d2q9_pf_velocity_autosym.md
+++ b/docs/models/multiphase/d2q9_pf_velocity/d2q9_pf_velocity_autosym.md
@@ -5,8 +5,8 @@ The 'd2q9_pf_velocity' model is a multiphase 2D(planar) lattice Boltzmann model
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_pf_velocity_autosym/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_pf_velocity_autosym/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_pf_velocity_autosym/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_pf_velocity_autosym/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/multiphase/d3q19_kuper.md b/docs/models/multiphase/d3q19_kuper.md
index 0d86ac554..fe5b4de41 100644
--- a/docs/models/multiphase/d3q19_kuper.md
+++ b/docs/models/multiphase/d3q19_kuper.md
@@ -5,8 +5,8 @@ d3q19_kuper
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q19_kuper/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q19_kuper/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q19_kuper/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q19_kuper/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/multiphase/d3q27_pf_velocity/d3q27_pf_velocity.md b/docs/models/multiphase/d3q27_pf_velocity/d3q27_pf_velocity.md
index 137e38541..9a8e5283a 100644
--- a/docs/models/multiphase/d3q27_pf_velocity/d3q27_pf_velocity.md
+++ b/docs/models/multiphase/d3q27_pf_velocity/d3q27_pf_velocity.md
@@ -5,8 +5,8 @@ The 'd3q27_pf_velocity' model is a multiphase 3D lattice Boltzmann model for the
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_pf_velocity/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_pf_velocity/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_pf_velocity/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_pf_velocity/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/multiphase/d3q27_pf_velocity/d3q27_pf_velocity_BGK.md b/docs/models/multiphase/d3q27_pf_velocity/d3q27_pf_velocity_BGK.md
index 2ce89dc1a..de943772d 100644
--- a/docs/models/multiphase/d3q27_pf_velocity/d3q27_pf_velocity_BGK.md
+++ b/docs/models/multiphase/d3q27_pf_velocity/d3q27_pf_velocity_BGK.md
@@ -5,8 +5,8 @@ The 'd3q27_pf_velocity' model is a multiphase 3D lattice Boltzmann model for the
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_pf_velocity_BGK/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_pf_velocity_BGK/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_pf_velocity_BGK/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_pf_velocity_BGK/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/multiphase/d3q27_pf_velocity/d3q27_pf_velocity_BGK_autosym.md b/docs/models/multiphase/d3q27_pf_velocity/d3q27_pf_velocity_BGK_autosym.md
index 996c7e9f7..12529e0ab 100644
--- a/docs/models/multiphase/d3q27_pf_velocity/d3q27_pf_velocity_BGK_autosym.md
+++ b/docs/models/multiphase/d3q27_pf_velocity/d3q27_pf_velocity_BGK_autosym.md
@@ -5,8 +5,8 @@ The 'd3q27_pf_velocity' model is a multiphase 3D lattice Boltzmann model for the
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_pf_velocity_BGK_autosym/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_pf_velocity_BGK_autosym/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_pf_velocity_BGK_autosym/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_pf_velocity_BGK_autosym/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/multiphase/d3q27_pf_velocity/d3q27_pf_velocity_OutFlow.md b/docs/models/multiphase/d3q27_pf_velocity/d3q27_pf_velocity_OutFlow.md
index ea063a72d..0e85213c3 100644
--- a/docs/models/multiphase/d3q27_pf_velocity/d3q27_pf_velocity_OutFlow.md
+++ b/docs/models/multiphase/d3q27_pf_velocity/d3q27_pf_velocity_OutFlow.md
@@ -5,8 +5,8 @@ The 'd3q27_pf_velocity' model is a multiphase 3D lattice Boltzmann model for the
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_pf_velocity_OutFlow/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_pf_velocity_OutFlow/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_pf_velocity_OutFlow/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_pf_velocity_OutFlow/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/multiphase/d3q27_pf_velocity/d3q27_pf_velocity_OutFlow_autosym.md b/docs/models/multiphase/d3q27_pf_velocity/d3q27_pf_velocity_OutFlow_autosym.md
index dc6d8f248..c94b0883d 100644
--- a/docs/models/multiphase/d3q27_pf_velocity/d3q27_pf_velocity_OutFlow_autosym.md
+++ b/docs/models/multiphase/d3q27_pf_velocity/d3q27_pf_velocity_OutFlow_autosym.md
@@ -5,8 +5,8 @@ The 'd3q27_pf_velocity' model is a multiphase 3D lattice Boltzmann model for the
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_pf_velocity_OutFlow_autosym/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_pf_velocity_OutFlow_autosym/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_pf_velocity_OutFlow_autosym/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_pf_velocity_OutFlow_autosym/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/multiphase/d3q27_pf_velocity/d3q27_pf_velocity_autosym.md b/docs/models/multiphase/d3q27_pf_velocity/d3q27_pf_velocity_autosym.md
index 4d4588d84..2cb314968 100644
--- a/docs/models/multiphase/d3q27_pf_velocity/d3q27_pf_velocity_autosym.md
+++ b/docs/models/multiphase/d3q27_pf_velocity/d3q27_pf_velocity_autosym.md
@@ -5,8 +5,8 @@ The 'd3q27_pf_velocity' model is a multiphase 3D lattice Boltzmann model for the
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_pf_velocity_autosym/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_pf_velocity_autosym/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_pf_velocity_autosym/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_pf_velocity_autosym/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/multiphase/experimental/d2q9_lee.md b/docs/models/multiphase/experimental/d2q9_lee.md
index e2f7aa208..813e035e0 100644
--- a/docs/models/multiphase/experimental/d2q9_lee.md
+++ b/docs/models/multiphase/experimental/d2q9_lee.md
@@ -5,8 +5,8 @@ d2q9_lee
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_lee/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_lee/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_lee/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_lee/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/multiphase/experimental/d2q9_pp_LBL.md b/docs/models/multiphase/experimental/d2q9_pp_LBL.md
index 855f6633b..355b2bddb 100644
--- a/docs/models/multiphase/experimental/d2q9_pp_LBL.md
+++ b/docs/models/multiphase/experimental/d2q9_pp_LBL.md
@@ -5,8 +5,8 @@ d2q9_pp_LBL
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_pp_LBL/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_pp_LBL/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_pp_LBL/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_pp_LBL/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/multiphase/experimental/d2q9_pp_MCMP.md b/docs/models/multiphase/experimental/d2q9_pp_MCMP.md
index 0dfea5f8d..e4ae59b6f 100644
--- a/docs/models/multiphase/experimental/d2q9_pp_MCMP.md
+++ b/docs/models/multiphase/experimental/d2q9_pp_MCMP.md
@@ -5,8 +5,8 @@ d2q9_pp_MCMP
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_pp_MCMP/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_pp_MCMP/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_pp_MCMP/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_pp_MCMP/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/multiphase/solidification/d2q9_solid.md b/docs/models/multiphase/solidification/d2q9_solid.md
index 600da8afa..31f6349f9 100644
--- a/docs/models/multiphase/solidification/d2q9_solid.md
+++ b/docs/models/multiphase/solidification/d2q9_solid.md
@@ -5,8 +5,8 @@ d2q9_solid
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_solid/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_solid/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_solid/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_solid/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/nonnewtonian/d3q27_viscoplastic.md b/docs/models/nonnewtonian/d3q27_viscoplastic.md
index 3a48c9741..eba7f00e5 100644
--- a/docs/models/nonnewtonian/d3q27_viscoplastic.md
+++ b/docs/models/nonnewtonian/d3q27_viscoplastic.md
@@ -5,8 +5,8 @@ d3q27_viscoplastic
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_viscoplastic/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q27_viscoplastic/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_viscoplastic/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q27_viscoplastic/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/optimization/d2q9_adj.md b/docs/models/optimization/d2q9_adj.md
index 87126ce84..f920671ba 100644
--- a/docs/models/optimization/d2q9_adj.md
+++ b/docs/models/optimization/d2q9_adj.md
@@ -5,8 +5,8 @@ d2q9_adj
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_adj/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_adj/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_adj/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_adj/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/optimization/d2q9_optimalMixing.md b/docs/models/optimization/d2q9_optimalMixing.md
index c3d7bfe51..18ac3cd67 100644
--- a/docs/models/optimization/d2q9_optimalMixing.md
+++ b/docs/models/optimization/d2q9_optimalMixing.md
@@ -5,8 +5,8 @@ d2q9_optimalMixing
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_optimalMixing/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_optimalMixing/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_optimalMixing/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_optimalMixing/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/optimization/d3q19_heat_adj.md b/docs/models/optimization/d3q19_heat_adj.md
index 83f366907..e9ed7833f 100644
--- a/docs/models/optimization/d3q19_heat_adj.md
+++ b/docs/models/optimization/d3q19_heat_adj.md
@@ -5,8 +5,8 @@ d3q19_heat_adj
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q19_heat_adj/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d3q19_heat_adj/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q19_heat_adj/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d3q19_heat_adj/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/optimization/experimental/d2q9_heat_adj.md b/docs/models/optimization/experimental/d2q9_heat_adj.md
index b6935d7fb..9cb6b05ea 100644
--- a/docs/models/optimization/experimental/d2q9_heat_adj.md
+++ b/docs/models/optimization/experimental/d2q9_heat_adj.md
@@ -5,8 +5,8 @@ d2q9_heat_adj
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_heat_adj/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/d2q9_heat_adj/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_heat_adj/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/d2q9_heat_adj/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/models/shallowwater/sw.md b/docs/models/shallowwater/sw.md
index a44ad112a..78b707684 100644
--- a/docs/models/shallowwater/sw.md
+++ b/docs/models/shallowwater/sw.md
@@ -10,8 +10,8 @@ Model has adjoint capabilities for unsteady optimization.
 
 ## Details
 [Model description files](Model description) files for this model:
-[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/sw/Dynamics.c.Rt)
-[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 7714130)/src/sw/Dynamics.R)
+[Dynamics.c](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/sw/Dynamics.c.Rt)
+[Dynamics.R](https://github.com/llaniewski/TCLB/blob/(HEAD detached at 21557e6)/src/sw/Dynamics.R)
 
 ### [Zonal Settings](Settings)
 
diff --git a/docs/tutorials/model-development/4.-D2Q9-HeatTransfer.md b/docs/tutorials/model-development/4.-D2Q9-HeatTransfer.md
index d438f7b73..b35c610a2 100644
--- a/docs/tutorials/model-development/4.-D2Q9-HeatTransfer.md
+++ b/docs/tutorials/model-development/4.-D2Q9-HeatTransfer.md
@@ -176,7 +176,7 @@ AddQuantity(name="U", vector=T)
 
 ### Settings
 
-Here, we define simulation settings (like properties of fluid). 
+Here, we define simulation settings (like properties of fluid).
 Once the simulation is launched, the solver will look for them in the `batch.xml` file.
 
 ```R
@@ -204,9 +204,8 @@ AddSetting(name="BoussinesqCoeff", default=0.0, comment='BoussinesqCoeff=rho_0*t
 
 ```R
 #	Boundary things
-AddNodeType(name="Solid", group="ADDITIONALS")
-AddNodeType(name="HeaterDirichletTemperature", group="ADDITIONALS")
-AddNodeType(name="HeaterNeumannHeatFlux", group="ADDITIONALS")
+AddNodeType(name="HeaterDirichletTemperature", group="ADDITIONALS_HEAT")
+AddNodeType(name="HeaterNeumannHeatFlux", group="ADDITIONALS_HEAT")
 ```
 
 ## Dynamics.c.Rt
@@ -295,7 +294,6 @@ CudaDeviceFunction void Run() {
 // This defines the dynamics that we run at each node in the domain.
 	switch (NodeType & NODE_BOUNDARY) {
 		case NODE_Wall:
-		case NODE_Solid:
 			BounceBack();
 			break;
 	}
@@ -305,10 +303,14 @@ CudaDeviceFunction void Run() {
 			CollisionBGK();
 			break;
 		case NODE_MRT:
-		case NODE_CM:
+		//case NODE_CM:
 			//CollisionCM();
 			break;
-			}
+		}
+
+	if ((NodeType & NODE_ADDITIONALS_HEAT) == NODE_HeaterDirichletTemperature) {
+		 HeatDirichletEquilibriumScheme();
+	}
 }
 ```
 
@@ -370,104 +372,33 @@ CudaDeviceFunction void CollisionBGK() {
 			// Calculate collision for temperature populations
 			// h[i] = ...;
 	}
-
-
-	// Below, you can add BCs which take place inside the fluid.
 }
 
 ```
 ### Boundary Conditions
 
-#### Wall Temperature: Dirichlet BC
+#### Temperature: Dirichlet BC
 
 ```c
-// at the wall
-CudaDeviceFunction void BounceBack()
+CudaDeviceFunction void HeatDirichletEquilibriumScheme()
 {
-	real_t tmp;
-	tmp = f[6];
-	f[6] = f[8];
-	f[8] = tmp;
-	tmp = f[2];
-	f[2] = f[4];
-	f[4] = tmp;
-	tmp = f[1];
-	f[1] = f[3];
-	f[3] = tmp;
-	tmp = f[5];
-	f[5] = f[7];
-	f[7] = tmp;
-	tmp = h[6];
-	h[6] = h[8];
-	h[8] = tmp;
-	tmp = h[2];
-	h[2] = h[4];
-	h[4] = tmp;
-	tmp = h[1];
-	h[1] = h[3];
-	h[3] = tmp;
-	tmp = h[5];
-	h[5] = h[7];
-	h[7] = tmp;
-
-	// Anti-Bounce-Back Scheme
-	// see chapter 8.5.2.1, eq 8.53, p318 from 'The Lattice Boltzmann Method: Principles and Practice'
-	// by T. Krüger, H. Kusumaatmaja, A. Kuzmin, O. Shardt, G. Silva, E.M. Viggen
-	if ((NodeType & NODE_ADDITIONALS) == NODE_HeaterDirichletTemperature)
-	{
-		for (int i = 0; i < 9; i++) {
-			h[i] = -h[i];} // reverse to obtain zero temperature
-		
-		real_t rho = getRho();
-		real_t rhoT = 2*rho*InitTemperature; // anti bounce-back - prescribe temperature
-		vector_t u = getU();
-		real_t h_source[9];
-		SetEquilibrium(h_source, rhoT, u);
-
-		for (int i = 0; i < 9; i++) {
-			h[i] += h_source[i];}
-	}
-}
-```
-
-#### Within Fluid
-
-```c
-CudaDeviceFunction void CollisionBGK() 
-{
-	// ... collision code ...
-	// apply BC to post-collision populations.
-
-	// equilibrium scheme for BC
+	// equilibrium scheme for BC - don't care and impose rho*Teq
 	// see chapter 5.3.4.2, eq 5.34, p191 from 'The Lattice Boltzmann Method: Principles and Practice'
 	// by T. Krüger, H. Kusumaatmaja, A. Kuzmin, O. Shardt, G. Silva, E.M. Viggen
 
-	// fluid Temperature: Dirichlet BC
-	if ((NodeType & NODE_ADDITIONALS) == NODE_HeaterDirichletTemperature)
-	{
-		// we are after collision in this if-scope {}
-		real_t d_rhoT = m00*(InitTemperature - localTemperature);
-		real_t h_source[9];
-		SetEquilibrium(h_source, d_rhoT, u);
-
-		for (int i = 0; i < 9; i++) {
-			h[i] += h_source[i];
-		}
-	}
+	real_t T = InitTemperature;
+	real_t m00 = getRho();
+	vector_t u = getU();
+	
+	real_t d_rhoT = m00*T;
+	real_t h_eq[9];
+	SetEquilibrium(h_eq, d_rhoT, u);
 
-	// heat source within fluid: Neumann BC
-	if ((NodeType & NODE_ADDITIONALS) == NODE_HeaterNeumannHeatFlux)
-	{
-		// TCLB is smart enough to distiguish values prescribed it batch.xml like
-		// InitTemperature-dirichlet_region="123"
-		// InitTemperature-neumann_heater_region="456"
-		real_t d_rhoT = m00*InitTemperature;  
-		real_t h_source[9];
-		SetEquilibrium(h_source, d_rhoT, u); 
-		for (int i = 0; i < 9; i++) {
-			h[i] += h_source[i];}
+	for (int i = 0; i < 9; i++) {
+		h[i] = h_eq[i];
 	}
 }
+
 ```
 
 ## Setting up a Simulation
@@ -476,33 +407,32 @@ Create a configuration file, `my_heat_cube.xml` in example/
 
 ```xml
 <?xml version="1.0"?>
-<?xml version="1.0"?>
 <CLBConfig version="2.0" output="output/">
       <Geometry nx="128" ny="128" >
-      <BGK><Box /></BGK>    
+      <BGK><Box /></BGK>
 
         <None name="heated_region">
             <!-- <Sphere dx="80" nx="10" dy="25" ny="10" />  -->
         </None>        
 
         <HeaterDirichletTemperature  name="dirichlet_hot"> 
-          <Box nx="1" dx="-3"/>
+          <Box nx="1" dx="-2"/>
         </HeaterDirichletTemperature>
 
         <HeaterDirichletTemperature  name="dirichlet_cold"> 
-          <Box dx="2" nx="1" />
+          <Box dx="1" nx="1" />
         </HeaterDirichletTemperature>
-      
+
         <Wall mask="ALL" name="border" >
-                <Box nx="2"/>
-                <Box dx="-2"/>
-                <Box ny="2"/>
-                <Box dy="-2"/>
+                <Box nx="1"/>
+                <Box dx="-1"/>
+                <Box ny="1"/>
+                <Box dy="-1"/>
         </Wall>
   </Geometry>
     <Model>
     <Params
-        k="0.1666666" 
+        k="0.1666666"
         cp="1.0"
 
         InitTemperature="0.0"
@@ -517,7 +447,7 @@ Create a configuration file, `my_heat_cube.xml` in example/
     <Params GravitationY="-1.0E-05" />
     <Params BoussinesqCoeff="10"/>
 
-    </Model> 
+    </Model>
 
 <Solve Iterations="10" output="output/"> <VTK Iterations="1"/> </Solve>
 <Failcheck Iterations="500" nx="128" ny="128" />
@@ -527,7 +457,8 @@ Create a configuration file, `my_heat_cube.xml` in example/
 ```
 
 !!! question
-    What is wrong with the following formulation for BC?
+		What is the difference between specifing the BC for temperature before/after the `<Wall mask="ALL" name="border" >... </Wall>` region?
+		Is the following description correct?
 		```
 				<HeaterDirichletTemperature  name="dirichlet_hot">
 					<Box dx="-1"/>
diff --git a/docs/tutorials/model-development/5.-SymbolicOperations_and_CLBM.md b/docs/tutorials/model-development/5.-SymbolicOperations_and_CLBM.md
index e553f6d93..0c1988509 100644
--- a/docs/tutorials/model-development/5.-SymbolicOperations_and_CLBM.md
+++ b/docs/tutorials/model-development/5.-SymbolicOperations_and_CLBM.md
@@ -229,33 +229,6 @@ CudaDeviceFunction void CollisionCM()
 
 	relax_and_collide_hydro_with_F(f, omega_nu, u, Force);
 	relax_and_collide_ADE(h, omega_k, u);
-
-	// equilibrium scheme for BC
-	// see chapter 5.3.4.2, eq 5.34, p191 from 'The Lattice Boltzmann Method: Principles and Practice'
-	// by T. Krüger, H. Kusumaatmaja, A. Kuzmin, O. Shardt, G. Silva, E.M. Viggen
-	if ((NodeType & NODE_ADDITIONALS) == NODE_HeaterDirichletTemperature) 
-	{
-		// we are after collision in this if-scope {}
-		real_t d_rhoT = m00*(InitTemperature - localTemperature);
-		real_t h_source[9];
-		set_eq(h_source, d_rhoT, u);
-
-		for (int i = 0; i < 9; i++) {
-			h[i] += h_source[i];
-		}
-	}
-
-	if ((NodeType & NODE_ADDITIONALS) == NODE_HeaterNeumannHeatFlux) 
-	{
-		// TCLB is smart enough to distiguish values prescribed it batch.xml like
-		// InitTemperature-dirichlet_region="123"
-		// InitTemperature-neumann_heater_region="456"
-		real_t d_rhoT = m00*InitTemperature;  
-		real_t h_source[9];
-		set_eq(h_source, d_rhoT, u); 
-		for (int i = 0; i < 9; i++) {
-			h[i] += h_source[i];}
-	}
 }
 ```
 
@@ -339,7 +312,7 @@ As it takes a while to compute the integrals, the result is hardcoded and used a
 
 ### Collision Kernel
 
-We are going to generate the missing functions:
+Now, it is time to generate the missing functions:
 
 ```c
 CudaDeviceFunction void SetEquilibrium(real_t x_in[9], real_t rho_Xeq, vector_t u){...}
@@ -356,70 +329,8 @@ then to density-probability functions - see the [streaming](#streaming) step.
 
 ## Setting up a Simulation
 
-To run the new collision kernel, change the configuration file, `my_heat_cube.xml` from [previous tutorial](/tutorials/model-development/4.-D2Q9-HeatTransfer) :
-`<BGK><Box /></BGK>` to `<CM><Box /></CM>`.  
-
-Finally, you may try a another geometry to see how the BC behavies
-
-```xml
-<?xml version="1.0"?>
-<CLBConfig version="2.0" output="output/">
-      <Geometry nx="500" ny="100" >
-      <CM><Box /></CM>    
-        <None name="heated_region">
-            <Sphere dx="175" nx="50" dy="25" ny="50" /> 
-        </None>        
-
-        <Solid name="solid_region">
-          <Box dx="100" dy="50" nx="20" ny="25"/>
-        </Solid>
-
-        <Wall mask="ALL" name="box" >     
-           <Box dx="100" dy="25" nx="20" ny="25"/>
-        </Wall>
-
-        <HeaterDirichletTemperature  name="dirichlet_region"> 
-          <Box dx="200" dy="50" nx="20" ny="25"/>
-        </HeaterDirichletTemperature>
-        
-        <HeaterNeumannHeatFlux  name="neumann_heater_region"> 
-          <Box dx="300" dy="25" nx="20" ny="25"/>
-        </HeaterNeumannHeatFlux>
-        <HeaterNeumannHeatFlux  name="neumann_cool_region"> 
-          <Box dx="350" dy="50" nx="20" ny="25"/>
-        </HeaterNeumannHeatFlux>
-
-      	<Solid name="solid_region">
-      	  <Box dx="400" dy="10" nx="20" ny="40"/>
-      	</Solid>
-
-  </Geometry>
-    <Model>
-    <Params
-        k="0.1666666" 
-        cp="1.0"
-
-        InitTemperature="0.0"
-        InitTemperature-heated_region="10"
-        InitTemperature-dirichlet_region="100.0"
-
-        InitTemperature-neumann_heater_region="1E0"
-        InitTemperature-neumann_cool_region="-1E0"
-        />
-    <Params  />
-    <Params GravitationX="-1.0E-5" />
-    <!-- <Params GravitationY="-1.0E-7" /> -->
-    <!-- <Params BoussinesqCoeff="1"/> -->
-    <!-- <Params VelocityX="-1.0E-1" /> -->
-
-    </Model> 
-
-<Solve Iterations="10" output="output/"> <VTK Iterations="1"/> </Solve>
-<Failcheck Iterations="10" nx="500" ny="100" />
-<Log Iterations="1000"/>
-<Solve Iterations="20000" output="output/"> <VTK Iterations="250"/> </Solve>
-</CLBConfig>
-```
+To run the new collision kernel, change the configuration file `my_heat_cube.xml` from [previous tutorial](/tutorials/model-development/4.-D2Q9-HeatTransfer) :
+`<BGK><Box /></BGK>` to `<CM><Box /></CM>` and uncomment the corresponding switch ` //case NODE_CM: ...` in Dynamics.c.Rt.  
 
 ## References
 
diff --git a/aux/boxes.R b/scripts/boxes.R
similarity index 100%
rename from aux/boxes.R
rename to scripts/boxes.R
diff --git a/aux/pipe.R b/scripts/pipe.R
similarity index 100%
rename from aux/pipe.R
rename to scripts/pipe.R