In [1]:
%matplotlib qt
sys.path.append('/home/gudni/Dropbox/code/python/MTview/')
import interactivePlotFunctions as iPf
In [ ]:
# Load data
In [2]:
In [14]:
ls ../ForwardModeling_noExtension_Coarse/*npy
In [8]:
gkrMTdata = np.load('../MTdataStArr_nsmesh_0.npy')
In [9]:
iPf.MTinteractiveMap([gkrMTdata])
Out[9]:
In [ ]:
gkrMTdata
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]:
In [ ]: