wassname/simpeg
1
0
Fork 0
mirror of https://github.com/wassname/simpeg.git synced 2026-07-13 17:45:30 +08:00
Code Issues Packages Projects Releases Wiki Activity

Updated DataMT while working 1D inversions.

Browse Source
This commit is contained in:
GudniRos
2015-08-05 10:44:45 -07:00
parent ba1a85e269
commit 295eeb6a25
2 changed files with 2199 additions and 182 deletions
Download Patch File Download Diff File Expand all files Collapse all files
notebooks/scipy2015/SeogiModelMT3Dfor21Dinv/MT3DforData1Dinv.ipynb
+2193 -176
View File
File diff suppressed because it is too large Load Diff
simpegMT/DataMT.py
+6 -6
View File
@@ -91,7 +91,7 @@ class DataMT(Survey.Data):
if srcType=='primary':
src = simpegMT.SurveyMT.srcMT_polxy_1Dprimary
elif srcType=='total':
src = simpegMT.SurveyMT.srcMT_polxy_1DhomotD
src = sdsimpegMT.SurveyMT.srcMT_polxy_1DhomotD
else:
raise NotImplementedError('{:s} is not a valid source type for MTdata')
@@ -105,14 +105,14 @@ class DataMT(Survey.Data):
# Find that data for freq
dFreq = recArray[recArray['freq'] == freq]
# Find the impedance rxTypes in the recArray.
rxTypes = [ comp for comp in recArray.dtype.names if len(comp)==4 and 'z' in comp and 'r' in comp or 'i' in comp]
rxTypes = [ comp for comp in recArray.dtype.names if len(comp)==4 and 'z' in comp and ('r' in comp or 'i' in comp)]
for rxType in rxTypes:
# Find index of not nan values in rxType
notNaNind = ~np.isnan(dFreq[rxType])
locs = rec2ndarr(dFreq[['x','y','z']][notNaNind].copy())
rxList.append(simpegMT.SurveyMT.RxMT(locs,rxType))
dataList.append(dFreq[rxType][notNaNind].data)
if np.any(notNaNind): # Make sure that there is any data to add.
locs = rec2ndarr(dFreq[['x','y','z']][notNaNind].copy())
rxList.append(simpegMT.SurveyMT.RxMT(locs,rxType))
dataList.append(dFreq[rxType][notNaNind])
srcList.append(src(rxList,freq))
# Make a survey