# In[87]:


pc = Sp.curtain_zenrad(meas)


# In[90]:


pcn = Sp.curtain_zenrad(meas,utc=False)


# In[89]:


pnc = Sp.plt_curtain_norm_zenrad(meas)


# In[91]:


pncn = Sp.plt_curtain_norm_zenrad(meas,utc=False)


# ## Now plot the measurements for only the selected sample airmass

# In[95]:


airmass = 1./np.cos(s['sza']*np.pi/180.0)
meas.airmass = 1.0/np.cos(meas.sza*np.pi/180.0)