nsidc / nsidc-data-access-notebook Goto Github PK

capability_url: https://n5eil02u.ecs.nsidc.org/egi/capabilities/ATL03.005.xml
response: <Response [400]> b'{"error":"invalid_request"}'

Traceback (most recent call last):

  File "/srv/conda/envs/notebook/lib/python3.7/site-packages/IPython/core/interactiveshell.py", line 3441, in run_code
    exec(code_obj, self.user_global_ns, self.user_ns)

  File "/tmp/ipykernel_168/938956396.py", line 13, in <module>
    root = ET.fromstring(response.content)

  File "/srv/conda/envs/notebook/lib/python3.7/xml/etree/ElementTree.py", line 1315, in XML
    parser.feed(text)

  File "<string>", line unknown
ParseError: not well-formed (invalid token): line 1, column 0

NameError                                 Traceback (most recent call last)
~\AppData\Local\Temp\ipykernel_15140\1286720228.py in <module>
     22                   'time': time_var,
     23                   'polygon': polygon,
---> 24                   'Boundingshape': Boundingshape,
     25                   'format': reformat,
     26                   'projection': projection,

NameError: name 'Boundingshape' is not defined

aoi = '2'
subdict = {
      'id': 'HEG',
       **'spatialSubsetting': 'true',**
       **'spatialSubsettingShapefile': 'false',**
       **'temporalSubsetting': 'false',**
       'type': 'both',
       'maxGransSyncRequest': '100',
       'maxGransAsyncRequest': '2000'
}

Input short name, e.g. ATL03, here:  MOD10_L2
Input start date in yyyy-MM-dd format:  2022-05-29
Input start time in HH:mm:ss format:  14:00:00
Input end date in yyyy-MM-dd format:  2022-06-22
Input end time in HH:mm:ss format:  16:00:00

# Here I just used the default shapefile that came with the notebook
Simplified polygon coordinates based on shapefile input: -51.1203803853486,69.21453533085385,-51.25214005333067,69.18690191925793,-50.04906889941042,69.08708570400971,-48.98867272012586,68.86377341849372,-48.47406719345924,68.95294818787565,-48.22490237564403,69.07628146713546,-48.64390544471899,69.42180752487434,-49.07439704242353,69.52418890858146,-50.41130570916474,69.24203271223362,-51.1203803853486,69.21453533085385

There are 168 granules of MOD10_L2 version 61 over my area and time of interest.
The average size of each granule is 10.55 MB and the total size of all 168 granules is 1772.72 MB

These reformatting options are available: ['GeoTIFF']
If you would like to reformat, copy and paste the reformatting option you would like (make sure to omit quotes, e.g. GeoTIFF), otherwise leave blank. GeoTIFF

These reprojection options are available with your requested format: ['GEOGRAPHIC', 'UNIVERSAL TRANSVERSE MERCATOR', 'POLAR STEREOGRAPHIC', 'SINUSOIDAL']
If you would like to reproject, copy and paste the reprojection option you would like (make sure to omit quotes), otherwise leave blank. UNIVERSAL TRANSVERSE MERCATOR

Variable subsetting is available. Would you like to subset a selection of variables? (y/n) n

There will be 1 total order(s) processed for our MOD10_L2 request.

---------------------------------------------------------------------------
NameError                                 Traceback (most recent call last)
~\AppData\Local\Temp\ipykernel_5656\1286720228.py in <module>
     22                   'time': time_var,
     23                   'polygon': polygon,
---> 24                   'Boundingshape': Boundingshape,
     25                   'format': reformat,
     26                   'projection': projection,

NameError: name 'Boundingshape' is not defined

ConnectionError: ('Connection aborted.', RemoteDisconnected('Remote end closed connection without response'))

# Query service capability URL 

from xml.etree import ElementTree as ET

capability_url = f'https://n5eil02u.ecs.nsidc.org/egi/capabilities/{short_name}.{latest_version}.xml'

print(capability_url)
# Create session to store cookie and pass credentials to capabilities url

session = requests.session()
s = session.get(capability_url)
response = session.get(s.url,auth=(uid,pswd))

root = ET.fromstring(response.content)