marcel/pylot
3
0
Fork 0
Code Issues 17 Pull Requests 1 Releases Wiki Activity

added minimum cushion for 2D case (for input file generation)

Browse Source
This commit is contained in:
Marcel Paffrath 2016-07-12 14:39:28 +02:00
parent 546dfad722
commit ed85d0ef12
2 changed files with 8 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
pylot/core/active/ActiveSeismoPick3D_GUI.py
View File

@ -424,6 +424,8 @@ class gui_control(object):
cwd = os.getcwd() cwd = os.getcwd()
interpolationMethod = 'linear' interpolationMethod = 'linear'
os.chdir(simuldir) os.chdir(simuldir)
if self.seisarray.twoDim:
interpolationMethod = 'nearest'
self.survey.seisarray.generateFMTOMOinputFromArray(propgrid, vgrid, (bbot, btop), cushionfactor, self.survey.seisarray.generateFMTOMOinputFromArray(propgrid, vgrid, (bbot, btop), cushionfactor,
interpolationMethod, customgrid = customgrid, interpolationMethod, customgrid = customgrid,
writeVTK = False) writeVTK = False)

6
pylot/core/active/seismicArrayPreparation.py
View File

@ -440,6 +440,7 @@ class SeisArray(object):
for phi in phiGrid: for phi in phiGrid:
xval = self._getDistance(phi) xval = self._getDistance(phi)
yval = self._getDistance(theta) yval = self._getDistance(theta)
z = griddata((measured_x, measured_y), measured_z, (xval, yval), method=method) z = griddata((measured_x, measured_y), measured_z, (xval, yval), method=method)
# in case the point lies outside, nan will be returned. Find nearest: # in case the point lies outside, nan will be returned. Find nearest:
if np.isnan(z) == True: if np.isnan(z) == True:
@ -601,6 +602,11 @@ class SeisArray(object):
theta_min, theta_max = (self._getAngle(min(y)), self._getAngle(max(y))) theta_min, theta_max = (self._getAngle(min(y)), self._getAngle(max(y)))
cushionPhi = abs(phi_max - phi_min) * cushionfactor cushionPhi = abs(phi_max - phi_min) * cushionfactor
cushionTheta = abs(theta_max - theta_min) * cushionfactor cushionTheta = abs(theta_max - theta_min) * cushionfactor
# 2D Case only:
if cushionPhi == 0.:
cushionPhi = 0.05
if cushionTheta == 0.:
cushionTheta = 0.05
phiWE = (phi_min - cushionPhi, phi_max + cushionPhi) phiWE = (phi_min - cushionPhi, phi_max + cushionPhi)
thetaSN = (theta_min - cushionTheta, theta_max + cushionTheta) thetaSN = (theta_min - cushionTheta, theta_max + cushionTheta)
return thetaSN, phiWE return thetaSN, phiWE