finetune flat area beautification

git-svn-id: https://svn.osgeo.org/grass/grass/trunk@42238 15284696-431f-4ddb-bdfa-cd5b030d7da7

raster/r.watershed/ram/do_flatarea.c

@@ -156,8 +156,6 @@ int do_flatarea(int index, CELL ele)
 
 	seg_index_rc(alt_seg, index_doer, &r, &c);
 
-	FLAG_SET(flat_done, r, c);
-
 	/* check all neighbours, breadth first search */
 	for (ct_dir = 0; ct_dir < sides; ct_dir++) {
 	    /* get r, c (upr, upc) for this neighbour */
@@ -354,7 +352,6 @@ int do_flatarea(int index, CELL ele)
 			    pq_add(index_up, up_pq);
 			    /* unset flag */
 			    nbr_order_found->flag = 0;
-			    //nbr_order_found->downhill = 0;
 			}
 		    }
 		}
@@ -375,6 +372,8 @@ int do_flatarea(int index, CELL ele)
 	seg_index_rc(alt_seg, index_doer, &r, &c);
 	this_in_list = FLAG_GET(in_list, r, c);
 
+	FLAG_SET(flat_done, r, c);
+
 	/* get uphill and downhill order for this point */
 	inc_order.index = index_doer;
 	if ((order_found = rbtree_find(order_tree, &inc_order)) == NULL)
@@ -394,7 +393,7 @@ int do_flatarea(int index, CELL ele)
 	    downhill_order = max_downhill_order;
 	    uphill_order = 0;
 	}
-	alt[index_doer] += uphill_order + max_downhill_order - downhill_order;
+	alt[index_doer] += uphill_order + (double)(max_downhill_order - downhill_order) / 2.0 + 0.5;
 
 	/* check all neighbours, breadth first search */
 	for (ct_dir = 0; ct_dir < sides; ct_dir++) {

raster/r.watershed/ram/init_vars.c

@@ -159,6 +159,8 @@ int init_vars(int argc, char *argv[])
 
     if (ele_map_type == FCELL_TYPE || ele_map_type == DCELL_TYPE)
 	ele_scale = 10000; 	/* should be enough to do the trick */
+    if (flat_flag)
+	ele_scale = 10000;
 
     /* read elevation input and mark NULL/masked cells */
     /* intialize accumulation and drainage direction */
@@ -186,6 +188,7 @@ int init_vars(int argc, char *argv[])
 	    else {
 		if (ele_map_type == CELL_TYPE) {
 		    alt_value = *((CELL *)ptr);
+		    alt_value *= ele_scale;
 		}
 		else if (ele_map_type == FCELL_TYPE) {
 		    dvalue = *((FCELL *)ptr);