diff --git a/models/electrokinetic/d2q9_poison_boltzmann/Dynamics.R b/models/electrokinetic/d2q9_poison_boltzmann/Dynamics.R
index ef44b30ae..f7da6873d 100644
--- a/models/electrokinetic/d2q9_poison_boltzmann/Dynamics.R
+++ b/models/electrokinetic/d2q9_poison_boltzmann/Dynamics.R
@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 
 AddDensity( name="g[0]", dx= 0, dy= 0, group="g")
 AddDensity( name="g[1]", dx= 1, dy= 0, group="g")
diff --git a/models/experimental/d3q27_csf/Dynamics.R b/models/experimental/d3q27_csf/Dynamics.R
index 34cd7cf65..89c0d41e5 100644
--- a/models/experimental/d3q27_csf/Dynamics.R
+++ b/models/experimental/d3q27_csf/Dynamics.R
@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 
 U = expand.grid(-1:1,-1:1,-1:1)
 
diff --git a/models/flow/auto/Dynamics.R b/models/flow/auto/Dynamics.R
index 2fc7e5a6b..50dd273ce 100644
--- a/models/flow/auto/Dynamics.R
+++ b/models/flow/auto/Dynamics.R
@@ -57,8 +57,9 @@ AddNodeType(name="Body", group="BODY")
 for (f in fname) AddField(f,dx=0,dy=0,dz=0) # Make f accessible also in present node (not only streamed)
 
 if (Options$part) {
-	AddStage("BaseIteration", "Run", save=Fields$group %in% c("f","Force"), load = DensityAll$group %in% c("f","Force"))
-	AddStage("CalcF", save=Fields$group == "Force", load = DensityAll$group %in% c("f","Force"), particle=TRUE)
+	AddStage("BaseIteration", "Run", save=Fields$group %in% c("f"), load = DensityAll$group %in% c("f","Force"))
+	AddStage("BaseInit", "Init", save=Fields$group %in% c("f"))
+	AddStage("CalcF", save=Fields$group == "Force", load = DensityAll$group %in% c("f"), particle=TRUE)
 	AddAction("Iteration", c("BaseIteration", "CalcF"))
 	AddAction("Init", c("BaseInit", "CalcF"))
 }
diff --git a/models/flow/d2q9_par/Dynamics.R b/models/flow/d2q9_par/Dynamics.R
index 3c26d0f02..c479e9e04 100644
--- a/models/flow/d2q9_par/Dynamics.R
+++ b/models/flow/d2q9_par/Dynamics.R
@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 # Density - table of variables of LB Node to stream
 #  name - variable name to stream
 #  dx,dy,dz - direction of streaming
diff --git a/models/flow/d3q27_cumulant_part/Dynamics.R b/models/flow/d3q27_cumulant_part/Dynamics.R
index e26c22cb1..66ecc70f9 100644
--- a/models/flow/d3q27_cumulant_part/Dynamics.R
+++ b/models/flow/d3q27_cumulant_part/Dynamics.R
@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 x = c(0,1,-1);
 P = expand.grid(x=0:2,y=0:2,z=0:2)
 U = expand.grid(x,x,x)
diff --git a/models/heat/d3q27_tePSM_per/Dynamics.R b/models/heat/d3q27_tePSM_per/Dynamics.R
index 0389fdecc..5b57833e4 100644
--- a/models/heat/d3q27_tePSM_per/Dynamics.R
+++ b/models/heat/d3q27_tePSM_per/Dynamics.R
@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 
 AddDensity( name="f[0]", dx= 0, dy= 0, dz= 0, group="f")
 AddDensity( name="f[1]", dx= 1, dy= 0, dz= 0, group="f")
diff --git a/models/multiphase/d2q9_ShanChen/Dynamics.R b/models/multiphase/d2q9_ShanChen/Dynamics.R
index 43f1b1a64..9aecb546f 100644
--- a/models/multiphase/d2q9_ShanChen/Dynamics.R
+++ b/models/multiphase/d2q9_ShanChen/Dynamics.R
@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 
 # Fluid Density Populations
 AddDensity( name="f[0]", dx= 0, dy= 0, group="f")
diff --git a/models/multiphase/d2q9_csf/Dynamics.R b/models/multiphase/d2q9_csf/Dynamics.R
index c4bc455f6..d1a2acdeb 100644
--- a/models/multiphase/d2q9_csf/Dynamics.R
+++ b/models/multiphase/d2q9_csf/Dynamics.R
@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 # Density - table of variables of LB Node to stream
 #  name - variable name to stream
 #  dx,dy,dz - direction of streaming
diff --git a/models/multiphase/d2q9_kuper/Dynamics.R b/models/multiphase/d2q9_kuper/Dynamics.R
index 3f75c35ad..3c8d41d5d 100644
--- a/models/multiphase/d2q9_kuper/Dynamics.R
+++ b/models/multiphase/d2q9_kuper/Dynamics.R
@@ -13,7 +13,7 @@ AddField("phi",stencil2d=1);
 
 AddStage("BaseIteration", "Run", save=Fields$group == "f", load=DensityAll$group == "f")
 AddStage("CalcPhi", save="phi",load=DensityAll$group == "f")
-AddStage("BaseInit", "Init", save=Fields$group == "f", load=DensityAll$group == "f")
+AddStage("BaseInit", "Init", save=Fields$group == "f")
 
 AddAction("Iteration", c("BaseIteration","CalcPhi"))
 AddAction("Init", c("BaseInit","CalcPhi"))
diff --git a/models/multiphase/d2q9_pf/Dynamics.R b/models/multiphase/d2q9_pf/Dynamics.R
index c09415efd..ddea29398 100644
--- a/models/multiphase/d2q9_pf/Dynamics.R
+++ b/models/multiphase/d2q9_pf/Dynamics.R
@@ -43,17 +43,16 @@ if (Options$fd) {
     AddField("phi"      ,stencil2d=1 );
 
     AddStage("BaseIteration", "Run", 
-             load=DensityAll$group == "h" | DensityAll$group == "f" | DensityAll$group == "BC",
-             save=Fields$group=="h" | Fields$group=="f"  
-         ) 
+        load=DensityAll$group == "h" | DensityAll$group == "f" | DensityAll$group == "BC",
+        save=Fields$group=="h" | Fields$group=="f" | Fields$group == "BC"
+    ) 
     AddStage("CalcPhi", 
-             save=Fields$name=="phi" ,  
-             load=DensityAll$group == "h"
-         )
+        save=Fields$name=="phi" ,  
+        load=DensityAll$group == "h"
+    )
     AddStage("BaseInit", "Init",  
-	    load=DensityAll$group == "BC",
-        save=Fields$group=="h" | Fields$group == "f" 
- 
+	    load=FALSE,
+        save=Fields$group=="h" | Fields$group == "f" | Fields$group == "BC"
     )
     AddAction("Iteration", c("BaseIteration","CalcPhi"))
     AddAction("Init", c("BaseInit","CalcPhi"))
diff --git a/models/multiphase/d2q9_pf_pressureEvolution/Dynamics.R b/models/multiphase/d2q9_pf_pressureEvolution/Dynamics.R
index 20555e7a4..1dc79e08d 100644
--- a/models/multiphase/d2q9_pf_pressureEvolution/Dynamics.R
+++ b/models/multiphase/d2q9_pf_pressureEvolution/Dynamics.R
@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 # Density - table of variables of LB Node to stream
 #  name - variable name to stream
 #  dx,dy,dz - direction of streaming
diff --git a/models/multiphase/d2q9_pf_velocity/Dynamics.R b/models/multiphase/d2q9_pf_velocity/Dynamics.R
index d2131991f..d2319bb53 100644
--- a/models/multiphase/d2q9_pf_velocity/Dynamics.R
+++ b/models/multiphase/d2q9_pf_velocity/Dynamics.R
@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 # 	Density - table of variables of LB Node to stream
 #	Pressure Evolution:
 AddDensity( name="g[0]", dx= 0, dy= 0, group="g")
diff --git a/models/multiphase/d2q9_scmp/Dynamics.R b/models/multiphase/d2q9_scmp/Dynamics.R
index 371d07e94..75e14c6f9 100644
--- a/models/multiphase/d2q9_scmp/Dynamics.R
+++ b/models/multiphase/d2q9_scmp/Dynamics.R
@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 
 AddDensity( name="f[0]", dx= 0, dy= 0, group="f")
 AddDensity( name="f[1]", dx= 1, dy= 0, group="f")
diff --git a/models/multiphase/d3q19_kuper/Dynamics.R b/models/multiphase/d3q19_kuper/Dynamics.R
index c7cd34df6..0947f4205 100644
--- a/models/multiphase/d3q19_kuper/Dynamics.R
+++ b/models/multiphase/d3q19_kuper/Dynamics.R
@@ -38,7 +38,7 @@ AddField("phi",stencil3d=1);
 
 AddStage("BaseIteration", "Run", save=Fields$group == "f", load=DensityAll$group == "f")
 AddStage("CalcPhi", save="phi",load=DensityAll$group == "f")
-AddStage("BaseInit", "Init", save=Fields$group == "f", load=DensityAll$group == "f")
+AddStage("BaseInit", "Init", save=Fields$group == "f")
 
 AddAction("Iteration", c("BaseIteration","CalcPhi"))
 AddAction("Init", c("BaseInit","CalcPhi"))
diff --git a/models/multiphase/d3q27_pf_velocity/Dynamics.R b/models/multiphase/d3q27_pf_velocity/Dynamics.R
index 98d4dc2c1..484d7069d 100644
--- a/models/multiphase/d3q27_pf_velocity/Dynamics.R
+++ b/models/multiphase/d3q27_pf_velocity/Dynamics.R
@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 # Density - table of variables of LB Node to stream
 #	Velocity-based Evolution d3q27:
 source("lattice.R")
diff --git a/models/multiphase/experimental/d2q9_pp_LBL/Dynamics.R b/models/multiphase/experimental/d2q9_pp_LBL/Dynamics.R
index dc9769532..6a33db7c4 100644
--- a/models/multiphase/experimental/d2q9_pp_LBL/Dynamics.R
+++ b/models/multiphase/experimental/d2q9_pp_LBL/Dynamics.R
@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 # Fluid Density Populations
 AddDensity( name="f[0]", dx= 0, dy= 0, group="f")
 AddDensity( name="f[1]", dx= 1, dy= 0, group="f")
diff --git a/models/multiphase/experimental/d2q9_pp_MCMP/Dynamics.R b/models/multiphase/experimental/d2q9_pp_MCMP/Dynamics.R
index 9e5aa24bb..cf4bfda23 100644
--- a/models/multiphase/experimental/d2q9_pp_MCMP/Dynamics.R
+++ b/models/multiphase/experimental/d2q9_pp_MCMP/Dynamics.R
@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 # Density - table of variables of LB Node to stream
 AddDensity( name="f[0]", dx= 0, dy= 0, group="f")
 AddDensity( name="f[1]", dx= 1, dy= 0, group="f")
diff --git a/models/reaction/d2q9_reaction_diffusion_system/Dynamics.R b/models/reaction/d2q9_reaction_diffusion_system/Dynamics.R
index 1cf0820ef..b6a6463d6 100644
--- a/models/reaction/d2q9_reaction_diffusion_system/Dynamics.R
+++ b/models/reaction/d2q9_reaction_diffusion_system/Dynamics.R
@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 if (Options$AllenCahn) {
  	Qname = 'Allen-Cahn'
 	DREs <- ('PHI')