Is it possible to model gas hydrate simulations in ice-bearing layers? Although there are different ice phases in TOUGH+HYDRATE, when we add ice into the system we got convergence failure problem.