I think the loops in glacier module go from `0` to `<= Nbands`. This is inconsistent, either in the definition of `Nbands` or in how  `GlAreaFract` is used (see earlier comment). This is a likely source of error and should be fixed.