area_wmean.tcl

This example shows how to build surface by approximating areas with constant real numbers and how to set additional condition (weighted mean value) on surface in area. In this example 6 different areas were used. This example looks like area_mean.tcl example, but for area6 weighted mean condition was set.

Here is the image of the resulting surface (green colour for areas 1,3 and 5, red colour for areas 2 and 4, magenta for area 6):

Pic 1. Surface with weighted mean value = 30, built from areas

This result was obtained with "examples/surfit/curvs/area_wmean.tcl" script.

#
# This script was generated with funner v2.1beta (www.gridding.info)

# 
# load plugins
# 
load libsurfit[info sharedlibextension]

# remove all previous data and gridding rules
clear_data 

# set name of surface
set map_name "map_area_wmean" 

# set solver
set_solver "cg" 

# set tolerance for solver
set tol 1e-005 

##
## load initial data 
##
 
# load area from text file 
area_read "area1.txt" "area1"  
 
# load area from text file 
area_read "area2.txt" "area2"  
 
# load area from text file 
area_read "area3.txt" "area3"  
 
# load area from text file 
area_read "area4.txt" "area4"  
 
# load area from text file 
area_read "area5.txt" "area5"  
 
# load area from text file 
area_read "area6.txt" "area6"  
 
# load surface from surfit datafile 
surf_load "func.func" "test_func"  
 
##
## construct grid 
##
grid_get -12 12 0.2 -12 12 0.2 
 
##
## create gridding rules
##

# resulting surface in area = value... 
area 10 "area1"  

# resulting surface in area = value... 
area 30 "area3"  

# resulting surface in area = value... 
area undef "area5"  

# resulting surface weighted mean value in area = value... 
area_wmean 30 "area6" "test_func"  

# resulting surface should tend to be constant or plane 
completer 1 1.6  

# add "resulting surface in area = value..." with weight 
area_add 20 0.07 "area2"  

# add "resulting surface in area = value..." with weight 
area_add 40 0.07 "area4"  

##
## run gridding algorithm
##
surfit 

##
## save results 
##

# unload grid from memory
grid_unload 

# save surface to surfit datafile 
surf_save "area_wmean.dat" "map_area_wmean" 

puts [concat "Surface weighted mean value in area6 : " [surf_area_wmean "area6" "test_func" "map_area_wmean"]]

---

surfit: open source gridding and contouring software.