In [1]:

    

import xray


    

In [2]:

    

url = 'http://geoport.whoi.edu/thredds/dodsC/usgs/data2/emontgomery/stellwagen/CF-1.6/RCNWR/9541aqd-cal.nc'


    

In [3]:

    

ds = xray.open_dataset(url)


    

In [4]:

    

def get_std_name_vars(ds,std_name):
    return {k: v for k, v in ds.data_vars.iteritems() if 'standard_name' in v.attrs.keys() and std_name in v.standard_name}


    

In [5]:

    

std_name = 'eastward_sea_water_velocity'
get_std_name_vars(ds,std_name)


    

    
Out[5]:





{u'u_1205': <xray.DataArray u'u_1205' (time: 17425, z: 8)>
 [139400 values with dtype=float64]
 Coordinates:
   * time       (time) datetime64[ns] 2013-03-27T01:31:00 2013-03-27T01:41:00 ...
     latitude   float64 ...
   * z          (z) float64 -2.16 -1.66 -1.16 -0.6596 -0.1596 0.3404 0.8404 1.34
     longitude  float64 ...
 Attributes:
     epic_code: 1205
     feneric_name: u
     name: u_1205
     standard_name: eastward_sea_water_velocity
     sensor_depth: -2.57962687162
     actual_max: 168.253
     initial_sensor_height: 0.25
     FORTRAN_format:  
     long_name: Eastward (u) velocity
     sensor_type: Nortek Aquadopp Profiler
     actual_min: -63.2215
     units: m/s
     serial_number: AQD 5379
     height_depth_units: m
     grid_mapping: crs
     _ChunkSize: [1000    8]}

In [5]: