trend.tcl

This example shows how to use surface as trend and how to modify trend with fault lines.

Here is the image of the resulting surface:

trend1

Pic 1. surface obtained from points with trend and 2 fault lines

trend2

Pic 2. surface obtained from points with trend and 2 fault lines

This result was obtained with "examples/surfit/surfs/surface_ineq.tcl" script.

00001 #
00002 # This script was generated with funner v2.1beta (www.gridding.info)
00003 
00004 # 
00005 # load plugins
00006 # 
00007 load libsurfit[info sharedlibextension]
00008 
00009 # remove all previous data and gridding rules
00010 clear_data 
00011 
00012 # set name of surface
00013 set map_name "map_trend" 
00014 
00015 # set solver
00016 set_solver "cg" 
00017 
00018 # set tolerance for solver
00019 set tol 1e-005 
00020 
00021 ##
00022 ## load initial data 
00023 ##
00024  
00025 # load points from text file 
00026 pnts_read "7points.txt" "7points"  
00027  
00028 # load curve from text file 
00029 curv_read "1fault2.txt" "1fault2"  
00030  
00031 # load curve from text file 
00032 curv_read "1fault.txt" "1fault"  
00033  
00034 # load surface from surfit datafile 
00035 surf_load "func.func" "test_func"  
00036  
00037 ##
00038 ## construct grid 
00039 ##
00040 grid_get -10 10 0.2 -10 10 0.2 
00041  
00042 ##
00043 ## create gridding rules
00044 ##
00045 
00046 # resulting surface at points = points values 
00047 points "7points" 
00048 
00049 # resulting surface looks like trend surface 
00050 trend 1 0.2 "test_func" 
00051 
00052 # fault line 
00053 fault "1fault"  
00054 
00055 # fault line 
00056 fault "1fault2"  
00057 
00058 ##
00059 ## run gridding algorithm
00060 ##
00061 surfit 
00062 
00063 ##
00064 ## save results 
00065 ##
00066 
00067 # unload grid from memory
00068 grid_unload 
00069 
00070 # save surface to surfit datafile 
00071 surf_save "trend.dat" "map_trend"



surfit: open source gridding and contouring software.