wassname/simpeg
1
0
Fork 0
mirror of https://github.com/wassname/simpeg.git synced 2026-07-11 04:59:44 +08:00
Code Issues Packages Projects Releases Wiki Activity

First commit for GRAV3D

Browse Source 
Add example for mag
Add example for grav
This commit is contained in:
D Fournier
2016-03-16 17:10:30 -07:00
parent c70256d028
commit 03dd9915b7
22 changed files with 200383 additions and 18 deletions
Download Patch File Download Diff File Expand all files Collapse all files
simpegPF/BaseGrav.py
+52
View File
@@ -0,0 +1,52 @@
from SimPEG import Maps, Survey, Utils, np, sp
from scipy.constants import mu_0
import re
class LinearSurvey(Survey.BaseSurvey):
"""Base Magnetics Survey"""
rxLoc = None #: receiver locations
rxType = None #: receiver type
def __init__(self, srcField, **kwargs):
self.srcField = srcField
Survey.BaseSurvey.__init__(self, **kwargs)
def eval(self, u):
return u
@property
def nD(self):
return self.prob.G.shape[0]
@property
def nRx(self):
return self.srcField.rxList[0].locs.shape[0]
# def setBackgroundField(self, SrcField):
# if getattr(self, 'B0', None) is None:
# self._B0 = SrcField.param[0] * dipazm_2_xyz(SrcField.param[1],SrcField.param[2])
# return self._B0
class SrcField(Survey.BaseSrc):
""" Define the inducing field """
param = None #: Inducing field param (Amp, Incl, Decl)
def __init__(self, rxList, **kwargs):
super(SrcField, self).__init__(rxList, **kwargs)
class RxObs(Survey.BaseRx):
"""A station location must have be located in 3-D"""
def __init__(self, locsXYZ, **kwargs):
locs = locsXYZ
assert locsXYZ.shape[1] == 3, 'locs must in 3-D (x,y,z).'
super(RxObs, self).__init__(locs, 'tmi', storeProjections=False, **kwargs)
@property
def nD(self):
"""Number of data in the receiver."""
return self.locs[0].shape[0]
simpegPF/Dev/GRAV/Block_sphere_model.den
+44180
View File
File diff suppressed because it is too large Load Diff
simpegPF/Dev/Gaussian.topo → simpegPF/Dev/GRAV/Gaussian.topo
View File
simpegPF/Dev/GRAV/Intgrl_GRAV_Inv_Driver.py
+247
View File
@@ -0,0 +1,247 @@
#%%
from SimPEG import *
import simpegPF as PF
import pylab as plt
import os
#home_dir = 'C:\Users\dominiquef.MIRAGEOSCIENCE\Documents\GIT\SimPEG\simpegpf\simpegPF\Dev'
#home_dir = 'C:\\Users\\dominiquef.MIRAGEOSCIENCE\\ownCloud\\Research\\Modelling\\Synthetic\\Nut_Cracker\\Induced_MAG3C'
#home_dir = 'C:\\Users\\dominiquef.MIRAGEOSCIENCE\\ownCloud\\Research\\Modelling\\Synthetic\\Block_Gaussian_topo\\GRAV'
#home_dir ='C:\\Users\\dominiquef.MIRAGEOSCIENCE\\ownCloud\\Research\\CraigModel'
home_dir = '.\\'
inpfile = 'PYGRAV3D_inv.inp'
dsep = '\\'
os.chdir(home_dir)
## New scripts to be added to basecode
#from fwr_MAG_data import fwr_MAG_data
#from read_MAGfwr_inp import read_MAGfwr_inp
#%%
# Read input file
[mshfile, obsfile, topofile, mstart, mref, wgtfile, chi, alphas, bounds, lpnorms] = PF.Gravity.read_GRAVinv_inp(home_dir + dsep + inpfile)
# Load mesh file
mesh = Mesh.TensorMesh.readUBC(mshfile)
#mesh = Utils.meshutils.readUBCTensorMesh(mshfile)
# Load in observation file
survey = PF.Gravity.readUBCgravObs(obsfile)
rxLoc = survey.srcField.rxList[0].locs
d = survey.dobs
wd = survey.std
ndata = survey.srcField.rxList[0].locs.shape[0]
beta_in = 1e+1
# Load in topofile or create flat surface
if topofile == 'null':
# All active
actv = np.asarray(range(mesh.nC))
else:
topo = np.genfromtxt(topofile,skip_header=1)
# Find the active cells
actv = PF.Magnetics.getActiveTopo(mesh,topo,'N')
nC = len(actv)
# Create active map to go from reduce set to full
actvMap = Maps.ActiveCells(mesh, actv, -100)
# Creat reduced identity map
idenMap = Maps.IdentityMap(nP = nC)
# Load starting model file
if isinstance(mstart, float):
mstart = np.ones(nC) * mstart
else:
mstart = Mesh.TensorMesh.readModelUBC(mesh,mstart)
mstart = mstart[actv]
# Load reference file
if isinstance(mref, float):
mref = np.ones(nC) * mref
else:
mref = Mesh.TensorMesh.readModelUBC(mesh,mref)
mref = mref[actv]
# Get index of the center
midx = int(mesh.nCx/2)
midy = int(mesh.nCy/2)
#%% Plot obs data
PF.Gravity.plot_obs_2D(survey,'Observed Data')
#%% Run inversion
prob = PF.Gravity.GravityIntegral(mesh, mapping = idenMap, actInd = actv)
prob.solverOpts['accuracyTol'] = 1e-4
#survey = Survey.LinearSurvey()
survey.pair(prob)
#survey.makeSyntheticData(data, std=0.01)
#survey.dobs=d
#survey.mtrue = model
# Write out the predicted
pred = prob.fields(mstart)
PF.Gravity.writeUBCobs(home_dir + dsep + 'Pred0.dat',survey,pred)
wr = np.sum(prob.G**2.,axis=0)**0.5
wr = ( wr/np.max(wr) )
wr_out = actvMap * wr
plt.figure()
ax = plt.subplot()
mesh.plotSlice(actvMap*mstart, ax = ax, normal = 'Y', ind=midx+1 ,clim = (-1e-1, mstart.max()))
plt.title('Distance weighting')
plt.xlabel('x');plt.ylabel('z')
plt.gca().set_aspect('equal', adjustable='box')
reg = Regularization.Simple(mesh, indActive = actv, mapping = idenMap)
reg.mref = mref
reg.wght = wr
#reg.alpha_s = 1.
# Create pre-conditioner
diagA = np.sum(prob.G**2.,axis=0) + beta_in*(reg.W.T*reg.W).diagonal()*wr
PC = Utils.sdiag(diagA**-1.)
dmis = DataMisfit.l2_DataMisfit(survey)
dmis.Wd = 1/wd
opt = Optimization.ProjectedGNCG(maxIter=10,lower=0.,upper=1., maxIterCG= 20, tolCG = 1e-3)
opt.approxHinv = PC
# opt = Optimization.InexactGaussNewton(maxIter=6)
invProb = InvProblem.BaseInvProblem(dmis, reg, opt, beta = beta_in)
beta = Directives.BetaSchedule(coolingFactor=2, coolingRate=1)
#betaest = Directives.BetaEstimate_ByEig()
target = Directives.TargetMisfit()
inv = Inversion.BaseInversion(invProb, directiveList=[beta,target])
m0 = mstart
# Run inversion
mrec = inv.run(m0)
m_out = actvMap*mrec
# Write result
Mesh.TensorMesh.writeModelUBC(mesh,'SimPEG_inv_l2l2.sus',m_out)
#Utils.meshutils.writeUBCTensorModel(home_dir+dsep+'wr.dat',mesh,wr_out)
# Plot predicted
pred = prob.fields(mrec)
#PF.Magnetics.plot_obs_2D(rxLoc,pred,wd,'Predicted Data')
#PF.Magnetics.plot_obs_2D(rxLoc,(d-pred),wd,'Residual Data')
print "Final misfit:" + str(np.sum( ((d-pred)/wd)**2. ) )
#%% Plot out a section of the model
yslice = midx+1
plt.figure()
ax = plt.subplot(221)
mesh.plotSlice(m_out, ax = ax, normal = 'Z', ind=-5, clim = (-1e-3, mrec.max()))
plt.plot(np.array([mesh.vectorCCx[0],mesh.vectorCCx[-1]]), np.array([mesh.vectorCCy[yslice],mesh.vectorCCy[yslice]]),c='w',linestyle = '--')
plt.title('Z: ' + str(mesh.vectorCCz[-5]) + ' m')
plt.xlabel('x');plt.ylabel('z')
plt.gca().set_aspect('equal', adjustable='box')
ax = plt.subplot(222)
mesh.plotSlice(m_out, ax = ax, normal = 'Z', ind=-8, clim = (-1e-3, mrec.max()))
plt.plot(np.array([mesh.vectorCCx[0],mesh.vectorCCx[-1]]), np.array([mesh.vectorCCy[yslice],mesh.vectorCCy[yslice]]),c='w',linestyle = '--')
plt.title('Z: ' + str(mesh.vectorCCz[-8]) + ' m')
plt.xlabel('x');plt.ylabel('z')
plt.gca().set_aspect('equal', adjustable='box')
ax = plt.subplot(212)
mesh.plotSlice(m_out, ax = ax, normal = 'Y', ind=yslice, clim = (-1e-3, mrec.max()))
plt.title('Cross Section')
plt.xlabel('x');plt.ylabel('z')
plt.gca().set_aspect('equal', adjustable='box')
#%% Run one more round for sparsity
phim = invProb.phi_m_last
phid = invProb.phi_d
reg = Regularization.Sparse(mesh, indActive = actv, mapping = idenMap)
reg.recModel = mrec
reg.mref = mref
reg.wght = wr
reg.eps = 1e-2
reg.p = lpnorms[0]
reg.qx = lpnorms[1]
reg.qz = lpnorms[2]
reg.qy = lpnorms[3]
diagA = np.sum(prob.G**2.,axis=0) + beta_in*(reg.W.T*reg.W).diagonal()*(wr)
PC = Utils.sdiag(diagA**-1.)
#reg.alpha_s = 1.
dmis = DataMisfit.l2_DataMisfit(survey)
dmis.Wd = 1/wd
opt = Optimization.ProjectedGNCG(maxIter=10 ,lower=0.,upper=1., maxIterCG= 10, tolCG = 1e-4)
opt.approxHinv = PC
#opt.phim_last = reg.eval(mrec)
# opt = Optimization.InexactGaussNewton(maxIter=6)
invProb = InvProblem.BaseInvProblem(dmis, reg, opt, beta = invProb.beta)
beta = Directives.BetaSchedule(coolingFactor=1, coolingRate=1)
#betaest = Directives.BetaEstimate_ByEig()
target = Directives.TargetMisfit()
IRLS =Directives.update_IRLS( phi_m_last = phim, phi_d_last = phid )
inv = Inversion.BaseInversion(invProb, directiveList=[beta,IRLS])
m0 = mrec
# Run inversion
mrec = inv.run(m0)
m_out = actvMap*mrec
Mesh.TensorMesh.writeModelUBC(mesh,'SimPEG_inv_l0l2.sus',m_out)
pred = prob.fields(mrec)
#%% Plot obs data
PF.Magnetics.plot_obs_2D(rxLoc,pred,wd,'Predicted Data')
PF.Magnetics.plot_obs_2D(rxLoc,d,wd,'Observed Data')
print "Final misfit:" + str(np.sum( ((d-pred)/wd)**2. ) )
#%% Plot out a section of the model
yslice = midx
plt.figure()
ax = plt.subplot(221)
mesh.plotSlice(m_out, ax = ax, normal = 'Z', ind=-5, clim = (-1e-2, mrec.max()))
plt.plot(np.array([mesh.vectorCCx[0],mesh.vectorCCx[-1]]), np.array([mesh.vectorCCy[yslice],mesh.vectorCCy[yslice]]),c='w',linestyle = '--')
plt.title('Z: ' + str(mesh.vectorCCz[-5]) + ' m')
plt.xlabel('x');plt.ylabel('z')
plt.gca().set_aspect('equal', adjustable='box')
ax = plt.subplot(222)
mesh.plotSlice(m_out, ax = ax, normal = 'Z', ind=-8, clim = (-1e-2, mrec.max()))
plt.plot(np.array([mesh.vectorCCx[0],mesh.vectorCCx[-1]]), np.array([mesh.vectorCCy[yslice],mesh.vectorCCy[yslice]]),c='w',linestyle = '--')
plt.title('Z: ' + str(mesh.vectorCCz[-8]) + ' m')
plt.xlabel('x');plt.ylabel('z')
plt.gca().set_aspect('equal', adjustable='box')
ax = plt.subplot(212)
mesh.plotSlice(m_out, ax = ax, normal = 'Y', ind=yslice, clim = (-1e-2, mrec.max()))
plt.title('Cross Section')
plt.xlabel('x');plt.ylabel('z')
plt.gca().set_aspect('equal', adjustable='box')
simpegPF/Dev/Mesh_10m.msh → simpegPF/Dev/GRAV/Mesh_10m.msh
View File
simpegPF/Dev/GRAV/PYGRAV3D_inv.inp
+10
View File
@@ -0,0 +1,10 @@
Mesh_10m.msh ! Mesh file
gzfor3d.dat ! Obsfile
Gaussian.topo ! Topofile | null
VALUE 1e-4 ! Starting model
VALUE 0 ! Reference model
DEFAULT ! Cell based weight file
1 ! target chi factor | DEFAULT=1
1 1 1 1 ! alpha s, x ,y ,z
VALUE 0 1 ! Lower and Upper Bounds for p-component
VALUE 0 2 2 2 1 ! lp-norm for amplitude inversion FILE pqxqyqzr.dat ! Norms VALUE p, qx, qy, qz, r | FILE m-by-5 matrix
simpegPF/Dev/GRAV/Pred0.dat
+617
View File
@@ -0,0 +1,617 @@
616
1.050000e+02 1.550000e+02 2.408312e+02 5.429956e-04 4.000000e-03
1.150000e+02 1.550000e+02 2.440469e+02 5.498174e-04 4.000000e-03
1.250000e+02 1.550000e+02 2.471986e+02 5.575820e-04 4.000000e-03
1.350000e+02 1.550000e+02 2.502569e+02 5.657694e-04 4.000000e-03
1.450000e+02 1.550000e+02 2.531923e+02 5.709679e-04 4.000000e-03
1.550000e+02 1.550000e+02 2.559753e+02 5.756866e-04 4.000000e-03
1.650000e+02 1.550000e+02 2.585770e+02 5.814230e-04 4.000000e-03
1.750000e+02 1.550000e+02 2.609695e+02 5.864902e-04 4.000000e-03
1.850000e+02 1.550000e+02 2.631267e+02 5.894365e-04 4.000000e-03
1.950000e+02 1.550000e+02 2.650244e+02 5.913321e-04 4.000000e-03
2.050000e+02 1.550000e+02 2.666411e+02 5.932002e-04 4.000000e-03
2.150000e+02 1.550000e+02 2.679578e+02 5.951067e-04 4.000000e-03
2.250000e+02 1.550000e+02 2.689593e+02 5.957423e-04 4.000000e-03
2.350000e+02 1.550000e+02 2.696337e+02 5.955739e-04 4.000000e-03
2.450000e+02 1.550000e+02 2.699729e+02 5.951335e-04 4.000000e-03
2.550000e+02 1.550000e+02 2.699729e+02 5.945539e-04 4.000000e-03
2.650000e+02 1.550000e+02 2.696337e+02 5.937014e-04 4.000000e-03
2.750000e+02 1.550000e+02 2.689593e+02 5.920908e-04 4.000000e-03
2.850000e+02 1.550000e+02 2.679578e+02 5.893730e-04 4.000000e-03
2.950000e+02 1.550000e+02 2.666411e+02 5.866030e-04 4.000000e-03
3.050000e+02 1.550000e+02 2.650244e+02 5.837291e-04 4.000000e-03
3.150000e+02 1.550000e+02 2.631267e+02 5.800023e-04 4.000000e-03
3.250000e+02 1.550000e+02 2.609695e+02 5.748185e-04 4.000000e-03
3.350000e+02 1.550000e+02 2.585770e+02 5.689802e-04 4.000000e-03
3.450000e+02 1.550000e+02 2.559753e+02 5.635854e-04 4.000000e-03
3.550000e+02 1.550000e+02 2.531923e+02 5.578471e-04 4.000000e-03
3.650000e+02 1.550000e+02 2.502569e+02 5.502964e-04 4.000000e-03
3.750000e+02 1.550000e+02 2.471986e+02 5.425413e-04 4.000000e-03
1.050000e+02 1.650000e+02 2.427601e+02 5.467004e-04 4.000000e-03
1.150000e+02 1.650000e+02 2.460869e+02 5.538664e-04 4.000000e-03
1.250000e+02 1.650000e+02 2.493475e+02 5.621375e-04 4.000000e-03
1.350000e+02 1.650000e+02 2.525115e+02 5.692413e-04 4.000000e-03
1.450000e+02 1.650000e+02 2.555483e+02 5.748121e-04 4.000000e-03
1.550000e+02 1.650000e+02 2.584275e+02 5.811887e-04 4.000000e-03
1.650000e+02 1.650000e+02 2.611190e+02 5.866881e-04 4.000000e-03
1.750000e+02 1.650000e+02 2.635942e+02 5.901621e-04 4.000000e-03
1.850000e+02 1.650000e+02 2.658260e+02 5.932671e-04 4.000000e-03
1.950000e+02 1.650000e+02 2.677892e+02 5.964597e-04 4.000000e-03
2.050000e+02 1.650000e+02 2.694617e+02 5.990576e-04 4.000000e-03
2.150000e+02 1.650000e+02 2.708240e+02 6.013334e-04 4.000000e-03
2.250000e+02 1.650000e+02 2.718600e+02 6.018728e-04 4.000000e-03
2.350000e+02 1.650000e+02 2.725577e+02 6.015521e-04 4.000000e-03
2.450000e+02 1.650000e+02 2.729086e+02 6.010120e-04 4.000000e-03
2.550000e+02 1.650000e+02 2.729086e+02 6.003679e-04 4.000000e-03
2.650000e+02 1.650000e+02 2.725577e+02 5.994891e-04 4.000000e-03
2.750000e+02 1.650000e+02 2.718600e+02 5.978165e-04 4.000000e-03
2.850000e+02 1.650000e+02 2.708240e+02 5.947749e-04 4.000000e-03
2.950000e+02 1.650000e+02 2.694617e+02 5.914120e-04 4.000000e-03
3.050000e+02 1.650000e+02 2.677892e+02 5.875223e-04 4.000000e-03
3.150000e+02 1.650000e+02 2.658260e+02 5.835428e-04 4.000000e-03
3.250000e+02 1.650000e+02 2.635942e+02 5.792868e-04 4.000000e-03
3.350000e+02 1.650000e+02 2.611190e+02 5.737387e-04 4.000000e-03
3.450000e+02 1.650000e+02 2.584275e+02 5.673629e-04 4.000000e-03
3.550000e+02 1.650000e+02 2.555483e+02 5.612567e-04 4.000000e-03
3.650000e+02 1.650000e+02 2.525115e+02 5.543355e-04 4.000000e-03
3.750000e+02 1.650000e+02 2.493475e+02 5.464295e-04 4.000000e-03
1.050000e+02 1.750000e+02 2.445174e+02 5.499201e-04 4.000000e-03
1.150000e+02 1.750000e+02 2.479454e+02 5.580225e-04 4.000000e-03
1.250000e+02 1.750000e+02 2.513052e+02 5.662716e-04 4.000000e-03
1.350000e+02 1.750000e+02 2.545655e+02 5.722191e-04 4.000000e-03
1.450000e+02 1.750000e+02 2.576947e+02 5.788615e-04 4.000000e-03
1.550000e+02 1.750000e+02 2.606614e+02 5.857153e-04 4.000000e-03
1.650000e+02 1.750000e+02 2.634348e+02 5.897775e-04 4.000000e-03
1.750000e+02 1.750000e+02 2.659854e+02 5.934387e-04 4.000000e-03
1.850000e+02 1.750000e+02 2.682850e+02 5.982697e-04 4.000000e-03
1.950000e+02 1.750000e+02 2.703080e+02 6.022473e-04 4.000000e-03
2.050000e+02 1.750000e+02 2.720313e+02 6.038864e-04 4.000000e-03
2.150000e+02 1.750000e+02 2.734350e+02 6.046617e-04 4.000000e-03
2.250000e+02 1.750000e+02 2.745027e+02 6.048116e-04 4.000000e-03
2.350000e+02 1.750000e+02 2.752216e+02 6.046296e-04 4.000000e-03
2.450000e+02 1.750000e+02 2.755831e+02 6.042147e-04 4.000000e-03
2.550000e+02 1.750000e+02 2.755831e+02 6.034214e-04 4.000000e-03
2.650000e+02 1.750000e+02 2.752216e+02 6.023779e-04 4.000000e-03
2.750000e+02 1.750000e+02 2.745027e+02 6.010023e-04 4.000000e-03
2.850000e+02 1.750000e+02 2.734350e+02 5.990618e-04 4.000000e-03
2.950000e+02 1.750000e+02 2.720313e+02 5.963515e-04 4.000000e-03
3.050000e+02 1.750000e+02 2.703080e+02 5.919350e-04 4.000000e-03
3.150000e+02 1.750000e+02 2.682850e+02 5.867276e-04 4.000000e-03
3.250000e+02 1.750000e+02 2.659854e+02 5.821963e-04 4.000000e-03
3.350000e+02 1.750000e+02 2.634348e+02 5.774408e-04 4.000000e-03
3.450000e+02 1.750000e+02 2.606614e+02 5.709620e-04 4.000000e-03
3.550000e+02 1.750000e+02 2.576947e+02 5.642088e-04 4.000000e-03
3.650000e+02 1.750000e+02 2.545655e+02 5.577908e-04 4.000000e-03
3.750000e+02 1.750000e+02 2.513052e+02 5.500977e-04 4.000000e-03
1.050000e+02 1.850000e+02 2.460835e+02 5.531692e-04 4.000000e-03
1.150000e+02 1.850000e+02 2.496018e+02 5.617843e-04 4.000000e-03
1.250000e+02 1.850000e+02 2.530500e+02 5.691570e-04 4.000000e-03
1.350000e+02 1.850000e+02 2.563960e+02 5.753053e-04 4.000000e-03
1.450000e+02 1.850000e+02 2.596076e+02 5.824057e-04 4.000000e-03
1.550000e+02 1.850000e+02 2.626524e+02 5.881723e-04 4.000000e-03
1.650000e+02 1.850000e+02 2.654988e+02 5.919221e-04 4.000000e-03
1.750000e+02 1.850000e+02 2.681164e+02 5.968115e-04 4.000000e-03
1.850000e+02 1.850000e+02 2.704766e+02 6.023523e-04 4.000000e-03
1.950000e+02 1.850000e+02 2.725528e+02 6.050219e-04 4.000000e-03
2.050000e+02 1.850000e+02 2.743215e+02 6.059367e-04 4.000000e-03
2.150000e+02 1.850000e+02 2.757621e+02 6.064240e-04 4.000000e-03
2.250000e+02 1.850000e+02 2.768578e+02 6.069147e-04 4.000000e-03
2.350000e+02 1.850000e+02 2.775956e+02 6.075274e-04 4.000000e-03
2.450000e+02 1.850000e+02 2.779668e+02 6.077914e-04 4.000000e-03
2.550000e+02 1.850000e+02 2.779668e+02 6.063119e-04 4.000000e-03
2.650000e+02 1.850000e+02 2.775956e+02 6.045111e-04 4.000000e-03
2.750000e+02 1.850000e+02 2.768578e+02 6.027972e-04 4.000000e-03
2.850000e+02 1.850000e+02 2.757621e+02 6.010481e-04 4.000000e-03
2.950000e+02 1.850000e+02 2.743215e+02 5.988891e-04 4.000000e-03
3.050000e+02 1.850000e+02 2.725528e+02 5.953109e-04 4.000000e-03
3.150000e+02 1.850000e+02 2.704766e+02 5.897272e-04 4.000000e-03
3.250000e+02 1.850000e+02 2.681164e+02 5.843633e-04 4.000000e-03
3.350000e+02 1.850000e+02 2.654988e+02 5.797596e-04 4.000000e-03
3.450000e+02 1.850000e+02 2.626524e+02 5.739973e-04 4.000000e-03
3.550000e+02 1.850000e+02 2.596076e+02 5.671189e-04 4.000000e-03
3.650000e+02 1.850000e+02 2.563960e+02 5.605218e-04 4.000000e-03
3.750000e+02 1.850000e+02 2.530500e+02 5.533262e-04 4.000000e-03
1.050000e+02 1.950000e+02 2.474409e+02 5.565420e-04 4.000000e-03
1.150000e+02 1.950000e+02 2.510374e+02 5.653731e-04 4.000000e-03
1.250000e+02 1.950000e+02 2.545621e+02 5.717217e-04 4.000000e-03
1.350000e+02 1.950000e+02 2.579825e+02 5.787385e-04 4.000000e-03
1.450000e+02 1.950000e+02 2.612655e+02 5.857708e-04 4.000000e-03
1.550000e+02 1.950000e+02 2.643779e+02 5.901023e-04 4.000000e-03
1.650000e+02 1.950000e+02 2.672876e+02 5.942135e-04 4.000000e-03
1.750000e+02 1.950000e+02 2.699634e+02 5.996769e-04 4.000000e-03
1.850000e+02 1.950000e+02 2.723760e+02 6.041994e-04 4.000000e-03
1.950000e+02 1.950000e+02 2.744984e+02 6.061458e-04 4.000000e-03
2.050000e+02 1.950000e+02 2.763063e+02 6.072375e-04 4.000000e-03
2.150000e+02 1.950000e+02 2.777790e+02 6.087786e-04 4.000000e-03
2.250000e+02 1.950000e+02 2.788990e+02 6.107456e-04 4.000000e-03
2.350000e+02 1.950000e+02 2.796533e+02 6.123173e-04 4.000000e-03
2.450000e+02 1.950000e+02 2.800326e+02 6.130981e-04 4.000000e-03
2.550000e+02 1.950000e+02 2.800326e+02 6.109806e-04 4.000000e-03
2.650000e+02 1.950000e+02 2.796533e+02 6.082019e-04 4.000000e-03
2.750000e+02 1.950000e+02 2.788990e+02 6.051279e-04 4.000000e-03
2.850000e+02 1.950000e+02 2.777790e+02 6.024151e-04 4.000000e-03
2.950000e+02 1.950000e+02 2.763063e+02 6.000594e-04 4.000000e-03
3.050000e+02 1.950000e+02 2.744984e+02 5.970238e-04 4.000000e-03
3.150000e+02 1.950000e+02 2.723760e+02 5.921762e-04 4.000000e-03
3.250000e+02 1.950000e+02 2.699634e+02 5.865202e-04 4.000000e-03
3.350000e+02 1.950000e+02 2.672876e+02 5.816262e-04 4.000000e-03
3.450000e+02 1.950000e+02 2.643779e+02 5.766444e-04 4.000000e-03
3.550000e+02 1.950000e+02 2.612655e+02 5.700371e-04 4.000000e-03
3.650000e+02 1.950000e+02 2.579825e+02 5.629817e-04 4.000000e-03
3.750000e+02 1.950000e+02 2.545621e+02 5.561657e-04 4.000000e-03
1.050000e+02 2.050000e+02 2.485739e+02 5.588184e-04 4.000000e-03
1.150000e+02 2.050000e+02 2.522356e+02 5.675440e-04 4.000000e-03
1.250000e+02 2.050000e+02 2.558243e+02 5.743326e-04 4.000000e-03
1.350000e+02 2.050000e+02 2.593068e+02 5.819145e-04 4.000000e-03
1.450000e+02 2.050000e+02 2.626493e+02 5.879312e-04 4.000000e-03
1.550000e+02 2.050000e+02 2.658182e+02 5.918407e-04 4.000000e-03
1.650000e+02 2.050000e+02 2.687807e+02 5.968891e-04 4.000000e-03
1.750000e+02 2.050000e+02 2.715050e+02 6.023919e-04 4.000000e-03
1.850000e+02 2.050000e+02 2.739614e+02 6.054107e-04 4.000000e-03
1.950000e+02 2.050000e+02 2.761223e+02 6.071607e-04 4.000000e-03
2.050000e+02 2.050000e+02 2.779631e+02 6.093873e-04 4.000000e-03
2.150000e+02 2.050000e+02 2.794624e+02 6.127621e-04 4.000000e-03
2.250000e+02 2.050000e+02 2.806028e+02 6.158952e-04 4.000000e-03
2.350000e+02 2.050000e+02 2.813707e+02 6.171642e-04 4.000000e-03
2.450000e+02 2.050000e+02 2.817570e+02 6.171369e-04 4.000000e-03
2.550000e+02 2.050000e+02 2.817570e+02 6.156456e-04 4.000000e-03
2.650000e+02 2.050000e+02 2.813707e+02 6.129636e-04 4.000000e-03
2.750000e+02 2.050000e+02 2.806028e+02 6.087952e-04 4.000000e-03
2.850000e+02 2.050000e+02 2.794624e+02 6.045031e-04 4.000000e-03
2.950000e+02 2.050000e+02 2.779631e+02 6.011401e-04 4.000000e-03
3.050000e+02 2.050000e+02 2.761223e+02 5.981751e-04 4.000000e-03
3.150000e+02 2.050000e+02 2.739614e+02 5.942954e-04 4.000000e-03
3.250000e+02 2.050000e+02 2.715050e+02 5.887768e-04 4.000000e-03
3.350000e+02 2.050000e+02 2.687807e+02 5.833256e-04 4.000000e-03
3.450000e+02 2.050000e+02 2.658182e+02 5.785639e-04 4.000000e-03
3.550000e+02 2.050000e+02 2.626493e+02 5.726284e-04 4.000000e-03
3.650000e+02 2.050000e+02 2.593068e+02 5.653867e-04 4.000000e-03
3.750000e+02 2.050000e+02 2.558243e+02 5.582196e-04 4.000000e-03
1.050000e+02 2.150000e+02 2.494692e+02 5.611836e-04 4.000000e-03
1.150000e+02 2.150000e+02 2.531825e+02 5.692806e-04 4.000000e-03
1.250000e+02 2.150000e+02 2.568218e+02 5.771409e-04 4.000000e-03
1.350000e+02 2.150000e+02 2.603533e+02 5.850442e-04 4.000000e-03
1.450000e+02 2.150000e+02 2.637429e+02 5.895888e-04 4.000000e-03
1.550000e+02 2.150000e+02 2.669565e+02 5.937853e-04 4.000000e-03
1.650000e+02 2.150000e+02 2.699606e+02 5.994047e-04 4.000000e-03
1.750000e+02 2.150000e+02 2.727234e+02 6.039962e-04 4.000000e-03
1.850000e+02 2.150000e+02 2.752143e+02 6.063251e-04 4.000000e-03
1.950000e+02 2.150000e+02 2.774056e+02 6.088183e-04 4.000000e-03
2.050000e+02 2.150000e+02 2.792724e+02 6.128416e-04 4.000000e-03
2.150000e+02 2.150000e+02 2.807928e+02 6.168619e-04 4.000000e-03
2.250000e+02 2.150000e+02 2.819493e+02 6.189064e-04 4.000000e-03
2.350000e+02 2.150000e+02 2.827280e+02 6.193568e-04 4.000000e-03
2.450000e+02 2.150000e+02 2.831197e+02 6.189143e-04 4.000000e-03
2.550000e+02 2.150000e+02 2.831197e+02 6.177902e-04 4.000000e-03
2.650000e+02 2.150000e+02 2.827280e+02 6.158064e-04 4.000000e-03
2.750000e+02 2.150000e+02 2.819493e+02 6.124216e-04 4.000000e-03
2.850000e+02 2.150000e+02 2.807928e+02 6.075425e-04 4.000000e-03
2.950000e+02 2.150000e+02 2.792724e+02 6.027165e-04 4.000000e-03
3.050000e+02 2.150000e+02 2.774056e+02 5.990841e-04 4.000000e-03
3.150000e+02 2.150000e+02 2.752143e+02 5.956646e-04 4.000000e-03
3.250000e+02 2.150000e+02 2.727234e+02 5.908023e-04 4.000000e-03
3.350000e+02 2.150000e+02 2.699606e+02 5.850802e-04 4.000000e-03
3.450000e+02 2.150000e+02 2.669565e+02 5.800663e-04 4.000000e-03
3.550000e+02 2.150000e+02 2.637429e+02 5.748501e-04 4.000000e-03
3.650000e+02 2.150000e+02 2.603533e+02 5.676490e-04 4.000000e-03
3.750000e+02 2.150000e+02 2.568218e+02 5.598991e-04 4.000000e-03
1.050000e+02 2.250000e+02 2.501165e+02 5.640607e-04 4.000000e-03
1.150000e+02 2.250000e+02 2.538671e+02 5.707407e-04 4.000000e-03
1.250000e+02 2.250000e+02 2.575429e+02 5.786016e-04 4.000000e-03
1.350000e+02 2.250000e+02 2.611099e+02 5.863602e-04 4.000000e-03
1.450000e+02 2.250000e+02 2.645335e+02 5.906318e-04 4.000000e-03
1.550000e+02 2.250000e+02 2.677793e+02 5.958751e-04 4.000000e-03
1.650000e+02 2.250000e+02 2.708137e+02 6.018782e-04 4.000000e-03
1.750000e+02 2.250000e+02 2.736041e+02 6.050130e-04 4.000000e-03
1.850000e+02 2.250000e+02 2.761201e+02 6.070007e-04 4.000000e-03
1.950000e+02 2.250000e+02 2.783334e+02 6.107352e-04 4.000000e-03
2.050000e+02 2.250000e+02 2.802189e+02 6.160540e-04 4.000000e-03
2.150000e+02 2.250000e+02 2.817546e+02 6.190299e-04 4.000000e-03
2.250000e+02 2.250000e+02 2.829227e+02 6.199593e-04 4.000000e-03
2.350000e+02 2.250000e+02 2.837092e+02 6.198498e-04 4.000000e-03
2.450000e+02 2.250000e+02 2.841048e+02 6.192416e-04 4.000000e-03
2.550000e+02 2.250000e+02 2.841048e+02 6.182828e-04 4.000000e-03
2.650000e+02 2.250000e+02 2.837092e+02 6.168222e-04 4.000000e-03
2.750000e+02 2.250000e+02 2.829227e+02 6.143722e-04 4.000000e-03
2.850000e+02 2.250000e+02 2.817546e+02 6.101772e-04 4.000000e-03
2.950000e+02 2.250000e+02 2.802189e+02 6.043558e-04 4.000000e-03
3.050000e+02 2.250000e+02 2.783334e+02 5.997473e-04 4.000000e-03
3.150000e+02 2.250000e+02 2.761201e+02 5.965115e-04 4.000000e-03
3.250000e+02 2.250000e+02 2.736041e+02 5.925221e-04 4.000000e-03
3.350000e+02 2.250000e+02 2.708137e+02 5.867084e-04 4.000000e-03
3.450000e+02 2.250000e+02 2.677793e+02 5.810503e-04 4.000000e-03
3.550000e+02 2.250000e+02 2.645335e+02 5.759395e-04 4.000000e-03
3.650000e+02 2.250000e+02 2.611099e+02 5.688848e-04 4.000000e-03
3.750000e+02 2.250000e+02 2.575429e+02 5.611484e-04 4.000000e-03
1.050000e+02 2.350000e+02 2.505080e+02 5.652222e-04 4.000000e-03
1.150000e+02 2.350000e+02 2.542811e+02 5.713655e-04 4.000000e-03
1.250000e+02 2.350000e+02 2.579790e+02 5.789862e-04 4.000000e-03
1.350000e+02 2.350000e+02 2.615674e+02 5.865473e-04 4.000000e-03
1.450000e+02 2.350000e+02 2.650117e+02 5.909287e-04 4.000000e-03
1.550000e+02 2.350000e+02 2.682770e+02 5.966606e-04 4.000000e-03
1.650000e+02 2.350000e+02 2.713296e+02 6.026911e-04 4.000000e-03
1.750000e+02 2.350000e+02 2.741368e+02 6.052751e-04 4.000000e-03
1.850000e+02 2.350000e+02 2.766679e+02 6.071619e-04 4.000000e-03
1.950000e+02 2.350000e+02 2.788945e+02 6.113895e-04 4.000000e-03
2.050000e+02 2.350000e+02 2.807913e+02 6.169624e-04 4.000000e-03
2.150000e+02 2.350000e+02 2.823363e+02 6.194572e-04 4.000000e-03
2.250000e+02 2.350000e+02 2.835114e+02 6.199276e-04 4.000000e-03
2.350000e+02 2.350000e+02 2.843026e+02 6.195898e-04 4.000000e-03
2.450000e+02 2.350000e+02 2.847007e+02 6.189162e-04 4.000000e-03
2.550000e+02 2.350000e+02 2.847007e+02 6.180257e-04 4.000000e-03
2.650000e+02 2.350000e+02 2.843026e+02 6.167881e-04 4.000000e-03
2.750000e+02 2.350000e+02 2.835114e+02 6.147503e-04 4.000000e-03
2.850000e+02 2.350000e+02 2.823363e+02 6.109355e-04 4.000000e-03
2.950000e+02 2.350000e+02 2.807913e+02 6.048990e-04 4.000000e-03
3.050000e+02 2.350000e+02 2.788945e+02 5.998911e-04 4.000000e-03
3.150000e+02 2.350000e+02 2.766679e+02 5.967026e-04 4.000000e-03
3.250000e+02 2.350000e+02 2.741368e+02 5.930923e-04 4.000000e-03
3.350000e+02 2.350000e+02 2.713296e+02 5.873126e-04 4.000000e-03
3.450000e+02 2.350000e+02 2.682770e+02 5.813427e-04 4.000000e-03
3.550000e+02 2.350000e+02 2.650117e+02 5.761443e-04 4.000000e-03
3.650000e+02 2.350000e+02 2.615674e+02 5.692119e-04 4.000000e-03
3.750000e+02 2.350000e+02 2.579791e+02 5.616526e-04 4.000000e-03
1.050000e+02 2.450000e+02 2.506390e+02 5.649581e-04 4.000000e-03
1.150000e+02 2.450000e+02 2.544197e+02 5.711114e-04 4.000000e-03
1.250000e+02 2.450000e+02 2.581250e+02 5.787055e-04 4.000000e-03
1.350000e+02 2.450000e+02 2.617206e+02 5.861887e-04 4.000000e-03
1.450000e+02 2.450000e+02 2.651717e+02 5.905714e-04 4.000000e-03
1.550000e+02 2.450000e+02 2.684435e+02 5.962678e-04 4.000000e-03
1.650000e+02 2.450000e+02 2.715022e+02 6.022046e-04 4.000000e-03
1.750000e+02 2.450000e+02 2.743151e+02 6.047792e-04 4.000000e-03
1.850000e+02 2.450000e+02 2.768512e+02 6.066780e-04 4.000000e-03
1.950000e+02 2.450000e+02 2.790823e+02 6.108633e-04 4.000000e-03
2.050000e+02 2.450000e+02 2.809829e+02 6.163223e-04 4.000000e-03
2.150000e+02 2.450000e+02 2.825310e+02 6.187874e-04 4.000000e-03
2.250000e+02 2.450000e+02 2.837084e+02 6.192588e-04 4.000000e-03
2.350000e+02 2.450000e+02 2.845013e+02 6.189239e-04 4.000000e-03
2.450000e+02 2.450000e+02 2.849001e+02 6.182517e-04 4.000000e-03
2.550000e+02 2.450000e+02 2.849001e+02 6.173600e-04 4.000000e-03
2.650000e+02 2.450000e+02 2.845013e+02 6.161196e-04 4.000000e-03
2.750000e+02 2.450000e+02 2.837084e+02 6.140847e-04 4.000000e-03
2.850000e+02 2.450000e+02 2.825310e+02 6.103024e-04 4.000000e-03
2.950000e+02 2.450000e+02 2.809829e+02 6.043597e-04 4.000000e-03
3.050000e+02 2.450000e+02 2.790823e+02 5.993960e-04 4.000000e-03
3.150000e+02 2.450000e+02 2.768512e+02 5.962072e-04 4.000000e-03
3.250000e+02 2.450000e+02 2.743151e+02 5.926122e-04 4.000000e-03
3.350000e+02 2.450000e+02 2.715022e+02 5.869006e-04 4.000000e-03
3.450000e+02 2.450000e+02 2.684435e+02 5.809731e-04 4.000000e-03
3.550000e+02 2.450000e+02 2.651717e+02 5.757873e-04 4.000000e-03
3.650000e+02 2.450000e+02 2.617206e+02 5.689081e-04 4.000000e-03
3.750000e+02 2.450000e+02 2.581250e+02 5.613858e-04 4.000000e-03
1.050000e+02 2.550000e+02 2.505080e+02 5.633083e-04 4.000000e-03
1.150000e+02 2.550000e+02 2.542811e+02 5.699913e-04 4.000000e-03
1.250000e+02 2.550000e+02 2.579790e+02 5.777635e-04 4.000000e-03
1.350000e+02 2.550000e+02 2.615674e+02 5.852805e-04 4.000000e-03
1.450000e+02 2.550000e+02 2.650117e+02 5.895554e-04 4.000000e-03
1.550000e+02 2.550000e+02 2.682770e+02 5.947215e-04 4.000000e-03
1.650000e+02 2.550000e+02 2.713296e+02 6.004685e-04 4.000000e-03
1.750000e+02 2.550000e+02 2.741368e+02 6.035409e-04 4.000000e-03
1.850000e+02 2.550000e+02 2.766679e+02 6.055483e-04 4.000000e-03
1.950000e+02 2.550000e+02 2.788945e+02 6.091852e-04 4.000000e-03
2.050000e+02 2.550000e+02 2.807913e+02 6.141943e-04 4.000000e-03
2.150000e+02 2.550000e+02 2.823363e+02 6.170474e-04 4.000000e-03
2.250000e+02 2.550000e+02 2.835114e+02 6.179555e-04 4.000000e-03
2.350000e+02 2.550000e+02 2.843026e+02 6.178487e-04 4.000000e-03
2.450000e+02 2.550000e+02 2.847007e+02 6.172422e-04 4.000000e-03
2.550000e+02 2.550000e+02 2.847007e+02 6.162818e-04 4.000000e-03
2.650000e+02 2.550000e+02 2.843026e+02 6.148231e-04 4.000000e-03
2.750000e+02 2.550000e+02 2.835114e+02 6.124040e-04 4.000000e-03
2.850000e+02 2.550000e+02 2.823363e+02 6.083360e-04 4.000000e-03
2.950000e+02 2.550000e+02 2.807913e+02 6.027706e-04 4.000000e-03
3.050000e+02 2.550000e+02 2.788945e+02 5.982683e-04 4.000000e-03
3.150000e+02 2.550000e+02 2.766679e+02 5.950419e-04 4.000000e-03
3.250000e+02 2.550000e+02 2.741368e+02 5.911251e-04 4.000000e-03
3.350000e+02 2.550000e+02 2.713296e+02 5.855016e-04 4.000000e-03
3.450000e+02 2.550000e+02 2.682770e+02 5.799453e-04 4.000000e-03
3.550000e+02 2.550000e+02 2.650117e+02 5.748704e-04 4.000000e-03
3.650000e+02 2.550000e+02 2.615674e+02 5.679810e-04 4.000000e-03
3.750000e+02 2.550000e+02 2.579791e+02 5.603613e-04 4.000000e-03
1.050000e+02 2.650000e+02 2.501165e+02 5.601579e-04 4.000000e-03
1.150000e+02 2.650000e+02 2.538671e+02 5.680712e-04 4.000000e-03
1.250000e+02 2.650000e+02 2.575429e+02 5.757842e-04 4.000000e-03
1.350000e+02 2.650000e+02 2.611099e+02 5.832999e-04 4.000000e-03
1.450000e+02 2.650000e+02 2.645335e+02 5.878014e-04 4.000000e-03
1.550000e+02 2.650000e+02 2.677793e+02 5.919907e-04 4.000000e-03
1.650000e+02 2.650000e+02 2.708137e+02 5.973301e-04 4.000000e-03
1.750000e+02 2.650000e+02 2.736041e+02 6.016067e-04 4.000000e-03
1.850000e+02 2.650000e+02 2.761201e+02 6.039017e-04 4.000000e-03
1.950000e+02 2.650000e+02 2.783334e+02 6.063698e-04 4.000000e-03
2.050000e+02 2.650000e+02 2.802189e+02 6.101235e-04 4.000000e-03
2.150000e+02 2.650000e+02 2.817546e+02 6.137454e-04 4.000000e-03
2.250000e+02 2.650000e+02 2.829227e+02 6.156175e-04 4.000000e-03
2.350000e+02 2.650000e+02 2.837092e+02 6.160280e-04 4.000000e-03
2.450000e+02 2.650000e+02 2.841048e+02 6.155805e-04 4.000000e-03
2.550000e+02 2.650000e+02 2.841048e+02 6.144664e-04 4.000000e-03
2.650000e+02 2.650000e+02 2.837092e+02 6.125318e-04 4.000000e-03
2.750000e+02 2.650000e+02 2.829227e+02 6.093205e-04 4.000000e-03
2.850000e+02 2.650000e+02 2.817546e+02 6.047823e-04 4.000000e-03
2.950000e+02 2.650000e+02 2.802189e+02 6.002061e-04 4.000000e-03
3.050000e+02 2.650000e+02 2.783334e+02 5.966362e-04 4.000000e-03
3.150000e+02 2.650000e+02 2.761201e+02 5.932768e-04 4.000000e-03
3.250000e+02 2.650000e+02 2.736041e+02 5.886569e-04 4.000000e-03
3.350000e+02 2.650000e+02 2.708137e+02 5.831866e-04 4.000000e-03
3.450000e+02 2.650000e+02 2.677793e+02 5.782525e-04 4.000000e-03
3.550000e+02 2.650000e+02 2.645335e+02 5.731215e-04 4.000000e-03
3.650000e+02 2.650000e+02 2.611099e+02 5.661918e-04 4.000000e-03
3.750000e+02 2.650000e+02 2.575429e+02 5.586299e-04 4.000000e-03
1.050000e+02 2.750000e+02 2.494692e+02 5.574935e-04 4.000000e-03
1.150000e+02 2.750000e+02 2.531825e+02 5.658536e-04 4.000000e-03
1.250000e+02 2.750000e+02 2.568218e+02 5.725907e-04 4.000000e-03
1.350000e+02 2.750000e+02 2.603533e+02 5.798230e-04 4.000000e-03
1.450000e+02 2.750000e+02 2.637429e+02 5.854812e-04 4.000000e-03
1.550000e+02 2.750000e+02 2.669565e+02 5.893845e-04 4.000000e-03
1.650000e+02 2.750000e+02 2.699606e+02 5.942035e-04 4.000000e-03
1.750000e+02 2.750000e+02 2.727234e+02 5.991691e-04 4.000000e-03
1.850000e+02 2.750000e+02 2.752143e+02 6.020174e-04 4.000000e-03
1.950000e+02 2.750000e+02 2.774056e+02 6.037791e-04 4.000000e-03
2.050000e+02 2.750000e+02 2.792724e+02 6.059423e-04 4.000000e-03
2.150000e+02 2.750000e+02 2.807928e+02 6.089528e-04 4.000000e-03
2.250000e+02 2.750000e+02 2.819493e+02 6.115935e-04 4.000000e-03
2.350000e+02 2.750000e+02 2.827280e+02 6.126836e-04 4.000000e-03
2.450000e+02 2.750000e+02 2.831197e+02 6.125799e-04 4.000000e-03
2.550000e+02 2.750000e+02 2.831197e+02 6.111748e-04 4.000000e-03
2.650000e+02 2.750000e+02 2.827280e+02 6.086782e-04 4.000000e-03
2.750000e+02 2.750000e+02 2.819493e+02 6.049410e-04 4.000000e-03
2.850000e+02 2.750000e+02 2.807928e+02 6.009787e-04 4.000000e-03
2.950000e+02 2.750000e+02 2.792724e+02 5.977084e-04 4.000000e-03
3.050000e+02 2.750000e+02 2.774056e+02 5.947750e-04 4.000000e-03
3.150000e+02 2.750000e+02 2.752143e+02 5.910739e-04 4.000000e-03
3.250000e+02 2.750000e+02 2.727234e+02 5.859694e-04 4.000000e-03
3.350000e+02 2.750000e+02 2.699606e+02 5.807728e-04 4.000000e-03
3.450000e+02 2.750000e+02 2.669565e+02 5.760928e-04 4.000000e-03
3.550000e+02 2.750000e+02 2.637429e+02 5.704210e-04 4.000000e-03
3.650000e+02 2.750000e+02 2.603533e+02 5.634999e-04 4.000000e-03
3.750000e+02 2.750000e+02 2.568218e+02 5.564841e-04 4.000000e-03
1.050000e+02 2.850000e+02 2.485739e+02 5.548906e-04 4.000000e-03
1.150000e+02 2.850000e+02 2.522356e+02 5.632369e-04 4.000000e-03
1.250000e+02 2.850000e+02 2.558243e+02 5.695373e-04 4.000000e-03
1.350000e+02 2.850000e+02 2.593068e+02 5.762976e-04 4.000000e-03
1.450000e+02 2.850000e+02 2.626493e+02 5.827335e-04 4.000000e-03
1.550000e+02 2.850000e+02 2.658182e+02 5.869501e-04 4.000000e-03
1.650000e+02 2.850000e+02 2.687807e+02 5.909947e-04 4.000000e-03
1.750000e+02 2.850000e+02 2.715050e+02 5.959483e-04 4.000000e-03
1.850000e+02 2.850000e+02 2.739614e+02 5.998862e-04 4.000000e-03
1.950000e+02 2.850000e+02 2.761223e+02 6.017368e-04 4.000000e-03
2.050000e+02 2.850000e+02 2.779631e+02 6.028749e-04 4.000000e-03
2.150000e+02 2.850000e+02 2.794624e+02 6.043483e-04 4.000000e-03
2.250000e+02 2.850000e+02 2.806028e+02 6.060550e-04 4.000000e-03
2.350000e+02 2.850000e+02 2.813707e+02 6.073290e-04 4.000000e-03
2.450000e+02 2.850000e+02 2.817570e+02 6.077357e-04 4.000000e-03
2.550000e+02 2.850000e+02 2.817570e+02 6.059738e-04 4.000000e-03
2.650000e+02 2.850000e+02 2.813707e+02 6.034653e-04 4.000000e-03
2.750000e+02 2.850000e+02 2.806028e+02 6.006274e-04 4.000000e-03
2.850000e+02 2.850000e+02 2.794624e+02 5.979992e-04 4.000000e-03
2.950000e+02 2.850000e+02 2.779631e+02 5.956369e-04 4.000000e-03
3.050000e+02 2.850000e+02 2.761223e+02 5.927212e-04 4.000000e-03
3.150000e+02 2.850000e+02 2.739614e+02 5.883355e-04 4.000000e-03
3.250000e+02 2.850000e+02 2.715050e+02 5.831390e-04 4.000000e-03
3.350000e+02 2.850000e+02 2.687807e+02 5.784084e-04 4.000000e-03
3.450000e+02 2.850000e+02 2.658182e+02 5.735911e-04 4.000000e-03
3.550000e+02 2.850000e+02 2.626493e+02 5.674213e-04 4.000000e-03
3.650000e+02 2.850000e+02 2.593068e+02 5.606748e-04 4.000000e-03
3.750000e+02 2.850000e+02 2.558243e+02 5.540031e-04 4.000000e-03
1.050000e+02 2.950000e+02 2.474409e+02 5.513553e-04 4.000000e-03
1.150000e+02 2.950000e+02 2.510374e+02 5.595548e-04 4.000000e-03
1.250000e+02 2.950000e+02 2.545621e+02 5.665269e-04 4.000000e-03
1.350000e+02 2.950000e+02 2.579825e+02 5.725830e-04 4.000000e-03
1.450000e+02 2.950000e+02 2.612655e+02 5.791413e-04 4.000000e-03
1.550000e+02 2.950000e+02 2.643779e+02 5.843586e-04 4.000000e-03
1.650000e+02 2.950000e+02 2.672876e+02 5.880684e-04 4.000000e-03
1.750000e+02 2.950000e+02 2.699634e+02 5.925860e-04 4.000000e-03
1.850000e+02 2.950000e+02 2.723760e+02 5.972158e-04 4.000000e-03
1.950000e+02 2.950000e+02 2.744984e+02 5.995951e-04 4.000000e-03
2.050000e+02 2.950000e+02 2.763063e+02 6.005279e-04 4.000000e-03
2.150000e+02 2.950000e+02 2.777790e+02 6.010715e-04 4.000000e-03
2.250000e+02 2.950000e+02 2.788990e+02 6.015764e-04 4.000000e-03
2.350000e+02 2.950000e+02 2.796533e+02 6.020888e-04 4.000000e-03
2.450000e+02 2.950000e+02 2.800326e+02 6.021598e-04 4.000000e-03
2.550000e+02 2.950000e+02 2.800326e+02 6.008473e-04 4.000000e-03
2.650000e+02 2.950000e+02 2.796533e+02 5.991307e-04 4.000000e-03
2.750000e+02 2.950000e+02 2.788990e+02 5.974100e-04 4.000000e-03
2.850000e+02 2.950000e+02 2.777790e+02 5.956308e-04 4.000000e-03
2.950000e+02 2.950000e+02 2.763063e+02 5.934962e-04 4.000000e-03
3.050000e+02 2.950000e+02 2.744984e+02 5.902128e-04 4.000000e-03
3.150000e+02 2.950000e+02 2.723760e+02 5.853408e-04 4.000000e-03
3.250000e+02 2.950000e+02 2.699634e+02 5.803802e-04 4.000000e-03
3.350000e+02 2.950000e+02 2.672876e+02 5.759133e-04 4.000000e-03
3.450000e+02 2.950000e+02 2.643779e+02 5.705687e-04 4.000000e-03
3.550000e+02 2.950000e+02 2.612655e+02 5.641798e-04 4.000000e-03
3.650000e+02 2.950000e+02 2.579825e+02 5.578032e-04 4.000000e-03
3.750000e+02 2.950000e+02 2.545621e+02 5.509221e-04 4.000000e-03
1.050000e+02 3.050000e+02 2.460835e+02 5.478507e-04 4.000000e-03
1.150000e+02 3.050000e+02 2.496018e+02 5.556758e-04 4.000000e-03
1.250000e+02 3.050000e+02 2.530500e+02 5.632535e-04 4.000000e-03
1.350000e+02 3.050000e+02 2.563960e+02 5.690940e-04 4.000000e-03
1.450000e+02 3.050000e+02 2.596076e+02 5.753502e-04 4.000000e-03
1.550000e+02 3.050000e+02 2.626524e+02 5.813195e-04 4.000000e-03
1.650000e+02 3.050000e+02 2.654988e+02 5.852051e-04 4.000000e-03
1.750000e+02 3.050000e+02 2.681164e+02 5.887946e-04 4.000000e-03
1.850000e+02 3.050000e+02 2.704766e+02 5.929786e-04 4.000000e-03
1.950000e+02 3.050000e+02 2.725528e+02 5.961642e-04 4.000000e-03
2.050000e+02 3.050000e+02 2.743215e+02 5.976224e-04 4.000000e-03
2.150000e+02 3.050000e+02 2.757621e+02 5.983041e-04 4.000000e-03
2.250000e+02 3.050000e+02 2.768578e+02 5.984595e-04 4.000000e-03
2.350000e+02 3.050000e+02 2.775956e+02 5.983125e-04 4.000000e-03
2.450000e+02 3.050000e+02 2.779668e+02 5.979012e-04 4.000000e-03
2.550000e+02 3.050000e+02 2.779668e+02 5.970956e-04 4.000000e-03
2.650000e+02 3.050000e+02 2.775956e+02 5.960208e-04 4.000000e-03
2.750000e+02 3.050000e+02 2.768578e+02 5.946522e-04 4.000000e-03
2.850000e+02 3.050000e+02 2.757621e+02 5.928166e-04 4.000000e-03
2.950000e+02 3.050000e+02 2.743215e+02 5.903114e-04 4.000000e-03
3.050000e+02 3.050000e+02 2.725528e+02 5.865195e-04 4.000000e-03
3.150000e+02 3.050000e+02 2.704766e+02 5.819050e-04 4.000000e-03
3.250000e+02 3.050000e+02 2.681164e+02 5.775720e-04 4.000000e-03
3.350000e+02 3.050000e+02 2.654988e+02 5.730517e-04 4.000000e-03
3.450000e+02 3.050000e+02 2.626524e+02 5.672268e-04 4.000000e-03
3.550000e+02 3.050000e+02 2.596076e+02 5.609182e-04 4.000000e-03
3.650000e+02 3.050000e+02 2.563960e+02 5.547160e-04 4.000000e-03
3.750000e+02 3.050000e+02 2.530500e+02 5.475155e-04 4.000000e-03
1.050000e+02 3.150000e+02 2.445174e+02 5.443478e-04 4.000000e-03
1.150000e+02 3.150000e+02 2.479454e+02 5.514585e-04 4.000000e-03
1.250000e+02 3.150000e+02 2.513052e+02 5.591445e-04 4.000000e-03
1.350000e+02 3.150000e+02 2.545655e+02 5.656576e-04 4.000000e-03
1.450000e+02 3.150000e+02 2.576947e+02 5.711102e-04 4.000000e-03
1.550000e+02 3.150000e+02 2.606614e+02 5.768536e-04 4.000000e-03
1.650000e+02 3.150000e+02 2.634348e+02 5.816035e-04 4.000000e-03
1.750000e+02 3.150000e+02 2.659854e+02 5.849171e-04 4.000000e-03
1.850000e+02 3.150000e+02 2.682850e+02 5.879063e-04 4.000000e-03
1.950000e+02 3.150000e+02 2.703080e+02 5.906821e-04 4.000000e-03
2.050000e+02 3.150000e+02 2.720313e+02 5.928071e-04 4.000000e-03
2.150000e+02 3.150000e+02 2.734350e+02 5.943868e-04 4.000000e-03
2.250000e+02 3.150000e+02 2.745027e+02 5.947804e-04 4.000000e-03
2.350000e+02 3.150000e+02 2.752216e+02 5.945003e-04 4.000000e-03
2.450000e+02 3.150000e+02 2.755831e+02 5.939834e-04 4.000000e-03
2.550000e+02 3.150000e+02 2.755831e+02 5.933191e-04 4.000000e-03
2.650000e+02 3.150000e+02 2.752216e+02 5.924121e-04 4.000000e-03
2.750000e+02 3.150000e+02 2.745027e+02 5.908904e-04 4.000000e-03
2.850000e+02 3.150000e+02 2.734350e+02 5.884179e-04 4.000000e-03
2.950000e+02 3.150000e+02 2.720313e+02 5.854754e-04 4.000000e-03
3.050000e+02 3.150000e+02 2.703080e+02 5.819845e-04 4.000000e-03
3.150000e+02 3.150000e+02 2.682850e+02 5.782192e-04 4.000000e-03
3.250000e+02 3.150000e+02 2.659854e+02 5.741842e-04 4.000000e-03
3.350000e+02 3.150000e+02 2.634348e+02 5.692098e-04 4.000000e-03
3.450000e+02 3.150000e+02 2.606614e+02 5.634212e-04 4.000000e-03
3.550000e+02 3.150000e+02 2.576947e+02 5.575839e-04 4.000000e-03
3.650000e+02 3.150000e+02 2.545655e+02 5.511218e-04 4.000000e-03
3.750000e+02 3.150000e+02 2.513052e+02 5.437527e-04 4.000000e-03
1.050000e+02 3.250000e+02 2.427601e+02 5.404234e-04 4.000000e-03
1.150000e+02 3.250000e+02 2.460869e+02 5.472008e-04 4.000000e-03
1.250000e+02 3.250000e+02 2.493475e+02 5.545913e-04 4.000000e-03
1.350000e+02 3.250000e+02 2.525115e+02 5.617819e-04 4.000000e-03
1.450000e+02 3.250000e+02 2.555483e+02 5.668409e-04 4.000000e-03
1.550000e+02 3.250000e+02 2.584275e+02 5.715018e-04 4.000000e-03
1.650000e+02 3.250000e+02 2.611190e+02 5.765357e-04 4.000000e-03
1.750000e+02 3.250000e+02 2.635942e+02 5.807427e-04 4.000000e-03
1.850000e+02 3.250000e+02 2.658260e+02 5.834651e-04 4.000000e-03
1.950000e+02 3.250000e+02 2.677892e+02 5.853640e-04 4.000000e-03
2.050000e+02 3.250000e+02 2.694617e+02 5.870820e-04 4.000000e-03
2.150000e+02 3.250000e+02 2.708240e+02 5.885576e-04 4.000000e-03
2.250000e+02 3.250000e+02 2.718600e+02 5.890350e-04 4.000000e-03
2.350000e+02 3.250000e+02 2.725577e+02 5.888417e-04 4.000000e-03
2.450000e+02 3.250000e+02 2.729086e+02 5.883872e-04 4.000000e-03
2.550000e+02 3.250000e+02 2.729086e+02 5.877853e-04 4.000000e-03
2.650000e+02 3.250000e+02 2.725577e+02 5.869246e-04 4.000000e-03
2.750000e+02 3.250000e+02 2.718600e+02 5.854428e-04 4.000000e-03
2.850000e+02 3.250000e+02 2.708240e+02 5.831063e-04 4.000000e-03
2.950000e+02 3.250000e+02 2.694617e+02 5.805170e-04 4.000000e-03
3.050000e+02 3.250000e+02 2.677892e+02 5.777134e-04 4.000000e-03
3.150000e+02 3.250000e+02 2.658260e+02 5.742408e-04 4.000000e-03
3.250000e+02 3.250000e+02 2.635942e+02 5.697216e-04 4.000000e-03
3.350000e+02 3.250000e+02 2.611190e+02 5.645182e-04 4.000000e-03
3.450000e+02 3.250000e+02 2.584275e+02 5.593523e-04 4.000000e-03
3.550000e+02 3.250000e+02 2.555483e+02 5.538546e-04 4.000000e-03
3.650000e+02 3.250000e+02 2.525115e+02 5.470277e-04 4.000000e-03
3.750000e+02 3.250000e+02 2.493475e+02 5.397383e-04 4.000000e-03
1.050000e+02 3.350000e+02 2.408312e+02 5.351498e-04 4.000000e-03
1.150000e+02 3.350000e+02 2.440469e+02 5.424034e-04 4.000000e-03
1.250000e+02 3.350000e+02 2.471986e+02 5.491569e-04 4.000000e-03
1.350000e+02 3.350000e+02 2.502569e+02 5.562992e-04 4.000000e-03
1.450000e+02 3.350000e+02 2.531923e+02 5.621954e-04 4.000000e-03
1.550000e+02 3.350000e+02 2.559753e+02 5.664371e-04 4.000000e-03
1.650000e+02 3.350000e+02 2.585770e+02 5.705807e-04 4.000000e-03
1.750000e+02 3.350000e+02 2.609695e+02 5.751452e-04 4.000000e-03
1.850000e+02 3.350000e+02 2.631267e+02 5.787395e-04 4.000000e-03
1.950000e+02 3.350000e+02 2.650244e+02 5.808021e-04 4.000000e-03
2.050000e+02 3.350000e+02 2.666411e+02 5.821589e-04 4.000000e-03
2.150000e+02 3.350000e+02 2.679578e+02 5.829232e-04 4.000000e-03
2.250000e+02 3.350000e+02 2.689593e+02 5.830906e-04 4.000000e-03
2.350000e+02 3.350000e+02 2.696337e+02 5.828630e-04 4.000000e-03
2.450000e+02 3.350000e+02 2.699729e+02 5.824423e-04 4.000000e-03
2.550000e+02 3.350000e+02 2.699729e+02 5.818963e-04 4.000000e-03
2.650000e+02 3.350000e+02 2.696337e+02 5.811537e-04 4.000000e-03
2.750000e+02 3.350000e+02 2.689593e+02 5.800181e-04 4.000000e-03
2.850000e+02 3.350000e+02 2.679578e+02 5.783115e-04 4.000000e-03
2.950000e+02 3.350000e+02 2.666411e+02 5.760619e-04 4.000000e-03
3.050000e+02 3.350000e+02 2.650244e+02 5.731967e-04 4.000000e-03
3.150000e+02 3.350000e+02 2.631267e+02 5.691950e-04 4.000000e-03
3.250000e+02 3.350000e+02 2.609695e+02 5.643932e-04 4.000000e-03
3.350000e+02 3.350000e+02 2.585770e+02 5.597115e-04 4.000000e-03
3.450000e+02 3.350000e+02 2.559753e+02 5.549189e-04 4.000000e-03
3.550000e+02 3.350000e+02 2.531923e+02 5.490422e-04 4.000000e-03
3.650000e+02 3.350000e+02 2.502569e+02 5.421354e-04 4.000000e-03
3.750000e+02 3.350000e+02 2.471986e+02 5.352196e-04 4.000000e-03
1.050000e+02 3.450000e+02 2.387514e+02 5.291055e-04 4.000000e-03
1.150000e+02 3.450000e+02 2.418474e+02 5.363853e-04 4.000000e-03
1.250000e+02 3.450000e+02 2.448817e+02 5.431556e-04 4.000000e-03
1.350000e+02 3.450000e+02 2.478261e+02 5.495701e-04 4.000000e-03
1.450000e+02 3.450000e+02 2.506521e+02 5.559364e-04 4.000000e-03
1.550000e+02 3.450000e+02 2.533315e+02 5.609401e-04 4.000000e-03
1.650000e+02 3.450000e+02 2.558362e+02 5.646659e-04 4.000000e-03
1.750000e+02 3.450000e+02 2.581396e+02 5.684769e-04 4.000000e-03
1.850000e+02 3.450000e+02 2.602165e+02 5.720889e-04 4.000000e-03
1.950000e+02 3.450000e+02 2.620435e+02 5.748991e-04 4.000000e-03
2.050000e+02 3.450000e+02 2.635999e+02 5.770222e-04 4.000000e-03
2.150000e+02 3.450000e+02 2.648676e+02 5.778473e-04 4.000000e-03
2.250000e+02 3.450000e+02 2.658318e+02 5.778804e-04 4.000000e-03
2.350000e+02 3.450000e+02 2.664810e+02 5.775891e-04 4.000000e-03
2.450000e+02 3.450000e+02 2.668076e+02 5.771647e-04 4.000000e-03
2.550000e+02 3.450000e+02 2.668076e+02 5.766570e-04 4.000000e-03
2.650000e+02 3.450000e+02 2.664810e+02 5.760080e-04 4.000000e-03
2.750000e+02 3.450000e+02 2.658318e+02 5.750274e-04 4.000000e-03
2.850000e+02 3.450000e+02 2.648676e+02 5.733201e-04 4.000000e-03
2.950000e+02 3.450000e+02 2.635999e+02 5.705590e-04 4.000000e-03
3.050000e+02 3.450000e+02 2.620435e+02 5.671843e-04 4.000000e-03
3.150000e+02 3.450000e+02 2.602165e+02 5.631316e-04 4.000000e-03
3.250000e+02 3.450000e+02 2.581396e+02 5.588117e-04 4.000000e-03
3.350000e+02 3.450000e+02 2.558362e+02 5.544935e-04 4.000000e-03
3.450000e+02 3.450000e+02 2.533315e+02 5.493405e-04 4.000000e-03
3.550000e+02 3.450000e+02 2.506521e+02 5.431094e-04 4.000000e-03
3.650000e+02 3.450000e+02 2.478261e+02 5.365783e-04 4.000000e-03
3.750000e+02 3.450000e+02 2.448817e+02 5.298225e-04 4.000000e-03
1.050000e+02 3.550000e+02 2.365428e+02 5.235083e-04 4.000000e-03
1.150000e+02 3.550000e+02 2.395115e+02 5.298688e-04 4.000000e-03
1.250000e+02 3.550000e+02 2.424211e+02 5.366340e-04 4.000000e-03
1.350000e+02 3.550000e+02 2.452446e+02 5.429506e-04 4.000000e-03
1.450000e+02 3.550000e+02 2.479545e+02 5.486717e-04 4.000000e-03
1.550000e+02 3.550000e+02 2.505237e+02 5.537523e-04 4.000000e-03
1.650000e+02 3.550000e+02 2.529255e+02 5.580486e-04 4.000000e-03
1.750000e+02 3.550000e+02 2.551343e+02 5.620033e-04 4.000000e-03
1.850000e+02 3.550000e+02 2.571258e+02 5.651344e-04 4.000000e-03
1.950000e+02 3.550000e+02 2.588778e+02 5.677294e-04 4.000000e-03
2.050000e+02 3.550000e+02 2.603702e+02 5.699353e-04 4.000000e-03
2.150000e+02 3.550000e+02 2.615858e+02 5.710071e-04 4.000000e-03
2.250000e+02 3.550000e+02 2.625104e+02 5.712583e-04 4.000000e-03
2.350000e+02 3.550000e+02 2.631330e+02 5.711114e-04 4.000000e-03
2.450000e+02 3.550000e+02 2.634461e+02 5.707644e-04 4.000000e-03
2.550000e+02 3.550000e+02 2.634461e+02 5.702735e-04 4.000000e-03
2.650000e+02 3.550000e+02 2.631330e+02 5.695823e-04 4.000000e-03
2.750000e+02 3.550000e+02 2.625104e+02 5.685039e-04 4.000000e-03
2.850000e+02 3.550000e+02 2.615858e+02 5.666697e-04 4.000000e-03
2.950000e+02 3.550000e+02 2.603702e+02 5.638591e-04 4.000000e-03
3.050000e+02 3.550000e+02 2.588778e+02 5.606492e-04 4.000000e-03
3.150000e+02 3.550000e+02 2.571258e+02 5.569690e-04 4.000000e-03
3.250000e+02 3.550000e+02 2.551343e+02 5.526733e-04 4.000000e-03
3.350000e+02 3.550000e+02 2.529255e+02 5.480608e-04 4.000000e-03
3.450000e+02 3.550000e+02 2.505237e+02 5.428038e-04 4.000000e-03
3.550000e+02 3.550000e+02 2.479545e+02 5.369193e-04 4.000000e-03
3.650000e+02 3.550000e+02 2.452446e+02 5.305918e-04 4.000000e-03
3.750000e+02 3.550000e+02 2.424211e+02 5.239117e-04 4.000000e-03
1.050000e+02 3.650000e+02 2.342278e+02 5.178718e-04 4.000000e-03
1.150000e+02 3.650000e+02 2.370632e+02 5.238494e-04 4.000000e-03
1.250000e+02 3.650000e+02 2.398421e+02 5.298155e-04 4.000000e-03
1.350000e+02 3.650000e+02 2.425387e+02 5.361508e-04 4.000000e-03
1.450000e+02 3.650000e+02 2.451270e+02 5.418885e-04 4.000000e-03
1.550000e+02 3.650000e+02 2.475808e+02 5.465390e-04 4.000000e-03
1.650000e+02 3.650000e+02 2.498748e+02 5.507862e-04 4.000000e-03
1.750000e+02 3.650000e+02 2.519843e+02 5.551481e-04 4.000000e-03
1.850000e+02 3.650000e+02 2.538864e+02 5.586489e-04 4.000000e-03
1.950000e+02 3.650000e+02 2.555596e+02 5.608019e-04 4.000000e-03
2.050000e+02 3.650000e+02 2.569850e+02 5.624156e-04 4.000000e-03
2.150000e+02 3.650000e+02 2.581461e+02 5.635730e-04 4.000000e-03
2.250000e+02 3.650000e+02 2.590291e+02 5.640903e-04 4.000000e-03
2.350000e+02 3.650000e+02 2.596237e+02 5.641136e-04 4.000000e-03
2.450000e+02 3.650000e+02 2.599228e+02 5.638430e-04 4.000000e-03
2.550000e+02 3.650000e+02 2.599228e+02 5.633496e-04 4.000000e-03
2.650000e+02 3.650000e+02 2.596237e+02 5.625662e-04 4.000000e-03
2.750000e+02 3.650000e+02 2.590291e+02 5.613112e-04 4.000000e-03
2.850000e+02 3.650000e+02 2.581461e+02 5.594587e-04 4.000000e-03
2.950000e+02 3.650000e+02 2.569850e+02 5.571517e-04 4.000000e-03
3.050000e+02 3.650000e+02 2.555596e+02 5.543435e-04 4.000000e-03
3.150000e+02 3.650000e+02 2.538864e+02 5.504999e-04 4.000000e-03
3.250000e+02 3.650000e+02 2.519843e+02 5.459411e-04 4.000000e-03
3.350000e+02 3.650000e+02 2.498748e+02 5.413349e-04 4.000000e-03
3.450000e+02 3.650000e+02 2.475808e+02 5.363574e-04 4.000000e-03
3.550000e+02 3.650000e+02 2.451270e+02 5.305824e-04 4.000000e-03
3.650000e+02 3.650000e+02 2.425387e+02 5.243051e-04 4.000000e-03
3.750000e+02 3.650000e+02 2.398421e+02 5.181239e-04 4.000000e-03
simpegPF/Dev/GRAV/SimPEG_inv_l0l2.sus
+44180
View File
File diff suppressed because it is too large Load Diff
simpegPF/Dev/GRAV/SimPEG_inv_l2l2.sus
+44180
View File
File diff suppressed because it is too large Load Diff
simpegPF/Dev/GRAV/gzfor3d.dat
+617
View File
@@ -0,0 +1,617 @@
616
105 155 240.8312 0.005799063 0.004
115 155 244.0469 0.007038728 0.004
125 155 247.1986 0.008397855 0.004
135 155 250.2569 0.009833759 0.004
145 155 253.1923 0.01128729 0.004
155 155 255.9753 0.01268985 0.004
165 155 258.577 0.01397494 0.004
175 155 260.9695 0.0150914 0.004
185 155 263.1267 0.01601382 0.004
195 155 265.0244 0.01674675 0.004
205 155 266.6411 0.01732118 0.004
215 155 267.9578 0.01778553 0.004
225 155 268.9593 0.01819414 0.004
235 155 269.6337 0.01859648 0.004
245 155 269.9729 0.01902865 0.004
255 155 269.9729 0.01950741 0.004
265 155 269.6337 0.02002644 0.004
275 155 268.9593 0.02055442 0.004
285 155 267.9578 0.02103531 0.004
295 155 266.6411 0.02139193 0.004
305 155 265.0244 0.02153416 0.004
315 155 263.1267 0.02137253 0.004
325 155 260.9695 0.02083605 0.004
335 155 258.577 0.01989047 0.004
345 155 255.9753 0.01855089 0.004
355 155 253.1923 0.0168833 0.004
365 155 250.2569 0.01499306 0.004
375 155 247.1986 0.01300364 0.004
105 165 242.7601 0.006972249 0.004
115 165 246.0869 0.008529513 0.004
125 165 249.3475 0.0102487 0.004
135 165 252.5115 0.0120693 0.004
145 165 255.5483 0.01390455 0.004
155 165 258.4275 0.01565232 0.004
165 165 261.119 0.0172141 0.004
175 165 263.5942 0.01851703 0.004
185 165 265.826 0.01953068 0.004
195 165 267.7892 0.02027232 0.004
205 165 269.4617 0.02079926 0.004
215 165 270.824 0.02119273 0.004
225 165 271.86 0.0215395 0.004
235 165 272.5577 0.02191636 0.004
245 165 272.9086 0.02237931 0.004
255 165 272.9086 0.02295668 0.004
265 165 272.5577 0.02364491 0.004
275 165 271.86 0.02440566 0.004
285 165 270.824 0.02516473 0.004
295 165 269.4617 0.02581447 0.004
305 165 267.7892 0.02622257 0.004
315 165 265.826 0.02624999 0.004
325 165 263.5942 0.0257775 0.004
335 165 261.119 0.02473599 0.004
345 165 258.4275 0.02313032 0.004
355 165 255.5483 0.02104608 0.004
365 165 252.5115 0.01863351 0.004
375 165 249.3475 0.01607355 0.004
105 175 244.5174 0.008325776 0.004
115 175 247.9454 0.01027856 0.004
125 175 251.3052 0.01245309 0.004
135 175 254.5655 0.01476452 0.004
145 175 257.6947 0.01708602 0.004
155 175 260.6614 0.01926537 0.004
165 175 263.4348 0.02115616 0.004
175 175 265.9854 0.02265431 0.004
185 175 268.285 0.02372502 0.004
195 175 270.308 0.02440855 0.004
205 175 272.0313 0.02480421 0.004
215 175 273.435 0.025042 0.004
225 175 274.5027 0.02525406 0.004
235 175 275.2216 0.02555334 0.004
245 175 275.5831 0.02602111 0.004
255 175 275.5831 0.02670055 0.004
265 175 275.2216 0.02759283 0.004
275 175 274.5027 0.02865335 0.004
285 175 273.435 0.02978808 0.004
295 175 272.0313 0.03085284 0.004
305 175 270.308 0.03166075 0.004
315 175 268.285 0.03200364 0.004
325 175 265.9854 0.03168994 0.004
335 175 263.4348 0.03059301 0.004
345 175 260.6614 0.02869458 0.004
355 175 257.6947 0.02610265 0.004
365 175 254.5655 0.0230307 0.004
375 175 251.3052 0.01974359 0.004
105 185 246.0835 0.009857299 0.004
115 185 249.6018 0.01229506 0.004
125 185 253.05 0.01503834 0.004
135 185 256.396 0.0179703 0.004
145 185 259.6076 0.02090728 0.004
155 185 262.6524 0.02362343 0.004
165 185 265.4988 0.02590153 0.004
175 185 268.1164 0.02759335 0.004
185 185 270.4766 0.02866242 0.004
195 185 272.5528 0.02918845 0.004
205 185 274.3215 0.02933548 0.004
215 185 275.7621 0.02930345 0.004
225 185 276.8578 0.02928489 0.004
235 185 277.5956 0.02943757 0.004
245 185 277.9668 0.02987237 0.004
255 185 277.9668 0.03064962 0.004
265 185 277.5956 0.03177723 0.004
275 185 276.8578 0.0332062 0.004
285 185 275.7621 0.03482317 0.004
295 185 274.3215 0.03644361 0.004
305 185 272.5528 0.0378143 0.004
315 185 270.4766 0.03863591 0.004
325 185 268.1164 0.03861361 0.004
335 185 265.4988 0.03753102 0.004
345 185 262.6524 0.0353252 0.004
355 185 259.6076 0.03212676 0.004
365 185 256.396 0.02823597 0.004
375 185 253.05 0.02403795 0.004
105 195 247.4409 0.01154312 0.004
115 195 251.0374 0.01456078 0.004
125 195 254.5621 0.01799904 0.004
135 195 257.9825 0.02170106 0.004
145 195 261.2655 0.0254053 0.004
155 195 264.3779 0.02878064 0.004
165 195 267.2876 0.03150758 0.004
175 195 269.9634 0.03337698 0.004
185 195 272.376 0.03435672 0.004
195 195 274.4984 0.03459087 0.004
205 195 276.3063 0.03433925 0.004
215 195 277.779 0.03389742 0.004
225 195 278.899 0.03353386 0.004
235 195 279.6533 0.03345782 0.004
245 195 280.0326 0.0338114 0.004
255 195 280.0326 0.03467241 0.004
265 195 279.6533 0.03605681 0.004
275 195 278.899 0.03791414 0.004
285 195 277.779 0.04011483 0.004
295 195 276.3063 0.04243402 0.004
305 195 274.4984 0.04454388 0.004
315 195 272.376 0.04603266 0.004
325 195 269.9634 0.04646737 0.004
335 195 267.2876 0.04550055 0.004
345 195 264.3779 0.04299266 0.004
355 195 261.2655 0.03909171 0.004
365 195 257.9825 0.03421217 0.004
375 195 254.5621 0.02890575 0.004
105 205 248.5739 0.01332811 0.004
115 205 252.2356 0.01701312 0.004
125 205 255.8243 0.02127089 0.004
135 205 259.3068 0.02589826 0.004
145 205 262.6493 0.03053213 0.004
155 205 265.8182 0.03469705 0.004
165 205 268.7807 0.03793137 0.004
175 205 271.505 0.03994489 0.004
185 205 273.9614 0.04072062 0.004
195 205 276.1223 0.040501 0.004
205 205 277.9631 0.03967994 0.004
215 205 279.4624 0.03867638 0.004
225 205 280.6028 0.03784785 0.004
235 205 281.3707 0.03745748 0.004
245 205 281.757 0.03767609 0.004
255 205 281.757 0.03859633 0.004
265 205 281.3707 0.0402421 0.004
275 205 280.6028 0.04256465 0.004
285 205 279.4624 0.04542375 0.004
295 205 277.9631 0.04855875 0.004
305 205 276.1223 0.05156521 0.004
315 205 273.9614 0.0539036 0.004
325 205 271.505 0.05497047 0.004
335 205 268.7807 0.05424284 0.004
345 205 265.8182 0.05146415 0.004
355 205 262.6493 0.04678684 0.004
365 205 259.3068 0.04076644 0.004
375 205 255.8243 0.03417279 0.004
105 215 249.4692 0.01511709 0.004
115 215 253.1825 0.01952782 0.004
125 215 256.8218 0.0247006 0.004
135 215 260.3533 0.03038382 0.004
145 215 263.7429 0.03609127 0.004
155 215 266.9565 0.04116193 0.004
165 215 269.9606 0.04494795 0.004
175 215 272.7234 0.04705669 0.004
185 215 275.2143 0.04750182 0.004
195 215 277.4056 0.04666402 0.004
205 215 279.2724 0.0451112 0.004
215 215 280.7928 0.0434099 0.004
225 215 281.9493 0.04201375 0.004
235 215 282.728 0.04123608 0.004
245 215 283.1197 0.04126942 0.004
255 215 283.1197 0.04221554 0.004
265 215 282.728 0.04410412 0.004
275 215 281.9493 0.04689043 0.004
285 215 280.7928 0.05042982 0.004
295 215 279.2724 0.05443395 0.004
305 215 277.4056 0.05842725 0.004
315 215 275.2143 0.06173844 0.004
325 215 272.7234 0.06357424 0.004
335 215 269.9606 0.06320119 0.004
345 215 266.9565 0.06020695 0.004
355 215 263.7429 0.05472921 0.004
365 215 260.3533 0.04748209 0.004
375 215 256.8218 0.03949725 0.004
105 225 250.1165 0.01677301 0.004
115 225 253.8671 0.02190986 0.004
125 225 257.5429 0.02802359 0.004
135 225 261.1099 0.03481869 0.004
145 225 264.5335 0.0416745 0.004
155 225 267.7793 0.04771334 0.004
165 225 270.8137 0.05206371 0.004
175 225 273.6041 0.05421411 0.004
185 225 276.1201 0.05422485 0.004
195 225 278.3334 0.0526493 0.004
205 225 280.2189 0.05025931 0.004
215 225 281.7546 0.04778167 0.004
225 225 282.9227 0.04576373 0.004
235 225 283.7092 0.04455936 0.004
245 225 284.1048 0.04437259 0.004
255 225 284.1048 0.0453067 0.004
265 225 283.7092 0.04739323 0.004
275 225 282.9227 0.05059164 0.004
285 225 281.7546 0.05475742 0.004
295 225 280.2189 0.05958364 0.004
305 225 278.3334 0.06453513 0.004
315 225 276.1201 0.0688198 0.004
325 225 273.6041 0.07145822 0.004
335 225 270.8137 0.07149845 0.004
345 225 267.7793 0.06835354 0.004
355 225 264.5335 0.06212873 0.004
365 225 261.1099 0.05369584 0.004
375 225 257.5429 0.04436213 0.004
105 235 250.508 0.01812822 0.004
115 235 254.2811 0.02390501 0.004
125 235 257.979 0.0308713 0.004
135 235 261.5674 0.0386997 0.004
145 235 265.0117 0.04664332 0.004
155 235 268.277 0.05360471 0.004
165 235 271.3296 0.05848048 0.004
175 235 274.1368 0.06063609 0.004
185 235 276.6679 0.060182 0.004
195 235 278.8945 0.05785612 0.004
205 235 280.7913 0.05463869 0.004
215 235 282.3363 0.05140971 0.004
225 235 283.5114 0.04879638 0.004
235 235 284.3026 0.04717992 0.004
245 235 284.7007 0.04676582 0.004
255 235 284.7007 0.04765184 0.004
265 235 284.3026 0.04986681 0.004
275 235 283.5114 0.05337203 0.004
285 235 282.3363 0.05802429 0.004
295 235 280.7913 0.06350478 0.004
305 235 278.8945 0.06923437 0.004
315 235 276.6679 0.07432524 0.004
325 235 274.1368 0.07764489 0.004
335 235 271.3296 0.07805449 0.004
345 235 268.277 0.07481179 0.004
355 235 265.0117 0.06798573 0.004
365 235 261.5674 0.05858048 0.004
375 235 257.9791 0.04814234 0.004
105 245 250.639 0.01901306 0.004
115 245 254.4197 0.02524244 0.004
125 245 258.125 0.03282988 0.004
135 245 261.7206 0.04143298 0.004
145 245 265.1717 0.05021275 0.004
155 245 268.4435 0.0578949 0.004
165 245 271.5022 0.06318747 0.004
175 245 274.3151 0.06534731 0.004
185 245 276.8512 0.06451759 0.004
195 245 279.0823 0.06159046 0.004
205 245 280.9829 0.0577174 0.004
215 245 282.531 0.05389807 0.004
225 245 283.7084 0.0508164 0.004
235 245 284.5013 0.04886841 0.004
245 245 284.9001 0.04825519 0.004
255 245 284.9001 0.04906546 0.004
265 245 284.5013 0.05132165 0.004
275 245 283.7084 0.05498284 0.004
285 245 282.531 0.05990451 0.004
295 245 280.9829 0.0657597 0.004
305 245 279.0823 0.07194213 0.004
315 245 276.8512 0.07750589 0.004
325 245 274.3151 0.08122502 0.004
335 245 271.5022 0.08184729 0.004
345 245 268.4435 0.07853589 0.004
355 245 265.1717 0.07133991 0.004
365 245 261.7206 0.06134745 0.004
375 245 258.125 0.05025301 0.004
105 255 250.508 0.0192972 0.004
115 255 254.2811 0.02570488 0.004
125 255 257.979 0.03355205 0.004
135 255 261.5674 0.04249819 0.004
145 255 265.0117 0.05166975 0.004
155 255 268.277 0.05971105 0.004
165 255 271.3296 0.06523315 0.004
175 255 274.1368 0.06742705 0.004
185 255 276.6679 0.06643808 0.004
195 255 278.8945 0.06322887 0.004
205 255 280.7913 0.05903525 0.004
215 255 282.3363 0.0549177 0.004
225 255 283.5114 0.05158873 0.004
235 255 284.3026 0.0494511 0.004
245 255 284.7007 0.04870165 0.004
255 255 284.7007 0.04942128 0.004
265 255 284.3026 0.05162471 0.004
275 255 283.5114 0.05526448 0.004
285 255 282.3363 0.06019031 0.004
295 255 280.7913 0.06606859 0.004
305 255 278.8945 0.07228396 0.004
315 255 276.6679 0.07787768 0.004
325 255 274.1368 0.08160897 0.004
335 255 271.3296 0.08221248 0.004
345 255 268.277 0.07884723 0.004
355 255 265.0117 0.07157118 0.004
365 255 261.5674 0.06149274 0.004
375 255 257.9791 0.05032607 0.004
105 265 250.1165 0.01892904 0.004
115 265 253.8671 0.02520065 0.004
125 265 257.5429 0.03288214 0.004
135 265 261.1099 0.04164867 0.004
145 265 264.5335 0.0506596 0.004
155 265 267.7793 0.05860138 0.004
165 265 270.8137 0.06411037 0.004
175 265 273.6041 0.06637393 0.004
185 265 276.1201 0.06550705 0.004
195 265 278.3334 0.0624285 0.004
205 265 280.2189 0.05834198 0.004
215 265 281.7546 0.05429351 0.004
225 265 282.9227 0.05099324 0.004
235 265 283.7092 0.04884508 0.004
245 265 284.1048 0.04804586 0.004
255 265 284.1048 0.04867341 0.004
265 265 283.7092 0.0507358 0.004
275 265 282.9227 0.05417622 0.004
285 265 281.7546 0.05883484 0.004
295 265 280.2189 0.06437394 0.004
305 265 278.3334 0.07018831 0.004
315 265 276.1201 0.07535485 0.004
325 265 273.6041 0.07870006 0.004
335 265 270.8137 0.07905139 0.004
345 265 267.7793 0.07565596 0.004
355 265 264.5335 0.06860885 0.004
365 265 261.1099 0.0589686 0.004
375 265 257.5429 0.04833401 0.004
105 275 249.4692 0.01795401 0.004
115 275 253.1825 0.02379951 0.004
125 275 256.8218 0.03092258 0.004
135 275 260.3533 0.03902478 0.004
145 275 263.7429 0.04735598 0.004
155 275 266.9565 0.05475529 0.004
165 275 269.9606 0.06000124 0.004
175 275 272.7234 0.06234236 0.004
185 275 275.2143 0.06183979 0.004
195 275 277.4056 0.05926803 0.004
205 275 279.2724 0.05568955 0.004
215 275 280.7928 0.05206215 0.004
225 275 281.9493 0.04905983 0.004
235 275 282.728 0.04707991 0.004
245 275 283.1197 0.04632227 0.004
255 275 283.1197 0.04686655 0.004
265 275 282.728 0.0487167 0.004
275 275 281.9493 0.05180627 0.004
285 275 280.7928 0.0559657 0.004
295 275 279.2724 0.06085887 0.004
305 275 277.4056 0.06591204 0.004
315 275 275.2143 0.07028226 0.004
325 275 272.7234 0.07293456 0.004
335 275 269.9606 0.07287328 0.004
345 275 266.9565 0.06950203 0.004
355 275 263.7429 0.06296467 0.004
365 275 260.3533 0.05420648 0.004
375 275 256.8218 0.04460298 0.004
105 285 248.5739 0.01650131 0.004
115 285 252.2356 0.02170977 0.004
125 285 255.8243 0.0279969 0.004
135 285 259.3068 0.03509919 0.004
145 285 262.6493 0.04239197 0.004
155 285 265.8182 0.04892891 0.004
165 285 268.7807 0.05370994 0.004
175 285 271.505 0.05609174 0.004
185 285 273.9614 0.05607118 0.004
195 285 276.1223 0.05423043 0.004
205 285 277.9631 0.05142543 0.004
215 285 279.4624 0.04846976 0.004
225 285 280.6028 0.045968 0.004
235 285 281.3707 0.04429695 0.004
245 285 281.757 0.04365721 0.004
255 285 281.757 0.04413205 0.004
265 285 281.3707 0.04572378 0.004
275 285 280.6028 0.04835886 0.004
285 285 279.4624 0.05186252 0.004
295 285 277.9631 0.05591103 0.004
305 285 276.1223 0.05998318 0.004
315 285 273.9614 0.06335168 0.004
325 285 271.505 0.06516555 0.004
335 285 268.7807 0.06465053 0.004
345 285 265.8182 0.06139148 0.004
355 285 262.6493 0.05557041 0.004
365 285 259.3068 0.04797648 0.004
375 285 255.8243 0.03970748 0.004
105 295 247.4409 0.01474585 0.004
115 295 251.0374 0.01920988 0.004
125 295 254.5621 0.02453186 0.004
135 295 257.9825 0.03048903 0.004
145 295 261.2655 0.03659167 0.004
155 295 264.3779 0.04211888 0.004
165 295 267.2876 0.04631082 0.004
175 295 269.9634 0.04865492 0.004
185 295 272.376 0.0490971 0.004
195 295 274.4984 0.04802939 0.004
205 295 276.3063 0.04608562 0.004
215 295 277.779 0.04390886 0.004
225 295 278.899 0.04200883 0.004
235 295 279.6533 0.04072504 0.004
245 295 280.0326 0.04025046 0.004
255 295 280.0326 0.04066963 0.004
265 295 279.6533 0.04198469 0.004
275 295 278.899 0.04411978 0.004
285 295 277.779 0.04690359 0.004
295 295 276.3063 0.05003831 0.004
305 295 274.4984 0.05307376 0.004
315 295 272.376 0.05541801 0.004
325 295 269.9634 0.05641811 0.004
335 295 267.2876 0.05552086 0.004
345 295 264.3779 0.05247445 0.004
355 295 261.2655 0.04747431 0.004
365 295 257.9825 0.04113925 0.004
375 295 254.5621 0.03429216 0.004
105 305 246.0835 0.01286437 0.004
115 305 249.6018 0.01657195 0.004
125 305 253.05 0.02093087 0.004
135 305 256.396 0.025761 0.004
145 305 259.6076 0.03069708 0.004
155 305 262.6524 0.03521911 0.004
165 305 265.4988 0.03878074 0.004
175 305 268.1164 0.04099702 0.004
185 305 270.4766 0.0417892 0.004
195 305 272.5528 0.04139455 0.004
205 305 274.3215 0.04024738 0.004
215 305 275.7621 0.03882269 0.004
225 305 276.8578 0.03752418 0.004
235 305 277.5956 0.03663987 0.004
245 305 277.9668 0.0363444 0.004
255 305 277.9668 0.03671838 0.004
265 305 277.5956 0.03776415 0.004
275 305 276.8578 0.03940897 0.004
285 305 275.7621 0.04149503 0.004
295 305 274.3215 0.04376318 0.004
305 305 272.5528 0.04584523 0.004
315 305 270.4766 0.04728653 0.004
325 305 268.1164 0.04761747 0.004
335 305 265.4988 0.04647246 0.004
345 305 262.6524 0.04372144 0.004
355 305 259.6076 0.03954818 0.004
365 305 256.396 0.03441396 0.004
375 305 253.05 0.02890733 0.004
105 315 244.5174 0.01100289 0.004
115 315 247.9454 0.01400962 0.004
125 315 251.3052 0.01749486 0.004
135 315 254.5655 0.02131889 0.004
145 315 257.6947 0.02521966 0.004
155 315 260.6614 0.0288366 0.004
165 315 263.4348 0.03179245 0.004
175 315 265.9854 0.03381035 0.004
185 315 268.285 0.03480784 0.004
195 315 270.308 0.03491417 0.004
205 315 272.0313 0.0344074 0.004
215 315 273.435 0.03361741 0.004
225 315 274.5027 0.03284509 0.004
235 315 275.2216 0.03231982 0.004
245 315 275.5831 0.03218912 0.004
255 315 275.5831 0.03252415 0.004
265 315 275.2216 0.03332707 0.004
275 315 274.5027 0.034533 0.004
285 315 273.435 0.03600577 0.004
295 315 272.0313 0.03753266 0.004
305 315 270.308 0.03882863 0.004
315 315 268.285 0.03956303 0.004
325 315 265.9854 0.03941708 0.004
335 315 263.4348 0.03816563 0.004
345 315 260.6614 0.03575625 0.004
355 315 257.6947 0.03234595 0.004
365 315 254.5655 0.02826607 0.004
375 315 251.3052 0.0239246 0.004
105 325 242.7601 0.009262438 0.004
115 325 246.0869 0.01165997 0.004
125 325 249.3475 0.01440206 0.004
135 325 252.5115 0.01738414 0.004
145 325 255.5483 0.0204234 0.004
155 325 258.4275 0.02327639 0.004
165 325 261.119 0.02568946 0.004
175 325 263.5942 0.02746935 0.004
185 325 265.826 0.0285419 0.004
195 325 267.7892 0.02896909 0.004
205 325 269.4617 0.0289186 0.004
215 325 270.824 0.02860734 0.004
225 325 271.86 0.02824721 0.004
235 325 272.5577 0.02800991 0.004
245 325 272.9086 0.02801165 0.004
255 325 272.9086 0.02831048 0.004
265 325 272.5577 0.02890786 0.004
275 325 271.86 0.02974958 0.004
285 325 270.824 0.03072531 0.004
295 325 269.4617 0.03167012 0.004
305 325 267.7892 0.03237453 0.004
315 325 265.826 0.03260964 0.004
325 325 263.5942 0.0321698 0.004
335 325 261.119 0.03092527 0.004
345 325 258.4275 0.02886682 0.004
355 325 255.5483 0.02612009 0.004
365 325 252.5115 0.02291743 0.004
375 325 249.3475 0.01953647 0.004
105 335 240.8312 0.007699974 0.004
115 335 244.0469 0.009590837 0.004
125 335 247.1986 0.01172752 0.004
135 335 250.2569 0.0140339 0.004
145 335 253.1923 0.01638488 0.004
155 335 255.9753 0.01861865 0.004
165 335 258.577 0.02056722 0.004
175 335 260.9695 0.02209812 0.004
185 335 263.1267 0.02315018 0.004
195 335 265.0244 0.02374723 0.004
205 335 266.6411 0.02398471 0.004
215 335 267.9578 0.02399836 0.004
225 335 268.9593 0.02392981 0.004
235 335 269.6337 0.02390022 0.004
245 335 269.9729 0.02399526 0.004
255 335 269.9729 0.02425872 0.004
265 335 269.6337 0.02469105 0.004
275 335 268.9593 0.02524948 0.004
285 335 267.9578 0.02584962 0.004
295 335 266.6411 0.02637017 0.004
305 335 265.0244 0.02666432 0.004
315 335 263.1267 0.0265805 0.004
325 335 260.9695 0.02599203 0.004
335 335 258.577 0.02482945 0.004
345 335 255.9753 0.02310427 0.004
355 335 253.1923 0.02091284 0.004
365 335 250.2569 0.0184159 0.004
375 335 247.1986 0.01580061 0.004
105 345 238.7514 0.00633773 0.004
115 345 241.8474 0.007819503 0.004
125 345 244.8817 0.009476531 0.004
135 345 247.8261 0.01125448 0.004
145 345 250.6521 0.01306869 0.004
155 345 253.3315 0.01481238 0.004
165 345 255.8362 0.01637497 0.004
175 345 258.1396 0.01766642 0.004
185 345 260.2165 0.01863876 0.004
195 345 262.0435 0.01929603 0.004
205 345 263.5999 0.01968909 0.004
215 345 264.8676 0.01989884 0.004
225 345 265.8318 0.02001533 0.004
235 345 266.481 0.02011947 0.004
245 345 266.8076 0.02027046 0.004
255 345 266.8076 0.02049899 0.004
265 345 266.481 0.02080442 0.004
275 345 265.8318 0.02115462 0.004
285 345 264.8676 0.02148826 0.004
295 345 263.5999 0.0217203 0.004
305 345 262.0435 0.02175223 0.004
315 345 260.2165 0.02148775 0.004
325 345 258.1396 0.02085237 0.004
335 345 255.8362 0.01981252 0.004
345 345 253.3315 0.0183876 0.004
355 345 250.6521 0.01664964 0.004
365 345 247.8261 0.01470919 0.004
375 345 244.8817 0.01269206 0.004
105 355 236.5428 0.005174561 0.004
115 355 239.5115 0.006332273 0.004
125 355 242.4211 0.007615574 0.004
135 355 245.2446 0.008986187 0.004
145 355 247.9545 0.01038705 0.004
155 355 250.5237 0.01174775 0.004
165 355 252.9255 0.01299539 0.004
175 355 255.1343 0.01406885 0.004
185 355 257.1258 0.01493152 0.004
195 355 258.8778 0.01557824 0.004
205 355 260.3702 0.01603358 0.004
215 355 261.5858 0.01634322 0.004
225 355 262.5104 0.01656157 0.004
235 355 263.133 0.01673973 0.004
245 355 263.4461 0.01691616 0.004
255 355 263.4461 0.01711076 0.004
265 355 263.133 0.01732203 0.004
275 355 262.5104 0.01752679 0.004
285 355 261.5858 0.01768243 0.004
295 355 260.3702 0.01773186 0.004
305 355 258.8778 0.01761165 0.004
315 355 257.1258 0.01726309 0.004
325 355 255.1343 0.01664471 0.004
335 355 252.9255 0.01574318 0.004
345 355 250.5237 0.01457928 0.004
355 355 247.9545 0.01320607 0.004
365 355 245.2446 0.01169955 0.004
375 355 242.4211 0.01014454 0.004
105 365 234.2278 0.00419585 0.004
115 365 237.0632 0.005099725 0.004
125 365 239.8421 0.006094392 0.004
135 365 242.5387 0.007153049 0.004
145 365 245.127 0.008237215 0.004
155 365 247.5808 0.00930021 0.004
165 365 249.8748 0.01029367 0.004
175 365 251.9843 0.01117584 0.004
185 365 253.8864 0.0119192 0.004
195 365 255.5596 0.01251499 0.004
205 365 256.985 0.012973 0.004
215 365 258.1461 0.01331717 0.004
225 365 259.0291 0.0135785 0.004
235 365 259.6237 0.0137874 0.004
245 365 259.9228 0.01396726 0.004
255 365 259.9228 0.01412991 0.004
265 365 259.6237 0.01427331 0.004
275 365 259.0291 0.01438129 0.004
285 365 258.1461 0.01442539 0.004
295 365 256.985 0.01436882 0.004
305 365 255.5596 0.01417239 0.004
315 365 253.8864 0.01380193 0.004
325 365 251.9843 0.01323605 0.004
335 365 249.8748 0.01247224 0.004
345 365 247.5808 0.01152964 0.004
355 365 245.127 0.01044711 0.004
365 365 242.5387 0.009277147 0.004
375 365 239.8421 0.008077218 0.004
simpegPF/Dev/MAG/Block_sphere_model.sus
+44180
View File
File diff suppressed because it is too large Load Diff
simpegPF/Dev/MAG/Gaussian.topo
+21610
View File
File diff suppressed because it is too large Load Diff
simpegPF/Dev/Intgrl_MAG_Inv_Driver.py → simpegPF/Dev/MAG/Intgrl_MAG_Inv_Driver.py
+1 -17
View File
@@ -84,26 +84,10 @@ else:
midx = int(mesh.nCx/2)
midy = int(mesh.nCy/2)
# Create forward operator
#F = PF.Magnetics.Intrgl_Fwr_Op(mesh,B,M_xyz,rxLoc,actv,'tmi')
# Get distance weighting function
#==============================================================================
# wr = PF.Magnetics.get_dist_wgt(mesh,rxLoc,actv,3.,np.min(mesh.hx)/4)
# #wrMap = PF.BaseMag.WeightMap(nC, wr)
#
#
# wr_out = actvMap * wr
# Mesh.TensorMesh.writeModelUBC(mesh,home_dir+dsep+'wr.dat',wr_out)
# #Utils.meshutils.writeUBCTensorModel(home_dir+dsep+'wr.dat',mesh,wr_out)
#
# #%%
# plt.figure()
# ax = plt.subplot()
# mesh.plotSlice(wr_out, ax = ax, normal = 'Y', ind=midx ,clim = (-1e-3, wr.max()))
# plt.title('Distance weighting')
# plt.xlabel('x');plt.ylabel('z')
# plt.gca().set_aspect('equal', adjustable='box')
#==============================================================================
#%% Plot obs data
@@ -144,7 +128,7 @@ PC = Utils.sdiag(diagA**-1.)
dmis = DataMisfit.l2_DataMisfit(survey)
dmis.Wd = wd
dmis.Wd = 1/wd
opt = Optimization.ProjectedGNCG(maxIter=10,lower=0.,upper=1., maxIterCG= 20, tolCG = 1e-3)
opt.approxHinv = PC
simpegPF/Dev/MAG/Mesh_10m.msh
+5
View File
@@ -0,0 +1,5 @@
47 47 20
-60.00 -60.00 280.00
50.00 30.00 20.00 41*10.00 20.00 30.00 50.00
50.00 30.00 20.00 41*10.00 20.00 30.00 50.00
17*10.00 20.00 30.00 50
simpegPF/Dev/Obs_loc_TMI.obs → simpegPF/Dev/MAG/Obs_loc_TMI.obs
View File
simpegPF/Dev/PYMAG3D_inv.inp → simpegPF/Dev/MAG/PYMAG3D_inv.inp
View File
simpegPF/Dev/Pred.dat → simpegPF/Dev/MAG/Pred.dat
View File
simpegPF/Dev/SimPEG_inv_l0l2.sus → simpegPF/Dev/MAG/SimPEG_inv_l0l2.sus
View File
simpegPF/Dev/SimPEG_inv_l2l2.sus → simpegPF/Dev/MAG/SimPEG_inv_l2l2.sus
View File
simpegPF/Gravity.py
+501
View File
@@ -0,0 +1,501 @@
from SimPEG import *
import BaseGrav as GRAV
import re
class GravityIntegral(Problem.BaseProblem):
#surveyPair = Survey.LinearSurvey
storeG = True #: Store the forward matrix by default, otherwise just compute d
actInd = None #: Active cell indices provided
def __init__(self, mesh, mapping=None, **kwargs):
Problem.BaseProblem.__init__(self, mesh, mapping=mapping, **kwargs)
def fwr_op(self):
# Add forward function
# kappa = self.curModel.kappa TODO
sus = self.mapping*self.curModel
return self.G.dot(sus)
def fields(self, m):
self.curModel = m
total = np.zeros(self.survey.nRx)
induced = self.fwr_op()
# rem = self.rem
if induced is not None:
total += induced
return total
# return self.G.dot(self.mapping*(m))
def Jvec(self, m, v, u=None):
dmudm = self.mapping.deriv(m)
return self.G.dot(dmudm*v)
def Jtvec(self, m, v, u=None):
dmudm = self.mapping.deriv(m)
return dmudm.T * (self.G.T.dot(v))
@property
def G(self):
if not self.ispaired:
raise Exception('Need to pair!')
if getattr(self,'_G', None) is None:
self._G = self.Intrgl_Fwr_Op( 'z' )
return self._G
def Intrgl_Fwr_Op(self, flag):
"""
Gravity forward operator in integral form
flag = 'z' | 'xyz'
Return
_G = Linear forward modeling operation
Created on March, 15th 2016
@author: dominiquef
"""
# Find non-zero cells
#inds = np.nonzero(actv)[0]
if getattr(self, 'actInd', None) is not None:
if self.actInd.dtype=='bool':
inds = np.asarray([inds for inds, elem in enumerate(self.actInd, 1) if elem], dtype = int) - 1
else:
inds = self.actInd
else:
inds = np.asarray(range(self.mesh.nC))
nC = len(inds)
# Create active cell projector
P = sp.csr_matrix((np.ones(nC),(inds, range(nC))),
shape=(self.mesh.nC, nC))
# Create vectors of nodal location (lower and upper coners for each cell)
xn = self.mesh.vectorNx;
yn = self.mesh.vectorNy;
zn = self.mesh.vectorNz;
yn2,xn2,zn2 = np.meshgrid(yn[1:], xn[1:], zn[1:])
yn1,xn1,zn1 = np.meshgrid(yn[0:-1], xn[0:-1], zn[0:-1])
Yn = P.T*np.c_[mkvc(yn1), mkvc(yn2)]
Xn = P.T*np.c_[mkvc(xn1), mkvc(xn2)]
Zn = P.T*np.c_[mkvc(zn1), mkvc(zn2)]
rxLoc = self.survey.srcField.rxList[0].locs
ndata = rxLoc.shape[0]
# Pre-allocate space and create magnetization matrix if required
# Pre-allocate space
if flag == 'z':
G = np.zeros((ndata, nC))
elif flag == 'xyz':
G = np.zeros((int(3*ndata), nC))
else:
print """Flag must be either 'z' | 'xyz', please revised"""
return
# Loop through all observations and create forward operator (ndata-by-nC)
print "Begin calculation of forward operator: " + flag
# Add counter to dsiplay progress. Good for large problems
count = -1;
for ii in range(ndata):
if flag=='z':
tt = get_T_mat(Xn,Yn,Zn,rxLoc[ii,:])
G[ii,:] = tt
elif flag == 'xyz':
print "Sorry 3-component not implemented yet"
# Display progress
count = progress(ii,count,ndata)
print "Done 100% ...forward operator completed!!\n"
return G
def get_T_mat(Xn,Yn,Zn,rxLoc):
"""
Load in the active nodes of a tensor mesh and computes the gravity tensor
for a given observation location rxLoc[obsx, obsy, obsz]
INPUT:
Xn, Yn, Zn: Node location matrix for the lower and upper most corners of
all cells in the mesh shape[nC,2]
M
OUTPUT:
Tx = [Txx Txy Txz]
Ty = [Tyx Tyy Tyz]
Tz = [Tzx Tzy Tzz]
where each elements have dimension 1-by-nC.
Only the upper half 5 elements have to be computed since symetric.
Currently done as for-loops but will eventually be changed to vector
indexing, once the topography has been figured out.
Created on Oct, 20th 2015
@author: dominiquef
"""
NewtG=6.6738e-3
eps = 1e-10 # add a small value to the locations to avoid /0
nC = Xn.shape[0]
# Pre-allocate space for 1D array
T = np.zeros((1,nC))
dz = rxLoc[2] - Zn + eps
dy = Yn - rxLoc[1] + eps
dx = Xn - rxLoc[0] + eps
# Compute contribution from each corners
for aa in range(2):
for bb in range(2):
for cc in range(2):
r = (dx[:,aa] ** 2 + dy[:,bb] ** 2 + dz[:,cc] ** 2) ** (0.50)
T = T - NewtG * (-1) ** aa * (-1) ** bb * (-1) ** cc * (dx[:,aa] * np.log ( dy[:,bb] + r ) + dy[:,bb] * np.log ( dx[:,aa] + r ) - dz[:,cc] * np.arctan ( dx[:,aa] * dy[:,bb] / ( dz[:,cc] * r )))
return T
def progress(iter,prog,final):
"""
progress(iter,prog,final)
Function measuring the progress of a process and print to screen the %.
Useful to estimate the remaining runtime of a large problem.
Created on Dec, 20th 2015
@author: dominiquef
"""
arg = np.floor(float(iter)/float(final)*10.);
if arg > prog:
strg = "Done " + str(arg*10) + " %"
print strg
prog = arg;
return prog
def writeUBCobs(filename,survey,d):
"""
writeUBCobs(filename,survey,d)
Function writing an observation file in UBC-GRAV3D format.
INPUT
filename : Name of out file including directory
survey
flag : dobs | dpred
OUTPUT
Obsfile
Created on Dec, 27th 2015
@author: dominiquef
"""
rxLoc = survey.srcField.rxList[0].locs
wd = survey.std
data = np.c_[rxLoc , d , wd]
with file(filename,'w') as fid:
fid.write('%i\n' %len(d) )
np.savetxt(fid, data, fmt='%e',delimiter=' ',newline='\n')
print "Observation file saved to: " + filename
def getActiveTopo(mesh,topo,flag):
"""
getActiveTopo(mesh,topo)
Function creates an active cell model from topography
INPUT
mesh : Mesh in SimPEG format
topo : Scatter points defining topography [x,y,z]
OUTPUT
actv : Active cell model
Created on Dec, 27th 2015
@founrdo
"""
import scipy.interpolate as interpolation
if (flag=='N'):
Zn = np.zeros((mesh.nNx,mesh.nNy))
# wght = np.zeros((mesh.nNx,mesh.nNy))
cx = mesh.vectorNx
cy = mesh.vectorNy
F = interpolation.NearestNDInterpolator(topo[:,0:2],topo[:,2])
[Y,X] = np.meshgrid(cy,cx)
Zn = F(X,Y)
actv = np.zeros((mesh.nCx, mesh.nCy, mesh.nCz))
if (flag=='N'):
Nz = mesh.vectorNz[1:]
for jj in range(mesh.nCy):
for ii in range(mesh.nCx):
temp = [kk for kk in range(len(Nz)) if np.all(Zn[ii:(ii+2),jj:(jj+2)] > Nz[kk]) ]
actv[ii,jj,temp] = 1
actv = mkvc(actv==1)
inds = np.asarray([inds for inds, elem in enumerate(actv, 1) if elem], dtype = int) - 1
return inds
def plot_obs_2D(survey,varstr):
""" Function plot_obs(rxLoc,d,wd)
Generate a 2d interpolated plot from scatter points of data
INPUT
rxLoc : Observation locations [x,y,z]
d : Data vector
wd : Uncertainty vector
OUTPUT
figure()
Created on Dec, 27th 2015
@author: dominiquef
"""
from scipy.interpolate import griddata
import pylab as plt
rxLoc = survey.srcField.rxList[0].locs
d = survey.dobs
wd = survey.std
# Create grid of points
x = np.linspace(rxLoc[:,0].min(), rxLoc[:,0].max(), 100)
y = np.linspace(rxLoc[:,1].min(), rxLoc[:,1].max(), 100)
X, Y = np.meshgrid(x,y)
# Interpolate
d_grid = griddata(rxLoc[:,0:2],d,(X,Y), method ='linear')
# Plot result
plt.figure()
plt.subplot()
plt.imshow(d_grid, extent=[x.min(), x.max(), y.min(), y.max()],origin = 'lower')
plt.colorbar(fraction=0.02)
plt.contour(X,Y, d_grid,10)
plt.scatter(rxLoc[:,0],rxLoc[:,1], c=d, s=20)
plt.title(varstr)
plt.gca().set_aspect('equal', adjustable='box')
def readUBCgravObs(obs_file):
"""
Read UBC grav file format
INPUT:
:param fileName, path to the UBC obs grav file
OUTPUT:
:param survey
"""
fid = open(obs_file,'r')
# First line has the number of rows
line = fid.readline()
ndat = np.array(line.split(),dtype=int)
# Pre-allocate space for obsx, obsy, obsz, data, uncert
line = fid.readline()
temp = np.array(line.split(),dtype=float)
d = np.zeros(ndat, dtype=float)
wd = np.zeros(ndat, dtype=float)
locXYZ = np.zeros( (ndat,3), dtype=float)
for ii in range(ndat):
temp = np.array(line.split(),dtype=float)
locXYZ[ii,:] = temp[:3]
d[ii] = temp[3]
wd[ii] = temp[4]
line = fid.readline()
rxLoc = GRAV.RxObs(locXYZ)
srcField = GRAV.SrcField([rxLoc])
survey = GRAV.LinearSurvey(srcField)
survey.dobs = d
survey.std = wd
return survey
def read_GRAVinv_inp(input_file):
"""Read input files for forward modeling MAG data with integral form
INPUT:
input_file: File name containing the forward parameter
OUTPUT:
mshfile
obsfile
topofile
start model
ref model
weightfile
chi_target
as, ax ,ay, az
upper, lower bounds
lp, lqx, lqy, lqz
# All files should be in the working directory, otherwise the path must
# be specified.
Created on Dec 21th, 2015
@author: dominiquef
"""
fid = open(input_file,'r')
# Line 1
line = fid.readline()
l_input = line.split('!')
mshfile = l_input[0].rstrip()
# Line 2
line = fid.readline()
l_input = line.split('!')
obsfile = l_input[0].rstrip()
# Line 3
line = fid.readline()
l_input = re.split('[!\s]',line)
if l_input=='null':
topofile = []
else:
topofile = l_input[0].rstrip()
# Line 4
line = fid.readline()
l_input = re.split('[!\s]',line)
if l_input[0]=='VALUE':
mstart = float(l_input[1])
else:
mstart = l_input[0].rstrip()
# Line 5
line = fid.readline()
l_input = re.split('[!\s]',line)
if l_input[0]=='VALUE':
mref = float(l_input[1])
else:
mref = l_input[0].rstrip()
# Line 7
line = fid.readline()
l_input = re.split('[!\s]',line)
if l_input=='DEFAULT':
wgtfile = []
else:
wgtfile = l_input[0].rstrip()
# Line 8
line = fid.readline()
l_input = re.split('[!\s]',line)
chi = float(l_input[0])
# Line 9
line = fid.readline()
l_input = re.split('[!\s]',line)
val = np.array(l_input[0:4])
alphas = val.astype(np.float)
# Line 10
line = fid.readline()
l_input = re.split('[!\s]',line)
if l_input[0]=='VALUE':
val = np.array(l_input[1:3])
bounds = val.astype(np.float)
else:
bounds = l_input[0].rstrip()
# Line 11
line = fid.readline()
l_input = re.split('[!\s]',line)
if l_input[0]=='VALUE':
val = np.array(l_input[1:6])
lpnorms = val.astype(np.float)
else:
lpnorms = l_input[0].rstrip()
return mshfile, obsfile, topofile, mstart, mref, wgtfile, chi, alphas, bounds, lpnorms
simpegPF/Magnetics.py
+1 -1
View File
@@ -169,7 +169,7 @@ class MagneticIntegral(Problem.BaseProblem):
else:
print """Flag must be either 'tmi' | 'xyz' | 'full', please revised"""
print """Flag must be either 'ind' | 'full', please revised"""
return
simpegPF/__init__.py
+2
View File
@@ -1,3 +1,5 @@
import MagAnalytics
import BaseMag
import Magnetics
import BaseGrav
import Gravity