Here is the image of the resulting surface:
Pic 1. surface obtained from points with trend and 2 fault lines
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.