2023-10-02 16:09:44,738 (natcap.invest.utils) utils.prepare_workspace(162) Level 100 Writing log messages to N:/PACIFIC/InVEST results\InVEST-Nutrient-Delivery-Ratio-log-2023-10-02--16_09_44.txt 2023-10-02 16:09:44,743 (natcap.invest.ui.model) model._logged_target(1643) Level 100 Starting model with parameters: Arguments for InVEST natcap.invest.ndr.ndr 3.10.2: biophysical_table_path N:/PACIFIC/InVEST workspace/biophysical_table_benta.csv calc_n False calc_p True dem_path N:/PACIFIC/InVEST workspace/dtm_8_buff.tif k_param 2 lulc_path N:/PACIFIC/InVEST workspace/lcm2015gb25m_1_reclass.tif n_workers -1 results_suffix _test2 runoff_proxy_path N:/PACIFIC/InVEST workspace/HadUK-Grid_MeanAnnualRainfall_1990_2017.tif subsurface_critical_length_n subsurface_eff_n threshold_flow_accumulation 100 watersheds_path N:/PACIFIC/InVEST workspace/catchment_sample2.shp workspace_dir N:/PACIFIC/InVEST results 2023-10-02 16:09:45,210 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(1741) INFO starting reprojection 2023-10-02 16:09:45,235 (osgeo) utils._log_gdal_errors(98) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer catchment_sample2 of geometry type POLYGON, which is not normally allowed by the GeoPackage specification, but the driver will however do it. To create a conformant GeoPackage, if using ogr2ogr, the -nlt option can be used to override the layer geometry type. This warning will no longer be emitted for this combination of layer and feature geometry type. 2023-10-02 16:09:45,293 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(1783) INFO reprojection 100.0% complete on watershed_results_ndr_test2.gpkg 2023-10-02 16:09:46,197 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-10-02 16:09:46,200 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 19896 2023-10-02 16:09:46,200 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-10-02 16:09:46,213 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-10-02 16:09:46,215 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-10-02 16:09:46,274 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 3 aligned: aligned_dem_test2.tif 2023-10-02 16:09:46,478 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-10-02 16:09:46,480 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 19896 2023-10-02 16:09:46,480 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-10-02 16:09:46,489 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-10-02 16:09:46,491 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-10-02 16:09:46,532 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 2 of 3 aligned: aligned_lulc_test2.tif 2023-10-02 16:09:46,728 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-10-02 16:09:46,731 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 19896 2023-10-02 16:09:46,731 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-10-02 16:09:46,743 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-10-02 16:09:46,745 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-10-02 16:09:46,801 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 3 of 3 aligned: aligned_runoff_proxy_test2.tif 2023-10-02 16:09:46,803 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 3 rasters. 2023-10-02 16:09:47,512 (pygeoprocessing.routing.routing) Task._call(1093) INFO (fill pits): complete 2023-10-02 16:09:48,369 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-10-02 16:09:48,458 (pygeoprocessing.routing.routing) Task._call(1093) DEBUG creating target flow accum raster layer 2023-10-02 16:09:48,569 (pygeoprocessing.routing.routing) Task._call(1093) DEBUG creating visited raster layer 2023-10-02 16:09:48,678 (pygeoprocessing.routing.routing) Task._call(1093) DEBUG starting search 2023-10-02 16:09:48,986 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-10-02 16:09:49,339 (pygeoprocessing.routing.routing) Task._call(1093) INFO filter out incomplete divergent streams 2023-10-02 16:09:49,389 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-10-02 16:09:49,804 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-10-02 16:09:49,807 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 19896 2023-10-02 16:09:49,807 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-10-02 16:09:49,834 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-10-02 16:09:49,835 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-10-02 16:09:50,087 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-10-02 16:09:50,090 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 19896 2023-10-02 16:09:50,090 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-10-02 16:09:50,106 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-10-02 16:09:50,108 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-10-02 16:09:50,271 (pygeoprocessing.routing.routing) Task._call(1093) DEBUG creating target flow accum raster layer 2023-10-02 16:09:50,359 (pygeoprocessing.routing.routing) Task._call(1093) DEBUG creating visited raster layer 2023-10-02 16:09:50,485 (pygeoprocessing.routing.routing) Task._call(1093) DEBUG starting search 2023-10-02 16:09:50,827 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-10-02 16:09:51,034 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-10-02 16:09:51,036 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 19896 2023-10-02 16:09:51,036 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-10-02 16:09:51,118 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-10-02 16:09:51,120 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-10-02 16:09:51,344 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-10-02 16:09:51,346 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 19896 2023-10-02 16:09:51,347 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-10-02 16:09:51,384 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-10-02 16:09:51,386 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-10-02 16:09:51,583 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-10-02 16:09:51,586 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 19896 2023-10-02 16:09:51,586 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-10-02 16:09:51,601 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-10-02 16:09:51,603 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-10-02 16:09:52,288 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-10-02 16:09:52,847 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-10-02 16:09:53,049 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-10-02 16:09:53,052 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 19896 2023-10-02 16:09:53,052 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-10-02 16:09:53,112 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-10-02 16:09:53,114 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-10-02 16:09:53,347 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-10-02 16:09:53,350 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 19896 2023-10-02 16:09:53,350 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-10-02 16:09:53,448 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-10-02 16:09:53,451 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-10-02 16:09:53,657 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-10-02 16:09:53,659 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 19896 2023-10-02 16:09:53,660 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-10-02 16:09:53,680 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-10-02 16:09:53,682 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-10-02 16:09:53,981 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-10-02 16:09:53,984 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 19896 2023-10-02 16:09:53,984 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-10-02 16:09:54,030 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-10-02 16:09:54,032 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-10-02 16:09:54,234 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-10-02 16:09:54,236 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 19896 2023-10-02 16:09:54,236 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-10-02 16:09:54,253 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-10-02 16:09:54,255 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-10-02 16:09:54,490 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-10-02 16:09:54,493 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 19896 2023-10-02 16:09:54,493 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-10-02 16:09:54,507 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-10-02 16:09:54,509 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-10-02 16:09:54,741 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-10-02 16:09:54,744 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 19896 2023-10-02 16:09:54,744 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-10-02 16:09:54,757 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-10-02 16:09:54,760 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-10-02 16:09:55,187 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-10-02 16:09:55,190 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 19896 2023-10-02 16:09:55,190 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-10-02 16:09:55,217 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-10-02 16:09:55,219 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-10-02 16:09:55,843 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-10-02 16:09:55,846 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 19896 2023-10-02 16:09:55,846 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-10-02 16:09:55,918 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-10-02 16:09:55,920 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-10-02 16:09:56,139 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-10-02 16:09:56,141 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 19896 2023-10-02 16:09:56,141 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-10-02 16:09:56,187 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-10-02 16:09:56,189 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-10-02 16:09:56,544 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2023-10-02 16:09:56,546 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2023-10-02 16:09:56,565 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2023-10-02 16:09:56,569 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2023-10-02 16:09:56,582 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2023-10-02 16:09:56,587 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2023-10-02 16:09:56,591 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_test2.gpkg 2023-10-02 16:09:56,592 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2023-10-02 16:09:56,595 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_test2.gpkg 2023-10-02 16:09:56,597 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_test2.gpkg 2023-10-02 16:09:56,654 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 1 disjoint polygon sets 2023-10-02 16:09:56,657 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 1 100.0% processed on watershed_results_ndr_test2.gpkg 2023-10-02 16:09:56,659 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 1 watershed_results_ndr_test2.gpkg 2023-10-02 16:09:56,699 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 1 watershed_results_ndr_test2.gpkg 2023-10-02 16:09:56,717 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 1 2023-10-02 16:09:56,724 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 1.3397298e+05 1.0000000e+02 0.0000000e+00 6.3613012e+04 0.0000000e+00 -1.0000000e+02] for ('N:/PACIFIC/InVEST results\\p_surface_export_test2.tif', 1) 2023-10-02 16:09:56,727 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 1 2023-10-02 16:09:56,729 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_test2.gpkg 2023-10-02 16:09:57,058 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2023-10-02 16:09:57,060 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2023-10-02 16:09:57,079 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2023-10-02 16:09:57,082 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2023-10-02 16:09:57,096 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2023-10-02 16:09:57,100 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2023-10-02 16:09:57,103 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_test2.gpkg 2023-10-02 16:09:57,105 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2023-10-02 16:09:57,107 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_test2.gpkg 2023-10-02 16:09:57,109 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_test2.gpkg 2023-10-02 16:09:57,165 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 1 disjoint polygon sets 2023-10-02 16:09:57,169 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 1 100.0% processed on watershed_results_ndr_test2.gpkg 2023-10-02 16:09:57,171 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 1 watershed_results_ndr_test2.gpkg 2023-10-02 16:09:57,210 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 1 watershed_results_ndr_test2.gpkg 2023-10-02 16:09:57,226 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 1 2023-10-02 16:09:57,233 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 1.3397298e+05 1.0000000e+02 0.0000000e+00 6.3613012e+04 0.0000000e+00 -1.0000000e+02] for ('N:/PACIFIC/InVEST results\\intermediate_outputs\\surface_load_p_test2.tif', 1) 2023-10-02 16:09:57,236 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 1 2023-10-02 16:09:57,238 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_test2.gpkg 2023-10-02 16:09:57,439 (natcap.invest.ndr.ndr) ndr.execute(730) INFO Writing summaries to output shapefile 2023-10-02 16:09:57,561 (osgeo) utils._log_gdal_errors(98) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer catchment_sample2 of geometry type POLYGON, which is not normally allowed by the GeoPackage specification, but the driver will however do it. To create a conformant GeoPackage, if using ogr2ogr, the -nlt option can be used to override the layer geometry type. This warning will no longer be emitted for this combination of layer and feature geometry type. 2023-10-02 16:09:57,717 (natcap.invest.ndr.ndr) ndr.execute(734) INFO NDR complete! 2023-10-02 16:09:57,719 (natcap.invest.ndr.ndr) ndr.execute(735) INFO _ _ ____ ____ 2023-10-02 16:09:57,721 (natcap.invest.ndr.ndr) ndr.execute(736) INFO | \ |"| | _"\U | _"\ u 2023-10-02 16:09:57,723 (natcap.invest.ndr.ndr) ndr.execute(737) INFO <| \| |>/| | | |\| |_) |/ 2023-10-02 16:09:57,725 (natcap.invest.ndr.ndr) ndr.execute(738) INFO U| |\ |uU| |_| |\| _ < 2023-10-02 16:09:57,727 (natcap.invest.ndr.ndr) ndr.execute(739) INFO |_| \_| |____/ u|_| \_\ 2023-10-02 16:09:57,729 (natcap.invest.ndr.ndr) ndr.execute(740) INFO || \\,-.|||_ // \\_ 2023-10-02 16:09:57,730 (natcap.invest.ndr.ndr) ndr.execute(741) INFO (_") (_/(__)_) (__) (__) 2023-10-02 16:09:57,733 (natcap.invest.utils) utils.prepare_workspace(170) INFO Elapsed time: 12.99s 2023-10-02 16:09:57,735 (natcap.invest.utils) utils.prepare_workspace(173) INFO Execution finished