import yt

data='/share/Part1/guido/cflows_gg/w_turbulence3/cflows_gg_hdf5_plt_cnt_0000'

ds      =yt.load(data)
# bs 256 pc
c_1=[0.25,0.5,0.5]
c_2=[0.75,0.5,0.5]

normal=[1,0,0]

radius=(32,'pc')
height=(32,'pc')


c_1=ds.disk(c_1,normal,radius,height)
c_2=ds.disk(c_2,normal,radius,height)


x_c1=c_1['x']
x_c2=c_2['x']