mask_wmean.tcl

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_mask_wmean" 
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 "points.txt" "points"  
00027  
00028 # load surface from surfit datafile 
00029 surf_load "func.func" "test_func"  
00030  
00031 # load mask from surfit datafile 
00032 mask_load "mask_not.mask" "test_mask"  
00033  
00034 ##
00035 ## construct grid 
00036 ##
00037 grid_get -10 10 0.2 -10 10 0.4 
00038  
00039 ##
00040 ## create gridding rules
00041 ##
00042 
00043 # resulting surface at points = points values 
00044 points "points" 
00045 
00046 # resulting surface weighted mean value = value where mask is true... 
00047 mask_wmean 2 "test_mask" "test_func"  
00048 
00049 # resulting surface should tend to be constant or plane 
00050 completer 
00051 
00052 ##
00053 ## run gridding algorithm
00054 ##
00055 surfit 
00056 
00057 ##
00058 ## save results 
00059 ##
00060 
00061 # unload grid from memory
00062 grid_unload 
00063 
00064 # save surface to surfit datafile 
00065 surf_save "mask_wmean.dat" "map_mask_wmean" 
00066 
00067 puts [concat "Surface weighted mean value for test_mask : " [surf_mask_wmean "test_mask" "test_func" "map_mask_wmean"]]



surfit: open source gridding and contouring software.