In [1]:

    

%matplotlib qt
sys.path.append('/home/gudni/Dropbox/code/python/MTview/')
import interactivePlotFunctions as iPf


    

    




/Applications/Anaconda/anaconda/lib/python2.7/site-packages/matplotlib/font_manager.py:273: UserWarning: Matplotlib is building the font cache using fc-list. This may take a moment.
  warnings.warn('Matplotlib is building the font cache using fc-list. This may take a moment.')






    




---------------------------------------------------------------------------
NameError                                 Traceback (most recent call last)
<ipython-input-1-ba62965f00ed> in <module>()
      1 get_ipython().magic(u'matplotlib qt')
----> 2 sys.path.append('/home/gudni/Dropbox/code/python/MTview/')
      3 import interactivePlotFunctions as iPf

NameError: name 'sys' is not defined

In [ ]:

    

# Load data


    

In [ ]:

    




    

In [ ]:

    

ls ../ForwardModeling_noExtension_Coarse/*npy


    

In [ ]:

    

gkrMTdata = np.load('../MTdataStArr_nsmesh_0.npy')
thiMTdata = np.load('../ForwardModeling_noExtension_Coarse/MTdataStArr_nsmesh_HKPK1Coarse_noExtension.npy')
gkrNewData = np.load('../ForwardModeling_noExtension_GKR/MTdataStArr_nsmesh_GKRHKPK1.npy')


    

In [2]:

    

iPf.MTinteractiveMap([thiMTdata,gkrNewData])


    

    




---------------------------------------------------------------------------
NameError                                 Traceback (most recent call last)
<ipython-input-2-73cd5d81df99> in <module>()
----> 1 iPf.MTinteractiveMap([thiMTdata,gkrNewData])

NameError: name 'iPf' is not defined

In [25]:

    

np.unique(gkrNewData['freq'])[12::2]


    

    
Out[25]:





array([    100.        ,     215.443469  ,     464.15888336,
          1000.        ,    2154.43469003,    4641.58883361,
         10000.        ,   21544.34690032,   46415.88833613,  100000.        ])

In [15]:

    

drecAll = np.load('../ForwardModeling_noExtension_Coarse/MTdataStArr_nsmesh_HKPK1Coarse_noExtension.npy')
# Select larger frequency band for the MT data
indMTFreq = np.sum( [drecAll['freq'] == val for val in  np.unique(drecAll['freq'])[5::]] ,axis=0,dtype=bool)
mtRecArr = drecAll[indMTFreq][['freq','x','y','z','zxx','zxy','zyx','zyy']]
dUse = NSEM.Data.fromRecArray(mtRecArr)


    

In [17]:

    

np.unique(drecAll['freq'])


    

    
Out[17]:





array([   1000.        ,    1359.35639088,    1847.84979742,
          2511.88643151,    3414.54887383,    4641.58883361,
          6309.5734448 ,    8576.95898591,   11659.1440118 ,
         15848.93192461,   21544.34690032,   29286.44564625,
         39810.71705535,   54116.95265465,   73564.22544596,  100000.        ])

In [26]:

    

4*475*10


    

    
Out[26]:





19000

In [28]:

    

np.unique(gkrNewData[['x','y','z']]).shape


    

    
Out[28]:





(475,)

In [ ]: