#!/usr/bin/env python
# coding: utf-8

# In[1]:


import threddsclient
cat = threddsclient.read_url('http://www.esrl.noaa.gov/psd/thredds/catalog.xml')
print cat.name, cat.url
print 'num refs =', len(cat.references)
print 'num datasets =', len(cat.datasets)


# In[2]:


for ref in cat.references:
    print ref.name


# In[3]:


cat2 = cat.references[0].follow()
print cat2.name, cat2.url
print 'num refs =', len(cat2.references)
print 'num datasets =', len(cat2.datasets)


# In[4]:


ds = cat2.datasets[0]
print ds.name, ds.url, ds.is_collection(), len(ds.datasets), len(ds.references)


# In[5]:


for i in range(0, len(ds.references)):
    print i, ds.references[i].name, ds.references[i].url


# In[6]:


cat3 = ds.references[33].follow()
print cat3.name, cat3.url, len(cat3.references), len(cat3.datasets)


# In[7]:


for ref in cat3.flat_references():
    print ref.name, ref.url


# In[8]:


cat4 = cat3.flat_references()[2].follow()
print cat4.name, cat4.url, len(cat4.flat_datasets())


# In[9]:


print '{0.name}, {0.url}, {0.modified}, {0.bytes}'.format(cat4.flat_datasets()[0])


# In[10]:


print cat4.flat_datasets()[0].download_url()


# In[11]:


print cat4.flat_datasets()[0].opendap_url()


# In[12]:


print cat4.download_urls()[0]


# In[13]:


print cat4.opendap_urls()[0]