Thermal Model iterative execution error

[jzmejia]

Error in ib.temperature (ThermalModel instance) when running model at the first thermal model timestep to recalculate ice block temperatures and refreezing

error is raised from
ice block.increment_time()

# Run temperature solver with updated geometry
if self.n % self.thermal_freq == 0:
    self.temperature.calc_temperature(self.crev_locs)

for num, crev in enumerate(self.crev_idx):
     # upstream of crevasse
     if num+1 < len(self.crev_idx) or max(crev) > rhs.size-nx:
         idx = [x-1 for x in crev if x < idx_surface]
         rhs[idx] = rhs[idx] + (self.Lf/self.heat_capacity_intercept * self.bluelayer_right[0])/self.ibg.dx

error there is pointing to the line with self.bluelayer_right which is updated by .calc_refreezing()