avoid tiny resolution differences

process/4_process_gebco.py

@@ -56,12 +56,20 @@ print('====================================================================')
 print(f'{"GENERAL BATHYMETRIC CHART OF THE OCEANS / GEBCO":^68}')
 
 # World
+
+# DEBUG
+if RESOLUTION_DEG == (1 / 60 / 60 * 15):
+       tr = ''
+else:
+       tr = f' -tr {RESOLUTION_DEG} {RESOLUTION_DEG}'
+#/DEBUG
+
 # ice surface
 msg = 'Resampling GEBCO ice surface...'
 cmd = (f'gdalwarp'
        f' -s_srs EPSG:{EPSG_WORLD}'
        f' -t_srs EPSG:{EPSG_WORLD}'
-       f' -tr {RESOLUTION_DEG} {RESOLUTION_DEG}'
+       f'{tr}'
        f' -r cubic'
        f' -of GTiff'
        f' {WARP_OPTIONS}'
@@ -100,7 +108,7 @@ msg = 'Resampling GEBCO sub ice...'
 cmd = (f'gdalwarp'
        f' -s_srs EPSG:{EPSG_WORLD}'
        f' -t_srs EPSG:{EPSG_WORLD}'
-       f' -tr {RESOLUTION_DEG} {RESOLUTION_DEG}'
+       f'{tr}'
        f' -r cubic'
        f' -of GTiff'
        f' {WARP_OPTIONS}'

process/config.py

@@ -194,8 +194,8 @@ GLIMS_IN = PATH_GLIMS / 'glims_polygons.shp'
 # warping in the lower latitudes.
 
 # Final AWI Basemap (World) resolution
-RESOLUTION_DEG = 0.004166666666666666  # 15 arc seconds (original GEBCO resolution)
-# RESOLUTION_DEG = 0.008333333333333333  # 30 arc seconds
+RESOLUTION_DEG =  1 / 60 / 60 * 15  # 15 arc seconds (original GEBCO resolution)
+# RESOLUTION_DEG =  1 / 60 / 60 * 30  # 30 arc seconds
 
 # Final AWI Basemap (Arctic and Antarctic) resolution
 RESOLUTION_M = 500.0  # approximate pendant to 15 arc seconds