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

improved text output and fixed some bugs

Browse Source
This commit is contained in:
Marcel Paffrath 2015-11-12 19:08:33 +01:00
parent d53e4b7c0c
commit db3429e780
1 changed files with 7 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -398,7 +398,8 @@ class SeisArray(object):
z.append(point[2]) z.append(point[2])
return min(z) return min(z)
self.generatePropgrid(nThetaP, nPhiP, nRP, Rbt, cushionpropgrid = 0.05) self.generatePropgrid(nThetaP, nPhiP, nRP, Rbt, cushionfactor = cushionfactor,
cushionpropgrid = 0.05)
surface = self.generateVgrid(nThetaI, nPhiI, nRI, Rbt, method = interpolationMethod, surface = self.generateVgrid(nThetaI, nPhiI, nRI, Rbt, method = interpolationMethod,
cushionfactor = cushionfactor, infilename = customgrid, cushionfactor = cushionfactor, infilename = customgrid,
returnTopo = True) returnTopo = True)
@ -486,7 +487,7 @@ class SeisArray(object):
thetaSN = (theta_min - cushionTheta, theta_max + cushionTheta) thetaSN = (theta_min - cushionTheta, theta_max + cushionTheta)
return thetaSN, phiWE return thetaSN, phiWE
def generatePropgrid(self, nTheta, nPhi, nR, Rbt, cushionpropgrid = 0.05, def generatePropgrid(self, nTheta, nPhi, nR, Rbt, cushionfactor, cushionpropgrid = 0.05,
refinement = (5, 5), outfilename = 'propgrid.in'): refinement = (5, 5), outfilename = 'propgrid.in'):
''' '''
Create a propergation grid file for FMTOMO using SeisArray boundaries Create a propergation grid file for FMTOMO using SeisArray boundaries
@ -503,6 +504,9 @@ class SeisArray(object):
:param: Rbt (bot, top) extensions of the model in m :param: Rbt (bot, top) extensions of the model in m
type: tuple type: tuple
:param: cushionfactor, add some extra space to the model (default: 0.1 = 10%)
type: float
:param: cushionpropogrid, cushionfactor for the propagationgrid (cushion direction :param: cushionpropogrid, cushionfactor for the propagationgrid (cushion direction
opposing to vgrids cushionfactor) opposing to vgrids cushionfactor)
type: float type: float
@ -519,7 +523,7 @@ class SeisArray(object):
print('Bottom of the grid: %s, top of the grid %s' print('Bottom of the grid: %s, top of the grid %s'
%(Rbt[0], Rbt[1])) %(Rbt[0], Rbt[1]))
thetaSN, phiWE = self.getThetaPhiFromArray(cushionfactor = 0) thetaSN, phiWE = self.getThetaPhiFromArray(cushionfactor)
thetaS = thetaSN[0] + cushionpropgrid thetaS = thetaSN[0] + cushionpropgrid
thetaN = thetaSN[1] - cushionpropgrid thetaN = thetaSN[1] - cushionpropgrid