TPT: A fix for ice with ctype = 0. e52f730969

This commit is contained in:
Simon Robertshaw 2012-04-18 21:01:54 +01:00
parent 1681429d7f
commit d85049fb93

View File

@ -2687,6 +2687,11 @@ void Simulation::update_particles_i(int start, int inc)
|| t==PT_WTRV) || t==PT_WTRV)
ctempl -= 2.0f*pv[y/CELL][x/CELL]; ctempl -= 2.0f*pv[y/CELL][x/CELL];
s = 1; 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) { if (ctemph>ptransitions[t].thv&&ptransitions[t].tht>-1) {
// particle type change due to high temperature // particle type change due to high temperature
#ifdef REALISTIC #ifdef REALISTIC