sub-grid interpolation now also affects basal melt in '-cold' mode; introduce option 'scale_bmr_gl' which scales bmr in sub-grid cells with provided factor