added file delta_tilde_rho.m which replaces epsilon.m
removed epsilon.m
removed grad_surf.m
removed scale_fac.m

substitute calculation of epsilon with \tilde{\Delta}\rho
stop wetting when number of added gridpoints stays constant or increases
removed grid cell lengths e1t and e2t from the code

fix bug in solve_lsqr() regarding the removal of northern boundary
fix but in find_regions() regarding northern and southern boundaries