negative values in subwatersheds output for water yield model ?

lyq021lyq021 Member
edited August 2015 in Freshwater Models
Dear all!
 The InVEST model version is 2.6.0(x 64). When I run the water  yield model, the process is sucessful. But when I load the subwatershed.shp file in  ArcGIS 10, there are some negative wield_mn values in several subbasins. I don't know how to solve the error. 
I hope someone can shed some light on this issue. All the thanks.
Below is the log.file:
Arguments:
biophysical_table_uri        D:/worksp/hydrology/data/biophysical_table.csv
depth_to_root_rest_layer_uri D:/worksp/hydrology/data/soil_ref/hdr.adf
eto_uri                      D:/worksp/hydrology/data/eto/hdr.adf
lulc_uri                     D:/worksp/hydrology/data/luc2011/hdr.adf
pawc_uri                     D:/worksp/hydrology/data/paw/hdr.adf
precipitation_uri            D:/worksp/hydrology/data/precp/hdr.adf
results_suffix               
seasonality_constant         6
sub_watersheds_uri           D:/worksp/hydrology/data/subbasin.shp
valuation_container          False
water_scarcity_container     False
watersheds_uri               D:/worksp/hydrology/data/water.shp
workspace_dir                D:\worksp\hydrology


08/22/2015 21:37:54  root               INFO     Logging will be saved to hydropower_water_yield-log-2015-08-22--21_37_54.txt
08/22/2015 21:37:54  root               DEBUG    Loaded the model from invest_natcap.hydropower.hydropower_water_yield
08/22/2015 21:37:54  root               INFO     Executing the loaded model
08/22/2015 21:37:54  root               INFO     Running InVEST version "2.6.0"
08/22/2015 21:37:54  root               INFO     Python architecture: ('64bit', 'WindowsPE')
08/22/2015 21:37:54  root               INFO     Disk space remaining for workspace: 39.4 GB
08/22/2015 21:37:54  root               INFO     Pointing temporary directory at the workspace at D:\worksp\hydrology
08/22/2015 21:37:54  root               INFO     Upddating os.environ["TMP"]=C:\Users\LIUYUQ~1\AppData\Local\Temp to D:\worksp\hydrology
08/22/2015 21:37:54  root               INFO     Upddating os.environ["TEMP"]=C:\Users\LIUYUQ~1\AppData\Local\Temp to D:\worksp\hydrology
08/22/2015 21:37:54  root               INFO     Setting os.environ["TMPDIR"]=D:\worksp\hydrology
08/22/2015 21:37:54  root               INFO     Starting hydropower_water_yield
08/22/2015 21:37:54  hydropower_water_yield INFO     Starting Water Yield Core Calculations
08/22/2015 21:37:54  raster_utils       INFO     clip_dataset nodata value is 255.0
masking out aoi
08/22/2015 21:37:54  hydropower_water_yield INFO     Reclassifying temp_Kc raster
08/22/2015 21:37:54  raster_utils       INFO     Reclassifying
08/22/2015 21:37:54  raster_utils       INFO     Creating lookup numpy array
08/22/2015 21:37:54  raster_utils       INFO     Looping through rows in the input data
08/22/2015 21:37:54  raster_utils       INFO     Flushing the cache and exiting reclassification
08/22/2015 21:37:54  hydropower_water_yield INFO     Reclassifying tmp_root raster
08/22/2015 21:37:54  raster_utils       INFO     Reclassifying
08/22/2015 21:37:54  raster_utils       INFO     Creating lookup numpy array
08/22/2015 21:37:54  raster_utils       INFO     Looping through rows in the input data
08/22/2015 21:37:55  raster_utils       INFO     Flushing the cache and exiting reclassification
08/22/2015 21:37:55  hydropower_water_yield INFO     Reclassifying tmp_veg raster
08/22/2015 21:37:55  raster_utils       INFO     Reclassifying
08/22/2015 21:37:55  raster_utils       INFO     Creating lookup numpy array
08/22/2015 21:37:55  raster_utils       INFO     Looping through rows in the input data
08/22/2015 21:37:55  raster_utils       INFO     Flushing the cache and exiting reclassification
08/22/2015 21:37:55  hydropower_water_yield DEBUG    Calculate PET from Ref Evap times Kc
masking out aoi
08/22/2015 21:37:55  hydropower_water_yield DEBUG    Performing fractp operation
08/22/2015 21:37:55  hydropower_water_yield DEBUG    {'Kc': -1.0, 'soil': -3.4028234663852886e+38, 'pawc': -3.4028234663852886e+38, 'veg': -1.0, 'precip': -3.4028234663852886e+38, 'root': -1.0, 'eto': -32768.0}
masking out aoi
08/22/2015 21:37:56  hydropower_water_yield DEBUG    Performing wyield operation
masking out aoi
08/22/2015 21:37:56  hydropower_water_yield DEBUG    Performing aet operation
masking out aoi
masking out aoi
08/22/2015 21:37:57  raster_utils       WARNING  couldn't remove file d:\worksp\hydrology\tmp\tmpwfjvkh
masking out aoi
08/22/2015 21:37:57  raster_utils       WARNING  couldn't remove file d:\worksp\hydrology\tmp\tmpsu3bwa
masking out aoi
08/22/2015 21:37:57  raster_utils       WARNING  couldn't remove file d:\worksp\hydrology\tmp\tmpg8awfw
masking out aoi
08/22/2015 21:37:58  raster_utils       WARNING  couldn't remove file d:\worksp\hydrology\tmp\tmpvb0o5y
08/22/2015 21:37:58  hydropower_water_yield DEBUG    wyield_value_dict_sws : {0: {'precip_mn': 94.68040903205545, 'AET_mn': 80.09244026572017, 'wyield_ha': 2.016919022e-06, 'PET_mn': 108.18394936141321, 'wyield_mn': 1.709253408141055, 'wyield_vol': 201691.9021606445, 'subws_id': 0, 'num_pixels': 118.0}, 1: {'precip_mn': 91.75475714206695, 'AET_mn': 228.1718223794647, 'wyield_ha': -0.000897001414948, 'PET_mn': 332.12037093928143, 'wyield_mn': -139.28593399806053, 'wyield_vol': -89700141.49475098, 'subws_id': 1, 'num_pixels': 644.0}, 2: {'precip_mn': 110.27273260250783, 'AET_mn': 79.97122088285279, 'wyield_ha': 9.2703211212e-05, 'PET_mn': 164.80644814692423, 'wyield_mn': 19.850794692110963, 'wyield_vol': 9270321.12121582, 'subws_id': 2, 'num_pixels': 467.0}, 3: {'precip_mn': 116.18838654245648, 'AET_mn': 7.134775772302047, 'wyield_ha': 0.00022797421875, 'PET_mn': 39.39676877415698, 'wyield_mn': 99.11922554347827, 'wyield_vol': 22797421.875, 'subws_id': 3, 'num_pixels': 230.0}, 4: {'precip_mn': 119.96211242675781, 'AET_mn': 144.74281734890408, 'wyield_ha': -2.0095611411e-05, 'PET_mn': 326.76690372419944, 'wyield_mn': -24.809396804114918, 'wyield_vol': -2009561.1411333084, 'subws_id': 4, 'num_pixels': 81.0}, 5: {'precip_mn': 117.58750069495476, 'AET_mn': 4.860923087792318, 'wyield_ha': 0.000332700286331, 'PET_mn': 6.746331310663067, 'wyield_mn': 109.08206109218911, 'wyield_vol': 33270028.633117676, 'subws_id': 5, 'num_pixels': 305.0}, 6: {'precip_mn': 116.28059213529757, 'AET_mn': 43.350967341158764, 'wyield_ha': 0.000112190448227, 'PET_mn': 144.28277084913597, 'wyield_mn': 45.056404910413136, 'wyield_vol': 11219044.822692871, 'subws_id': 6, 'num_pixels': 249.0}, 7: {'precip_mn': 114.08815143075407, 'AET_mn': 46.16606831244219, 'wyield_ha': 0.000128144076087, 'PET_mn': 236.18082420096218, 'wyield_mn': 59.88040938667048, 'wyield_vol': 12814407.608747482, 'subws_id': 7, 'num_pixels': 214.0}, 8: {'precip_mn': 117.6692830939042, 'AET_mn': 144.0590070132856, 'wyield_ha': -4.6946582565e-05, 'PET_mn': 560.7017549635872, 'wyield_mn': -24.839461674765936, 'wyield_vol': -4694658.256530762, 'subws_id': 8, 'num_pixels': 189.0}, 9: {'precip_mn': 115.09147755522706, 'AET_mn': 77.04945406978486, 'wyield_ha': 7.9249727303e-05, 'PET_mn': 398.27324930898743, 'wyield_mn': 35.85960511465418, 'wyield_vol': 7924972.7303385725, 'subws_id': 9, 'num_pixels': 221.0}, 10: {'precip_mn': 115.84700024317182, 'AET_mn': 92.98984793693789, 'wyield_ha': 1.4123274048e-05, 'PET_mn': 226.9400642918002, 'wyield_mn': 22.779474270680257, 'wyield_vol': 1412327.4047821758, 'subws_id': 10, 'num_pixels': 62.0}, 11: {'precip_mn': 122.98930208126228, 'AET_mn': 1.818269688658919, 'wyield_ha': 0.000197399085541, 'PET_mn': 7.900521759606578, 'wyield_mn': 121.103733460596, 'wyield_vol': 19739908.55407715, 'subws_id': 11, 'num_pixels': 163.0}, 12: {'precip_mn': 116.20279437807488, 'AET_mn': 10.396053988213414, 'wyield_ha': 0.000688003247528, 'PET_mn': 50.747819513869075, 'wyield_mn': 100.8802415730317, 'wyield_vol': 68800324.75280762, 'subws_id': 12, 'num_pixels': 682.0}, 13: {'precip_mn': 120.07355535986589, 'AET_mn': 5.361086511487112, 'wyield_ha': 0.000215321357465, 'PET_mn': 30.595203548825847, 'wyield_mn': 112.73369500774363, 'wyield_vol': 21532135.746479034, 'subws_id': 13, 'num_pixels': 191.0}, 14: {'precip_mn': 116.12877772404597, 'AET_mn': 56.47128796694326, 'wyield_ha': 6.0724433966e-05, 'PET_mn': 290.1352654964316, 'wyield_mn': 59.533758789885276, 'wyield_vol': 6072443.396568298, 'subws_id': 14, 'num_pixels': 102.0}, 15: {'precip_mn': 117.83882500143612, 'AET_mn': 29.66981084346771, 'wyield_ha': 1.7582225131e-05, 'PET_mn': 210.75364952087403, 'wyield_mn': 87.91112565547228, 'wyield_vol': 1758222.5131094456, 'subws_id': 15, 'num_pixels': 20.0}, 16: {'precip_mn': 122.14829092965999, 'AET_mn': 76.29377536433084, 'wyield_ha': 3.2045594025e-05, 'PET_mn': 346.1887512786048, 'wyield_mn': 45.779420035226, 'wyield_vol': 3204559.4024658203, 'subws_id': 16, 'num_pixels': 70.0}, 17: {'precip_mn': 120.76739959716797, 'AET_mn': 94.28053406552152, 'wyield_ha': 2.9235493511e-05, 'PET_mn': 458.3909723737218, 'wyield_mn': 26.338282442307687, 'wyield_vol': 2923549.3510961533, 'subws_id': 17, 'num_pixels': 111.0}, 18: {'precip_mn': 121.50240071614583, 'AET_mn': 358.61179787772045, 'wyield_ha': -1.6603122976e-05, 'PET_mn': 872.1714564732143, 'wyield_mn': -237.1874710874898, 'wyield_vol': -1660312.2976124284, 'subws_id': 18, 'num_pixels': 7.0}, 19: {'precip_mn': 121.58491473727756, 'AET_mn': 107.75005700704935, 'wyield_ha': 1.9882028053e-05, 'PET_mn': 500.6987369477332, 'wyield_mn': 13.903516120843953, 'wyield_vol': 1988202.8052806852, 'subws_id': 19, 'num_pixels': 143.0}, 20: {'precip_mn': 121.05517578125, 'AET_mn': 120.72503852844238, 'wyield_ha': 1.4885545e-08, 'PET_mn': 838.5, 'wyield_mn': 0.372138634324074, 'wyield_vol': 1488.554537296296, 'subws_id': 20, 'num_pixels': 4.0}, 21: {'precip_mn': 126.12533150010437, 'AET_mn': 39.374562714451045, 'wyield_ha': 0.000223511512146, 'PET_mn': 108.01824212628742, 'wyield_mn': 86.6323690488357, 'wyield_vol': 22351151.21459961, 'subws_id': 21, 'num_pixels': 258.0}, 22: {'precip_mn': 116.7535880055147, 'AET_mn': 62.753835160391674, 'wyield_ha': 0.000230075553856, 'PET_mn': 320.994959738141, 'wyield_mn': 54.77989377535525, 'wyield_vol': 23007555.385649204, 'subws_id': 22, 'num_pixels': 420.0}, 23: {'precip_mn': 116.69160829023882, 'AET_mn': 111.11112003809555, 'wyield_ha': 1.1051798284e-05, 'PET_mn': 458.2039068532935, 'wyield_mn': 4.868633605310045, 'wyield_vol': 1105179.8284053802, 'subws_id': 23, 'num_pixels': 227.0}, 24: {'precip_mn': 115.84453999627497, 'AET_mn': 115.7869582051382, 'wyield_ha': -7.58244606e-07, 'PET_mn': 513.8520931873022, 'wyield_mn': -0.396986704846327, 'wyield_vol': -75824.46062564845, 'subws_id': 24, 'num_pixels': 191.0}, 25: {'precip_mn': 111.32856623331706, 'AET_mn': 100.44997762917448, 'wyield_ha': 2.6674297726e-05, 'PET_mn': 490.30975996680525, 'wyield_mn': 11.448196448993274, 'wyield_vol': 2667429.7726154327, 'subws_id': 25, 'num_pixels': 233.0}, 26: {'precip_mn': 120.70830049415083, 'AET_mn': 86.23941883151814, 'wyield_ha': 0.000101033854554, 'PET_mn': 446.4670914989407, 'wyield_mn': 34.24876425549135, 'wyield_vol': 10103385.45536995, 'subws_id': 26, 'num_pixels': 295.0}, 27: {'precip_mn': 119.70776633060339, 'AET_mn': 102.20001030696257, 'wyield_ha': 4.6236936467e-05, 'PET_mn': 534.5571010895358, 'wyield_mn': 17.64768567431064, 'wyield_vol': 4623693.646669388, 'subws_id': 27, 'num_pixels': 262.0}, 28: {'precip_mn': 128.8796287457423, 'AET_mn': 119.87039529227874, 'wyield_ha': 2.5613790634e-05, 'PET_mn': 517.7974502729697, 'wyield_mn': 9.01894036397128, 'wyield_vol': 2561379.0633678436, 'subws_id': 28, 'num_pixels': 284.0}, 29: {'precip_mn': 135.6611581802368, 'AET_mn': 113.08556454937633, 'wyield_ha': 3.6551443346e-05, 'PET_mn': 458.0802193225884, 'wyield_mn': 22.28746545460166, 'wyield_vol': 3655144.3345546722, 'subws_id': 29, 'num_pixels': 164.0}, 30: {'precip_mn': 132.23133850097656, 'AET_mn': 265.4214274088542, 'wyield_ha': -2.4082667878e-05, 'PET_mn': 894.5166829427084, 'wyield_mn': -133.79259932041168, 'wyield_vol': -2408266.7877674103, 'subws_id': 30, 'num_pixels': 18.0}, 31: {'precip_mn': 133.91978149414064, 'AET_mn': 44.99355355898539, 'wyield_ha': 2.657294312e-06, 'PET_mn': 279.9933348099391, 'wyield_mn': 88.57647705078125, 'wyield_vol': 265729.43115234375, 'subws_id': 31, 'num_pixels': 3.0}, 32: {'precip_mn': 128.64646006266275, 'AET_mn': 91.19583073616027, 'wyield_ha': 0.000109156419015, 'PET_mn': 363.54674392382304, 'wyield_mn': 36.38547300497691, 'wyield_vol': 10915641.901493073, 'subws_id': 32, 'num_pixels': 300.0}, 33: {'precip_mn': 135.75130806717218, 'AET_mn': 49.942105917497116, 'wyield_ha': 4.7116961498e-05, 'PET_mn': 280.808301301436, 'wyield_mn': 85.66720272411, 'wyield_vol': 4711696.14982605, 'subws_id': 33, 'num_pixels': 55.0}, 34: {'precip_mn': 131.7352280970982, 'AET_mn': 117.29601548138787, 'wyield_ha': 1.7322630849e-05, 'PET_mn': 511.8771067899816, 'wyield_mn': 10.189782852285049, 'wyield_vol': 1732263.0848884583, 'subws_id': 34, 'num_pixels': 170.0}, 35: {'precip_mn': 129.91161957383156, 'AET_mn': 103.88904232801524, 'wyield_ha': 0.00012154373076, 'PET_mn': 470.6769936238716, 'wyield_mn': 23.78546590213673, 'wyield_vol': 12154373.075991869, 'subws_id': 35, 'num_pixels': 511.0}, 36: {'precip_mn': 137.4862803674233, 'AET_mn': 225.3450376247162, 'wyield_ha': -0.000383299798794, 'PET_mn': 622.0765672786051, 'wyield_mn': -87.5113695876239, 'wyield_vol': -38329979.87937927, 'subws_id': 36, 'num_pixels': 438.0}, 37: {'precip_mn': 139.29782470985043, 'AET_mn': 170.98465461465432, 'wyield_ha': -0.000148081276399, 'PET_mn': 515.2076809699463, 'wyield_mn': -34.35760473292278, 'wyield_vol': -14808127.639889717, 'subws_id': 37, 'num_pixels': 431.0}, 38: {'precip_mn': 132.28338991940677, 'AET_mn': 114.0924626103154, 'wyield_ha': 5.8728864365e-05, 'PET_mn': 423.7644431785301, 'wyield_mn': 21.751431246157047, 'wyield_vol': 5872886.436462402, 'subws_id': 38, 'num_pixels': 270.0}}
masking out aoi
08/22/2015 21:37:58  raster_utils       WARNING  couldn't remove file d:\worksp\hydrology\tmp\tmppe2d9u
masking out aoi
08/22/2015 21:37:58  raster_utils       WARNING  couldn't remove file d:\worksp\hydrology\tmp\tmp3thybv
masking out aoi
08/22/2015 21:37:59  raster_utils       WARNING  couldn't remove file d:\worksp\hydrology\tmp\tmpwxrk2v
masking out aoi
08/22/2015 21:37:59  raster_utils       WARNING  couldn't remove file d:\worksp\hydrology\tmp\tmplbrikd
08/22/2015 21:37:59  hydropower_water_yield DEBUG    wyield_value_dict_ws : {0: {'precip_mn': 119.97197446123387, 'AET_mn': 97.71238500460495, 'wyield_ha': 0.001769716872678, 'PET_mn': 334.87214666570117, 'wyield_mn': 20.306561935488272, 'wyield_vol': 176971687.2677803, 'num_pixels': 8715.0, 'ws_id': 0}}
08/22/2015 21:37:59  hydropower_water_yield DEBUG    Water Scarcity Not Selected
08/22/2015 21:37:59  root               INFO     Opening file explorer to workspace directory
08/22/2015 21:37:59  root               INFO     Using windows explorer to view files
08/22/2015 21:37:59  root               INFO     Disk space free: 39.4 GB
08/22/2015 21:37:59  root               INFO     Elapsed time: 5.05s
08/22/2015 21:37:59  root               INFO     Finished.
08/22/2015 21:37:59  root               INFO     Elapsed time: 5.47s
08/22/2015 21:37:59  root               INFO     Operations completed successfully
1164 x 494 - 426K
Post edited by lyq021 on

Comments

This discussion has been closed.