make sure that sub-grid interp. of basal melt also is applied in cold mode
rewrote sub-grid gl interpolation --> now only direct assignment of array elements, e.g. it is written to gl_mask(i,j) but not gl_mask(i+1,j), works
avoid seldomly occuring negative values, clean up
sub gl interpolation: Fix number of vWork2d vectors