PyPSA / pypsa-usa

PyPSA-USA: An Open-Source Energy System Optimization Model for the United States
https://pypsa-usa.readthedocs.io
MIT License
35 stars 15 forks source link

Error in Calculating Eastern Interconnect Renewable Profiles #236

Closed ktehranchi closed 2 months ago

ktehranchi commented 3 months ago

Checklist

The Issue

I'm integrating the eastern interconnect and am running into this issue. I've seen it before when I first was setting up the repo, and am logging this issue for future references if we encounter it again.

Steps To Reproduce

No response

Expected Behavior

No response

Error Message

wildcards: interconnect=eastern, technology=solar
    resources: tmpdir=/var/folders/00/vqry7y9s78q3rvf_m_jmk_bc0000gn/T, mem_mb=40000, mem_mib=38147

INFO:__main__:correction_factor is set as 0.854337
INFO:__main__:Calculate landuse availabilities...
multiprocessing.pool.RemoteTraceback:
"""
Traceback (most recent call last):
  File "/Users/kamrantehranchi/miniforge3/envs/pypsa-usa/lib/python3.11/multiprocessing/pool.py", line 125, in worker
    result = (True, func(*args, **kwds))
                    ^^^^^^^^^^^^^^^^^^^
  File "/Users/kamrantehranchi/miniforge3/envs/pypsa-usa/lib/python3.11/multiprocessing/pool.py", line 48, in mapstar
    return list(map(*args))
           ^^^^^^^^^^^^^^^^
  File "/Users/kamrantehranchi/miniforge3/envs/pypsa-usa/lib/python3.11/site-packages/atlite/gis.py", line 673, in _process_func
    return shape_availability_reprojected(shapes.loc[[i]], *args)[0]
           ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/Users/kamrantehranchi/miniforge3/envs/pypsa-usa/lib/python3.11/site-packages/atlite/gis.py", line 373, in shape_availability_reprojected
    masked, transform = shape_availability(geometry, excluder)
                        ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/Users/kamrantehranchi/miniforge3/envs/pypsa-usa/lib/python3.11/site-packages/atlite/gis.py", line 302, in shape_availability
    transform, shape = padded_transform_and_shape(bounds, res=excluder.res)
                       ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/Users/kamrantehranchi/miniforge3/envs/pypsa-usa/lib/python3.11/site-packages/atlite/gis.py", line 206, in padded_transform_and_shape
    shape = int((top - bottom) / res), int((right - left) / res)
            ^^^^^^^^^^^^^^^^^^^^^^^^^
ValueError: cannot convert float NaN to integer
"""

The above exception was the direct cause of the following exception:

Traceback (most recent call last):
  File "/Users/kamrantehranchi/Local_Documents/pypsa-usa/workflow/.snakemake/scripts/tmprs6ifr7i.build_renewable_profiles.py", line 304, in <module>
    availability = cutout.availabilitymatrix(regions, excluder, **kwargs)
                   ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/Users/kamrantehranchi/miniforge3/envs/pypsa-usa/lib/python3.11/site-packages/atlite/gis.py", line 748, in compute_availabilitymatrix
    availability = list(pool.map(_process_func, shapes.index))
                        ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/Users/kamrantehranchi/miniforge3/envs/pypsa-usa/lib/python3.11/multiprocessing/pool.py", line 367, in map
    return self._map_async(func, iterable, mapstar, chunksize).get()
           ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/Users/kamrantehranchi/miniforge3/envs/pypsa-usa/lib/python3.11/multiprocessing/pool.py", line 774, in get
    raise self._value
  File "/Users/kamrantehranchi/miniforge3/envs/pypsa-usa/lib/python3.11/multiprocessing/pool.py", line 125, in worker
    result = (True, func(*args, **kwds))
^^^^^^^^^^^^^^^
  File "/Users/kamrantehranchi/miniforge3/envs/pypsa-usa/lib/python3.11/multiprocessing/pool.py", line 48, in mapstar
    return list(map(*args))
      ^^^^^^^^^^^^^^^^^
  File "/Users/kamrantehranchi/miniforge3/envs/pypsa-usa/lib/python3.11/site-packages/atlite/gis.py", line 673, in _process_func
    return shape_availability_reprojected(shapes.loc[[i]], *args)[0]
  ^^^^^^^^^^^^^^^^^
  File "/Users/kamrantehranchi/miniforge3/envs/pypsa-usa/lib/python3.11/site-packages/atlite/gis.py", line 373, in shape_availability_reprojected
    masked, transform = shape_availability(geometry, excluder)
      ^^^^^^^^^^^^^^^^^
  File "/Users/kamrantehranchi/miniforge3/envs/pypsa-usa/lib/python3.11/site-packages/atlite/gis.py", line 302, in shape_availability
    transform, shape = padded_transform_and_shape(bounds, res=excluder.res)
      ^^^^^^^^^^^^^^^^^
  File "/Users/kamrantehranchi/miniforge3/envs/pypsa-usa/lib/python3.11/site-packages/atlite/gis.py", line 206, in padded_transform_and_shape
    shape = int((top - bottom) / res), int((right - left) / res)
      ^^^^^^^^^^^^^^^^^
ValueError: cannot convert float NaN to integer

Anything else?

No response