Solve 3D Heat transfer equation in matlab (Energy balance method - enthalpy method)

3D Heat transfer equation with phase change (enthalpy method). Initial temperature from 0.5 to -3.2 [degC], temperature of well 60 [degC], heat exc...