diff --git a/Applications/FileIO/GocadIO/CoordinateSystem.h b/Applications/FileIO/GocadIO/CoordinateSystem.h
index 4328887cc17ecb8f606615ad9db83a92834b6cfe..cd500b6eb45be66437fc71c00e7ca093cc6f069d 100644
--- a/Applications/FileIO/GocadIO/CoordinateSystem.h
+++ b/Applications/FileIO/GocadIO/CoordinateSystem.h
@@ -7,6 +7,8 @@
  *              http://www.opengeosys.org/project/license
  */
 
+#pragma once
+
 #include <fstream>
 
 namespace FileIO
diff --git a/MaterialLib/Adsorption/Density100MPa.cpp b/MaterialLib/Adsorption/Density100MPa.cpp
index f5efe3ce2f0a303c8ece0cd8251517960d200f60..23e1dd81fc75a0437fe5df97ab6691f2a6af2813 100644
--- a/MaterialLib/Adsorption/Density100MPa.cpp
+++ b/MaterialLib/Adsorption/Density100MPa.cpp
@@ -14,7 +14,7 @@ namespace
 
 // NaX_HighP_polyfrac_CC.pickle
 // date extracted 2015-06-23 15:38:35 file mtime 2015-06-23 15:19:57
-const double c[] = {
+const double Density100MPa_c[] = {
     0.3490302932983226,        /* a0 */
     -0.0014061345691831226,    /* a1 */
     -0.0007399303393402753,    /* a2 */
@@ -46,7 +46,7 @@ double Density100MPa::getAlphaT(const double T_Ads) const
 // Characteristic curve. Return W (A)
 double Density100MPa::characteristicCurve(const double A) const
 {
-    double W = curvePolyfrac(c, A); // cm^3/g
+    double W = curvePolyfrac(Density100MPa_c, A); // cm^3/g
 
     if (W < 0.0) {
         W = 0.0; // TODO [CL] debug output
@@ -57,7 +57,7 @@ double Density100MPa::characteristicCurve(const double A) const
 
 double Density100MPa::dCharacteristicCurve(const double A) const
 {
-    return dCurvePolyfrac(c, A);
+    return dCurvePolyfrac(Density100MPa_c, A);
 }
 
 }  // namespace Adsorption
diff --git a/MaterialLib/Adsorption/DensityConst.cpp b/MaterialLib/Adsorption/DensityConst.cpp
index b3048a2044a92aea01548a2a71eb9dd856965404..07a87c3e49a8e8026a370a0ebcd3e6780dece3f1 100644
--- a/MaterialLib/Adsorption/DensityConst.cpp
+++ b/MaterialLib/Adsorption/DensityConst.cpp
@@ -15,7 +15,7 @@ namespace
 
 // NaX_Constant_polyfrac_CC.pickle
 // date extracted 2015-06-23 15:38:35 file mtime 2015-06-23 15:20:05
-const double c[] = {
+const double DensityConst_c[] = {
     0.3824098506898007,            /* a0 */
     -0.001316857559708455,        /* a1 */
     -0.0007935756090263691,        /* a2 */
@@ -43,7 +43,7 @@ double DensityConst::getAlphaT(const double /*T_Ads*/) const
 // Characteristic curve. Return W (A)
 double DensityConst::characteristicCurve(const double A) const
 {
-    double W = curvePolyfrac(c, A); //cm^3/g
+    double W = curvePolyfrac(DensityConst_c, A); //cm^3/g
 
     if (W < 0.0) {
         W = 0.0; // TODO [CL] debug output
@@ -54,7 +54,7 @@ double DensityConst::characteristicCurve(const double A) const
 
 double DensityConst::dCharacteristicCurve(const double A) const
 {
-    return dCurvePolyfrac(c, A);
+    return dCurvePolyfrac(DensityConst_c, A);
 }
 
 }  // namespace Adsorption
diff --git a/MaterialLib/Adsorption/DensityCook.cpp b/MaterialLib/Adsorption/DensityCook.cpp
index d78f6d578112773c82ce58a928af658635af727b..ab71403c0cc835aa245915499fdde52715894491 100644
--- a/MaterialLib/Adsorption/DensityCook.cpp
+++ b/MaterialLib/Adsorption/DensityCook.cpp
@@ -14,7 +14,7 @@ namespace
 
 // NaX_Dean_polyfrac_CC.pickle
 // date extracted 2015-06-23 15:38:35 file mtime 2015-06-23 15:19:42
-const double c[] = {
+const double DensityCook_c[] = {
     0.3632627555646154,        /* a0 */
     -0.0014090624975800715,    /* a1 */
     -0.0007717609035743321,    /* a2 */
@@ -42,7 +42,7 @@ double DensityCook::getAlphaT(const double T_Ads) const
 // Characteristic curve. Return W (A)
 double DensityCook::characteristicCurve(const double A) const
 {
-    double W = curvePolyfrac(c, A); //cm^3/g
+    double W = curvePolyfrac(DensityCook_c, A); //cm^3/g
 
     if (W < 0.0) {
         W = 0.0; // TODO [CL] debug output
@@ -53,7 +53,7 @@ double DensityCook::characteristicCurve(const double A) const
 
 double DensityCook::dCharacteristicCurve(const double A) const
 {
-    return dCurvePolyfrac(c, A);
+    return dCurvePolyfrac(DensityCook_c, A);
 }
 
 }  // namespace Adsorption
diff --git a/MaterialLib/Adsorption/DensityDubinin.cpp b/MaterialLib/Adsorption/DensityDubinin.cpp
index 7587fd75f9d82680860b157002e0b9a773d85186..c55e65ccf658cb2ce74f6eed48ed28f6fc321402 100644
--- a/MaterialLib/Adsorption/DensityDubinin.cpp
+++ b/MaterialLib/Adsorption/DensityDubinin.cpp
@@ -17,7 +17,7 @@ namespace
 
 // NaX_Dubinin_polyfrac_CC.pickle
 // date extracted 2015-06-23 16:47:50 file mtime 2015-06-23 16:47:23
-const double c[] = {
+const double DensityDubinin_c[] = {
     0.3635538371322433,        /* a0 */
     -0.0014521033261199435,    /* a1 */
     -0.0007855160157616825,    /* a2 */
@@ -77,7 +77,7 @@ double DensityDubinin::getAlphaT(const double T_Ads) const
 // Characteristic curve. Return W (A)
 double DensityDubinin::characteristicCurve(const double A) const
 {
-    double W = curvePolyfrac(c, A); // cm^3/g
+    double W = curvePolyfrac(DensityDubinin_c, A); // cm^3/g
 
     if (W < 0.0) {
         W = 0.0; // TODO [CL] debug output
@@ -88,7 +88,7 @@ double DensityDubinin::characteristicCurve(const double A) const
 
 double DensityDubinin::dCharacteristicCurve(const double A) const
 {
-    return dCurvePolyfrac(c, A);
+    return dCurvePolyfrac(DensityDubinin_c, A);
 }
 
 }  // namespace Adsorption
diff --git a/MaterialLib/Adsorption/DensityHauer.cpp b/MaterialLib/Adsorption/DensityHauer.cpp
index b0777ecd5971b0f9df5b0272049740e0aacfec78..e40af09893ff87ec74b24cb2689b86ab44e96f9a 100644
--- a/MaterialLib/Adsorption/DensityHauer.cpp
+++ b/MaterialLib/Adsorption/DensityHauer.cpp
@@ -14,7 +14,7 @@ namespace
 
 // NaX_Hauer_polyfrac_CC.pickle
 // date extracted 2015-06-23 15:38:35 file mtime 2015-06-23 15:19:19
-const double c[] = {
+const double DensityHauer_c[] = {
     0.36490158988356747,    /* a0 */
     -0.0013723270478333963,    /* a1 */
     -0.0007655780628099964,    /* a2 */
@@ -47,7 +47,7 @@ double DensityHauer::getAlphaT(const double T_Ads) const
 // Characteristic curve. Return W (A)
 double DensityHauer::characteristicCurve(const double A) const
 {
-    double W = curvePolyfrac(c, A); // cm^3/g
+    double W = curvePolyfrac(DensityHauer_c, A); // cm^3/g
 
     if (W < 0.0) {
         W = 0.0; // TODO [CL] debug output
@@ -58,7 +58,7 @@ double DensityHauer::characteristicCurve(const double A) const
 
 double DensityHauer::dCharacteristicCurve(const double A) const
 {
-    return dCurvePolyfrac(c, A);
+    return dCurvePolyfrac(DensityHauer_c, A);
 }
 
 }  // namespace Adsorption
diff --git a/MaterialLib/Adsorption/DensityLegacy.cpp b/MaterialLib/Adsorption/DensityLegacy.cpp
index 38f612676aef4da59e5037db520a64765fd9eb99..cabf257d5bf0bdcd68dddaeea4b1bae1c6cf5fe5 100644
--- a/MaterialLib/Adsorption/DensityLegacy.cpp
+++ b/MaterialLib/Adsorption/DensityLegacy.cpp
@@ -13,7 +13,7 @@ namespace
 {
 
 //parameters from least squares fit (experimental data)
-const double c[] = { 0.34102920966608297,
+const double DensityLegacy_c[] = { 0.34102920966608297,
                      -0.0013106032830951296,
                      -0.00060754147575378876,
                      3.7843404172683339e-07,
@@ -50,7 +50,7 @@ double DensityLegacy::getAlphaT(const double T_Ads) const
 // Characteristic curve. Return W (A)
 double DensityLegacy::characteristicCurve(const double A) const
 {
-    double W = curvePolyfrac(c, A); // cm^3/g
+    double W = curvePolyfrac(DensityLegacy_c, A); // cm^3/g
 
     if (W < 0.0) {
         W = 0.0; // TODO [CL] debug output
@@ -61,7 +61,7 @@ double DensityLegacy::characteristicCurve(const double A) const
 
 double DensityLegacy::dCharacteristicCurve(const double A) const
 {
-    return dCurvePolyfrac(c, A);
+    return dCurvePolyfrac(DensityLegacy_c, A);
 }
 
 }  // namespace Adsorption
diff --git a/MaterialLib/Adsorption/DensityMette.cpp b/MaterialLib/Adsorption/DensityMette.cpp
index 325d53f38215dc0f5736082546471f6823613701..c2b26e7f9b75cf7ab84f0f9d49f57dfdd26bbf15 100644
--- a/MaterialLib/Adsorption/DensityMette.cpp
+++ b/MaterialLib/Adsorption/DensityMette.cpp
@@ -14,7 +14,7 @@ namespace
 {
 // NaX_Mette_polyfrac_CC.pickle
 // date extracted 2015-06-23 15:38:35 file mtime 2015-06-23 15:19:26
-const double c[] = {
+const double DensityMette_c[] = {
     0.36340572890087813,    /* a0 */
     -0.0013449597038375108,    /* a1 */
     -0.0007581210111121073,    /* a2 */
@@ -49,7 +49,7 @@ double DensityMette::getAlphaT(const double T_Ads) const
 // Characteristic curve. Return W (A)
 double DensityMette::characteristicCurve(const double A) const
 {
-    double W = curvePolyfrac(c, A); // cm^3/g
+    double W = curvePolyfrac(DensityMette_c, A); // cm^3/g
 
     if (W < 0.0) {
         W = 0.0; // TODO [CL] debug output
@@ -60,7 +60,7 @@ double DensityMette::characteristicCurve(const double A) const
 
 double DensityMette::dCharacteristicCurve(const double A) const
 {
-    return dCurvePolyfrac(c, A);
+    return dCurvePolyfrac(DensityMette_c, A);
 }
 
 }  // namespace Adsorption
diff --git a/MaterialLib/Adsorption/DensityNunez.cpp b/MaterialLib/Adsorption/DensityNunez.cpp
index 23ae5dcbc59d6346095eb5120ae43dd733ccbc1f..a516fba6df7b63bfbfa23fb9925f3b94b582423c 100644
--- a/MaterialLib/Adsorption/DensityNunez.cpp
+++ b/MaterialLib/Adsorption/DensityNunez.cpp
@@ -14,7 +14,7 @@ namespace
 
 // NaX_Nunez_polyfrac_CC.pickle
 // date extracted 2015-06-23 15:38:35 file mtime 2015-06-23 15:19:34
-const double c[] = {
+const double DensityNunez_c[] = {
     0.3631900485031771,        /* a0 */
     -0.0014242280940080726,    /* a1 */
     -0.0007751726942386291,    /* a2 */
@@ -57,7 +57,7 @@ double DensityNunez::getAlphaT(const double T_Ads) const
 // Characteristic curve. Return W (A)
 double DensityNunez::characteristicCurve(const double A) const
 {
-    double W = curvePolyfrac(c, A); // cm^3/g
+    double W = curvePolyfrac(DensityNunez_c, A); // cm^3/g
 
     if (W < 0.0) {
         W = 0.0; // TODO [CL] debug output
@@ -68,7 +68,7 @@ double DensityNunez::characteristicCurve(const double A) const
 
 double DensityNunez::dCharacteristicCurve(const double A) const
 {
-    return dCurvePolyfrac(c, A);
+    return dCurvePolyfrac(DensityNunez_c, A);
 }
 
 }  // namespace Adsorption
diff --git a/MaterialLib/FractureModels/LogPenalty.h b/MaterialLib/FractureModels/LogPenalty.h
index 68aaffa3ce9202ee791e13824741ee9e05626d4e..af37f2405117e36b5a3670c65b74a78c59ef9b0a 100644
--- a/MaterialLib/FractureModels/LogPenalty.h
+++ b/MaterialLib/FractureModels/LogPenalty.h
@@ -7,6 +7,8 @@
  *              http://www.opengeosys.org/project/license
  */
 
+#pragma once
+
 #include "LinearElasticIsotropic.h"
 
 namespace MaterialLib
diff --git a/ProcessLib/SourceTerms/LineSourceTermFEM.h b/ProcessLib/SourceTerms/LineSourceTermFEM.h
index 6face47302041ab1d447a539a4072c90e5cb99ff..1e394307a347325a22510ba051040e0a6817c59d 100644
--- a/ProcessLib/SourceTerms/LineSourceTermFEM.h
+++ b/ProcessLib/SourceTerms/LineSourceTermFEM.h
@@ -31,13 +31,13 @@ public:
     virtual ~LineSourceTermLocalAssemblerInterface() = default;
 };
 
-const unsigned NUM_NODAL_DOF = 1;
-
 template <typename ShapeFunction, typename IntegrationMethod,
           unsigned GlobalDim>
 class LineSourceTermLocalAssembler final
     : public LineSourceTermLocalAssemblerInterface
 {
+    static const unsigned NUM_NODAL_DOF = 1;
+
     using ShapeMatricesType = ShapeMatrixPolicyType<ShapeFunction, GlobalDim>;
 
     using LocalAssemblerTraits = ProcessLib::LocalAssemblerTraits<