canyon_trend.tcl

This script is a part of a grand canyon reconstruction example. See globe library description for details.

00001 #
00002 # This script was generated with funner v2.1beta (www.gridding.info)
00003 
00004 # 
00005 # load plugins
00006 # 
00007 load libglobe[info sharedlibextension]
00008 load libsurfit[info sharedlibextension]
00009 
00010 # remove all previous data and gridding rules
00011 clear_data 
00012 
00013 # set name of surface
00014 set map_name "canyon_trend" 
00015 
00016 # set solver
00017 set_solver "cg" 
00018 
00019 # set tolerance for solver
00020 set tol 1e-005 
00021 
00022 ##
00023 ## load initial data 
00024 ##
00025  
00026 # load surface from surfit datafile 
00027 surf_load "N35W112.dat" "N35W112"  
00028  
00029 # load surface from surfit datafile 
00030 surf_load "N35W113.dat" "N35W113"  
00031  
00032 # load surface from surfit datafile 
00033 surf_load "N35W114.dat" "N35W114"  
00034  
00035 # load surface from surfit datafile 
00036 surf_load "N36W112.dat" "N36W112"  
00037  
00038 # load surface from surfit datafile 
00039 surf_load "N36W113.dat" "N36W113"  
00040  
00041 # load surface from surfit datafile 
00042 surf_load "N36W114.dat" "N36W114"  
00043  
00044 # load surface from surfit datafile 
00045 surf_load "N37W112.dat" "N37W112"  
00046  
00047 # load surface from surfit datafile 
00048 surf_load "N37W113.dat" "N37W113"  
00049  
00050 # load surface from surfit datafile 
00051 surf_load "N37W114.dat" "N37W114"  
00052  
00053 ##
00054 ## construct grid 
00055 ##
00056 grid_get -114 -111.6 0.0006 35.6 37 0.0007 
00057  
00058 ##
00059 ## create gridding rules
00060 ##
00061 
00062 # resulting surface = surface values 
00063 surface "N35W112" 
00064 
00065 # resulting surface = surface values 
00066 surface "N35W113" 
00067 
00068 # resulting surface = surface values 
00069 surface "N35W114" 
00070 
00071 # resulting surface = surface values 
00072 surface "N36W112" 
00073 
00074 # resulting surface = surface values 
00075 surface "N36W113" 
00076 
00077 # resulting surface = surface values 
00078 surface "N36W114" 
00079 
00080 # resulting surface = surface values 
00081 surface "N37W112" 
00082 
00083 # resulting surface = surface values 
00084 surface "N37W113" 
00085 
00086 # resulting surface = surface values 
00087 surface "N37W114" 
00088 
00089 ##
00090 ## run gridding algorithm
00091 ##
00092 surfit 
00093 
00094 ##
00095 ## save results 
00096 ##
00097 
00098 # unload grid from memory
00099 grid_unload 
00100 
00101 # save surface to surfit datafile 
00102 surf_save "canyon_trend.dat" "canyon_trend" 



surfit: open source gridding and contouring software.