2023-09-05 16:12:00,960 (natcap.invest.utils) utils.prepare_workspace(161) Level 100 Writing log messages to [/Users/kevinli/Documents/InVEST_PdS/SeasonalWaterYield/InVEST-natcap.invest.seasonal_water_yield.seasonal_water_yield-log-2023-09-05--16_12_00.txt] 2023-09-05 16:12:00,961 (__main__) cli.main(458) Level 100 Starting model with parameters: Arguments for InVEST natcap.invest.seasonal_water_yield.seasonal_water_yield 3.13.0: alpha_m 0.08333333 aoi_path /Users/kevinli/Box Sync/InVEST/cleaned_input/watersheds_PdS_ws.gpkg beta_i 1 biophysical_table_path /Users/kevinli/Box Sync/InVEST/Data requirements/BiophysicalTables for PdS/biophys_swy.csv climate_zone_raster_path climate_zone_table_path dem_raster_path /Users/kevinli/Box Sync/InVEST/cleaned_input/DEM_pdsws.tif et0_dir /Users/kevinli/Box Sync/InVEST/cleaned_input/evap_reproj gamma 1 l_path lulc_raster_path /Users/kevinli/Box Sync/InVEST/Scenarios/flrstrp/flrstrp_fld0.tif monthly_alpha False monthly_alpha_path n_workers -1 precip_dir /Users/kevinli/Box Sync/InVEST/cleaned_input/prec_reproj rain_events_table_path /Users/kevinli/Box Sync/InVEST/cleaned_input/rain_events/PdS_prec_events-mly-prcp-avgnds-ge001hi.csv results_suffix PdS_swy soil_group_path /Users/kevinli/Box Sync/InVEST/cleaned_input/hydrogrp.tif threshold_flow_accumulation 1000 user_defined_climate_zones False user_defined_local_recharge False workspace_dir /Users/kevinli/Documents/InVEST_PdS/SeasonalWaterYield 2023-09-05 16:12:00,961 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(561) INFO prepare and test inputs for common errors 2023-09-05 16:12:01,031 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(627) INFO Building file registry 2023-09-05 16:12:01,031 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(633) INFO Checking that the AOI is not the output aggregate vector 2023-09-05 16:12:01,031 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(642) INFO Aligning and clipping dataset list 2023-09-05 16:12:01,178 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 1 of 27 aligned: prcp_a0_PdS_swy.tif 2023-09-05 16:12:01,188 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 2 of 27 aligned: prcp_a1_PdS_swy.tif 2023-09-05 16:12:01,199 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 3 of 27 aligned: prcp_a2_PdS_swy.tif 2023-09-05 16:12:01,211 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 4 of 27 aligned: prcp_a3_PdS_swy.tif 2023-09-05 16:12:01,220 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 5 of 27 aligned: prcp_a4_PdS_swy.tif 2023-09-05 16:12:01,232 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 6 of 27 aligned: prcp_a5_PdS_swy.tif 2023-09-05 16:12:01,249 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 7 of 27 aligned: prcp_a6_PdS_swy.tif 2023-09-05 16:12:01,260 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 8 of 27 aligned: prcp_a7_PdS_swy.tif 2023-09-05 16:12:01,285 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 9 of 27 aligned: prcp_a8_PdS_swy.tif 2023-09-05 16:12:01,296 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 10 of 27 aligned: prcp_a9_PdS_swy.tif 2023-09-05 16:12:01,306 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 11 of 27 aligned: prcp_a10_PdS_swy.tif 2023-09-05 16:12:01,317 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 12 of 27 aligned: prcp_a11_PdS_swy.tif 2023-09-05 16:12:01,326 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 13 of 27 aligned: soil_group_aligned_PdS_swy.tif 2023-09-05 16:12:01,335 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 14 of 27 aligned: et0_a0_PdS_swy.tif 2023-09-05 16:12:01,344 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 15 of 27 aligned: et0_a1_PdS_swy.tif 2023-09-05 16:12:01,353 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 16 of 27 aligned: et0_a2_PdS_swy.tif 2023-09-05 16:12:01,362 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 17 of 27 aligned: et0_a3_PdS_swy.tif 2023-09-05 16:12:01,372 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 18 of 27 aligned: et0_a4_PdS_swy.tif 2023-09-05 16:12:01,381 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 19 of 27 aligned: et0_a5_PdS_swy.tif 2023-09-05 16:12:01,390 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 20 of 27 aligned: et0_a6_PdS_swy.tif 2023-09-05 16:12:01,399 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 21 of 27 aligned: et0_a7_PdS_swy.tif 2023-09-05 16:12:01,408 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 22 of 27 aligned: et0_a8_PdS_swy.tif 2023-09-05 16:12:01,417 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 23 of 27 aligned: et0_a9_PdS_swy.tif 2023-09-05 16:12:01,428 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 24 of 27 aligned: et0_a10_PdS_swy.tif 2023-09-05 16:12:01,441 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 25 of 27 aligned: et0_a11_PdS_swy.tif 2023-09-05 16:12:01,458 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 26 of 27 aligned: lulc_aligned_PdS_swy.tif 2023-09-05 16:12:01,475 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 27 of 27 aligned: dem_aligned_PdS_swy.tif 2023-09-05 16:12:01,475 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1000) INFO aligned all 27 rasters. 2023-09-05 16:12:01,537 (pygeoprocessing.routing.routing) Task._call(1093) INFO (fill pits): complete 2023-09-05 16:12:01,612 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-09-05 16:12:01,673 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-09-05 16:12:01,698 (pygeoprocessing.routing.routing) Task._call(1093) INFO filter out incomplete divergent streams 2023-09-05 16:12:01,700 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-09-05 16:12:01,701 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(744) INFO quick flow 2023-09-05 16:12:01,701 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(758) INFO loading number of monthly events 2023-09-05 16:12:01,821 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:01,822 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:01,826 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:01,827 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:01,848 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:01,848 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:01,850 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:01,850 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:01,853 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(824) INFO calculate quick flow for month 1 2023-09-05 16:12:01,875 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:01,876 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:01,885 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:01,885 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:01,888 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(824) INFO calculate quick flow for month 2 2023-09-05 16:12:01,918 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:01,918 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:01,925 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:01,925 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:01,929 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(824) INFO calculate quick flow for month 3 2023-09-05 16:12:01,991 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:01,991 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:02,001 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:02,001 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:02,005 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(824) INFO calculate quick flow for month 4 2023-09-05 16:12:02,028 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:02,029 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:02,037 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:02,037 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:02,059 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(824) INFO calculate quick flow for month 5 2023-09-05 16:12:02,142 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:02,142 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:02,151 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:02,151 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:02,157 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(824) INFO calculate quick flow for month 6 2023-09-05 16:12:02,192 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:02,193 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:02,201 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:02,202 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:02,206 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(824) INFO calculate quick flow for month 7 2023-09-05 16:12:02,231 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:02,231 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:02,241 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:02,241 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:02,246 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(824) INFO calculate quick flow for month 8 2023-09-05 16:12:02,267 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:02,267 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:02,277 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:02,278 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:02,281 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(824) INFO calculate quick flow for month 9 2023-09-05 16:12:02,302 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:02,302 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:02,311 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:02,312 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:02,315 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(824) INFO calculate quick flow for month 10 2023-09-05 16:12:02,337 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:02,337 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:02,345 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:02,345 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:02,348 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(824) INFO calculate quick flow for month 11 2023-09-05 16:12:02,369 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:02,369 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:02,376 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:02,377 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:02,380 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(824) INFO calculate quick flow for month 12 2023-09-05 16:12:02,403 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:02,403 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:02,410 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:02,410 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:02,447 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:02,447 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:02,453 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:02,454 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:02,472 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(851) INFO calculate local recharge 2023-09-05 16:12:02,484 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:02,485 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:02,488 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:02,488 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:02,504 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:02,504 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:02,506 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:02,506 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:02,510 (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 (classify kc month 1 (33)) as a dependent task to this one (classify kc month 2 (34)). 2023-09-05 16:12:02,511 (taskgraph.Task) Task.is_precalculated(1252) INFO not precalculated (classify kc month 2 (34)), Task hash exists, but there are these mismatches: Recorded path not in target path list /Users/kevinli/Documents/InVEST_PdS/SeasonalWaterYield/cache_dir/kc_1_PdS_swy.tif 2023-09-05 16:12:02,524 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:02,525 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:02,527 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:02,527 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:02,531 (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 (classify kc month 2 (34)) as a dependent task to this one (classify kc month 3 (34)). 2023-09-05 16:12:02,532 (taskgraph.Task) Task.is_precalculated(1252) INFO not precalculated (classify kc month 3 (34)), Task hash exists, but there are these mismatches: Recorded path not in target path list /Users/kevinli/Documents/InVEST_PdS/SeasonalWaterYield/cache_dir/kc_2_PdS_swy.tif 2023-09-05 16:12:02,544 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:02,544 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:02,546 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:02,547 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:02,562 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:02,563 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:02,565 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:02,565 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:02,586 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:02,587 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:02,589 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:02,589 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:02,606 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:02,606 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:02,609 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:02,609 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:02,624 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:02,624 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:02,626 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:02,626 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:02,642 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:02,642 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:02,645 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:02,645 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:02,666 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:02,666 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:02,668 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:02,668 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:02,692 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:02,692 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:02,694 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:02,694 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:02,698 (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 (classify kc month 0 (32)) as a dependent task to this one (classify kc month 11 (41)). 2023-09-05 16:12:02,699 (taskgraph.Task) Task.is_precalculated(1252) INFO not precalculated (classify kc month 11 (41)), Task hash exists, but there are these mismatches: Recorded path not in target path list /Users/kevinli/Documents/InVEST_PdS/SeasonalWaterYield/cache_dir/kc_0_PdS_swy.tif 2023-09-05 16:12:02,709 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:02,709 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:02,711 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:02,712 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:03,075 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-05 16:12:03,075 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-05 16:12:03,076 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-05 16:12:03,076 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-05 16:12:03,082 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._aggregate_recharge(1318) WARNING /Users/kevinli/Documents/InVEST_PdS/SeasonalWaterYield/aggregated_results_swy_PdS_swy.shp exists, deleting and writing new output 2023-09-05 16:12:03,137 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 1 of 1 aligned: clipped_raster.tif 2023-09-05 16:12:03,137 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1000) INFO aligned all 1 rasters. 2023-09-05 16:12:03,141 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1524) INFO build a lookup of aggregate field value to FID 2023-09-05 16:12:03,142 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1531) INFO creating disjoint polygon set 2023-09-05 16:12:03,143 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2554) INFO build shapely polygon list 2023-09-05 16:12:03,149 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2580) INFO build shapely rtree index 2023-09-05 16:12:03,154 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2593) INFO poly feature lookup 100.0% complete on aggregated_results_swy_PdS_swy.shp 2023-09-05 16:12:03,154 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2597) INFO build poly intersection lookup 2023-09-05 16:12:03,154 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2623) INFO poly intersection feature lookup 100.0% complete on aggregated_results_swy_PdS_swy.shp 2023-09-05 16:12:03,154 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2655) INFO maximal subset build 100.0% complete on aggregated_results_swy_PdS_swy.shp 2023-09-05 16:12:03,159 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1551) INFO processing 1 disjoint polygon sets 2023-09-05 16:12:03,160 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 1 of 1 100.0% processed on aggregated_results_swy_PdS_swy.shp 2023-09-05 16:12:03,160 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 1 of 1 aggregated_results_swy_PdS_swy.shp 2023-09-05 16:12:03,165 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 1 of 1 aggregated_results_swy_PdS_swy.shp 2023-09-05 16:12:03,167 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1756) INFO all done processing polygon sets for aggregated_results_swy_PdS_swy.shp 2023-09-05 16:12:03,191 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 1 of 1 aligned: clipped_raster.tif 2023-09-05 16:12:03,191 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1000) INFO aligned all 1 rasters. 2023-09-05 16:12:03,196 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1524) INFO build a lookup of aggregate field value to FID 2023-09-05 16:12:03,196 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1531) INFO creating disjoint polygon set 2023-09-05 16:12:03,196 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2554) INFO build shapely polygon list 2023-09-05 16:12:03,200 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2580) INFO build shapely rtree index 2023-09-05 16:12:03,201 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2593) INFO poly feature lookup 100.0% complete on aggregated_results_swy_PdS_swy.shp 2023-09-05 16:12:03,202 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2597) INFO build poly intersection lookup 2023-09-05 16:12:03,202 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2623) INFO poly intersection feature lookup 100.0% complete on aggregated_results_swy_PdS_swy.shp 2023-09-05 16:12:03,202 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2655) INFO maximal subset build 100.0% complete on aggregated_results_swy_PdS_swy.shp 2023-09-05 16:12:03,217 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1551) INFO processing 1 disjoint polygon sets 2023-09-05 16:12:03,217 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 1 of 1 100.0% processed on aggregated_results_swy_PdS_swy.shp 2023-09-05 16:12:03,217 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 1 of 1 aggregated_results_swy_PdS_swy.shp 2023-09-05 16:12:03,220 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 1 of 1 aggregated_results_swy_PdS_swy.shp 2023-09-05 16:12:03,222 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1756) INFO all done processing polygon sets for aggregated_results_swy_PdS_swy.shp 2023-09-05 16:12:03,227 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(924) INFO calculate L_sum 2023-09-05 16:12:03,290 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-09-05 16:12:03,395 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(960) INFO (\w/) SWY Complete! 2023-09-05 16:12:03,395 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(961) INFO (.. \ 2023-09-05 16:12:03,395 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(962) INFO _/ ) \______ 2023-09-05 16:12:03,395 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(963) INFO (oo /'\ )`, 2023-09-05 16:12:03,395 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(964) INFO `--' (v __( / || 2023-09-05 16:12:03,395 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(965) INFO ||| ||| || 2023-09-05 16:12:03,395 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(966) INFO //_| //_| 2023-09-05 16:12:03,396 (natcap.invest.utils) utils.prepare_workspace(169) INFO Elapsed time: 2.44s 2023-09-05 16:12:03,396 (natcap.invest.utils) utils.prepare_workspace(172) INFO Execution finished