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