GenericMappingTools / gmt

The Generic Mapping Tools
https://www.generic-mapping-tools.org
Other
865 stars 360 forks source link

Surface crashes in 6.4.0 #7132

Closed mkbarker closed 1 year ago

mkbarker commented 2 years ago

GMT 6.4.0 'surface' routine crashes with this:

gmt blockmedian Site01_scatInterp_v1_final_adj.xyzi -bi4 -R-19/-3/-20/-4 -I0.005 -C -bo3 -r | gmt surface -r -bi3 -R-19/-3/-20/-4 -I0.005 -T0.75 -Gtmp.grd -C0.001 -Ll-11 -Lu11 -Vd

If I leave out the -L flags, then it doesn't crash, but the resulting grid file (tmp.grd) is not correct (it looks like two big triangles).

If I leave out the -r flags (from both blockmedian and surface), then it doesn't crash, and the resulting grid file is correct, but it is gridline registered, which is not what I want.

If I decrease the resolution to 0.01 instead of 0.005 (and keep the -r and -L flags), then it produces the correct pixel-registered grid file.

Full script that generated the error

Link to input file. 1

gmt blockmedian Site01_scatInterp_v1_final_adj.xyzi -bi4 -R-19/-3/-20/-4 -I0.005 -C -bo3 -r | gmt surface -r -bi3 -R-19/-3/-20/-4 -I0.005 -T0.75 -Gtmp.grd -C0.001 -Ll-11 -Lu11 -Vd

Full error message

gmt [DEBUG]: GMT_Create_Session: Terminal width = 195
gmt [DEBUG]: Obtained the ppid from parent: 1674
gmt [DEBUG]: Enter: gmtinit_new_GMT_ctrl
gmt [DEBUG]: GMT->session.SHAREDIR = /panfs/ccds02/nobackup/people/sberton2/.spack_repo/opt/spack/linux-centos7-haswell/gcc-11.1.0/gmt-6.4.0-xdmacvrgjo3w4eprh25qlo4fjsxfjx5u/share
gmt [DEBUG]: GMT->session.HOMEDIR = /home/mkbarker
gmt [DEBUG]: GMT->session.USERDIR = /home/mkbarker/.gmt [created]
gmt [DEBUG]: GMT->session.CACHEDIR = /home/mkbarker/.gmt/cache [created]
gmt [DEBUG]: GMT: 0. Will try to find subdir=postscriptlight stem = PSL_custom_fonts suffix=.txt
gmt [DEBUG]: GMT: 1. gmt_getsharepath trying current dir
gmt [DEBUG]: GMT: 2. gmt_getsharepath trying USERDIR /home/mkbarker/.gmt
gmt [DEBUG]: GMT: 3. gmt_getsharepath trying USERDIR subdir /home/mkbarker/.gmt/postscriptlight
gmt [DEBUG]: GMT: 4. gmt_getsharepath trying SHAREDIR subdir /panfs/ccds02/nobackup/people/sberton2/.spack_repo/opt/spack/linux-centos7-haswell/gcc-11.1.0/gmt-6.4.0-xdmacvrgjo3w4eprh25qlo4fjsxfjx5u/share/postscriptlight
gmt [DEBUG]: GMT: 5. gmt_getsharepath trying SHAREDIR /panfs/ccds02/nobackup/people/sberton2/.spack_repo/opt/spack/linux-centos7-haswell/gcc-11.1.0/gmt-6.4.0-xdmacvrgjo3w4eprh25qlo4fjsxfjx5u/share
gmt [DEBUG]: GMT: 6. gmt_getsharepath failed
gmt [DEBUG]: Map distance calculation will be Cartesian
gmt [DEBUG]: Exit:  gmtinit_new_GMT_ctrl
gmt [DEBUG]: Enter: New_PSL_Ctrl
gmt [DEBUG]: Exit:  New_PSL_Ctrl
gmt [DEBUG]: Enter: gmt_manage_workflow
gmt [DEBUG]: Exit : gmt_manage_workflow
gmt [DEBUG]: Enter: PSL_beginsession
gmt [DEBUG]: Exit : PSL_beginsession
gmt [DEBUG]: Enter: PSL_setdefaults
gmt [DEBUG]: Exit : PSL_setdefaults
gmt [DEBUG]: Enter: gmtlib_io_init
gmt [DEBUG]: Exit : gmtlib_io_init
gmt [DEBUG]: Enter: gmt_hash_init
gmt [DEBUG]: Exit:  gmt_hash_init
gmt [DEBUG]: Enter: gmt_hash_init
gmt [DEBUG]: Exit:  gmt_hash_init
gmt [DEBUG]: Enter: gmt_reload_settings
gmt [DEBUG]: The PROJ_GEODESIC set to Vincenty
gmt [DEBUG]: Look for file /home/mkbarker/gmt.conf
gmt [DEBUG]: Look for file /home/mkbarker/.gmt/gmt.conf
gmt [DEBUG]: Look for file /home/mkbarker/.gmt/server/gmt.conf
gmt [DEBUG]: Look for file /home/mkbarker/.gmt/cache/gmt.conf
gmt [DEBUG]: Could not find file gmt.conf
gmt [DEBUG]: Exit:  gmt_reload_settings
gmt [DEBUG]: Enter: gmtlib_plot_C_format
gmt [DEBUG]: Exit:  gmtlib_plot_C_format
gmt [DEBUG]: Enter: gmtinit_get_history
gmt [DEBUG]: Enter: gmt_hash_init
gmt [DEBUG]: Exit:  gmt_hash_init
gmt [DEBUG]: Exit:  gmtinit_get_history
gmt [DEBUG]: Initialize FFTW with 12 threads.
gmt [DEBUG]: GMT_Create_Session initialized GMT structure
gmt [DEBUG]: Loading core GMT shared library: libgmt.so
gmt [DEBUG]: Shared Library # 0 (core). Path = libgmt.so
gmt [DEBUG]: Loading GMT plugins from: /panfs/ccds02/nobackup/people/sberton2/.spack_repo/opt/spack/linux-centos7-haswell/gcc-11.1.0/gmt-6.4.0-xdmacvrgjo3w4eprh25qlo4fjsxfjx5u/lib64/gmt/plugins
gmt [DEBUG]: Shared Library # 1 (supplements). Path = /panfs/ccds02/nobackup/people/sberton2/.spack_repo/opt/spack/linux-centos7-haswell/gcc-11.1.0/gmt-6.4.0-xdmacvrgjo3w4eprh25qlo4fjsxfjx5u/lib64/gmt/plugins/supplements.so
gmt [DEBUG]: Revised options: -r -bi3 -R-19/-3/-20/-4 -I0.005 -T0.75 -Gtmp.grd -C0.001 -Ll-11 -Lu11 -Vd
surface [DEBUG]: History: Process -R-19/-3/-20/-4
surface [DEBUG]: Look for file -19/-3/-20/-4 in /home/mkbarker/.gmt
surface [DEBUG]: Look for file -19/-3/-20/-4 in /home/mkbarker/.gmt/cache
surface [DEBUG]: Look for file -19/-3/-20/-4 in /home/mkbarker/.gmt/server
surface [DEBUG]: Got regular w/e/s/n for region (-19/-3/-20/-4)
surface [DEBUG]: Look for file -11 in /home/mkbarker/.gmt
surface [DEBUG]: Look for file -11 in /home/mkbarker/.gmt/cache
surface [DEBUG]: Look for file -11 in /home/mkbarker/.gmt/server
surface [DEBUG]: Look for file 11 in /home/mkbarker/.gmt
surface [DEBUG]: Look for file 11 in /home/mkbarker/.gmt/cache
surface [DEBUG]: Look for file 11 in /home/mkbarker/.gmt/server
surface [INFORMATION]: Provides 3, expects 3-column binary data
surface [INFORMATION]: Cartesian input grid
surface [DEBUG]: Chosen boundary condition for all edges: natural
surface [INFORMATION]: Cartesian input grid
surface [DEBUG]: Object ID 0 : Registered Grid Memory Reference a03cd0 as an Input resource with geometry Surface [n_objects = 1]
surface [DEBUG]: Successfully created a new Grid container
surface [INFORMATION]: Cannot suggest any n_columns,n_rows better than your current -R -I settings.
surface [INFORMATION]: Cartesian input grid
surface [DEBUG]: Chosen boundary condition for all edges: natural
surface [INFORMATION]: Cartesian input grid
surface [DEBUG]: Object ID 1 : Registered Grid Memory Reference a09610 as an Input resource with geometry Surface [n_objects = 2]
surface [DEBUG]: Successfully created a new Grid container
surface [INFORMATION]: Cartesian input grid
surface [DEBUG]: Chosen boundary condition for all edges: natural
surface [INFORMATION]: Cartesian input grid
surface [DEBUG]: Object ID 2 : Registered Grid Memory Reference a094d0 as an Input resource with geometry Surface [n_objects = 3]
surface [DEBUG]: Successfully created a new Grid container
surface [INFORMATION]: Grid domain: W: -19 E: -3 S: -20 N: -4 n_columns: 3200 n_rows: 3200 [pixel registration]
surface [INFORMATION]: Processing input table data
surface [INFORMATION]: Input 3 columns via binary records using format ddd
surface [DEBUG]: gmtapi_init_import: Passed family = Data Table and geometry = Point
surface [DEBUG]: gmtapi_init_import: Added 0 new sources
surface [DEBUG]: Object ID 3 : Registered Data Table Stream 2b0dbcd47640 as an Input resource with geometry Point [n_objects = 4]
surface [DEBUG]: gmtapi_init_import: Added stdin to registered sources
surface [DEBUG]: GMT_Init_IO: Returned first Input object ID = 3
surface [DEBUG]: GMT_Begin_IO: Mode value 1 not considered (ignored)
surface [DEBUG]: GMT_Begin_IO: Initialize record-by-record access for Input
surface [DEBUG]: gmtapi_next_io_source: Selected object 3
surface [INFORMATION]: Reading Data Table from Standard Input stream
surface [DEBUG]: GMT_Begin_IO: Input resource access is now enabled [record-by-record]
surface [DEBUG]: GMT_End_IO: Input resource access is now disabled
surface [INFORMATION]: Minimum value of your dataset x,y,z at: -18.9974994659 -19.9424991608 -0.523202717304
surface [INFORMATION]: Maximum value of your dataset x,y,z at: -11.5725002289 -12.4724998474 1.9595233202
surface [INFORMATION]: Eliminate data points that are not nearest a node.
surface [INFORMATION]: Plane fit z = 1.05921 + (8.01847e-05 * col) + (-6.01073e-06 * row)
surface [INFORMATION]: Normalize detrended data constraints by z rms = 0.391272
surface [DEBUG]: Load any data constraint limit grids
surface [DEBUG]: Object ID 4 : Registered Grid Memory Reference a09390 as an Input resource with geometry Surface [n_objects = 5]
surface [DEBUG]: Successfully duplicated a Grid
surface [DEBUG]: Object ID 5 : Registered Grid Memory Reference a09250 as an Input resource with geometry Surface [n_objects = 6]
surface [DEBUG]: Successfully duplicated a Grid
surface [INFORMATION]: Recompute data index for next iteration [stride = 640]
surface [INFORMATION]: ------------------------------------------
surface [INFORMATION]: Memory for data array          :   22.2 Mb
surface [INFORMATION]: Memory for final grid          :   39.2 Mb
surface [INFORMATION]: Memory for constraint grid     :   39.2 Mb
surface [INFORMATION]: Memory for constraint grid     :   39.2 Mb
surface [INFORMATION]: Memory for Briggs coefficients :   22.2 Mb
surface [INFORMATION]: Memory for node status         :    9.8 Mb
surface [INFORMATION]: ------------------------------------------
surface [INFORMATION]: Total memory use               :  171.8 Mb
surface [INFORMATION]: ==========================================
surface [DEBUG]: Successfully added data array to previously registered Grid container
surface [INFORMATION]: Grid Mode    Iteration   Max Change  Conv Limit  Total Iterations
surface [INFORMATION]: Set finite-difference coefficients [stride = 640]
surface [INFORMATION]: Determine nearest point and set Briggs coefficients [stride = 640]
surface [DEBUG]: Starting iterations, mode = data Max iterations = 320000 [stride = 640]
surface [DEBUG]: Apply all boundary conditions [stride = 640]
surface [DEBUG]: Iteration 0
surface [DEBUG]:  640   D          1    -0.391271693085 1.5625e-06           1
surface [INFORMATION]:  640 D          1    -0.391271693085 1.5625e-06           1
surface [INFORMATION]: Recompute data index for next iteration [stride = 128]
surface [INFORMATION]: Expand grid by factor of 5 when going from stride = 640 to 128
surface [INFORMATION]: Fill in expanded grid by bilinear interpolation [stride = 128]
surface [DEBUG]: Starting iterations, mode = node Max iterations = 64000 [stride = 128]
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 0
surface [DEBUG]:  128   I          1    0.865824664794  7.8125e-06           2
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 1
surface [DEBUG]:  128   I          2    0.262993189222  7.8125e-06           3
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 2
surface [DEBUG]:  128   I          3    0.169091881958  7.8125e-06           4
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 3
surface [DEBUG]:  128   I          4    0.113299972914  7.8125e-06           5
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 4
surface [DEBUG]:  128   I          5    0.0890424228621 7.8125e-06           6
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 5
surface [DEBUG]:  128   I          6    0.0785462813039 7.8125e-06           7
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 6
surface [DEBUG]:  128   I          7    0.0626517169683 7.8125e-06           8
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 7
surface [DEBUG]:  128   I          8    0.0524155797165 7.8125e-06           9
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 8
surface [DEBUG]:  128   I          9    0.0449849287096 7.8125e-06          10
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 9
surface [DEBUG]:  128   I         10    0.0384314512745 7.8125e-06          11
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 10
surface [DEBUG]:  128   I         11    0.0335729883286 7.8125e-06          12
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 11
surface [DEBUG]:  128   I         12    0.029087566783  7.8125e-06          13
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 12
surface [DEBUG]:  128   I         13    0.0253638440422 7.8125e-06          14
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 13
surface [DEBUG]:  128   I         14    0.0220909274451 7.8125e-06          15
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 14
surface [DEBUG]:  128   I         15    0.0192735012827 7.8125e-06          16
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 15
surface [DEBUG]:  128   I         16    0.0167797024012 7.8125e-06          17
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 16
surface [DEBUG]:  128   I         17    0.0145976267219 7.8125e-06          18
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 17
surface [DEBUG]:  128   I         18    0.0127165333314 7.8125e-06          19
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 18
surface [DEBUG]:  128   I         19    0.0115074168671 7.8125e-06          20
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 19
surface [DEBUG]:  128   I         20    0.0104645412136 7.8125e-06          21
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 20
surface [DEBUG]:  128   I         21    0.00949032088252    7.8125e-06          22
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 21
surface [DEBUG]:  128   I         22    0.00868005176322    7.8125e-06          23
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 22
surface [DEBUG]:  128   I         23    0.00793927054783    7.8125e-06          24
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 23
surface [DEBUG]:  128   I         24    0.00724042654355    7.8125e-06          25
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 24
surface [DEBUG]:  128   I         25    0.00658572160193    7.8125e-06          26
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 25
surface [DEBUG]:  128   I         26    0.00597606176751    7.8125e-06          27
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 26
surface [DEBUG]:  128   I         27    0.00541238165384    7.8125e-06          28
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 27
surface [DEBUG]:  128   I         28    0.00489469146411    7.8125e-06          29
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 28
surface [DEBUG]:  128   I         29    0.00442171171649    7.8125e-06          30
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 29
surface [DEBUG]:  128   I         30    0.00399069162703    7.8125e-06          31
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 30
surface [DEBUG]:  128   I         31    0.00359855594991    7.8125e-06          32
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 31
surface [DEBUG]:  128   I         32    0.00324241309697    7.8125e-06          33
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 32
surface [DEBUG]:  128   I         33    0.00291905518072    7.8125e-06          34
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 33
surface [DEBUG]:  128   I         34    0.00262604567592    7.8125e-06          35
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 34
surface [DEBUG]:  128   I         35    0.00236053584788    7.8125e-06          36
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 35
surface [DEBUG]:  128   I         36    0.00212018304085    7.8125e-06          37
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 36
surface [DEBUG]:  128   I         37    0.00190292008559    7.8125e-06          38
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 37
surface [DEBUG]:  128   I         38    0.00171011013022    7.8125e-06          39
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 38
surface [DEBUG]:  128   I         39    0.0015365515811 7.8125e-06          40
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 39
surface [DEBUG]:  128   I         40    0.00137952426391    7.8125e-06          41
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 40
surface [DEBUG]:  128   I         41    0.00123775685935    7.8125e-06          42
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 41
surface [DEBUG]:  128   I         42    0.00110975969959    7.8125e-06          43
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 42
surface [DEBUG]:  128   I         43    0.000994412472762   7.8125e-06          44
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 43
surface [DEBUG]:  128   I         44    0.000890390802923   7.8125e-06          45
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 44
surface [DEBUG]:  128   I         45    0.000796835580256   7.8125e-06          46
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 45
surface [DEBUG]:  128   I         46    0.000712761175208   7.8125e-06          47
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 46
surface [DEBUG]:  128   I         47    0.000637145226687   7.8125e-06          48
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 47
surface [DEBUG]:  128   I         48    0.000569408192678   7.8125e-06          49
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 48
surface [DEBUG]:  128   I         49    0.00050852771209    7.8125e-06          50
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 49
surface [DEBUG]:  128   I         50    0.000454034437516   7.8125e-06          51
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 50
surface [DEBUG]:  128   I         51    0.000405236591693   7.8125e-06          52
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 51
surface [DEBUG]:  128   I         52    0.000361473006968   7.8125e-06          53
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 52
surface [DEBUG]:  128   I         53    0.000322462074897   7.8125e-06          54
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 53
surface [DEBUG]:  128   I         54    0.000287665066284   7.8125e-06          55
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 54
surface [DEBUG]:  128   I         55    0.000256451423089   7.8125e-06          56
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 55
surface [DEBUG]:  128   I         56    0.00022860891866    7.8125e-06          57
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 56
surface [DEBUG]:  128   I         57    0.000203733506098   7.8125e-06          58
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 57
surface [DEBUG]:  128   I         58    0.000181466032605   7.8125e-06          59
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 58
surface [DEBUG]:  128   I         59    0.000161651409474   7.8125e-06          60
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 59
surface [DEBUG]:  128   I         60    0.000143799882887   7.8125e-06          61
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 60
surface [DEBUG]:  128   I         61    0.000128058378991   7.8125e-06          62
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 61
surface [DEBUG]:  128   I         62    0.000113953469095   7.8125e-06          63
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 62
surface [DEBUG]:  128   I         63    0.000101313739365   7.8125e-06          64
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 63
surface [DEBUG]:  128   I         64    9.00085887834e-05   7.8125e-06          65
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 64
surface [DEBUG]:  128   I         65    8.00829114497e-05   7.8125e-06          66
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 65
surface [DEBUG]:  128   I         66    7.11489855927e-05   7.8125e-06          67
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 66
surface [DEBUG]:  128   I         67    6.31619171131e-05   7.8125e-06          68
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 67
surface [DEBUG]:  128   I         68    5.61910878013e-05   7.8125e-06          69
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 68
surface [DEBUG]:  128   I         69    4.98773448591e-05   7.8125e-06          70
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 69
surface [DEBUG]:  128   I         70    4.41839567502e-05   7.8125e-06          71
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 70
surface [DEBUG]:  128   I         71    3.9306825001e-05    7.8125e-06          72
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 71
surface [DEBUG]:  128   I         72    3.48296588681e-05   7.8125e-06          73
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 72
surface [DEBUG]:  128   I         73    3.09728475687e-05   7.8125e-06          74
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 73
surface [DEBUG]:  128   I         74    2.73405067679e-05   7.8125e-06          75
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 74
surface [DEBUG]:  128   I         75    2.43856587459e-05   7.8125e-06          76
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 75
surface [DEBUG]:  128   I         76    2.16022245593e-05   7.8125e-06          77
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 76
surface [DEBUG]:  128   I         77    1.91003988169e-05   7.8125e-06          78
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 77
surface [DEBUG]:  128   I         78    1.6949563309e-05    7.8125e-06          79
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 78
surface [DEBUG]:  128   I         79    1.49374913824e-05   7.8125e-06          80
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 79
surface [DEBUG]:  128   I         80    1.3386604299e-05    7.8125e-06          81
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 80
surface [DEBUG]:  128   I         81    1.17704167067e-05   7.8125e-06          82
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 81
surface [DEBUG]:  128   I         82    1.04154311497e-05   7.8125e-06          83
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 82
surface [DEBUG]:  128   I         83    9.19920917383e-06   7.8125e-06          84
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 83
surface [DEBUG]:  128   I         84    8.21153897856e-06   7.8125e-06          85
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 84
surface [DEBUG]:  128   I         85    7.28100672837e-06   7.8125e-06          86
surface [INFORMATION]:  128 I         85    7.28100672837e-06   7.8125e-06          86
surface [INFORMATION]: Determine nearest point and set Briggs coefficients [stride = 128]
surface [DEBUG]: Starting iterations, mode = data Max iterations = 64000 [stride = 128]
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 0
surface [DEBUG]:  128   D          1    0.118818296963  7.8125e-06          87
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 1
surface [DEBUG]:  128   D          2    0.0476370290872 7.8125e-06          88
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 2
surface [DEBUG]:  128   D          3    0.0190988285752 7.8125e-06          89
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 3
surface [DEBUG]:  128   D          4    0.00765714654773    7.8125e-06          90
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 4
surface [DEBUG]:  128   D          5    0.00306991729833    7.8125e-06          91
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 5
surface [DEBUG]:  128   D          6    0.00123080115504    7.8125e-06          92
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 6
surface [DEBUG]:  128   D          7    0.000493485158766   7.8125e-06          93
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 7
surface [DEBUG]:  128   D          8    0.00019784993928    7.8125e-06          94
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 8
surface [DEBUG]:  128   D          9    7.93376649789e-05   7.8125e-06          95
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 9
surface [DEBUG]:  128   D         10    3.18119403741e-05   7.8125e-06          96
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 10
surface [DEBUG]:  128   D         11    1.2752511801e-05    7.8125e-06          97
surface [DEBUG]: Apply all boundary conditions [stride = 128]
surface [DEBUG]: Iteration 11
surface [DEBUG]:  128   D         12    5.08637224413e-06   7.8125e-06          98
surface [INFORMATION]:  128 D         12    5.08637224413e-06   7.8125e-06          98
surface [INFORMATION]: Recompute data index for next iteration [stride = 64]
surface [INFORMATION]: Expand grid by factor of 2 when going from stride = 128 to 64
surface [INFORMATION]: Fill in expanded grid by bilinear interpolation [stride = 64]
surface [DEBUG]: Starting iterations, mode = node Max iterations = 32000 [stride = 64]
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 0
surface [DEBUG]:   64   I          1    0.1577321464    1.5625e-05          99
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 1
surface [DEBUG]:   64   I          2    0.0693129781223 1.5625e-05         100
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 2
surface [DEBUG]:   64   I          3    0.0255644839367 1.5625e-05         101
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 3
surface [DEBUG]:   64   I          4    0.0121759685775 1.5625e-05         102
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 4
surface [DEBUG]:   64   I          5    0.00788582106108    1.5625e-05         103
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 5
surface [DEBUG]:   64   I          6    0.00492819266576    1.5625e-05         104
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 6
surface [DEBUG]:   64   I          7    0.00312282337771    1.5625e-05         105
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 7
surface [DEBUG]:   64   I          8    0.00183320490531    1.5625e-05         106
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 8
surface [DEBUG]:   64   I          9    0.00123509433314    1.5625e-05         107
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 9
surface [DEBUG]:   64   I         10    0.000835588973488   1.5625e-05         108
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 10
surface [DEBUG]:   64   I         11    0.000600285578507   1.5625e-05         109
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 11
surface [DEBUG]:   64   I         12    0.000413092166375   1.5625e-05         110
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 12
surface [DEBUG]:   64   I         13    0.00027395004634    1.5625e-05         111
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 13
surface [DEBUG]:   64   I         14    0.000182971770509   1.5625e-05         112
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 14
surface [DEBUG]:   64   I         15    0.000120712199293   1.5625e-05         113
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 15
surface [DEBUG]:   64   I         16    8.74526860573e-05   1.5625e-05         114
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 16
surface [DEBUG]:   64   I         17    5.85990440624e-05   1.5625e-05         115
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 17
surface [DEBUG]:   64   I         18    4.23907435601e-05   1.5625e-05         116
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 18
surface [DEBUG]:   64   I         19    2.84391567402e-05   1.5625e-05         117
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 19
surface [DEBUG]:   64   I         20    1.97521763316e-05   1.5625e-05         118
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 20
surface [DEBUG]:   64   I         21    1.34892102742e-05   1.5625e-05         119
surface [INFORMATION]:   64 I         21    1.34892102742e-05   1.5625e-05         119
surface [INFORMATION]: Determine nearest point and set Briggs coefficients [stride = 64]
surface [DEBUG]: Starting iterations, mode = data Max iterations = 32000 [stride = 64]
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 0
surface [DEBUG]:   64   D          1    0.0642454035003 1.5625e-05         120
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 1
surface [DEBUG]:   64   D          2    0.0256981623016 1.5625e-05         121
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 2
surface [DEBUG]:   64   D          3    0.0102792672842 1.5625e-05         122
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 3
surface [DEBUG]:   64   D          4    0.00411170781516    1.5625e-05         123
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 4
surface [DEBUG]:   64   D          5    0.00164467895955    1.5625e-05         124
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 5
surface [DEBUG]:   64   D          6    0.000657865955255   1.5625e-05         125
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 6
surface [DEBUG]:   64   D          7    0.000263148745643   1.5625e-05         126
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 7
surface [DEBUG]:   64   D          8    0.000105260399742   1.5625e-05         127
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 8
surface [DEBUG]:   64   D          9    4.21065234374e-05   1.5625e-05         128
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 9
surface [DEBUG]:   64   D         10    1.68435108597e-05   1.5625e-05         129
surface [DEBUG]: Apply all boundary conditions [stride = 64]
surface [DEBUG]: Iteration 10
surface [DEBUG]:   64   D         11    6.72997280825e-06   1.5625e-05         130
surface [INFORMATION]:   64 D         11    6.72997280825e-06   1.5625e-05         130
surface [INFORMATION]: Recompute data index for next iteration [stride = 32]
surface [INFORMATION]: Expand grid by factor of 2 when going from stride = 64 to 32
surface [INFORMATION]: Fill in expanded grid by bilinear interpolation [stride = 32]
surface [DEBUG]: Starting iterations, mode = node Max iterations = 16000 [stride = 32]
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 0
surface [DEBUG]:   32   I          1    0.0936112617264 3.125e-05          131
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 1
surface [DEBUG]:   32   I          2    0.0368228303788 3.125e-05          132
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 2
surface [DEBUG]:   32   I          3    0.0119288793598 3.125e-05          133
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 3
surface [DEBUG]:   32   I          4    0.00720486123364    3.125e-05          134
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 4
surface [DEBUG]:   32   I          5    0.00405907401373    3.125e-05          135
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 5
surface [DEBUG]:   32   I          6    0.00263565556408    3.125e-05          136
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 6
surface [DEBUG]:   32   I          7    0.00163061288271    3.125e-05          137
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 7
surface [DEBUG]:   32   I          8    0.00108315178291    3.125e-05          138
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 8
surface [DEBUG]:   32   I          9    0.000738765062187   3.125e-05          139
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 9
surface [DEBUG]:   32   I         10    0.000485443982238   3.125e-05          140
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 10
surface [DEBUG]:   32   I         11    0.000309095877007   3.125e-05          141
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 11
surface [DEBUG]:   32   I         12    0.000220479005317   3.125e-05          142
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 12
surface [DEBUG]:   32   I         13    0.000133006933171   3.125e-05          143
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 13
surface [DEBUG]:   32   I         14    9.71794009035e-05   3.125e-05          144
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 14
surface [DEBUG]:   32   I         15    6.14253317079e-05   3.125e-05          145
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 15
surface [DEBUG]:   32   I         16    3.84273087735e-05   3.125e-05          146
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 16
surface [DEBUG]:   32   I         17    2.5906446376e-05    3.125e-05          147
surface [INFORMATION]:   32 I         17    2.5906446376e-05    3.125e-05          147
surface [INFORMATION]: Determine nearest point and set Briggs coefficients [stride = 32]
surface [DEBUG]: Starting iterations, mode = data Max iterations = 16000 [stride = 32]
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 0
surface [DEBUG]:   32   D          1    0.060841571938  3.125e-05          148
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 1
surface [DEBUG]:   32   D          2    0.0247257232407 3.125e-05          149
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 2
surface [DEBUG]:   32   D          3    0.0100883622323 3.125e-05          150
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 3
surface [DEBUG]:   32   D          4    0.00457551467217    3.125e-05          151
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 4
surface [DEBUG]:   32   D          5    0.00225391627865    3.125e-05          152
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 5
surface [DEBUG]:   32   D          6    0.00103353375153    3.125e-05          153
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 6
surface [DEBUG]:   32   D          7    0.000449056989959   3.125e-05          154
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 7
surface [DEBUG]:   32   D          8    0.000186447801501   3.125e-05          155
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 8
surface [DEBUG]:   32   D          9    7.42512304588e-05   3.125e-05          156
surface [DEBUG]: Apply all boundary conditions [stride = 32]
surface [DEBUG]: Iteration 9
surface [DEBUG]:   32   D         10    2.83904991878e-05   3.125e-05          157
surface [INFORMATION]:   32 D         10    2.83904991878e-05   3.125e-05          157
surface [INFORMATION]: Recompute data index for next iteration [stride = 16]
surface [INFORMATION]: Expand grid by factor of 2 when going from stride = 32 to 16
surface [INFORMATION]: Fill in expanded grid by bilinear interpolation [stride = 16]
surface [DEBUG]: Starting iterations, mode = node Max iterations = 8000 [stride = 16]
surface [DEBUG]: Apply all boundary conditions [stride = 16]
surface [DEBUG]: Iteration 0
surface [DEBUG]:   16   I          1    0.051711525956  6.25e-05           158
surface [DEBUG]: Apply all boundary conditions [stride = 16]
surface [DEBUG]: Iteration 1
surface [DEBUG]:   16   I          2    0.0206482596298 6.25e-05           159
surface [DEBUG]: Apply all boundary conditions [stride = 16]
surface [DEBUG]: Iteration 2
surface [DEBUG]:   16   I          3    0.0065597412142 6.25e-05           160
surface [DEBUG]: Apply all boundary conditions [stride = 16]
surface [DEBUG]: Iteration 3
surface [DEBUG]:   16   I          4    0.00376810872095    6.25e-05           161
surface [DEBUG]: Apply all boundary conditions [stride = 16]
surface [DEBUG]: Iteration 4
surface [DEBUG]:   16   I          5    0.00212375212178    6.25e-05           162
surface [DEBUG]: Apply all boundary conditions [stride = 16]
surface [DEBUG]: Iteration 5
surface [DEBUG]:   16   I          6    0.00134731886923    6.25e-05           163
surface [DEBUG]: Apply all boundary conditions [stride = 16]
surface [DEBUG]: Iteration 6
surface [DEBUG]:   16   I          7    0.000961566316847   6.25e-05           164
surface [DEBUG]: Apply all boundary conditions [stride = 16]
surface [DEBUG]: Iteration 7
surface [DEBUG]:   16   I          8    0.000578417622248   6.25e-05           165
surface [DEBUG]: Apply all boundary conditions [stride = 16]
surface [DEBUG]: Iteration 8
surface [DEBUG]:   16   I          9    0.000399841137156   6.25e-05           166
surface [DEBUG]: Apply all boundary conditions [stride = 16]
surface [DEBUG]: Iteration 9
surface [DEBUG]:   16   I         10    0.000283485833722   6.25e-05           167
surface [DEBUG]: Apply all boundary conditions [stride = 16]
surface [DEBUG]: Iteration 10
surface [DEBUG]:   16   I         11    0.000184159678593   6.25e-05           168
surface [DEBUG]: Apply all boundary conditions [stride = 16]
surface [DEBUG]: Iteration 11
surface [DEBUG]:   16   I         12    0.000131874377472   6.25e-05           169
surface [DEBUG]: Apply all boundary conditions [stride = 16]
surface [DEBUG]: Iteration 12
surface [DEBUG]:   16   I         13    8.07583635873e-05   6.25e-05           170
surface [DEBUG]: Apply all boundary conditions [stride = 16]
surface [DEBUG]: Iteration 13
surface [DEBUG]:   16   I         14    6.23272949854e-05   6.25e-05           171
surface [INFORMATION]:   16 I         14    6.23272949854e-05   6.25e-05           171
surface [INFORMATION]: Determine nearest point and set Briggs coefficients [stride = 16]
surface [DEBUG]: Starting iterations, mode = data Max iterations = 8000 [stride = 16]
surface [DEBUG]: Apply all boundary conditions [stride = 16]
surface [DEBUG]: Iteration 0
surface [DEBUG]:   16   D          1    0.0554780186788 6.25e-05           172
surface [DEBUG]: Apply all boundary conditions [stride = 16]
surface [DEBUG]: Iteration 1
surface [DEBUG]:   16   D          2    0.0201522055679 6.25e-05           173
surface [DEBUG]: Apply all boundary conditions [stride = 16]
surface [DEBUG]: Iteration 2
surface [DEBUG]:   16   D          3    0.00693798993274    6.25e-05           174
surface [DEBUG]: Apply all boundary conditions [stride = 16]
surface [DEBUG]: Iteration 3
surface [DEBUG]:   16   D          4    0.00322166908744    6.25e-05           175
surface [DEBUG]: Apply all boundary conditions [stride = 16]
surface [DEBUG]: Iteration 4
surface [DEBUG]:   16   D          5    0.00168155630072    6.25e-05           176
surface [DEBUG]: Apply all boundary conditions [stride = 16]
surface [DEBUG]: Iteration 5
surface [DEBUG]:   16   D          6    0.000789646560534   6.25e-05           177
surface [DEBUG]: Apply all boundary conditions [stride = 16]
surface [DEBUG]: Iteration 6
surface [DEBUG]:   16   D          7    0.000416565901638   6.25e-05           178
surface [DEBUG]: Apply all boundary conditions [stride = 16]
surface [DEBUG]: Iteration 7
surface [DEBUG]:   16   D          8    0.000210901259038   6.25e-05           179
surface [DEBUG]: Apply all boundary conditions [stride = 16]
surface [DEBUG]: Iteration 8
surface [DEBUG]:   16   D          9    9.86207673503e-05   6.25e-05           180
surface [DEBUG]: Apply all boundary conditions [stride = 16]
surface [DEBUG]: Iteration 9
surface [DEBUG]:   16   D         10    4.48079757044e-05   6.25e-05           181
surface [INFORMATION]:   16 D         10    4.48079757044e-05   6.25e-05           181
surface [INFORMATION]: Recompute data index for next iteration [stride = 8]
surface [INFORMATION]: Expand grid by factor of 2 when going from stride = 16 to 8
surface [INFORMATION]: Fill in expanded grid by bilinear interpolation [stride = 8]
surface [DEBUG]: Starting iterations, mode = node Max iterations = 4000 [stride = 8]
surface [DEBUG]: Apply all boundary conditions [stride = 8]
surface [DEBUG]: Iteration 0
surface [DEBUG]:    8   I          1    0.022222551844  0.000125           182
surface [DEBUG]: Apply all boundary conditions [stride = 8]
surface [DEBUG]: Iteration 1
surface [DEBUG]:    8   I          2    0.00889093465471    0.000125           183
surface [DEBUG]: Apply all boundary conditions [stride = 8]
surface [DEBUG]: Iteration 2
surface [DEBUG]:    8   I          3    0.0033690062958 0.000125           184
surface [DEBUG]: Apply all boundary conditions [stride = 8]
surface [DEBUG]: Iteration 3
surface [DEBUG]:    8   I          4    0.00182076184234    0.000125           185
surface [DEBUG]: Apply all boundary conditions [stride = 8]
surface [DEBUG]: Iteration 4
surface [DEBUG]:    8   I          5    0.000990429141726   0.000125           186
surface [DEBUG]: Apply all boundary conditions [stride = 8]
surface [DEBUG]: Iteration 5
surface [DEBUG]:    8   I          6    0.000636256527575   0.000125           187
surface [DEBUG]: Apply all boundary conditions [stride = 8]
surface [DEBUG]: Iteration 6
surface [DEBUG]:    8   I          7    0.000422173911157   0.000125           188
surface [DEBUG]: Apply all boundary conditions [stride = 8]
surface [DEBUG]: Iteration 7
surface [DEBUG]:    8   I          8    0.000274260096217   0.000125           189
surface [DEBUG]: Apply all boundary conditions [stride = 8]
surface [DEBUG]: Iteration 8
surface [DEBUG]:    8   I          9    0.000189812253883   0.000125           190
surface [DEBUG]: Apply all boundary conditions [stride = 8]
surface [DEBUG]: Iteration 9
surface [DEBUG]:    8   I         10    0.000126219761541   0.000125           191
surface [DEBUG]: Apply all boundary conditions [stride = 8]
surface [DEBUG]: Iteration 10
surface [DEBUG]:    8   I         11    8.64639955415e-05   0.000125           192
surface [INFORMATION]:    8 I         11    8.64639955415e-05   0.000125           192
surface [INFORMATION]: Determine nearest point and set Briggs coefficients [stride = 8]
surface [DEBUG]: Starting iterations, mode = data Max iterations = 4000 [stride = 8]
surface [DEBUG]: Apply all boundary conditions [stride = 8]
surface [DEBUG]: Iteration 0
surface [DEBUG]:    8   D          1    0.0222586199914 0.000125           193
surface [DEBUG]: Apply all boundary conditions [stride = 8]
surface [DEBUG]: Iteration 1
surface [DEBUG]:    8   D          2    0.00968459238908    0.000125           194
surface [DEBUG]: Apply all boundary conditions [stride = 8]
surface [DEBUG]: Iteration 2
surface [DEBUG]:    8   D          3    0.004279767711  0.000125           195
surface [DEBUG]: Apply all boundary conditions [stride = 8]
surface [DEBUG]: Iteration 3
surface [DEBUG]:    8   D          4    0.00208018263738    0.000125           196
surface [DEBUG]: Apply all boundary conditions [stride = 8]
surface [DEBUG]: Iteration 4
surface [DEBUG]:    8   D          5    0.00107431067263    0.000125           197
surface [DEBUG]: Apply all boundary conditions [stride = 8]
surface [DEBUG]: Iteration 5
surface [DEBUG]:    8   D          6    0.00056394335635    0.000125           198
surface [DEBUG]: Apply all boundary conditions [stride = 8]
surface [DEBUG]: Iteration 6
surface [DEBUG]:    8   D          7    0.000371539488524   0.000125           199
surface [DEBUG]: Apply all boundary conditions [stride = 8]
surface [DEBUG]: Iteration 7
surface [DEBUG]:    8   D          8    0.000298441112514   0.000125           200
surface [DEBUG]: Apply all boundary conditions [stride = 8]
surface [DEBUG]: Iteration 8
surface [DEBUG]:    8   D          9    0.000207440380621   0.000125           201
surface [DEBUG]: Apply all boundary conditions [stride = 8]
surface [DEBUG]: Iteration 9
surface [DEBUG]:    8   D         10    0.000122391005128   0.000125           202
surface [INFORMATION]:    8 D         10    0.000122391005128   0.000125           202
surface [INFORMATION]: Recompute data index for next iteration [stride = 4]
surface [INFORMATION]: Expand grid by factor of 2 when going from stride = 8 to 4
surface [INFORMATION]: Fill in expanded grid by bilinear interpolation [stride = 4]
surface [DEBUG]: Starting iterations, mode = node Max iterations = 2000 [stride = 4]
surface [DEBUG]: Apply all boundary conditions [stride = 4]
surface [DEBUG]: Iteration 0
surface [DEBUG]:    4   I          1    0.0112391215642 0.00025        203
surface [DEBUG]: Apply all boundary conditions [stride = 4]
surface [DEBUG]: Iteration 1
surface [DEBUG]:    4   I          2    0.00393448319378    0.00025        204
surface [DEBUG]: Apply all boundary conditions [stride = 4]
surface [DEBUG]: Iteration 2
surface [DEBUG]:    4   I          3    0.00159852176393    0.00025        205
surface [DEBUG]: Apply all boundary conditions [stride = 4]
surface [DEBUG]: Iteration 3
surface [DEBUG]:    4   I          4    0.000916129406585   0.00025        206
surface [DEBUG]: Apply all boundary conditions [stride = 4]
surface [DEBUG]: Iteration 4
surface [DEBUG]:    4   I          5    0.000496802189483   0.00025        207
surface [DEBUG]: Apply all boundary conditions [stride = 4]
surface [DEBUG]: Iteration 5
surface [DEBUG]:    4   I          6    0.000298278439601   0.00025        208
surface [DEBUG]: Apply all boundary conditions [stride = 4]
surface [DEBUG]: Iteration 6
surface [DEBUG]:    4   I          7    0.000200962315757   0.00025        209
surface [INFORMATION]:    4 I          7    0.000200962315757   0.00025        209
surface [INFORMATION]: Determine nearest point and set Briggs coefficients [stride = 4]
surface [DEBUG]: Starting iterations, mode = data Max iterations = 2000 [stride = 4]
surface [DEBUG]: Apply all boundary conditions [stride = 4]
surface [DEBUG]: Iteration 0
surface [DEBUG]:    4   D          1    0.0107981694766 0.00025        210
surface [DEBUG]: Apply all boundary conditions [stride = 4]
surface [DEBUG]: Iteration 1
surface [DEBUG]:    4   D          2    0.00725090217361    0.00025        211
surface [DEBUG]: Apply all boundary conditions [stride = 4]
surface [DEBUG]: Iteration 2
surface [DEBUG]:    4   D          3    0.00375568529915    0.00025        212
surface [DEBUG]: Apply all boundary conditions [stride = 4]
surface [DEBUG]: Iteration 3
surface [DEBUG]:    4   D          4    0.00180171245955    0.00025        213
surface [DEBUG]: Apply all boundary conditions [stride = 4]
surface [DEBUG]: Iteration 4
surface [DEBUG]:    4   D          5    0.00120869405167    0.00025        214
surface [DEBUG]: Apply all boundary conditions [stride = 4]
surface [DEBUG]: Iteration 5
surface [DEBUG]:    4   D          6    0.000761121303559   0.00025        215
surface [DEBUG]: Apply all boundary conditions [stride = 4]
surface [DEBUG]: Iteration 6
surface [DEBUG]:    4   D          7    0.000396464896788   0.00025        216
surface [DEBUG]: Apply all boundary conditions [stride = 4]
surface [DEBUG]: Iteration 7
surface [DEBUG]:    4   D          8    0.000270606328686   0.00025        217
surface [DEBUG]: Apply all boundary conditions [stride = 4]
surface [DEBUG]: Iteration 8
surface [DEBUG]:    4   D          9    0.000198040117989   0.00025        218
surface [INFORMATION]:    4 D          9    0.000198040117989   0.00025        218
surface [INFORMATION]: Recompute data index for next iteration [stride = 2]
surface [INFORMATION]: Expand grid by factor of 2 when going from stride = 4 to 2
surface [INFORMATION]: Fill in expanded grid by bilinear interpolation [stride = 2]
surface [DEBUG]: Starting iterations, mode = node Max iterations = 1000 [stride = 2]
surface [DEBUG]: Apply all boundary conditions [stride = 2]
surface [DEBUG]: Iteration 0
surface [DEBUG]:    2   I          1    0.00422898338676    0.0005         219
surface [DEBUG]: Apply all boundary conditions [stride = 2]
surface [DEBUG]: Iteration 1
surface [DEBUG]:    2   I          2    0.00160943307082    0.0005         220
surface [DEBUG]: Apply all boundary conditions [stride = 2]
surface [DEBUG]: Iteration 2
surface [DEBUG]:    2   I          3    0.000568059329043   0.0005         221
surface [DEBUG]: Apply all boundary conditions [stride = 2]
surface [DEBUG]: Iteration 3
surface [DEBUG]:    2   I          4    0.000384485314389   0.0005         222
surface [INFORMATION]:    2 I          4    0.000384485314389   0.0005         222
surface [INFORMATION]: Determine nearest point and set Briggs coefficients [stride = 2]
surface [DEBUG]: Starting iterations, mode = data Max iterations = 1000 [stride = 2]
surface [DEBUG]: Apply all boundary conditions [stride = 2]
surface [DEBUG]: Iteration 0
surface [DEBUG]:    2   D          1    0.00645081760696    0.0005         223
surface [DEBUG]: Apply all boundary conditions [stride = 2]
surface [DEBUG]: Iteration 1
surface [DEBUG]:    2   D          2    0.0045083787558 0.0005         224
surface [DEBUG]: Apply all boundary conditions [stride = 2]
surface [DEBUG]: Iteration 2
surface [DEBUG]:    2   D          3    0.00235320178449    0.0005         225
surface [DEBUG]: Apply all boundary conditions [stride = 2]
surface [DEBUG]: Iteration 3
surface [DEBUG]:    2   D          4    0.00123288814525    0.0005         226
surface [DEBUG]: Apply all boundary conditions [stride = 2]
surface [DEBUG]: Iteration 4
surface [DEBUG]:    2   D          5    0.000803301733566   0.0005         227
surface [DEBUG]: Apply all boundary conditions [stride = 2]
surface [DEBUG]: Iteration 5
surface [DEBUG]:    2   D          6    0.000539614835549   0.0005         228
surface [DEBUG]: Apply all boundary conditions [stride = 2]
surface [DEBUG]: Iteration 6
surface [DEBUG]:    2   D          7    0.000423879886948   0.0005         229
surface [INFORMATION]:    2 D          7    0.000423879886948   0.0005         229
surface [INFORMATION]: Recompute data index for next iteration [stride = 1]
surface [INFORMATION]: Expand grid by factor of 2 when going from stride = 2 to 1
surface [INFORMATION]: Fill in expanded grid by bilinear interpolation [stride = 1]
surface [DEBUG]: Starting iterations, mode = node Max iterations = 500 [stride = 1]
surface [DEBUG]: Apply all boundary conditions [stride = 1]
surface [DEBUG]: Iteration 0
surface [DEBUG]:    1   I          1    0.00075495244612    0.001          230
surface [INFORMATION]:    1 I          1    0.00075495244612    0.001          230
surface [INFORMATION]: Determine nearest point and set Briggs coefficients [stride = 1]
surface [DEBUG]: Starting iterations, mode = data Max iterations = 500 [stride = 1]
surface [DEBUG]: Apply all boundary conditions [stride = 1]
surface [DEBUG]: Iteration 0
surface [DEBUG]:    1   D          1    0.00405275567934    0.001          231
surface [DEBUG]: Apply all boundary conditions [stride = 1]
surface [DEBUG]: Iteration 1
surface [DEBUG]:    1   D          2    0.00252914431234    0.001          232
surface [DEBUG]: Apply all boundary conditions [stride = 1]
surface [DEBUG]: Iteration 2
surface [DEBUG]:    1   D          3    0.00137233610634    0.001          233
surface [DEBUG]: Apply all boundary conditions [stride = 1]
surface [DEBUG]: Iteration 3
surface [DEBUG]:    1   D          4    0.000891460098754   0.001          234
surface [INFORMATION]:    1 D          4    0.000891460098754   0.001          234
surface [INFORMATION]: Compute rms misfit and curvature.
surface [DEBUG]: Apply all boundary conditions [stride = 1]
surface [INFORMATION]: Fit info: N data points  N nodes mean error  rms error   curvature
surface [INFORMATION]:    972012    10246401    0   0   711.859772643
surface [DEBUG]: GMT_Destroy_Data: freed memory for a Grid for object 4
surface [DEBUG]: gmtlib_unregister_io: Unregistering object no 4 [n_objects = 5]
surface [DEBUG]: gmtlib_unregister_io: Object no 4 has non-NULL resource pointer
surface [DEBUG]: GMT_Destroy_Data: freed memory for a Grid for object 5
surface [DEBUG]: gmtlib_unregister_io: Unregistering object no 5 [n_objects = 4]
surface [DEBUG]: gmtlib_unregister_io: Object no 5 has non-NULL resource pointer
surface [INFORMATION]: Increase pad by -1 0 -1 0
surface [DEBUG]: Prepare final output grid [stride = 1]
surface [DEBUG]: Load any data constraint limit grids
surface [DEBUG]: Object ID 6 : Registered Grid Memory Reference a09250 as an Input resource with geometry Surface [n_objects = 5]
surface [DEBUG]: Successfully duplicated a Grid
ERROR: Caught signal number 11 (Segmentation fault) at
/panfs/ccds02/nobackup/people/sberton2/.spack_repo/opt/spack/linux-centos7-haswell/gcc-11.1.0/gmt-6.4.0-xdmacvrgjo3w4eprh25qlo4fjsxfjx5u/lib64/libgmt.so.6(+0x381d48)[0x2b0db896bd48]
[0x2b0dc5f81000]
Stack backtrace:
/panfs/ccds02/nobackup/people/sberton2/.spack_repo/opt/spack/linux-centos7-haswell/gcc-11.1.0/gmt-6.4.0-xdmacvrgjo3w4eprh25qlo4fjsxfjx5u/lib64/libgmt.so.6(sig_handler_unix+0xc5)[0x2b0db8719f55]
/lib64/libpthread.so.0(+0xf630)[0x2b0dbbe0b630]
/panfs/ccds02/nobackup/people/sberton2/.spack_repo/opt/spack/linux-centos7-haswell/gcc-11.1.0/gmt-6.4.0-xdmacvrgjo3w4eprh25qlo4fjsxfjx5u/lib64/libgmt.so.6(+0x381d48)[0x2b0db896bd48]
/panfs/ccds02/nobackup/people/sberton2/.spack_repo/opt/spack/linux-centos7-haswell/gcc-11.1.0/gmt-6.4.0-xdmacvrgjo3w4eprh25qlo4fjsxfjx5u/lib64/libgmt.so.6(+0x3828c7)[0x2b0db896c8c7]
/panfs/ccds02/nobackup/people/sberton2/.spack_repo/opt/spack/linux-centos7-haswell/gcc-11.1.0/gmt-6.4.0-xdmacvrgjo3w4eprh25qlo4fjsxfjx5u/lib64/libgmt.so.6(GMT_surface+0x31b4)[0x2b0db8970354]
/panfs/ccds02/nobackup/people/sberton2/.spack_repo/opt/spack/linux-centos7-haswell/gcc-11.1.0/gmt-6.4.0-xdmacvrgjo3w4eprh25qlo4fjsxfjx5u/lib64/libgmt.so.6(GMT_Call_Module+0x33e)[0x2b0db863f2de]
gmt(main+0x412)[0x404ea2]
/lib64/libc.so.6(__libc_start_main+0xf5)[0x2b0dbc9a1555]
gmt[0x405bb3]

Actual outcome

Incorrect grid file with two triangle-like features:

Screen Shot 2022-12-04 at 3 21 23 PM

Expected outcome

Correct grid file with topography:

Screen Shot 2022-12-04 at 3 21 37 PM

System information

welcome[bot] commented 2 years ago

👋 Thanks for opening your first issue here! Please make sure you filled out the template with as much detail as possible. We appreciate that you took the time to contribute!

Please make sure you read our Contributing Guide and abide by our Code of Conduct.

PaulWessel commented 1 year ago

The triangular feature is most likely due to a dimension (say length of a column) that is off by 1 and hence points are placed successively off by one as we go across. This also explains why the triangle might go the other way (the dimension is 1 less or 1 more).