In [1]:

    

import pysac
pysac.__version__


    

    
Out[1]:





'0.2.dev202'

In [2]:

    

import numpy as np


    

In [3]:

    

import pysac.yt
import pysac.analysis.tube3D.process_utils as tpu
import pysac.analysis.tube3D.tvtk_tube_functions as ttf


    

    




WARNING:traits.has_traits:DEPRECATED: traits.has_traits.wrapped_class, 'the 'implements' class advisor has been deprecated. Use the 'provides' class decorator.

In [4]:

    

ds = pysac.yt.SACGDFDataset('../thesis/Chapter2/Data/Slog_p180-0_A20-r3_B015_00001.gdf')
cg = ds.index.grids[0]
xmax, ymax, zmax = np.array(cg['mag_field_x'].shape) - 1
domain = {'xmax':xmax, 'ymax':ymax, 'zmax':zmax}


    

    




yt : [INFO     ] 2015-02-09 15:02:09,419 Parameters: current_time              = 1.05812632277
yt : [INFO     ] 2015-02-09 15:02:09,419 Parameters: domain_dimensions         = [128 128 128]
yt : [INFO     ] 2015-02-09 15:02:09,420 Parameters: domain_left_edge          = [  781250.    781250.   3664122.1]
yt : [INFO     ] 2015-02-09 15:02:09,420 Parameters: domain_right_edge         = [  1.99218750e+08   1.99218750e+08   1.58778630e+08]
yt : [INFO     ] 2015-02-09 15:02:09,421 Parameters: cosmological_simulation   = 0.0
yt : [WARNING  ] 2015-02-09 15:02:09,428 Cannot interpret units: 0.001 * , setting to dimensionless.
yt : [WARNING  ] 2015-02-09 15:02:09,428 Cannot interpret units: 0.001 * , setting to dimensionless.
yt : [WARNING  ] 2015-02-09 15:02:09,429 Cannot interpret units: 10 * , setting to dimensionless.
yt : [WARNING  ] 2015-02-09 15:02:09,429 Cannot interpret units: 10 * , setting to dimensionless.
yt : [WARNING  ] 2015-02-09 15:02:09,429 Cannot interpret units: 10000.0 * , setting to dimensionless.
yt : [WARNING  ] 2015-02-09 15:02:09,429 Cannot interpret units: 10000.0 * , setting to dimensionless.
yt : [WARNING  ] 2015-02-09 15:02:09,430 Cannot interpret units: 10000.0 * , setting to dimensionless.
yt : [WARNING  ] 2015-02-09 15:02:09,430 Cannot interpret units: 10000.0 * , setting to dimensionless.
yt : [WARNING  ] 2015-02-09 15:02:09,430 Cannot interpret units: 10000.0 * , setting to dimensionless.
yt : [WARNING  ] 2015-02-09 15:02:09,430 Cannot interpret units: 10000.0 * , setting to dimensionless.
yt : [INFO     ] 2015-02-09 15:02:09,453 Loading field plugins.
yt : [INFO     ] 2015-02-09 15:02:09,454 Loaded angular_momentum (8 new fields)
yt : [INFO     ] 2015-02-09 15:02:09,454 Loaded astro (15 new fields)
yt : [INFO     ] 2015-02-09 15:02:09,454 Loaded cosmology (22 new fields)
yt : [INFO     ] 2015-02-09 15:02:09,455 Loaded fluid (64 new fields)
yt : [INFO     ] 2015-02-09 15:02:09,456 Loaded fluid_vector (96 new fields)
yt : [INFO     ] 2015-02-09 15:02:09,456 Loaded geometric (112 new fields)
yt : [INFO     ] 2015-02-09 15:02:09,457 Loaded local (112 new fields)
yt : [INFO     ] 2015-02-09 15:02:09,457 Loaded magnetic_field (120 new fields)
yt : [INFO     ] 2015-02-09 15:02:09,458 Loaded my_plugins (120 new fields)
yt : [INFO     ] 2015-02-09 15:02:09,458 Loaded species (122 new fields)

In [5]:

    

bfield, vfield = tpu.get_yt_mlab(ds, np.s_[:,:,:], flux=False)


    

In [6]:

    

pd_seeds = ttf.make_circle_seeds(100, 60, **domain)
surf_field_lines, surface = ttf.create_flux_surface(bfield, pd_seeds)


    

In [7]:

    

bfield


    

    
Out[7]:





<mayavi.sources.array_source.ArraySource at 0x7fd16b30e710>

In [9]:

    

surf_field_lines, surface = ttf.create_flux_surface(bfield, pd_seeds)


    

In [ ]:

    

ttf.make_poly_norms()