fault_approx.tcl

This example shows how one goes from randomly spaced data points to an evenly sampled surface with taking fault curves into account. In this sample we have 11 scattered data points to approximate (file examples/surfit/curvs/11points.txt) and 4 fault curves, which are divide whole surface domain into 3 independent subdomains. The main idea of fault curves is that parts of resulting surface between fault curve are independent of each other. Important! You should place fault rule only after fault modifiable rules!

See also:

fault

Here is the image of the resulting surface:

Pic 1. Surface with faults

This result was obtained with "examples/surfit/curvs/fault_approx.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_faults_approx" 

# set solver
set_solver "cg" 

# set tolerance for solver
set tol 1e-005 

##
## load initial data 
##
 
# load points from text file 
pnts_read "11points.txt" "11points"  
 
# load curve from text file 
curv_read "flt_line1.txt" "flt_line1"  
 
# load curve from text file 
curv_read "flt_line2.txt" "flt_line2"  
 
# load curve from text file 
curv_read "flt_line3.txt" "flt_line3"  
 
# load curve from text file 
curv_read "flt_line4.txt" "flt_line4"  
 
##
## construct grid 
##
grid_get 1 99 1 -50 50 1 
 
##
## create gridding rules
##

# resulting surface should tend to be constant or plane 
completer 0.5 2  

# fault line 
fault "flt_line1"  

# fault line 
fault "flt_line2"  

# fault line 
fault "flt_line3"  

# fault line 
fault "flt_line4"  

# add "resulting surface at points = points values" with weight 
points_add 1 "11points" 

##
## run gridding algorithm
##
surfit 

##
## save results 
##

# unload grid from memory
grid_unload 

# save surface to surfit datafile 
surf_save "fault_approx.dat" "map_faults_approx" 

---

surfit: open source gridding and contouring software.