In [ ]:

    




    

In [ ]:

    




    

In [1]:

    

plt.style.use('hagen_default')


    

In [ ]:

    




    

In [18]:

    

from atmPy.data_archives import arm


    

In [125]:

    

reload(arm)
reload(arm._read_data)
reload(arm.aosaps)
reload(arm._tools)


    

    
Out[125]:





<module 'atmPy.data_archives.arm._tools' from '/Users/htelg/prog/atm-py/atmPy/data_archives/arm/_tools.py'>

In [126]:

    

window = ('2016-11-15', '2016-11-26')
fname = '/Volumes/HTelg_4TB_Backup/arm_data/SGP/aosapsE13/'
aosaps = arm.aosaps.open_path(fname, window = window, average=(60, 's'))


    

    




Opening 12 files.
{'suffix': '.nc', 'product': 'aosaps', 'facility': 'E13', 'site': 'sgp', 'timestamp': Timestamp('2016-11-15 21:29:39'), 'qc_level': 'a1'}

In [117]:

    

%matplotlib nbagg


    

In [127]:

    

aosaps.convert2dSdlogDp().plot(norm = 'log')


    

    















    












    
Out[127]:





(<Figure size 800x640 with 2 Axes>,
 <matplotlib.axes._subplots.AxesSubplot at 0x115de9160>,
 <matplotlib.collections.QuadMesh at 0x115e0c978>,
 <matplotlib.colorbar.Colorbar at 0x115e17128>)

In [122]:

    

fname = '/Volumes/HTelg_4TB_Backup/arm_data/SGP/aosapsE13/sgpaosapsE13.a1.20161115.212939.nc'
path = Path(fname)
path.name.split('.')[1]


    

    
Out[122]:





'a1'

In [ ]: