puts "area by geometry: $g_area"
expr ($t_area - $g_area) / $g_area * 100
}
+
+# Check if list of xdistcs-command is valid
+proc checkList {List Tolerance D_good} {
+ set L1 [llength ${List}]
+ set L2 10
+ set L3 5
+ set N [expr (${L1} - ${L2})/${L3} + 1]
+
+ for {set i 1} {${i} <= ${N}} {incr i} {
+ set j1 [expr ${L2} + (${i}-1)*${L3}]
+ set j2 [expr ${j1} + 2]
+ set T [lindex ${List} ${j1}]
+ set D [lindex ${List} ${j2}]
+ #puts "i=${i} j1=${j1} j2=${j2} T=${T} D=${D}"
+ if { [expr abs(${D} - ${D_good})] > ${Tolerance} } {
+ puts "Error: i=${i} T=${T} D=${D}"
+ }
+ }
+}