3 and 1 second topo data producing error messages about resampling

[smalleyr]

Description of the problem

I get the following error message below I try to use either the 3 or 1 second dem data. It still makes the map, and seems to be using the correct resolution data.

Full script that generated the error

# PASTE CODE HERE

cat earth-relief-ocean.sh
#!/usr/bin/env bash
#
# Pseudo-color plot of the Earth relief data

gmt begin earth-relief-ocean png E100
#    gmt basemap -R-72/-71/-30.5/-29.5 -JM-71.5/15c -Bf
    gmt basemap -R-71.7/-71.45/-30.5/-30.25 -JM-71.5/60c -Bf

    # Plot the GMT Earth relief data using the default colormap/CPT and with
    # automatic hill shading (-I+d)

# at 30s, 15 s does not complain
#    gmt grdimage @earth_relief_30s -I+d
#    gmt grdimage @earth_relief_15s -I+d

#below - complains that output grid and input grid not same and has to resample output grid
#still produces map
#grdblend [WARNING]: File @S40W080.earth_relief_15s_p.nc has different increments (0.00416666666667/0.00416666666667) than the output grid (0.000833333333333/0.000833333333333) - must resample
#but seems to be reading wrong file, 15s!!
#looks a bit out of focus
#    gmt grdimage @earth_relief_03s -I+d

#same problem here, still produces map
#grdblend [WARNING]: File @S40W080.earth_relief_15s_p.nc has different increments (0.00416666666667/0.00416666666667) than the output grid (0.000277777777778/0.000277777777778) - must resample
#seems to be reading wrong file - 15s.
#good and sharp
    gmt grdimage @earth_relief_01s -I+d
gmt end show

**Short error message

earth-relief-ocean.sh grdblend [WARNING]: File @S40W080.earth_relief_15s_p.nc has different increments (0.00416666666667/0.00416666666667) than the output grid (0.000277777777778/0.000277777777778) - must resample grdblend [WARNING]: File @S40W080.earth_relief_15s_p.nc coordinates are phase-shifted w.r.t. the output grid - must resample

Full error message

earth-relief-ocean.sh gmt [DEBUG]: Obtained the ppid from parent: 45734 gmt [DEBUG]: Enter: gmtinit_new_GMT_ctrl gmt [DEBUG]: GMT->session.SHAREDIR = /Users/robertsmalley/anaconda3/envs/gmtsar/share/gmt gmt [DEBUG]: GMT->session.HOMEDIR = /Users/robertsmalley gmt [DEBUG]: GMT->session.USERDIR = /Users/robertsmalley/.gmt [created] gmt [DEBUG]: GMT->session.CACHEDIR = /Users/robertsmalley/.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 /Users/robertsmalley/.gmt gmt [DEBUG]: GMT: 3. gmt_getsharepath trying USERDIR subdir /Users/robertsmalley/.gmt/postscriptlight gmt [DEBUG]: GMT: 4. gmt_getsharepath trying SHAREDIR subdir /Users/robertsmalley/anaconda3/envs/gmtsar/share/gmt/postscriptlight gmt [DEBUG]: GMT: 5. gmt_getsharepath trying SHAREDIR /Users/robertsmalley/anaconda3/envs/gmtsar/share/gmt 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]: GMT now running in modern mode [Session ID = 45734] 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]: gmtlib_get_graphics_item: Fig: 0 Subplot: 2 Panel: () Inset: 0 gmt [DEBUG]: Reading GMT Default parameters from file: /Users/robertsmalley/.gmt/sessions/gmt_session.45734/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]: gmtlib_get_graphics_item: Fig: 0 Subplot: 2 Panel: () Inset: 0 gmt [DEBUG]: Enter: gmt_hash_init gmt [DEBUG]: Exit: gmt_hash_init gmt [DEBUG]: Exit: gmtinit_get_history gmt [DEBUG]: Initialize FFTW with 8 threads. gmt [DEBUG]: Local file /Users/robertsmalley/.gmt/server/gmt_data_server.txt found gmt [DEBUG]: File /Users/robertsmalley/.gmt/server/gmt_data_server.txt less than 24 hours old, refresh is premature. gmt [DEBUG]: Load contents from /Users/robertsmalley/.gmt/server/gmt_data_server.txt gmt [DEBUG]: Local file /Users/robertsmalley/.gmt/server/gmt_hash_server.txt found gmt [DEBUG]: File /Users/robertsmalley/.gmt/server/gmt_hash_server.txt less than 24 hours old, refresh is premature. gmt [DEBUG]: GMT_Create_Session initialized GMT structure gmt [DEBUG]: Loading core GMT shared library: libgmt.dylib gmt [DEBUG]: Shared Library # 0 (core). Path = libgmt.dylib gmt [DEBUG]: Loading GMT plugins from: /Users/robertsmalley/anaconda3/envs/gmtsar/lib/gmt/plugins gmt [DEBUG]: Shared Library # 1 (supplements). Path = /Users/robertsmalley/anaconda3/envs/gmtsar/lib/gmt/plugins/supplements.so gmt [DEBUG]: GMT now running in modern mode [Session ID = 45734] gmt [DEBUG]: Map distance calculation will be using great circle approximation with authalic auxiliary latitudes and authalic (R_2) radius = 6371007.1809 m, in meter. gmt [DEBUG]: Use PS filename /Users/robertsmalley/.gmt/sessions/gmt_session.45734/gmt_0.ps- gmt [DEBUG]: gmtinit_get_current_panel: No current panel selected so not in subplot mode gmt [DEBUG]: Modern mode: Added -R to options since history is available. gmt [DEBUG]: Modern: Adding -J to options since there is history available. gmt [DEBUG]: Map distance calculation will be using great circle approximation with authalic auxiliary latitudes and authalic (R_2) radius = 6371007.1809 m, in meter. gmt [DEBUG]: gmt_get_filename: In: -71.7/-71.45/-30.5/-30.25 Out: -71.7/-71.45/-30.5/-30.25 gmt [DEBUG]: Look for file -71.7/-71.45/-30.5/-30.25 in /Users/robertsmalley/.gmt gmt [DEBUG]: Look for file -71.7/-71.45/-30.5/-30.25 in /Users/robertsmalley/.gmt/cache gmt [DEBUG]: Look for file -71.7/-71.45/-30.5/-30.25 in /Users/robertsmalley/.gmt/server gmt [DEBUG]: Got regular w/e/s/n for region (-71.7/-71.45/-30.5/-30.25) gmt [DEBUG]: Remote file @srtm_tiles.nc exists locally as /Users/robertsmalley/.gmt/cache/srtm_tiles.nc gmt [DEBUG]: gmt_get_filename: In: srtm_tiles.nc Out: srtm_tiles.nc gmt [DEBUG]: Look for file srtm_tiles.nc in /Users/robertsmalley/.gmt gmt [DEBUG]: Look for file srtm_tiles.nc in /Users/robertsmalley/.gmt/cache gmt [DEBUG]: Found readable file /Users/robertsmalley/.gmt/cache/srtm_tiles.nc gmt [DEBUG]: Replace file srtm_tiles.nc with path /Users/robertsmalley/.gmt/cache/srtm_tiles.nc gmt [DEBUG]: gmt_get_filename: In: srtm_tiles.nc Out: srtm_tiles.nc gmt [DEBUG]: Look for file srtm_tiles.nc in /Users/robertsmalley/.gmt gmt [DEBUG]: Look for file srtm_tiles.nc in /Users/robertsmalley/.gmt/cache gmt [DEBUG]: gmt_get_filename: In: srtm_tiles.nc Out: srtm_tiles.nc gmt [DEBUG]: Look for file srtm_tiles.nc in /Users/robertsmalley/.gmt gmt [DEBUG]: Look for file srtm_tiles.nc in /Users/robertsmalley/.gmt/cache gmt [DEBUG]: Found readable file /Users/robertsmalley/.gmt/cache/srtm_tiles.nc gmt [DEBUG]: Object ID 0 : Registered Grid File srtm_tiles.nc as an Input resource with geometry Surface [n_objects = 1] gmt [DEBUG]: gmtapi_begin_io: Input resource access is now enabled [container] gmt [DEBUG]: gmtapi_import_grid: Passed ID = 0 and mode = 0 gmt [DEBUG]: Look for file srtm_tiles.nc in /Users/robertsmalley/.gmt gmt [DEBUG]: Look for file srtm_tiles.nc in /Users/robertsmalley/.gmt/cache gmt [DEBUG]: Found readable file /Users/robertsmalley/.gmt/cache/srtm_tiles.nc gmt [DEBUG]: gmt_get_filename: In: /Users/robertsmalley/.gmt/cache/srtm_tiles.nc Out: /Users/robertsmalley/.gmt/cache/srtm_tiles.nc gmt [DEBUG]: Call gmtgrdio_doctor_geo_increments on a geographic grid gmt [DEBUG]: Geographic input grid, longitudes span exactly 360 gmt [INFORMATION]: Reading grid from file srtm_tiles.nc gmt [DEBUG]: packed z-range: [0,1] gmt [DEBUG]: Geographic input grid, longitudes span exactly 360 gmt [DEBUG]: Grid is considered to have a 360-degree longitude range. gmt [DEBUG]: Chosen boundary condition for left and right edges: geographic gmt [DEBUG]: Chosen boundary condition for bottom and top edges: natural gmt [INFORMATION]: Set boundary condition for top edge: natural gmt [INFORMATION]: Set boundary condition for bottom edge: natural gmt [DEBUG]: GMT_End_IO: Input resource access is now disabled gmt [DEBUG]: GMT_Destroy_Data: freed memory for a Grid for object 0 gmt [DEBUG]: gmtlib_unregister_io: Unregistering object no 0 [n_objects = 0] gmt [DEBUG]: gmtlib_unregister_io: Object no 0 has non-NULL resource pointer gmt [DEBUG]: Revised options: @earth_relief_01s_g/ -I+d -Vd -R-71.7/-71.45/-30.5/-30.25 -JM-71.5/60c grdimage [DEBUG]: History: Process -R-71.7/-71.45/-30.5/-30.25 grdimage [DEBUG]: History: Process -JM-71.5/60c grdimage [DEBUG]: Map distance calculation will be using great circle approximation with authalic auxiliary latitudes and authalic (R_2) radius = 6371007.1809 m, in meter. grdimage [DEBUG]: gmt_get_filename: In: -71.7/-71.45/-30.5/-30.25 Out: -71.7/-71.45/-30.5/-30.25 grdimage [DEBUG]: Look for file -71.7/-71.45/-30.5/-30.25 in /Users/robertsmalley/.gmt grdimage [DEBUG]: Look for file -71.7/-71.45/-30.5/-30.25 in /Users/robertsmalley/.gmt/cache grdimage [DEBUG]: Look for file -71.7/-71.45/-30.5/-30.25 in /Users/robertsmalley/.gmt/server grdimage [DEBUG]: Got regular w/e/s/n for region (-71.7/-71.45/-30.5/-30.25) grdimage [DEBUG]: gmt_get_filename: In: /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 Out: /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdimage [DEBUG]: Given full path to file /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdimage [DEBUG]: Replace file /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 with /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdimage [DEBUG]: Projected grid is non-orthogonal, nonlinear, or dpi was changed grdimage [DEBUG]: gmt_get_filename: In: /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 Out: /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdimage [DEBUG]: Given full path to file /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdimage [DEBUG]: gmt_get_filename: In: /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 Out: /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdimage [DEBUG]: Given full path to file /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdimage [DEBUG]: gmt_get_filename: In: /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 Out: /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdimage [DEBUG]: gmt_get_filename: In: /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 Out: /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdimage [DEBUG]: Found readable file /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdimage [DEBUG]: Object ID 1 : Registered Grid File /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 as an Input resource with geometry Surface [n_objects = 1] grdimage [DEBUG]: gmtapi_begin_io: Input resource access is now enabled [container] grdimage [DEBUG]: gmtapi_import_grid: Passed ID = 1 and mode = 0 grdimage [DEBUG]: Object ID 2 : Registered Grid Memory Copy 7fe977607630 as an Output resource with geometry Surface [n_objects = 2] grdimage [DEBUG]: Successfully created a new Grid container grdimage [DEBUG]: VirtualFile name created: @GMTAPI@-S-O-G-G-G-Y-000002 grdimage [DEBUG]: GMT now running in modern mode [Session ID = 45734] grdimage [DEBUG]: Revised options: @earth_relief_01s_g/ -R-71.7/-71.45/-30.5/-30.25 -I01s -rg -G@GMTAPI@-S-O-G-G-G-Y-000002 -Co+n -N0 grdimage (gmtlib_free_tmp_arrays): tried to free unallocated memory grdblend [DEBUG]: History: Process -R-71.7/-71.45/-30.5/-30.25 grdblend [DEBUG]: gmt_get_filename: In: -71.7/-71.45/-30.5/-30.25 Out: -71.7/-71.45/-30.5/-30.25 grdblend [DEBUG]: Look for file -71.7/-71.45/-30.5/-30.25 in /Users/robertsmalley/.gmt grdblend [DEBUG]: Look for file -71.7/-71.45/-30.5/-30.25 in /Users/robertsmalley/.gmt/cache grdblend [DEBUG]: Look for file -71.7/-71.45/-30.5/-30.25 in /Users/robertsmalley/.gmt/server grdblend [DEBUG]: Got regular w/e/s/n for region (-71.7/-71.45/-30.5/-30.25) grdblend [DEBUG]: gmt_get_filename: In: /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 Out: /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdblend [DEBUG]: Given full path to file /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdblend [DEBUG]: Replace file /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 with /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdblend [DEBUG]: Map distance calculation will be using great circle approximation with authalic auxiliary latitudes and authalic (R_2) radius = 6371007.1809 m, in meter. grdblend [DEBUG]: gmtapi_init_import: Passed family = Data Table and geometry = Text grdblend [DEBUG]: gmt_get_filename: In: /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 Out: /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdblend [DEBUG]: Given full path to file /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdblend [DEBUG]: gmt_get_filename: In: /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 Out: /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdblend [DEBUG]: gmt_get_filename: In: /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 Out: /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdblend [DEBUG]: Found readable file /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdblend [DEBUG]: Object ID 3 : Registered Data Table File /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 as an Input resource with geometry Text [n_objects = 3] grdblend (gmtapi_init_import): tried to free unallocated memory grdblend [DEBUG]: gmtapi_init_import: Added 1 new sources grdblend [DEBUG]: GMT_Init_IO: Returned first Input object ID = 3 grdblend [DEBUG]: GMT_Begin_IO: Mode value 1 not considered (ignored) grdblend [DEBUG]: GMT_Begin_IO: Initialize record-by-record access for Input grdblend [DEBUG]: gmtapi_next_io_source: Selected object 3 grdblend [DEBUG]: gmt_get_filename: In: /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 Out: /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdblend [DEBUG]: Given full path to file /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdblend [DEBUG]: Found readable file /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdblend [DEBUG]: Found readable file /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdblend [INFORMATION]: Reading Data Table from file /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdblend [DEBUG]: GMT_Begin_IO: Input resource access is now enabled [record-by-record] grdblend [DEBUG]: Got tiled list determined from a plot region grdblend [DEBUG]: Geographic input grid, longitudes span less than 360 grdblend [DEBUG]: Chosen boundary condition for all edges: geographic grdblend [DEBUG]: Geographic input grid, longitudes span less than 360 grdblend [DEBUG]: Object ID 4 : Registered Grid Memory Reference 7fe977500ad0 as an Input resource with geometry Surface [n_objects = 4] grdblend [DEBUG]: Successfully created a new Grid container grdblend [DEBUG]: gmt_get_filename: In: S31W072.earth_relief_01s_g.nc Out: S31W072.earth_relief_01s_g.nc grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt/cache grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt/server grdblend [DEBUG]: gmt_get_filename: In: S40W080.earth_relief_15s_p.nc Out: S40W080.earth_relief_15s_p.nc grdblend [DEBUG]: Look for file S40W080.earth_relief_15s_p.nc in /Users/robertsmalley/.gmt grdblend [DEBUG]: Look for file S40W080.earth_relief_15s_p.nc in /Users/robertsmalley/.gmt/cache grdblend [DEBUG]: Look for file S40W080.earth_relief_15s_p.nc in /Users/robertsmalley/.gmt/server grdblend [DEBUG]: Remote file @S31W072.earth_relief_01s_g.nc exists locally as /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_01s_g/S31W072.earth_relief_01s_g.nc grdblend [DEBUG]: gmt_get_filename: In: S31W072.earth_relief_01s_g.nc Out: S31W072.earth_relief_01s_g.nc grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt/cache grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt/server grdblend [DEBUG]: Found readable file /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_01s_g/S31W072.earth_relief_01s_g.nc grdblend [DEBUG]: Replace file S31W072.earth_relief_01s_g.nc with path /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_01s_g/S31W072.earth_relief_01s_g.nc grdblend [DEBUG]: gmt_get_filename: In: S31W072.earth_relief_01s_g.nc Out: S31W072.earth_relief_01s_g.nc grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt/cache grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt/server grdblend [DEBUG]: gmt_get_filename: In: S31W072.earth_relief_01s_g.nc Out: S31W072.earth_relief_01s_g.nc grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt/cache grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt/server grdblend [DEBUG]: Found readable file /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_01s_g/S31W072.earth_relief_01s_g.nc grdblend [DEBUG]: Object ID 5 : Registered Grid File S31W072.earth_relief_01s_g.nc as an Input resource with geometry Surface [n_objects = 5] grdblend [DEBUG]: gmtapi_import_grid: Passed ID = 5 and mode = 33 grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt/cache grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt/server grdblend [DEBUG]: Found readable file /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_01s_g/S31W072.earth_relief_01s_g.nc grdblend [DEBUG]: gmt_get_filename: In: /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_01s_g/S31W072.earth_relief_01s_g.nc Out: /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_01s_g/S31W072.earth_relief_01s_g.nc grdblend [DEBUG]: Call gmtgrdio_doctor_geo_increments on a geographic grid grdblend [DEBUG]: Geographic input grid, longitudes span less than 360 grdblend [DEBUG]: Remote file @S40W080.earth_relief_15s_p.nc exists locally as /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc grdblend [DEBUG]: gmt_get_filename: In: S40W080.earth_relief_15s_p.nc Out: S40W080.earth_relief_15s_p.nc grdblend [DEBUG]: Look for file S40W080.earth_relief_15s_p.nc in /Users/robertsmalley/.gmt grdblend [DEBUG]: Look for file S40W080.earth_relief_15s_p.nc in /Users/robertsmalley/.gmt/cache grdblend [DEBUG]: Look for file S40W080.earth_relief_15s_p.nc in /Users/robertsmalley/.gmt/server grdblend [DEBUG]: Found readable file /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc grdblend [DEBUG]: Replace file S40W080.earth_relief_15s_p.nc with path /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc grdblend [DEBUG]: gmt_get_filename: In: S40W080.earth_relief_15s_p.nc Out: S40W080.earth_relief_15s_p.nc grdblend [DEBUG]: Look for file S40W080.earth_relief_15s_p.nc in /Users/robertsmalley/.gmt grdblend [DEBUG]: Look for file S40W080.earth_relief_15s_p.nc in /Users/robertsmalley/.gmt/cache grdblend [DEBUG]: Look for file S40W080.earth_relief_15s_p.nc in /Users/robertsmalley/.gmt/server grdblend [DEBUG]: gmt_get_filename: In: S40W080.earth_relief_15s_p.nc Out: S40W080.earth_relief_15s_p.nc grdblend [DEBUG]: Look for file S40W080.earth_relief_15s_p.nc in /Users/robertsmalley/.gmt grdblend [DEBUG]: Look for file S40W080.earth_relief_15s_p.nc in /Users/robertsmalley/.gmt/cache grdblend [DEBUG]: Look for file S40W080.earth_relief_15s_p.nc in /Users/robertsmalley/.gmt/server grdblend [DEBUG]: Found readable file /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc grdblend [DEBUG]: Object ID 6 : Registered Grid File S40W080.earth_relief_15s_p.nc as an Input resource with geometry Surface [n_objects = 6] grdblend [DEBUG]: gmtapi_import_grid: Passed ID = 6 and mode = 33 grdblend [DEBUG]: Look for file S40W080.earth_relief_15s_p.nc in /Users/robertsmalley/.gmt grdblend [DEBUG]: Look for file S40W080.earth_relief_15s_p.nc in /Users/robertsmalley/.gmt/cache grdblend [DEBUG]: Look for file S40W080.earth_relief_15s_p.nc in /Users/robertsmalley/.gmt/server grdblend [DEBUG]: Found readable file /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc grdblend [DEBUG]: gmt_get_filename: In: /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc Out: /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc grdblend [DEBUG]: Call gmtgrdio_doctor_geo_increments on a geographic grid grdblend [INFORMATION]: Round-off patrol changed geographic grid increment for longitude from 0.00416666666666666834 to 0.00416666666666666661 grdblend [DEBUG]: Geographic input grid, longitudes span less than 360 grdblend [DEBUG]: Grid @S31W072.earth_relief_01s_g.nc: out: -1080/2520/2700/-900 in: -1081/2521/2701/-901 skip: 0 offset: -900 grdblend [INFORMATION]: Blend file /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_01s_g/S31W072.earth_relief_01s_g.nc in -72/-71/-31/-30 with normal weight 1 [-900-2700] grdblend (gmtlib_close_grd): tried to free unallocated memory grdblend [DEBUG]: GMT_Destroy_Data: freed memory for a Grid for object 5 grdblend [DEBUG]: gmtlib_unregister_io: Unregistering object no 5 [n_objects = 5] grdblend [DEBUG]: gmtlib_unregister_io: Object no 5 has non-NULL resource pointer grdblend [WARNING]: File @S40W080.earth_relief_15s_p.nc has different increments (0.00416666666667/0.00416666666667) than the output grid (0.000277777777778/0.000277777777778) - must resample grdblend [WARNING]: File @S40W080.earth_relief_15s_p.nc coordinates are phase-shifted w.r.t. the output grid - must resample grdblend [INFORMATION]: Resample @S40W080.earth_relief_15s_p.nc via grdsample @S40W080.earth_relief_15s_p.nc -rg -I0.000277777777778/0.000277777777778 -R-71.7/-71.45/-30.5/-30.25 -G/var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc -Vd -fg --GMT_HISTORY=false grdblend [DEBUG]: GMT now running in modern mode [Session ID = 45734] grdblend [DEBUG]: Revised options: @S40W080.earth_relief_15s_p.nc -rg -I0.000277777777778/0.000277777777778 -R-71.7/-71.45/-30.5/-30.25 -G/var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc -Vd -fg --GMT_HISTORY=false grdblend (gmtlib_free_tmp_arrays): tried to free unallocated memory grdsample [DEBUG]: History: Process -R-71.7/-71.45/-30.5/-30.25 grdsample [DEBUG]: Map distance calculation will be using great circle approximation with authalic auxiliary latitudes and authalic (R_2) radius = 6371007.1809 m, in meter. grdsample [DEBUG]: gmt_get_filename: In: -71.7/-71.45/-30.5/-30.25 Out: -71.7/-71.45/-30.5/-30.25 grdsample [DEBUG]: Look for file -71.7/-71.45/-30.5/-30.25 in /Users/robertsmalley/.gmt grdsample [DEBUG]: Look for file -71.7/-71.45/-30.5/-30.25 in /Users/robertsmalley/.gmt/cache grdsample [DEBUG]: Look for file -71.7/-71.45/-30.5/-30.25 in /Users/robertsmalley/.gmt/server grdsample [DEBUG]: Got regular w/e/s/n for region (-71.7/-71.45/-30.5/-30.25) grdsample [DEBUG]: Remote file @S40W080.earth_relief_15s_p.nc exists locally as /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc grdsample [DEBUG]: Replace file @S40W080.earth_relief_15s_p.nc with /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc grdsample [INFORMATION]: Processing input grid grdsample [DEBUG]: gmt_get_filename: In: /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc Out: /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc grdsample [DEBUG]: Given full path to file /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc grdsample [DEBUG]: gmt_get_filename: In: /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc Out: /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc grdsample [DEBUG]: Given full path to file /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc grdsample [DEBUG]: gmt_get_filename: In: /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc Out: /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc grdsample [DEBUG]: gmt_get_filename: In: /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc Out: /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc grdsample [DEBUG]: Found readable file /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc grdsample [DEBUG]: Object ID 7 : Registered Grid File /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc as an Input resource with geometry Surface [n_objects = 6] grdsample [DEBUG]: gmtapi_import_grid: Passed ID = 7 and mode = 1 grdsample [DEBUG]: Found readable file /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc grdsample [DEBUG]: gmt_get_filename: In: /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc Out: /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc grdsample [DEBUG]: Call gmtgrdio_doctor_geo_increments on a geographic grid grdsample [INFORMATION]: Round-off patrol changed geographic grid increment for longitude from 0.00416666666666666834 to 0.00416666666666666661 grdsample [DEBUG]: Geographic input grid, longitudes span less than 360 grdsample [DEBUG]: Grid/Image dimensions imply w/e/s/n = -71.7/-71.45/-30.5/-30.25, inc = 0.000277778/0.000277778, gridline registration, n_layers = 1 grdsample [DEBUG]: Geographic input grid, longitudes span less than 360 grdsample [INFORMATION]: Given domain implies x_inc = 0.000277778 grdsample [INFORMATION]: Given domain implies y_inc = 0.000277778 grdsample [DEBUG]: Chosen boundary condition for all edges: geographic grdsample [DEBUG]: Geographic input grid, longitudes span less than 360 grdsample [DEBUG]: Object ID 8 : Registered Grid Memory Reference 7fe979980490 as an Input resource with geometry Surface [n_objects = 7] grdsample [DEBUG]: Successfully created a new Grid container grdsample [INFORMATION]: Input grid (-80/-70/-40/-30) n_columns = 2400 n_rows = 2400 dx = 0.00416666666667 dy = 0.00416666666667 registration = 1 grdsample [INFORMATION]: Output grid (-71.7/-71.45/-30.5/-30.25) n_columns = 901 n_rows = 901 dx = 0.000277777777778 dy = 0.000277777777778 registration = 0 grdsample [DEBUG]: gmtapi_import_grid: Passed ID = 7 and mode = 2 grdsample [INFORMATION]: Reading grid from file /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc grdsample [DEBUG]: packed z-range: [-12895,10882] grdsample [DEBUG]: Geographic input grid, longitudes span less than 360 grdsample [DEBUG]: Chosen boundary condition for left and right edges: geographic grdsample [DEBUG]: Chosen boundary condition for bottom edge: extended data grdsample [DEBUG]: Chosen boundary condition for top edge: geographic grdsample [INFORMATION]: Set boundary condition for left edge: natural grdsample [INFORMATION]: Set boundary condition for right edge: natural grdsample [INFORMATION]: Set boundary condition for top edge: natural grdsample [DEBUG]: Object ID 9 : Registered Grid File /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc as an Output resource with geometry Surface [n_objects = 8] grdsample [DEBUG]: gmtapi_begin_io: Output resource access is now enabled [container] grdsample [DEBUG]: gmtapi_export_grid: Passed ID = 9 and mode = 0 grdsample [INFORMATION]: Writing grid to file /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc grdsample [DEBUG]: region: 0 0, grid: -71.7 -71.45 grdsample [DEBUG]: -> region: -71.7 -71.45, grid: -71.7 -71.45 grdsample [DEBUG]: row: 0 900, col: 0 900 grdsample [DEBUG]: netCDF Library version: 4 grdsample [DEBUG]: packed z-range: [-643.765,655.608] grdsample [DEBUG]: GMT_End_IO: Output resource access is now disabled grdsample [DEBUG]: gmtlib_garbage_collection: Destroying object: C=0 A=0 ID=7 W=Input F=Grid M=File S=Used P=7fe977400020 N=/Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_15s_p/S40W080.earth_relief_15s_p.nc grdsample [DEBUG]: gmtlib_garbage_collection: Destroying object: C=0 A=0 ID=8 W=Input F=Grid M=Memory Reference S=Unused P=7fe979980490 N=(null) grdsample [DEBUG]: GMTAPI_Garbage_Collection freed 2 memory objects grdsample [DEBUG]: gmtlib_unregister_io: Unregistering object no 7 [n_objects = 7] grdsample [DEBUG]: gmtlib_unregister_io: Unregistering object no 8 [n_objects = 6] grdsample [DEBUG]: gmtlib_unregister_io: Unregistering object no 9 [n_objects = 5] grdsample (gmtlib_free_tmp_arrays): tried to free unallocated memory grdblend (gmtlib_close_grd): tried to free unallocated memory grdblend [DEBUG]: GMT_Destroy_Data: freed memory for a Grid for object 6 grdblend [DEBUG]: gmtlib_unregister_io: Unregistering object no 6 [n_objects = 4] grdblend [DEBUG]: gmtlib_unregister_io: Object no 6 has non-NULL resource pointer grdblend [DEBUG]: gmt_get_filename: In: /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc Out: /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc grdblend [DEBUG]: Given full path to file /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc grdblend [DEBUG]: gmt_get_filename: In: /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc Out: /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc grdblend [DEBUG]: Given full path to file /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc grdblend [DEBUG]: gmt_get_filename: In: /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc Out: /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc grdblend [DEBUG]: gmt_get_filename: In: /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc Out: /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc grdblend [DEBUG]: Found readable file /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc grdblend [DEBUG]: Object ID 10 : Registered Grid File /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc as an Input resource with geometry Surface [n_objects = 5] grdblend [DEBUG]: gmtapi_import_grid: Passed ID = 10 and mode = 33 grdblend [DEBUG]: Found readable file /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc grdblend [DEBUG]: gmt_get_filename: In: /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc Out: /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc grdblend [DEBUG]: Call gmtgrdio_doctor_geo_increments on a geographic grid grdblend [DEBUG]: Geographic input grid, longitudes span less than 360 grdblend [DEBUG]: Grid /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc: out: 0/900/900/0 in: -29881/6121/35101/-901 skip: 0 offset: 0 grdblend [INFORMATION]: Blend file /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc in -80/-70/-40/-30 with normal weight 1 [0-900] grdblend (gmtlib_close_grd): tried to free unallocated memory grdblend [DEBUG]: GMT_Destroy_Data: freed memory for a Grid for object 10 grdblend [DEBUG]: gmtlib_unregister_io: Unregistering object no 10 [n_objects = 4] grdblend [DEBUG]: gmtlib_unregister_io: Object no 10 has non-NULL resource pointer grdblend [DEBUG]: GMT_End_IO: Input resource access is now disabled grdblend [INFORMATION]: Processing input grids grdblend [DEBUG]: Successfully added data array to previously registered Grid container grdblend [DEBUG]: Remote file @S31W072.earth_relief_01s_g.nc exists locally as /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_01s_g/S31W072.earth_relief_01s_g.nc grdblend [DEBUG]: gmt_get_filename: In: S31W072.earth_relief_01s_g.nc Out: S31W072.earth_relief_01s_g.nc grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt/cache grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt/server grdblend [DEBUG]: Found readable file /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_01s_g/S31W072.earth_relief_01s_g.nc grdblend [DEBUG]: Replace file S31W072.earth_relief_01s_g.nc with path /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_01s_g/S31W072.earth_relief_01s_g.nc grdblend [DEBUG]: gmt_get_filename: In: S31W072.earth_relief_01s_g.nc Out: S31W072.earth_relief_01s_g.nc grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt/cache grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt/server grdblend [DEBUG]: gmt_get_filename: In: S31W072.earth_relief_01s_g.nc Out: S31W072.earth_relief_01s_g.nc grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt/cache grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt/server grdblend [DEBUG]: Found readable file /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_01s_g/S31W072.earth_relief_01s_g.nc grdblend [DEBUG]: Object ID 11 : Registered Grid File S31W072.earth_relief_01s_g.nc as an Input resource with geometry Surface [n_objects = 5] grdblend [DEBUG]: gmtapi_begin_io: Input resource access is now enabled [container] grdblend [DEBUG]: gmtapi_import_grid: Passed ID = 11 and mode = 33 grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt/cache grdblend [DEBUG]: Look for file S31W072.earth_relief_01s_g.nc in /Users/robertsmalley/.gmt/server grdblend [DEBUG]: Found readable file /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_01s_g/S31W072.earth_relief_01s_g.nc grdblend [DEBUG]: gmt_get_filename: In: /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_01s_g/S31W072.earth_relief_01s_g.nc Out: /Users/robertsmalley/.gmt/server/earth/earth_relief/earth_relief_01s_g/S31W072.earth_relief_01s_g.nc grdblend [DEBUG]: Call gmtgrdio_doctor_geo_increments on a geographic grid grdblend [DEBUG]: Geographic input grid, longitudes span less than 360 grdblend [DEBUG]: GMT_End_IO: Input resource access is now disabled grdblend [DEBUG]: gmt_get_filename: In: /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc Out: /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc grdblend [DEBUG]: Given full path to file /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc grdblend [DEBUG]: gmt_get_filename: In: /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc Out: /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc grdblend [DEBUG]: Given full path to file /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc grdblend [DEBUG]: gmt_get_filename: In: /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc Out: /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc grdblend [DEBUG]: gmt_get_filename: In: /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc Out: /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc grdblend [DEBUG]: Found readable file /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc grdblend [DEBUG]: Object ID 12 : Registered Grid File /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc as an Input resource with geometry Surface [n_objects = 6] grdblend [DEBUG]: gmtapi_begin_io: Input resource access is now enabled [container] grdblend [DEBUG]: gmtapi_import_grid: Passed ID = 12 and mode = 33 grdblend [DEBUG]: Found readable file /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc grdblend [DEBUG]: gmt_get_filename: In: /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc Out: /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc grdblend [DEBUG]: Call gmtgrdio_doctor_geo_increments on a geographic grid grdblend [DEBUG]: Geographic input grid, longitudes span less than 360 grdblend [DEBUG]: GMT_End_IO: Input resource access is now disabled grdblend [INFORMATION]: Processed row 901 of 901 grdblend [DEBUG]: GMT_Write_Data: Writing Grid to memory object 2 from object 4 which transfers ownership grdblend [DEBUG]: gmtapi_begin_io: Output resource access is now enabled [container] grdblend [DEBUG]: gmtapi_export_data: Messenger dummy output container for object 2 [item 1] freed and set resource=data=NULL grdblend [DEBUG]: gmtapi_export_grid: Passed ID = 2 and mode = 0 grdblend [INFORMATION]: Referencing grid data to GMT_GRID memory location grdblend [DEBUG]: Chosen boundary condition for all edges: geographic grdblend [INFORMATION]: Set boundary condition for all edges: natural grdblend [INFORMATION]: Set boundary condition for left edge: natural grdblend [INFORMATION]: Set boundary condition for right edge: natural grdblend [INFORMATION]: Set boundary condition for bottom edge: natural grdblend [INFORMATION]: Set boundary condition for top edge: natural grdblend [DEBUG]: GMT_End_IO: Output resource access is now disabled grdblend (gmtlib_close_grd): tried to free unallocated memory grdblend [DEBUG]: GMT_Destroy_Data: freed memory for a Grid for object 11 grdblend [DEBUG]: gmtlib_unregister_io: Unregistering object no 11 [n_objects = 5] grdblend [DEBUG]: gmtlib_unregister_io: Object no 11 has non-NULL resource pointer grdblend [DEBUG]: Delete /var/folders/xx/rznqg73d2kz2_0z38wvr_dzw0000gn/T/grdblend_resampled_45737_1.nc grdblend (gmtlib_close_grd): tried to free unallocated memory grdblend [DEBUG]: GMT_Destroy_Data: freed memory for a Grid for object 12 grdblend [DEBUG]: gmtlib_unregister_io: Unregistering object no 12 [n_objects = 4] grdblend [DEBUG]: gmtlib_unregister_io: Object no 12 has non-NULL resource pointer grdblend [INFORMATION]: Blended grid size of @GMTAPI@-S-O-G-G-G-Y-000002 is 901 x 901 grdblend [INFORMATION]: All nodes assigned values grdblend [DEBUG]: gmtlib_unregister_io: Unregistering object no 3 [n_objects = 3] grdblend [DEBUG]: gmtlib_unregister_io: Unregistering object no 4 [n_objects = 2] grdblend (gmtlib_free_tmp_arrays): tried to free unallocated memory grdimage [INFORMATION]: Set boundary condition for all edges: natural grdimage [INFORMATION]: Set boundary condition for left edge: natural grdimage [INFORMATION]: Set boundary condition for right edge: natural grdimage [INFORMATION]: Set boundary condition for bottom edge: natural grdimage [INFORMATION]: Set boundary condition for top edge: natural grdimage [DEBUG]: GMT_End_IO: Input resource access is now disabled grdimage [DEBUG]: VirtualFile name created: @GMTAPI@-S-I-G-G-G-N-000001 grdimage [DEBUG]: Found readable file /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdimage [DEBUG]: /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 is most likely a grid. Open in GMT as grid grdimage [DEBUG]: Projected values in meters: -22263.9 5565.97 -3.5463e+06 -3.5142e+06 grdimage [INFORMATION]: Map scale is 0.463831 km per cm or 1:46383.1. grdimage [INFORMATION]: Derive intensity grid from data grid grdimage [DEBUG]: Object ID 13 : Registered Grid Memory Copy 7fe977400020 as an Output resource with geometry Surface [n_objects = 3] grdimage [DEBUG]: Successfully created a new Grid container grdimage [DEBUG]: VirtualFile name created: @GMTAPI@-S-O-G-G-G-Y-000013 grdimage [INFORMATION]: Calling grdgradient with args -G@GMTAPI@-S-O-G-G-G-Y-000013 -A-45.0 -Nt1+a0 -R-71.7/-71.45/-30.5/-30.25 --GMT_HISTORY=false @GMTAPI@-S-I-G-G-G-N-000001 grdimage [DEBUG]: GMT now running in modern mode [Session ID = 45734] grdimage [DEBUG]: Revised options: -G@GMTAPI@-S-O-G-G-G-Y-000013 -A-45.0 -Nt1+a0 -R-71.7/-71.45/-30.5/-30.25 --GMT_HISTORY=false @GMTAPI@-S-I-G-G-G-N-000001 grdimage (gmtlib_free_tmp_arrays): tried to free unallocated memory grdgradient [DEBUG]: History: Process -R-71.7/-71.45/-30.5/-30.25 grdgradient [DEBUG]: gmt_get_filename: In: -71.7/-71.45/-30.5/-30.25 Out: -71.7/-71.45/-30.5/-30.25 grdgradient [DEBUG]: Look for file -71.7/-71.45/-30.5/-30.25 in /Users/robertsmalley/.gmt grdgradient [DEBUG]: Look for file -71.7/-71.45/-30.5/-30.25 in /Users/robertsmalley/.gmt/cache grdgradient [DEBUG]: Look for file -71.7/-71.45/-30.5/-30.25 in /Users/robertsmalley/.gmt/server grdgradient [DEBUG]: Got regular w/e/s/n for region (-71.7/-71.45/-30.5/-30.25) grdgradient [DEBUG]: gmt_get_filename: In: -45.0 Out: -45.0 grdgradient [DEBUG]: Look for file -45.0 in /Users/robertsmalley/.gmt grdgradient [DEBUG]: Look for file -45.0 in /Users/robertsmalley/.gmt/cache grdgradient [DEBUG]: Look for file -45.0 in /Users/robertsmalley/.gmt/server grdgradient [DEBUG]: amplitude_set = Y offset_set = N sigma_set = N grdgradient [INFORMATION]: Processing input grid grdgradient [DEBUG]: Geographic input grid, longitudes span less than 360 grdgradient [DEBUG]: Geographic input grid, longitudes span less than 360 grdgradient [DEBUG]: gmtlib_expand_headerpad: No pad adjustment needed grdgradient [DEBUG]: Chosen boundary condition for all edges: natural grdgradient [INFORMATION]: Set boundary condition for all edges: natural grdgradient [INFORMATION]: Set boundary condition for left edge: natural grdgradient [INFORMATION]: Set boundary condition for right edge: natural grdgradient [INFORMATION]: Set boundary condition for bottom edge: natural grdgradient [INFORMATION]: Set boundary condition for top edge: natural grdgradient [DEBUG]: Object ID 14 : Registered Grid Memory Reference 7fe977607c10 as an Input resource with geometry Surface [n_objects = 4] grdgradient [DEBUG]: Successfully duplicated a Grid grdgradient [INFORMATION]: Min Mean Max sigma intensities:grdgradient [INFORMATION]: -0.891533229346 0.0199924687431 0.945773560813 0.135881893739 grdgradient [DEBUG]: GMT_Write_Data: Writing Grid to memory object 13 from object 14 which transfers ownership grdgradient [DEBUG]: gmtapi_begin_io: Output resource access is now enabled [container] grdgradient [DEBUG]: gmtapi_export_data: Messenger dummy output container for object 13 [item 2] freed and set resource=data=NULL grdgradient [DEBUG]: gmtapi_export_grid: Passed ID = 13 and mode = 0 grdgradient [INFORMATION]: Referencing grid data to GMT_GRID memory location grdgradient [DEBUG]: Chosen boundary condition for all edges: natural grdgradient [INFORMATION]: Set boundary condition for all edges: natural grdgradient [INFORMATION]: Set boundary condition for left edge: natural grdgradient [INFORMATION]: Set boundary condition for right edge: natural grdgradient [INFORMATION]: Set boundary condition for bottom edge: natural grdgradient [INFORMATION]: Set boundary condition for top edge: natural grdgradient [DEBUG]: GMT_End_IO: Output resource access is now disabled grdgradient [DEBUG]: gmtlib_unregister_io: Unregistering object no 14 [n_objects = 3] grdgradient (gmtlib_free_tmp_arrays): tried to free unallocated memory grdimage [DEBUG]: Running in PS mode modern grdimage [DEBUG]: Use PS filename /Users/robertsmalley/.gmt/sessions/gmt_session.45734/gmt_0.ps- grdimage [DEBUG]: Append to hidden PS file /Users/robertsmalley/.gmt/sessions/gmt_session.45734/gmt_0.ps- grdimage [DEBUG]: No figure file /Users/robertsmalley/.gmt/sessions/gmt_session.45734/gmt.figures - nothing to do grdimage [DEBUG]: Got session name as earth-relief-ocean and default graphics formats as png grdimage [INFORMATION]: Allocate and read data from file /Users/robertsmalley/.gmt/sessions/gmt_session.45734/=tiled_75_PO.000000 grdimage [INFORMATION]: Allocates memory and read intensity file grdimage [INFORMATION]: project grid files grdimage [DEBUG]: Object ID 15 : Registered Grid Memory Reference 7fe9776075d0 as an Input resource with geometry Surface [n_objects = 4] grdimage [DEBUG]: Successfully duplicated a Grid grdimage [DEBUG]: gmt_project_init: IN: Inc [0/0] n_columns/n_rows [901/901] dpi = 0 offset = 0 grdimage [DEBUG]: gmt_project_init: OUT: Inc [0/0] n_columns/n_rows [901/901] dpi = 0 offset = 0 grdimage [INFORMATION]: Grid projection from size 901x901 to 901x901 grdimage [DEBUG]: Successfully added data array to previously registered Grid container grdimage [DEBUG]: gmt_grd_project: In [-71.7/-71.45/-30.5/-30.25] and out [0/23.6220472441/0/27.2438145383] grdimage [DEBUG]: GMT_Destroy_Data: freed memory for a Grid for object 1 grdimage [DEBUG]: gmtlib_unregister_io: Unregistering object no 1 [n_objects = 3] grdimage [DEBUG]: gmtlib_unregister_io: Object no 1 has non-NULL resource pointer grdimage [DEBUG]: Object ID 16 : Registered Grid Memory Reference 7fe9774032f0 as an Input resource with geometry Surface [n_objects = 4] grdimage [DEBUG]: Successfully duplicated a Grid grdimage [DEBUG]: gmt_project_init: IN: Inc [0/0] n_columns/n_rows [901/901] dpi = 0 offset = 0 grdimage [DEBUG]: gmt_project_init: OUT: Inc [0/0] n_columns/n_rows [901/901] dpi = 0 offset = 0 grdimage [INFORMATION]: Grid projection from size 901x901 to 901x901 grdimage [DEBUG]: Successfully added data array to previously registered Grid container grdimage [DEBUG]: gmt_grd_project: In [-71.7/-71.45/-30.5/-30.25] and out [0/23.6220472441/0/27.2438145383] grdimage [DEBUG]: GMT_Destroy_Data: freed memory for a Grid for object 13 grdimage [DEBUG]: gmtlib_unregister_io: Unregistering object no 13 [n_objects = 3] grdimage [DEBUG]: gmtlib_unregister_io: Object no 13 has non-NULL resource pointer grdimage [DEBUG]: gmtlib_get_graphics_item: Fig: 0 Subplot: 2 Panel: () Inset: 0 grdimage [DEBUG]: No current cpt file found grdimage [DEBUG]: gmt_get_filename: In: geo Out: geo grdimage [DEBUG]: Look for file geo in /Users/robertsmalley/.gmt grdimage [DEBUG]: Look for file geo in /Users/robertsmalley/.gmt/cache grdimage [DEBUG]: Look for file geo in /Users/robertsmalley/.gmt/server grdimage [DEBUG]: GMT: 0. Will try to find subdir=cpt stem = geo suffix=.cpt grdimage [DEBUG]: GMT: 1. gmt_getsharepath trying current dir grdimage [DEBUG]: GMT: 2. gmt_getsharepath trying USERDIR /Users/robertsmalley/.gmt grdimage [DEBUG]: GMT: 3. gmt_getsharepath trying USERDIR subdir /Users/robertsmalley/.gmt/cpt grdimage [DEBUG]: GMT: 4. gmt_getsharepath trying SHAREDIR subdir /Users/robertsmalley/anaconda3/envs/gmtsar/share/gmt/cpt grdimage [DEBUG]: gmt_get_filename: In: /Users/robertsmalley/anaconda3/envs/gmtsar/share/gmt/cpt/geo.cpt Out: /Users/robertsmalley/anaconda3/envs/gmtsar/share/gmt/cpt/geo.cpt grdimage [DEBUG]: Given full path to file /Users/robertsmalley/anaconda3/envs/gmtsar/share/gmt/cpt/geo.cpt grdimage [DEBUG]: gmt_get_filename: In: /Users/robertsmalley/anaconda3/envs/gmtsar/share/gmt/cpt/geo.cpt Out: /Users/robertsmalley/anaconda3/envs/gmtsar/share/gmt/cpt/geo.cpt grdimage [DEBUG]: gmt_get_filename: In: /Users/robertsmalley/anaconda3/envs/gmtsar/share/gmt/cpt/geo.cpt Out: /Users/robertsmalley/anaconda3/envs/gmtsar/share/gmt/cpt/geo.cpt grdimage [DEBUG]: Found readable file /Users/robertsmalley/anaconda3/envs/gmtsar/share/gmt/cpt/geo.cpt grdimage [DEBUG]: Object ID 17 : Registered CPT File /Users/robertsmalley/anaconda3/envs/gmtsar/share/gmt/cpt/geo.cpt as an Input resource with geometry Non-Geographical [n_objects = 4] grdimage [DEBUG]: gmtapi_begin_io: Input resource access is now enabled [container] grdimage [DEBUG]: gmtapi_import_palette: Passed ID = 17 and mode = 8 grdimage [INFORMATION]: Reading CPT from File /Users/robertsmalley/anaconda3/envs/gmtsar/share/gmt/cpt/geo.cpt grdimage [DEBUG]: Reading CPT from /Users/robertsmalley/anaconda3/envs/gmtsar/share/gmt/cpt/geo.cpt grdimage [DEBUG]: GMT_End_IO: Input resource access is now disabled grdimage [INFORMATION]: Auto-stretching CPT file geo to fit data range -643.765 to 741 grdimage [DEBUG]: Found CPT hinge at z' = 0 for slice k = 8! grdimage [DEBUG]: gmtlib_get_graphics_item: Fig: 0 Subplot: 2 Panel: () Inset: 0 grdimage [DEBUG]: Object ID 18 : Registered CPT File /Users/robertsmalley/.gmt/sessions/gmt_session.45734/gmt.cpt as an Output resource with geometry Non-Geographical [n_objects = 5] grdimage [DEBUG]: gmtapi_begin_io: Output resource access is now enabled [container] grdimage [DEBUG]: gmtapi_export_palette: Passed ID = 18 and mode = 0 grdimage [INFORMATION]: Write CPT to File /Users/robertsmalley/.gmt/sessions/gmt_session.45734/gmt.cpt grdimage [DEBUG]: Writing CPT to /Users/robertsmalley/.gmt/sessions/gmt_session.45734/gmt.cpt grdimage [DEBUG]: GMT_End_IO: Output resource access is now disabled grdimage [INFORMATION]: Save current CPT file to /Users/robertsmalley/.gmt/sessions/gmt_session.45734/gmt.cpt ! grdimage [INFORMATION]: Evaluate image pixel colors grdimage [INFORMATION]: Basic z(x,y) with optional illumination and no PostScript colormasking. grdimage [DEBUG]: GMT_Destroy_Data: freed memory for a Grid for object 16 grdimage [DEBUG]: gmtlib_unregister_io: Unregistering object no 16 [n_objects = 4] grdimage [DEBUG]: gmtlib_unregister_io: Object no 16 has non-NULL resource pointer grdimage [INFORMATION]: Creating 24-bit color image PSL: Too many colors to make colormap - using 24-bit direct color instead. PSL: DEFLATE compressed 2435403 to 1987634 bytes (18.4% savings at compression level 5) grdimage [DEBUG]: Current size of half-baked PS file /Users/robertsmalley/.gmt/sessions/gmt_session.45734/gmt_0.ps- = 2538555. grdimage (GMT_grdimage): tried to free unallocated memory grdimage [DEBUG]: GMT_Destroy_Data: freed memory for a Grid for object 15 grdimage [DEBUG]: gmtlib_unregister_io: Unregistering object no 15 [n_objects = 3] grdimage [DEBUG]: gmtlib_unregister_io: Object no 15 has non-NULL resource pointer grdimage [DEBUG]: gmtlib_garbage_collection: Destroying object: C=0 A=0 ID=17 W=Input F=CPT M=File S=Used P=7fe9776168f0 N=/Users/robertsmalley/anaconda3/envs/gmtsar/share/gmt/cpt/geo.cpt grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage (gmtsupport_free_range): tried to free unallocated memory grdimage [DEBUG]: GMTAPI_Garbage_Collection freed 1 memory objects grdimage [DEBUG]: gmtlib_unregister_io: Unregistering object no 2 [n_objects = 2] grdimage [DEBUG]: gmtlib_unregister_io: Unregistering object no 17 [n_objects = 1] grdimage [DEBUG]: gmtlib_unregister_io: Unregistering object no 18 [n_objects = 0] grdimage (gmtlib_free_tmp_arrays): tried to free unallocated memory gmt [DEBUG]: Entering GMT_Destroy_Session gmt [DEBUG]: gmtlib_get_graphics_item: Fig: 0 Subplot: 2 Panel: () Inset: 0

Actual outcome

Expected outcome

System information

• Operating system:
• GMT version (gmt --version):

[welcome[bot]]

👋 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.

[joa-quim]

Please, use the forum to ask questions. This place is to report GMT bugs.

[PaulWessel]

I see no errors, just warnings. As for mixing 15s and 1/3s over ocean, all explained in the docs: file:///Users/pwessel/GMTdev/gmt-dev/dbuild/gmt6/share/doc/html/datasets/remote-data.html#technical-information