diff --git a/src/simulation/Simulation.cpp b/src/simulation/Simulation.cpp
index 836047932..0f0cd9cbc 100644
--- a/src/simulation/Simulation.cpp
+++ b/src/simulation/Simulation.cpp
@@ -2687,6 +2687,11 @@ void Simulation::update_particles_i(int start, int inc)
 					         || t==PT_WTRV)
 						ctempl -= 2.0f*pv[y/CELL][x/CELL];
 					s = 1;
+
+					//A fix for ice with ctype = 0
+					if (t==PT_ICEI && parts[i].ctype==0)
+						parts[i].ctype = PT_WATR;
+
 					if (ctemph>ptransitions[t].thv&&ptransitions[t].tht>-1) {
 						// particle type change due to high temperature
 #ifdef REALISTIC