2022-12-20 12:19:53,333 (natcap.invest.utils) utils.prepare_workspace(161) Level 100 Writing log messages to [C:\invest_test666\InVEST-natcap.invest.seasonal_water_yield.seasonal_water_yield-log-2022-12-20--12_19_53.txt] 2022-12-20 12:19:53,334 (__main__) cli.main(475) Level 100 Starting model with parameters: Arguments for InVEST natcap.invest.seasonal_water_yield.seasonal_water_yield 3.12.0: alpha_m 1/12 aoi_path C:\invest_test666\area_29N.shp beta_i 1 biophysical_table_path C:\invest_test666\biophysical_table_99.csv climate_zone_raster_path climate_zone_table_path dem_raster_path C:\invest_test666\DEM_cut_ocs.tif et0_dir C:\invest_test666\ET0_29N gamma 1 l_path lulc_raster_path C:\invest_test666\ocs_nodata7_sans0raster.tif monthly_alpha False monthly_alpha_path n_workers -1 precip_dir C:\invest_test666\monthly_precipitation_mm_32229 rain_events_table_path C:\invest_test666\rain_events_good.csv results_suffix soil_group_path C:\invest_test666\HYdro_sans0_nodata9999.tif threshold_flow_accumulation 1000 user_defined_climate_zones False user_defined_local_recharge False workspace_dir C:\invest_test666 2022-12-20 12:19:53,335 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(423) INFO prepare and test inputs for common errors 2022-12-20 12:19:53,456 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(489) INFO Building file registry 2022-12-20 12:19:53,457 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(495) INFO Checking that the AOI is not the output aggregate vector 2022-12-20 12:19:53,458 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(504) INFO Aligning and clipping dataset list 2022-12-20 12:19:55,024 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 27 aligned: prcp_a0.tif 2022-12-20 12:19:56,247 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 2 of 27 aligned: prcp_a1.tif 2022-12-20 12:19:57,248 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 3 of 27 aligned: prcp_a2.tif 2022-12-20 12:19:58,471 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 4 of 27 aligned: prcp_a3.tif 2022-12-20 12:19:59,424 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 5 of 27 aligned: prcp_a4.tif 2022-12-20 12:20:00,570 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 6 of 27 aligned: prcp_a5.tif 2022-12-20 12:20:01,774 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 7 of 27 aligned: prcp_a6.tif 2022-12-20 12:20:03,096 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 8 of 27 aligned: prcp_a7.tif 2022-12-20 12:20:04,235 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 9 of 27 aligned: prcp_a8.tif 2022-12-20 12:20:05,409 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 10 of 27 aligned: prcp_a9.tif 2022-12-20 12:20:06,536 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 11 of 27 aligned: prcp_a10.tif 2022-12-20 12:20:07,854 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 12 of 27 aligned: prcp_a11.tif 2022-12-20 12:20:08,577 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 13 of 27 aligned: soil_group_aligned.tif 2022-12-20 12:20:09,055 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 14 of 27 aligned: et0_a0.tif 2022-12-20 12:20:09,495 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 15 of 27 aligned: et0_a1.tif 2022-12-20 12:20:09,930 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 16 of 27 aligned: et0_a2.tif 2022-12-20 12:20:10,436 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 17 of 27 aligned: et0_a3.tif 2022-12-20 12:20:10,976 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 18 of 27 aligned: et0_a4.tif 2022-12-20 12:20:11,448 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 19 of 27 aligned: et0_a5.tif 2022-12-20 12:20:12,006 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 20 of 27 aligned: et0_a6.tif 2022-12-20 12:20:12,404 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 21 of 27 aligned: et0_a7.tif 2022-12-20 12:20:12,858 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 22 of 27 aligned: et0_a8.tif 2022-12-20 12:20:13,388 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 23 of 27 aligned: et0_a9.tif 2022-12-20 12:20:14,472 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 24 of 27 aligned: et0_a10.tif 2022-12-20 12:20:15,621 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 25 of 27 aligned: et0_a11.tif 2022-12-20 12:20:16,448 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 26 of 27 aligned: lulc_aligned.tif 2022-12-20 12:20:18,152 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 27 of 27 aligned: dem_aligned.tif 2022-12-20 12:20:18,153 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 27 rasters. 2022-12-20 12:20:29,036 (pygeoprocessing.routing.routing) Task._call(1093) INFO (fill pits): 9093120 of 15862920 pixels complete 2022-12-20 12:20:37,352 (pygeoprocessing.routing.routing) Task._call(1093) INFO (fill pits): complete 2022-12-20 12:20:48,089 (pygeoprocessing.routing.routing) Task._call(1093) INFO 31.9% complete 2022-12-20 12:20:59,062 (pygeoprocessing.routing.routing) Task._call(1093) INFO 70.1% complete 2022-12-20 12:21:07,959 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2022-12-20 12:21:18,012 (pygeoprocessing.routing.routing) Task._call(1093) INFO mfd flow accum 50.8% complete 2022-12-20 12:21:29,459 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2022-12-20 12:21:38,396 (pygeoprocessing.routing.routing) Task._call(1093) INFO filter out incomplete divergent streams 2022-12-20 12:21:38,631 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2022-12-20 12:21:38,651 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(606) INFO quick flow 2022-12-20 12:21:38,652 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(620) INFO loading number of monthly events 2022-12-20 12:21:39,066 (taskgraph.Task) Task.add_task(653) INFO A task was created that had an identical args signature sans target paths, but a different target_path_list of the same length. To avoid recomputation, dynamically adding previous Task (n_events as a constant raster month 0 (5)) as a dependent task to this one (n_events as a constant raster month 1 (6)). 2022-12-20 12:21:39,071 (taskgraph.Task) Task.is_precalculated(1252) INFO not precalculated (n_events as a constant raster month 1 (6)), Task hash exists, but there are these mismatches: Recorded path not in target path list c:\invest_test666\cache_dir\n_events0.tif 2022-12-20 12:21:39,483 (taskgraph.Task) Task.add_task(653) INFO A task was created that had an identical args signature sans target paths, but a different target_path_list of the same length. To avoid recomputation, dynamically adding previous Task (n_events as a constant raster month 1 (6)) as a dependent task to this one (n_events as a constant raster month 2 (6)). 2022-12-20 12:21:39,486 (taskgraph.Task) Task.is_precalculated(1252) INFO not precalculated (n_events as a constant raster month 2 (6)), Task hash exists, but there are these mismatches: Recorded path not in target path list c:\invest_test666\cache_dir\n_events1.tif 2022-12-20 12:21:43,096 (taskgraph.Task) Task.add_task(653) INFO A task was created that had an identical args signature sans target paths, but a different target_path_list of the same length. To avoid recomputation, dynamically adding previous Task (n_events as a constant raster month 2 (6)) as a dependent task to this one (n_events as a constant raster month 10 (13)). 2022-12-20 12:21:43,099 (taskgraph.Task) Task.is_precalculated(1252) INFO not precalculated (n_events as a constant raster month 10 (13)), Task hash exists, but there are these mismatches: Recorded path not in target path list c:\invest_test666\cache_dir\n_events2.tif 2022-12-20 12:21:43,561 (taskgraph.Task) Task.add_task(653) INFO A task was created that had an identical args signature sans target paths, but a different target_path_list of the same length. To avoid recomputation, dynamically adding previous Task (n_events as a constant raster month 10 (13)) as a dependent task to this one (n_events as a constant raster month 11 (13)). 2022-12-20 12:21:43,564 (taskgraph.Task) Task.is_precalculated(1252) INFO not precalculated (n_events as a constant raster month 11 (13)), Task hash exists, but there are these mismatches: Recorded path not in target path list c:\invest_test666\cache_dir\n_events10.tif 2022-12-20 12:21:44,033 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-12-20 12:21:44,035 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-12-20 12:21:44,045 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(498) ERROR exception encountered in raster_calculator Traceback (most recent call last): File "pygeoprocessing\geoprocessing.py", line 443, in raster_calculator File "natcap\invest\seasonal_water_yield\seasonal_water_yield.py", line 1093, in cn_op ValueError: Values in the LULC raster were found that are not represented under the 'lucode' key column of the Biophysical table. The missing values found in the LULC raster but not the table are: [0.0]. 2022-12-20 12:21:44,050 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(513) INFO Waiting for raster stats worker result. 2022-12-20 12:21:44,051 (pygeoprocessing.geoprocessing_core) threading.run(953) WARNING No valid pixels were received, sending None. 2022-12-20 12:21:44,051 (taskgraph.Task) Task.add_task(706) ERROR Something went wrong when adding task calculate curve number (13), terminating taskgraph. Traceback (most recent call last): File "taskgraph\Task.py", line 674, in add_task File "taskgraph\Task.py", line 1093, in _call File "natcap\invest\seasonal_water_yield\seasonal_water_yield.py", line 1109, in _calculate_curve_number_raster File "pygeoprocessing\geoprocessing.py", line 443, in raster_calculator File "natcap\invest\seasonal_water_yield\seasonal_water_yield.py", line 1093, in cn_op ValueError: Values in the LULC raster were found that are not represented under the 'lucode' key column of the Biophysical table. The missing values found in the LULC raster but not the table are: [0.0]. 2022-12-20 12:21:44,053 (natcap.invest.utils) utils.prepare_workspace(166) ERROR Exception while executing natcap.invest.seasonal_water_yield.seasonal_water_yield Traceback (most recent call last): File "natcap\invest\utils.py", line 164, in prepare_workspace File "invest\cli.py", line 486, in main File "natcap\invest\seasonal_water_yield\seasonal_water_yield.py", line 408, in execute File "natcap\invest\seasonal_water_yield\seasonal_water_yield.py", line 665, in _execute File "taskgraph\Task.py", line 674, in add_task File "taskgraph\Task.py", line 1093, in _call File "natcap\invest\seasonal_water_yield\seasonal_water_yield.py", line 1109, in _calculate_curve_number_raster File "pygeoprocessing\geoprocessing.py", line 443, in raster_calculator File "natcap\invest\seasonal_water_yield\seasonal_water_yield.py", line 1093, in cn_op ValueError: Values in the LULC raster were found that are not represented under the 'lucode' key column of the Biophysical table. The missing values found in the LULC raster but not the table are: [0.0]. 2022-12-20 12:21:44,058 (natcap.invest.utils) utils.prepare_workspace(169) INFO Elapsed time: 1m 50.72s 2022-12-20 12:21:44,058 (natcap.invest.utils) utils.prepare_workspace(172) INFO Execution finished