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

Merge branch 'Dom' of https://github.com/simpeg/simpegPF into ref/merge

Browse Source 
# Conflicts:
#	simpegPF/Dev/MAG/Intgrl_MAG_Inv_Driver.py
This commit is contained in:
Rowan Cockett
2016-04-29 17:12:57 -07:00
parent b01b14efb9 6a3d0e4ffe
commit 64610d3a0a
5 changed files with 691 additions and 285 deletions
Download Patch File Download Diff File Expand all files Collapse all files
simpegPF/Dev/GRAV/Intgrl_GRAV_Inv_Driver.py
+25 -30
View File
@@ -3,16 +3,14 @@ from SimPEG import *
import simpegPF as PF
import pylab as plt
import os
home_dir = 'C:\\Users\\dominiquef.MIRAGEOSCIENCE\\ownCloud\\Research\\CraigModel\\GRAV\\checkerboard_tests'
home_dir = '.'
inpfile = 'PYGRAV3D_inv.inp'
dsep = '\\'
os.chdir(home_dir)
dsep = os.path.sep
plt.close('all')
#%% User input
# Initial beta
beta_in = 1e-2
# Treshold values for compact norm
eps_p = 1e-3 # Small model values
eps_q = 1e-3 # Small model gradient
@@ -122,22 +120,20 @@ reg.mref = mref
reg.wght = wr
# Create pre-conditioner
diagA = np.sum(prob.G**2.,axis=0) + beta_in*(reg.W.T*reg.W).diagonal()
PC = Utils.sdiag(diagA**-1.)
# Data misfit function
dmis = DataMisfit.l2_DataMisfit(survey)
dmis.Wd = 1./wd
opt = Optimization.ProjectedGNCG(maxIter=20,lower=bounds[0],upper=bounds[1], maxIterCG= 50, tolCG = 1e-4)
opt.approxHinv = PC
invProb = InvProblem.BaseInvProblem(dmis, reg, opt)
invProb = InvProblem.BaseInvProblem(dmis, reg, opt, beta = beta_in)
# Add directives to the inversion
beta = Directives.BetaSchedule(coolingFactor=2, coolingRate=1)
beta_init = Directives.BetaEstimate_ByEig()
target = Directives.TargetMisfit()
update_Jacobi = Directives.Update_lin_PreCond(onlyOnStart=True)
inv = Inversion.BaseInversion(invProb, directiveList=[beta,target])
inv = Inversion.BaseInversion(invProb, directiveList=[beta,target,beta_init,update_Jacobi])
m0 = mstart
@@ -148,7 +144,6 @@ 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)
@@ -204,25 +199,23 @@ reg.eps_p = eps_p
reg.eps_q = eps_q
reg.norms = lpnorms
diagA = np.sum(prob.G**2.,axis=0) + beta_in*(reg.W.T*reg.W).diagonal()
PC = Utils.sdiag(diagA**-1.)
#reg.alpha_s = 1.
dmis = DataMisfit.l2_DataMisfit(survey)
dmis.Wd = 1./wd
opt = Optimization.ProjectedGNCG(maxIter=20 , maxIterLS = 20,lower=bounds[0],upper=bounds[1], maxIterCG= 50, tolCG = 1e-4)
opt.approxHinv = PC
#opt.phim_last = reg.eval(mrec)
# opt = Optimization.InexactGaussNewton(maxIter=6)
opt = Optimization.ProjectedGNCG(maxIter=20 , maxIterLS = 20,lower=bounds[0],upper=bounds[1], maxIterCG= 50, tolCG = 1e-4)
invProb = InvProblem.BaseInvProblem(dmis, reg, opt, beta = invProb.beta)
# Add directives to the inversion
beta = Directives.BetaSchedule(coolingFactor=1, coolingRate=1)
#betaest = Directives.BetaEstimate_ByEig()
update_beta = Directives.Scale_Beta(tol = 0.05)
target = Directives.TargetMisfit()
IRLS =Directives.Update_IRLS( phi_m_last = phim, phi_d_last = phid )
update_Jacobi = Directives.Update_lin_PreCond(onlyOnStart=False)
save_log = Directives.SaveOutputEveryIteration()
save_log.fileName = 'LogName_blabla'
inv = Inversion.BaseInversion(invProb, directiveList=[beta,IRLS])
inv = Inversion.BaseInversion(invProb, directiveList=[beta,IRLS,update_beta,update_Jacobi,save_log])
m0 = mrec
@@ -242,24 +235,26 @@ print "Final misfit:" + str(np.sum( ((d-pred)/wd)**2. ) )
#%% Plot out a section of the model
yslice = midx
m_out[m_out==-100] = np.nan
plt.figure()
ax = plt.subplot(221)
mesh.plotSlice(m_out, ax = ax, normal = 'Z', ind=-18, clim = (vmin,vmax) , pcolorOpts = {'cmap':'bwr'})
mesh.plotSlice(m_out, ax = ax, normal = 'Z', ind=-5, clim = (mrec.min(), mrec.max()), pcolorOpts = {'cmap':'bwr'})
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[-18]) + ' m')
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=-23, clim = (vmin,vmax), pcolorOpts = {'cmap':'bwr'})
mesh.plotSlice(m_out, ax = ax, normal = 'Z', ind=-8, clim = (mrec.min(), mrec.max()), pcolorOpts = {'cmap':'bwr'})
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[-23]) + ' m')
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 = (vmin,vmax), pcolorOpts = {'cmap':'bwr'})
mesh.plotSlice(m_out, ax = ax, normal = 'Y', ind=yslice, clim = (mrec.min(), mrec.max()), pcolorOpts = {'cmap':'bwr'})
plt.title('Cross Section')
plt.xlabel('x');plt.ylabel('z')
plt.gca().set_aspect('equal', adjustable='box')
simpegPF/Dev/GRAV/PYGRAV3D_inv.inp
+1 -1
View File
@@ -1,7 +1,7 @@
Mesh_10m.msh ! Mesh file
Grav.dat ! Obsfile
Gaussian.topo ! Topofile | null
VALUE 1e-4 ! Starting model
VALUE 0.0 ! Starting model
VALUE 0.0 ! Reference model
DEFAULT ! Cell based weight file
1 ! target chi factor | DEFAULT=1
simpegPF/Dev/GRAV/Pred0.dat
+617 -192
View File
@@ -1,192 +1,617 @@
191
3.624724e+05 6.008600e+06 2.185513e+03 2.950277e-02 5.000000e-02
3.630474e+05 6.007742e+06 2.182354e+03 2.958855e-02 5.000000e-02
3.658851e+05 6.010287e+06 2.181814e+03 2.909199e-02 5.000000e-02
3.634466e+05 6.005839e+06 2.190788e+03 2.947114e-02 5.000000e-02
3.641340e+05 6.005157e+06 2.173695e+03 2.936905e-02 5.000000e-02
3.644170e+05 6.003660e+06 2.189776e+03 2.896865e-02 5.000000e-02
3.673445e+05 6.005950e+06 2.194062e+03 2.885706e-02 5.000000e-02
3.605751e+05 6.011702e+06 2.267910e+03 2.855588e-02 5.000000e-02
3.687754e+05 6.014567e+06 2.230211e+03 2.611820e-02 5.000000e-02
3.600941e+05 6.011233e+06 2.205465e+03 2.849530e-02 5.000000e-02
3.662178e+05 6.013153e+06 2.310513e+03 2.835436e-02 5.000000e-02
3.623110e+05 6.012305e+06 2.287369e+03 2.878394e-02 5.000000e-02
3.666170e+05 6.013924e+06 2.328716e+03 2.779231e-02 5.000000e-02
3.633072e+05 6.013055e+06 2.350074e+03 2.838244e-02 5.000000e-02
3.593416e+05 6.011481e+06 2.188002e+03 2.819680e-02 5.000000e-02
3.592337e+05 6.010404e+06 2.202372e+03 2.842105e-02 5.000000e-02
3.670648e+05 6.010916e+06 2.176537e+03 2.865192e-02 5.000000e-02
3.594689e+05 6.013023e+06 2.187709e+03 2.760522e-02 5.000000e-02
3.606836e+05 6.009231e+06 2.185043e+03 2.913644e-02 5.000000e-02
3.634058e+05 6.008465e+06 2.171398e+03 2.959731e-02 5.000000e-02
3.630400e+05 6.007708e+06 2.181436e+03 2.959071e-02 5.000000e-02
3.634243e+05 6.005860e+06 2.188994e+03 2.947753e-02 5.000000e-02
3.619294e+05 6.007571e+06 2.196007e+03 2.945406e-02 5.000000e-02
3.613222e+05 6.010390e+06 2.161111e+03 2.912337e-02 5.000000e-02
3.609065e+05 6.005685e+06 2.192528e+03 2.913821e-02 5.000000e-02
3.600812e+05 6.004739e+06 2.201676e+03 2.870265e-02 5.000000e-02
3.582056e+05 6.003746e+06 2.328486e+03 2.813477e-02 5.000000e-02
3.644598e+05 6.010922e+06 2.168886e+03 2.914594e-02 5.000000e-02
3.660609e+05 6.004171e+06 2.184777e+03 2.884816e-02 5.000000e-02
3.662338e+05 6.004030e+06 2.186003e+03 2.876446e-02 5.000000e-02
3.594967e+05 6.014888e+06 2.161400e+03 2.641712e-02 5.000000e-02
3.596059e+05 6.013750e+06 2.154978e+03 2.730693e-02 5.000000e-02
3.598438e+05 6.013214e+06 2.185253e+03 2.765375e-02 5.000000e-02
3.599208e+05 6.012727e+06 2.203143e+03 2.790024e-02 5.000000e-02
3.600763e+05 6.012250e+06 2.220439e+03 2.817421e-02 5.000000e-02
3.611366e+05 6.011684e+06 2.297049e+03 2.871383e-02 5.000000e-02
3.615742e+05 6.011880e+06 2.309140e+03 2.912090e-02 5.000000e-02
3.620552e+05 6.012062e+06 2.303044e+03 2.912972e-02 5.000000e-02
3.626264e+05 6.012780e+06 2.289942e+03 2.863952e-02 5.000000e-02
3.630185e+05 6.012978e+06 2.325430e+03 2.841604e-02 5.000000e-02
3.634810e+05 6.013343e+06 2.355786e+03 2.840198e-02 5.000000e-02
3.639229e+05 6.013535e+06 2.357851e+03 2.835332e-02 5.000000e-02
3.644825e+05 6.012914e+06 2.363058e+03 2.837733e-02 5.000000e-02
3.646570e+05 6.012452e+06 2.351115e+03 2.838988e-02 5.000000e-02
3.649237e+05 6.012810e+06 2.359124e+03 2.843872e-02 5.000000e-02
3.651010e+05 6.013278e+06 2.358159e+03 2.837036e-02 5.000000e-02
3.656116e+05 6.013387e+06 2.346855e+03 2.825637e-02 5.000000e-02
3.661738e+05 6.013631e+06 2.317519e+03 2.820977e-02 5.000000e-02
3.668565e+05 6.014278e+06 2.312226e+03 2.754386e-02 5.000000e-02
3.674190e+05 6.014362e+06 2.298255e+03 2.724274e-02 5.000000e-02
3.679217e+05 6.014372e+06 2.261726e+03 2.672379e-02 5.000000e-02
3.684490e+05 6.014348e+06 2.244935e+03 2.643617e-02 5.000000e-02
3.691523e+05 6.014957e+06 2.201091e+03 2.568247e-02 5.000000e-02
3.696587e+05 6.014932e+06 2.201995e+03 2.542427e-02 5.000000e-02
3.590002e+05 6.012880e+06 2.180407e+03 2.753383e-02 5.000000e-02
3.588834e+05 6.012462e+06 2.177161e+03 2.767358e-02 5.000000e-02
3.585091e+05 6.012213e+06 2.167542e+03 2.764672e-02 5.000000e-02
3.589367e+05 6.011716e+06 2.184280e+03 2.797252e-02 5.000000e-02
3.588939e+05 6.011250e+06 2.215462e+03 2.804477e-02 5.000000e-02
3.585927e+05 6.010922e+06 2.213598e+03 2.802561e-02 5.000000e-02
3.582362e+05 6.010625e+06 2.228342e+03 2.800097e-02 5.000000e-02
3.577469e+05 6.010694e+06 2.268556e+03 2.831717e-02 5.000000e-02
3.574252e+05 6.011061e+06 2.307562e+03 2.796307e-02 5.000000e-02
3.584220e+05 6.011412e+06 2.303356e+03 2.808073e-02 5.000000e-02
3.593649e+05 6.011057e+06 2.167662e+03 2.838074e-02 5.000000e-02
3.591553e+05 6.009881e+06 2.166452e+03 2.857786e-02 5.000000e-02
3.595398e+05 6.009477e+06 2.160818e+03 2.880609e-02 5.000000e-02
3.600429e+05 6.009469e+06 2.165836e+03 2.896599e-02 5.000000e-02
3.609785e+05 6.009612e+06 2.160491e+03 2.920621e-02 5.000000e-02
3.612110e+05 6.009171e+06 2.161913e+03 2.931579e-02 5.000000e-02
3.615351e+05 6.008806e+06 2.162689e+03 2.941500e-02 5.000000e-02
3.616756e+05 6.008330e+06 2.161829e+03 2.947662e-02 5.000000e-02
3.586542e+05 6.009829e+06 2.201714e+03 2.833175e-02 5.000000e-02
3.583539e+05 6.009427e+06 2.297539e+03 2.866337e-02 5.000000e-02
3.579539e+05 6.009058e+06 2.333486e+03 2.859184e-02 5.000000e-02
3.574419e+05 6.008757e+06 2.449519e+03 2.803439e-02 5.000000e-02
3.571929e+05 6.008326e+06 2.470826e+03 2.787830e-02 5.000000e-02
3.568360e+05 6.008077e+06 2.500811e+03 2.762602e-02 5.000000e-02
3.613670e+05 6.006269e+06 2.186354e+03 2.931716e-02 5.000000e-02
3.609547e+05 6.005975e+06 2.201936e+03 2.916378e-02 5.000000e-02
3.605453e+05 6.005187e+06 2.192131e+03 2.896361e-02 5.000000e-02
3.596760e+05 6.004399e+06 2.216005e+03 2.846017e-02 5.000000e-02
3.591839e+05 6.004206e+06 2.223183e+03 2.822088e-02 5.000000e-02
3.586900e+05 6.004318e+06 2.240749e+03 2.800929e-02 5.000000e-02
3.582144e+05 6.004122e+06 2.267646e+03 2.815678e-02 5.000000e-02
3.574795e+05 6.004361e+06 2.273342e+03 2.785160e-02 5.000000e-02
3.573075e+05 6.005061e+06 2.257032e+03 2.768883e-02 5.000000e-02
3.577691e+05 6.005261e+06 2.257991e+03 2.782176e-02 5.000000e-02
3.582753e+05 6.005377e+06 2.253481e+03 2.824865e-02 5.000000e-02
3.589901e+05 6.005582e+06 2.225082e+03 2.843000e-02 5.000000e-02
3.594262e+05 6.005850e+06 2.208702e+03 2.867859e-02 5.000000e-02
3.599277e+05 6.005832e+06 2.206934e+03 2.884661e-02 5.000000e-02
3.616705e+05 6.006416e+06 2.159826e+03 2.944395e-02 5.000000e-02
3.689559e+05 6.013437e+06 2.160188e+03 2.686265e-02 5.000000e-02
3.689809e+05 6.012858e+06 2.160165e+03 2.713001e-02 5.000000e-02
3.689575e+05 6.012164e+06 2.166898e+03 2.742852e-02 5.000000e-02
3.687815e+05 6.011610e+06 2.169930e+03 2.772659e-02 5.000000e-02
3.685491e+05 6.011021e+06 2.161041e+03 2.805256e-02 5.000000e-02
3.683284e+05 6.010552e+06 2.162389e+03 2.829110e-02 5.000000e-02
3.681385e+05 6.009930e+06 2.160917e+03 2.852168e-02 5.000000e-02
3.683213e+05 6.009170e+06 2.161018e+03 2.855719e-02 5.000000e-02
3.681582e+05 6.008637e+06 2.160607e+03 2.868224e-02 5.000000e-02
3.678554e+05 6.008233e+06 2.161901e+03 2.883354e-02 5.000000e-02
3.676810e+05 6.007757e+06 2.161971e+03 2.891626e-02 5.000000e-02
3.675351e+05 6.007272e+06 2.161757e+03 2.896734e-02 5.000000e-02
3.674284e+05 6.006771e+06 2.161644e+03 2.898216e-02 5.000000e-02
3.672650e+05 6.006275e+06 2.161193e+03 2.899572e-02 5.000000e-02
3.670451e+05 6.005801e+06 2.161079e+03 2.900593e-02 5.000000e-02
3.668871e+05 6.005298e+06 2.168517e+03 2.895119e-02 5.000000e-02
3.669246e+05 6.004782e+06 2.160798e+03 2.884703e-02 5.000000e-02
3.674886e+05 6.004466e+06 2.160969e+03 2.858275e-02 5.000000e-02
3.680010e+05 6.007243e+06 2.186598e+03 2.873252e-02 5.000000e-02
3.685053e+05 6.007066e+06 2.247928e+03 2.844525e-02 5.000000e-02
3.689689e+05 6.006886e+06 2.294719e+03 2.876087e-02 5.000000e-02
3.694464e+05 6.006798e+06 2.340813e+03 2.847357e-02 5.000000e-02
3.699384e+05 6.006629e+06 2.452104e+03 2.790729e-02 5.000000e-02
3.702331e+05 6.006504e+06 2.478823e+03 2.768065e-02 5.000000e-02
3.648724e+05 6.003990e+06 2.162273e+03 2.903588e-02 5.000000e-02
3.615819e+05 6.006914e+06 2.166497e+03 2.945036e-02 5.000000e-02
3.619111e+05 6.005892e+06 2.161213e+03 2.943030e-02 5.000000e-02
3.624320e+05 6.006129e+06 2.161200e+03 2.951822e-02 5.000000e-02
3.628240e+05 6.006286e+06 2.162131e+03 2.955949e-02 5.000000e-02
3.632946e+05 6.006160e+06 2.161080e+03 2.956698e-02 5.000000e-02
3.676218e+05 6.013497e+06 2.169281e+03 2.740393e-02 5.000000e-02
3.679297e+05 6.013094e+06 2.168942e+03 2.751932e-02 5.000000e-02
3.677949e+05 6.012264e+06 2.167650e+03 2.796881e-02 5.000000e-02
3.672235e+05 6.011537e+06 2.168687e+03 2.843018e-02 5.000000e-02
3.666300e+05 6.012199e+06 2.168907e+03 2.836406e-02 5.000000e-02
3.659040e+05 6.011713e+06 2.165875e+03 2.871230e-02 5.000000e-02
3.659362e+05 6.011209e+06 2.167131e+03 2.887116e-02 5.000000e-02
3.657312e+05 6.010915e+06 2.166870e+03 2.899740e-02 5.000000e-02
3.655597e+05 6.011327e+06 2.167944e+03 2.891109e-02 5.000000e-02
3.650121e+05 6.011050e+06 2.167057e+03 2.906459e-02 5.000000e-02
3.644476e+05 6.011501e+06 2.165267e+03 2.897208e-02 5.000000e-02
3.642410e+05 6.011986e+06 2.164931e+03 2.880292e-02 5.000000e-02
3.640583e+05 6.012445e+06 2.165348e+03 2.861521e-02 5.000000e-02
3.637110e+05 6.012686e+06 2.165028e+03 2.851203e-02 5.000000e-02
3.632662e+05 6.012413e+06 2.165095e+03 2.863726e-02 5.000000e-02
3.628353e+05 6.012198e+06 2.165437e+03 2.872215e-02 5.000000e-02
3.626366e+05 6.011643e+06 2.167954e+03 2.892334e-02 5.000000e-02
3.619752e+05 6.011006e+06 2.169750e+03 2.905700e-02 5.000000e-02
3.608859e+05 6.012099e+06 2.332313e+03 2.884199e-02 5.000000e-02
3.607787e+05 6.012623e+06 2.334431e+03 2.864554e-02 5.000000e-02
3.604558e+05 6.012994e+06 2.340932e+03 2.825683e-02 5.000000e-02
3.604541e+05 6.013514e+06 2.417598e+03 2.775451e-02 5.000000e-02
3.605505e+05 6.014052e+06 2.428877e+03 2.743902e-02 5.000000e-02
3.605665e+05 6.014532e+06 2.431259e+03 2.703843e-02 5.000000e-02
3.609731e+05 6.014173e+06 2.431723e+03 2.750775e-02 5.000000e-02
3.614862e+05 6.014207e+06 2.463004e+03 2.751270e-02 5.000000e-02
3.619207e+05 6.014382e+06 2.454613e+03 2.747475e-02 5.000000e-02
3.624145e+05 6.014385e+06 2.450971e+03 2.754979e-02 5.000000e-02
3.640306e+05 6.003363e+06 2.230716e+03 2.912848e-02 5.000000e-02
3.636639e+05 6.003018e+06 2.254717e+03 2.944835e-02 5.000000e-02
3.632341e+05 6.002875e+06 2.270788e+03 2.935057e-02 5.000000e-02
3.630920e+05 6.002515e+06 2.342944e+03 2.897828e-02 5.000000e-02
3.627738e+05 6.002067e+06 2.449163e+03 2.845754e-02 5.000000e-02
3.625173e+05 6.007599e+06 2.169206e+03 2.958475e-02 5.000000e-02
3.635370e+05 6.008139e+06 2.164473e+03 2.963331e-02 5.000000e-02
3.627917e+05 6.008156e+06 2.169742e+03 2.959324e-02 5.000000e-02
3.633881e+05 6.009013e+06 2.164473e+03 2.956369e-02 5.000000e-02
3.663215e+05 6.003654e+06 2.202612e+03 2.858425e-02 5.000000e-02
3.664925e+05 6.003153e+06 2.185419e+03 2.839115e-02 5.000000e-02
3.664874e+05 6.002625e+06 2.175218e+03 2.818385e-02 5.000000e-02
3.664629e+05 6.002138e+06 2.182963e+03 2.793174e-02 5.000000e-02
3.665085e+05 6.001658e+06 2.183867e+03 2.764954e-02 5.000000e-02
3.667288e+05 6.001228e+06 2.186336e+03 2.730574e-02 5.000000e-02
3.669117e+05 6.000761e+06 2.194478e+03 2.690493e-02 5.000000e-02
3.672005e+05 6.000532e+06 2.193351e+03 2.667164e-02 5.000000e-02
3.615060e+05 6.010861e+06 2.164254e+03 2.903521e-02 5.000000e-02
3.608114e+05 6.010798e+06 2.166447e+03 2.890123e-02 5.000000e-02
3.678002e+05 6.003356e+06 2.180167e+03 2.808437e-02 5.000000e-02
3.682322e+05 6.003116e+06 2.219280e+03 2.777212e-02 5.000000e-02
3.685762e+05 6.002753e+06 2.243722e+03 2.751404e-02 5.000000e-02
3.690196e+05 6.002528e+06 2.262227e+03 2.771937e-02 5.000000e-02
3.604802e+05 6.004043e+06 2.252471e+03 2.865938e-02 5.000000e-02
3.609087e+05 6.003672e+06 2.320061e+03 2.916318e-02 5.000000e-02
3.609935e+05 6.003061e+06 2.341764e+03 2.896186e-02 5.000000e-02
3.610941e+05 6.002503e+06 2.353045e+03 2.870830e-02 5.000000e-02
3.613520e+05 6.002090e+06 2.402026e+03 2.841739e-02 5.000000e-02
3.616648e+05 6.001616e+06 2.460820e+03 2.804598e-02 5.000000e-02
3.621238e+05 6.001462e+06 2.490509e+03 2.797261e-02 5.000000e-02
3.657493e+05 6.004426e+06 2.169680e+03 2.901485e-02 5.000000e-02
3.636824e+05 6.000022e+06 2.865274e+03 2.708391e-02 5.000000e-02
3.571072e+05 6.007140e+06 2.504079e+03 2.777789e-02 5.000000e-02
3.687904e+05 6.004471e+06 2.394426e+03 2.822128e-02 5.000000e-02
3.660711e+05 6.001991e+06 2.247601e+03 2.794515e-02 5.000000e-02
3.565341e+05 6.013529e+06 2.495766e+03 2.596333e-02 5.000000e-02
3.572975e+05 6.002759e+06 2.586419e+03 2.739037e-02 5.000000e-02
3.602847e+05 6.000220e+06 2.741396e+03 2.681741e-02 5.000000e-02
3.592838e+05 6.002004e+06 2.769861e+03 2.752123e-02 5.000000e-02
3.624843e+05 5.999891e+06 2.693787e+03 2.696451e-02 5.000000e-02
616
1.050000e+02 1.550000e+02 2.408312e+02 0.000000e+00 2.000000e-03
1.150000e+02 1.550000e+02 2.440469e+02 0.000000e+00 2.000000e-03
1.250000e+02 1.550000e+02 2.471986e+02 0.000000e+00 2.000000e-03
1.350000e+02 1.550000e+02 2.502569e+02 0.000000e+00 2.000000e-03
1.450000e+02 1.550000e+02 2.531923e+02 0.000000e+00 2.000000e-03
1.550000e+02 1.550000e+02 2.559753e+02 0.000000e+00 2.000000e-03
1.650000e+02 1.550000e+02 2.585770e+02 0.000000e+00 2.000000e-03
1.750000e+02 1.550000e+02 2.609695e+02 0.000000e+00 2.000000e-03
1.850000e+02 1.550000e+02 2.631267e+02 0.000000e+00 2.000000e-03
1.950000e+02 1.550000e+02 2.650244e+02 0.000000e+00 2.000000e-03
2.050000e+02 1.550000e+02 2.666411e+02 0.000000e+00 2.000000e-03
2.150000e+02 1.550000e+02 2.679578e+02 0.000000e+00 2.000000e-03
2.250000e+02 1.550000e+02 2.689593e+02 0.000000e+00 2.000000e-03
2.350000e+02 1.550000e+02 2.696337e+02 0.000000e+00 2.000000e-03
2.450000e+02 1.550000e+02 2.699729e+02 0.000000e+00 2.000000e-03
2.550000e+02 1.550000e+02 2.699729e+02 0.000000e+00 2.000000e-03
2.650000e+02 1.550000e+02 2.696337e+02 0.000000e+00 2.000000e-03
2.750000e+02 1.550000e+02 2.689593e+02 0.000000e+00 2.000000e-03
2.850000e+02 1.550000e+02 2.679578e+02 0.000000e+00 2.000000e-03
2.950000e+02 1.550000e+02 2.666411e+02 0.000000e+00 2.000000e-03
3.050000e+02 1.550000e+02 2.650244e+02 0.000000e+00 2.000000e-03
3.150000e+02 1.550000e+02 2.631267e+02 0.000000e+00 2.000000e-03
3.250000e+02 1.550000e+02 2.609695e+02 0.000000e+00 2.000000e-03
3.350000e+02 1.550000e+02 2.585770e+02 0.000000e+00 2.000000e-03
3.450000e+02 1.550000e+02 2.559753e+02 0.000000e+00 2.000000e-03
3.550000e+02 1.550000e+02 2.531923e+02 0.000000e+00 2.000000e-03
3.650000e+02 1.550000e+02 2.502569e+02 0.000000e+00 2.000000e-03
3.750000e+02 1.550000e+02 2.471986e+02 0.000000e+00 2.000000e-03
1.050000e+02 1.650000e+02 2.427601e+02 0.000000e+00 2.000000e-03
1.150000e+02 1.650000e+02 2.460869e+02 0.000000e+00 2.000000e-03
1.250000e+02 1.650000e+02 2.493475e+02 0.000000e+00 2.000000e-03
1.350000e+02 1.650000e+02 2.525115e+02 0.000000e+00 2.000000e-03
1.450000e+02 1.650000e+02 2.555483e+02 0.000000e+00 2.000000e-03
1.550000e+02 1.650000e+02 2.584275e+02 0.000000e+00 2.000000e-03
1.650000e+02 1.650000e+02 2.611190e+02 0.000000e+00 2.000000e-03
1.750000e+02 1.650000e+02 2.635942e+02 0.000000e+00 2.000000e-03
1.850000e+02 1.650000e+02 2.658260e+02 0.000000e+00 2.000000e-03
1.950000e+02 1.650000e+02 2.677892e+02 0.000000e+00 2.000000e-03
2.050000e+02 1.650000e+02 2.694617e+02 0.000000e+00 2.000000e-03
2.150000e+02 1.650000e+02 2.708240e+02 0.000000e+00 2.000000e-03
2.250000e+02 1.650000e+02 2.718600e+02 0.000000e+00 2.000000e-03
2.350000e+02 1.650000e+02 2.725577e+02 0.000000e+00 2.000000e-03
2.450000e+02 1.650000e+02 2.729086e+02 0.000000e+00 2.000000e-03
2.550000e+02 1.650000e+02 2.729086e+02 0.000000e+00 2.000000e-03
2.650000e+02 1.650000e+02 2.725577e+02 0.000000e+00 2.000000e-03
2.750000e+02 1.650000e+02 2.718600e+02 0.000000e+00 2.000000e-03
2.850000e+02 1.650000e+02 2.708240e+02 0.000000e+00 2.000000e-03
2.950000e+02 1.650000e+02 2.694617e+02 0.000000e+00 2.000000e-03
3.050000e+02 1.650000e+02 2.677892e+02 0.000000e+00 2.000000e-03
3.150000e+02 1.650000e+02 2.658260e+02 0.000000e+00 2.000000e-03
3.250000e+02 1.650000e+02 2.635942e+02 0.000000e+00 2.000000e-03
3.350000e+02 1.650000e+02 2.611190e+02 0.000000e+00 2.000000e-03
3.450000e+02 1.650000e+02 2.584275e+02 0.000000e+00 2.000000e-03
3.550000e+02 1.650000e+02 2.555483e+02 0.000000e+00 2.000000e-03
3.650000e+02 1.650000e+02 2.525115e+02 0.000000e+00 2.000000e-03
3.750000e+02 1.650000e+02 2.493475e+02 0.000000e+00 2.000000e-03
1.050000e+02 1.750000e+02 2.445174e+02 0.000000e+00 2.000000e-03
1.150000e+02 1.750000e+02 2.479454e+02 0.000000e+00 2.000000e-03
1.250000e+02 1.750000e+02 2.513052e+02 0.000000e+00 2.000000e-03
1.350000e+02 1.750000e+02 2.545655e+02 0.000000e+00 2.000000e-03
1.450000e+02 1.750000e+02 2.576947e+02 0.000000e+00 2.000000e-03
1.550000e+02 1.750000e+02 2.606614e+02 0.000000e+00 2.000000e-03
1.650000e+02 1.750000e+02 2.634348e+02 0.000000e+00 2.000000e-03
1.750000e+02 1.750000e+02 2.659854e+02 0.000000e+00 2.000000e-03
1.850000e+02 1.750000e+02 2.682850e+02 0.000000e+00 2.000000e-03
1.950000e+02 1.750000e+02 2.703080e+02 0.000000e+00 2.000000e-03
2.050000e+02 1.750000e+02 2.720313e+02 0.000000e+00 2.000000e-03
2.150000e+02 1.750000e+02 2.734350e+02 0.000000e+00 2.000000e-03
2.250000e+02 1.750000e+02 2.745027e+02 0.000000e+00 2.000000e-03
2.350000e+02 1.750000e+02 2.752216e+02 0.000000e+00 2.000000e-03
2.450000e+02 1.750000e+02 2.755831e+02 0.000000e+00 2.000000e-03
2.550000e+02 1.750000e+02 2.755831e+02 0.000000e+00 2.000000e-03
2.650000e+02 1.750000e+02 2.752216e+02 0.000000e+00 2.000000e-03
2.750000e+02 1.750000e+02 2.745027e+02 0.000000e+00 2.000000e-03
2.850000e+02 1.750000e+02 2.734350e+02 0.000000e+00 2.000000e-03
2.950000e+02 1.750000e+02 2.720313e+02 0.000000e+00 2.000000e-03
3.050000e+02 1.750000e+02 2.703080e+02 0.000000e+00 2.000000e-03
3.150000e+02 1.750000e+02 2.682850e+02 0.000000e+00 2.000000e-03
3.250000e+02 1.750000e+02 2.659854e+02 0.000000e+00 2.000000e-03
3.350000e+02 1.750000e+02 2.634348e+02 0.000000e+00 2.000000e-03
3.450000e+02 1.750000e+02 2.606614e+02 0.000000e+00 2.000000e-03
3.550000e+02 1.750000e+02 2.576947e+02 0.000000e+00 2.000000e-03
3.650000e+02 1.750000e+02 2.545655e+02 0.000000e+00 2.000000e-03
3.750000e+02 1.750000e+02 2.513052e+02 0.000000e+00 2.000000e-03
1.050000e+02 1.850000e+02 2.460835e+02 0.000000e+00 2.000000e-03
1.150000e+02 1.850000e+02 2.496018e+02 0.000000e+00 2.000000e-03
1.250000e+02 1.850000e+02 2.530500e+02 0.000000e+00 2.000000e-03
1.350000e+02 1.850000e+02 2.563960e+02 0.000000e+00 2.000000e-03
1.450000e+02 1.850000e+02 2.596076e+02 0.000000e+00 2.000000e-03
1.550000e+02 1.850000e+02 2.626524e+02 0.000000e+00 2.000000e-03
1.650000e+02 1.850000e+02 2.654988e+02 0.000000e+00 2.000000e-03
1.750000e+02 1.850000e+02 2.681164e+02 0.000000e+00 2.000000e-03
1.850000e+02 1.850000e+02 2.704766e+02 0.000000e+00 2.000000e-03
1.950000e+02 1.850000e+02 2.725528e+02 0.000000e+00 2.000000e-03
2.050000e+02 1.850000e+02 2.743215e+02 0.000000e+00 2.000000e-03
2.150000e+02 1.850000e+02 2.757621e+02 0.000000e+00 2.000000e-03
2.250000e+02 1.850000e+02 2.768578e+02 0.000000e+00 2.000000e-03
2.350000e+02 1.850000e+02 2.775956e+02 0.000000e+00 2.000000e-03
2.450000e+02 1.850000e+02 2.779668e+02 0.000000e+00 2.000000e-03
2.550000e+02 1.850000e+02 2.779668e+02 0.000000e+00 2.000000e-03
2.650000e+02 1.850000e+02 2.775956e+02 0.000000e+00 2.000000e-03
2.750000e+02 1.850000e+02 2.768578e+02 0.000000e+00 2.000000e-03
2.850000e+02 1.850000e+02 2.757621e+02 0.000000e+00 2.000000e-03
2.950000e+02 1.850000e+02 2.743215e+02 0.000000e+00 2.000000e-03
3.050000e+02 1.850000e+02 2.725528e+02 0.000000e+00 2.000000e-03
3.150000e+02 1.850000e+02 2.704766e+02 0.000000e+00 2.000000e-03
3.250000e+02 1.850000e+02 2.681164e+02 0.000000e+00 2.000000e-03
3.350000e+02 1.850000e+02 2.654988e+02 0.000000e+00 2.000000e-03
3.450000e+02 1.850000e+02 2.626524e+02 0.000000e+00 2.000000e-03
3.550000e+02 1.850000e+02 2.596076e+02 0.000000e+00 2.000000e-03
3.650000e+02 1.850000e+02 2.563960e+02 0.000000e+00 2.000000e-03
3.750000e+02 1.850000e+02 2.530500e+02 0.000000e+00 2.000000e-03
1.050000e+02 1.950000e+02 2.474409e+02 0.000000e+00 2.000000e-03
1.150000e+02 1.950000e+02 2.510374e+02 0.000000e+00 2.000000e-03
1.250000e+02 1.950000e+02 2.545621e+02 0.000000e+00 2.000000e-03
1.350000e+02 1.950000e+02 2.579825e+02 0.000000e+00 2.000000e-03
1.450000e+02 1.950000e+02 2.612655e+02 0.000000e+00 2.000000e-03
1.550000e+02 1.950000e+02 2.643779e+02 0.000000e+00 2.000000e-03
1.650000e+02 1.950000e+02 2.672876e+02 0.000000e+00 2.000000e-03
1.750000e+02 1.950000e+02 2.699634e+02 0.000000e+00 2.000000e-03
1.850000e+02 1.950000e+02 2.723760e+02 0.000000e+00 2.000000e-03
1.950000e+02 1.950000e+02 2.744984e+02 0.000000e+00 2.000000e-03
2.050000e+02 1.950000e+02 2.763063e+02 0.000000e+00 2.000000e-03
2.150000e+02 1.950000e+02 2.777790e+02 0.000000e+00 2.000000e-03
2.250000e+02 1.950000e+02 2.788990e+02 0.000000e+00 2.000000e-03
2.350000e+02 1.950000e+02 2.796533e+02 0.000000e+00 2.000000e-03
2.450000e+02 1.950000e+02 2.800326e+02 0.000000e+00 2.000000e-03
2.550000e+02 1.950000e+02 2.800326e+02 0.000000e+00 2.000000e-03
2.650000e+02 1.950000e+02 2.796533e+02 0.000000e+00 2.000000e-03
2.750000e+02 1.950000e+02 2.788990e+02 0.000000e+00 2.000000e-03
2.850000e+02 1.950000e+02 2.777790e+02 0.000000e+00 2.000000e-03
2.950000e+02 1.950000e+02 2.763063e+02 0.000000e+00 2.000000e-03
3.050000e+02 1.950000e+02 2.744984e+02 0.000000e+00 2.000000e-03
3.150000e+02 1.950000e+02 2.723760e+02 0.000000e+00 2.000000e-03
3.250000e+02 1.950000e+02 2.699634e+02 0.000000e+00 2.000000e-03
3.350000e+02 1.950000e+02 2.672876e+02 0.000000e+00 2.000000e-03
3.450000e+02 1.950000e+02 2.643779e+02 0.000000e+00 2.000000e-03
3.550000e+02 1.950000e+02 2.612655e+02 0.000000e+00 2.000000e-03
3.650000e+02 1.950000e+02 2.579825e+02 0.000000e+00 2.000000e-03
3.750000e+02 1.950000e+02 2.545621e+02 0.000000e+00 2.000000e-03
1.050000e+02 2.050000e+02 2.485739e+02 0.000000e+00 2.000000e-03
1.150000e+02 2.050000e+02 2.522356e+02 0.000000e+00 2.000000e-03
1.250000e+02 2.050000e+02 2.558243e+02 0.000000e+00 2.000000e-03
1.350000e+02 2.050000e+02 2.593068e+02 0.000000e+00 2.000000e-03
1.450000e+02 2.050000e+02 2.626493e+02 0.000000e+00 2.000000e-03
1.550000e+02 2.050000e+02 2.658182e+02 0.000000e+00 2.000000e-03
1.650000e+02 2.050000e+02 2.687807e+02 0.000000e+00 2.000000e-03
1.750000e+02 2.050000e+02 2.715050e+02 0.000000e+00 2.000000e-03
1.850000e+02 2.050000e+02 2.739614e+02 0.000000e+00 2.000000e-03
1.950000e+02 2.050000e+02 2.761223e+02 0.000000e+00 2.000000e-03
2.050000e+02 2.050000e+02 2.779631e+02 0.000000e+00 2.000000e-03
2.150000e+02 2.050000e+02 2.794624e+02 0.000000e+00 2.000000e-03
2.250000e+02 2.050000e+02 2.806028e+02 0.000000e+00 2.000000e-03
2.350000e+02 2.050000e+02 2.813707e+02 0.000000e+00 2.000000e-03
2.450000e+02 2.050000e+02 2.817570e+02 0.000000e+00 2.000000e-03
2.550000e+02 2.050000e+02 2.817570e+02 0.000000e+00 2.000000e-03
2.650000e+02 2.050000e+02 2.813707e+02 0.000000e+00 2.000000e-03
2.750000e+02 2.050000e+02 2.806028e+02 0.000000e+00 2.000000e-03
2.850000e+02 2.050000e+02 2.794624e+02 0.000000e+00 2.000000e-03
2.950000e+02 2.050000e+02 2.779631e+02 0.000000e+00 2.000000e-03
3.050000e+02 2.050000e+02 2.761223e+02 0.000000e+00 2.000000e-03
3.150000e+02 2.050000e+02 2.739614e+02 0.000000e+00 2.000000e-03
3.250000e+02 2.050000e+02 2.715050e+02 0.000000e+00 2.000000e-03
3.350000e+02 2.050000e+02 2.687807e+02 0.000000e+00 2.000000e-03
3.450000e+02 2.050000e+02 2.658182e+02 0.000000e+00 2.000000e-03
3.550000e+02 2.050000e+02 2.626493e+02 0.000000e+00 2.000000e-03
3.650000e+02 2.050000e+02 2.593068e+02 0.000000e+00 2.000000e-03
3.750000e+02 2.050000e+02 2.558243e+02 0.000000e+00 2.000000e-03
1.050000e+02 2.150000e+02 2.494692e+02 0.000000e+00 2.000000e-03
1.150000e+02 2.150000e+02 2.531825e+02 0.000000e+00 2.000000e-03
1.250000e+02 2.150000e+02 2.568218e+02 0.000000e+00 2.000000e-03
1.350000e+02 2.150000e+02 2.603533e+02 0.000000e+00 2.000000e-03
1.450000e+02 2.150000e+02 2.637429e+02 0.000000e+00 2.000000e-03
1.550000e+02 2.150000e+02 2.669565e+02 0.000000e+00 2.000000e-03
1.650000e+02 2.150000e+02 2.699606e+02 0.000000e+00 2.000000e-03
1.750000e+02 2.150000e+02 2.727234e+02 0.000000e+00 2.000000e-03
1.850000e+02 2.150000e+02 2.752143e+02 0.000000e+00 2.000000e-03
1.950000e+02 2.150000e+02 2.774056e+02 0.000000e+00 2.000000e-03
2.050000e+02 2.150000e+02 2.792724e+02 0.000000e+00 2.000000e-03
2.150000e+02 2.150000e+02 2.807928e+02 0.000000e+00 2.000000e-03
2.250000e+02 2.150000e+02 2.819493e+02 0.000000e+00 2.000000e-03
2.350000e+02 2.150000e+02 2.827280e+02 0.000000e+00 2.000000e-03
2.450000e+02 2.150000e+02 2.831197e+02 0.000000e+00 2.000000e-03
2.550000e+02 2.150000e+02 2.831197e+02 0.000000e+00 2.000000e-03
2.650000e+02 2.150000e+02 2.827280e+02 0.000000e+00 2.000000e-03
2.750000e+02 2.150000e+02 2.819493e+02 0.000000e+00 2.000000e-03
2.850000e+02 2.150000e+02 2.807928e+02 0.000000e+00 2.000000e-03
2.950000e+02 2.150000e+02 2.792724e+02 0.000000e+00 2.000000e-03
3.050000e+02 2.150000e+02 2.774056e+02 0.000000e+00 2.000000e-03
3.150000e+02 2.150000e+02 2.752143e+02 0.000000e+00 2.000000e-03
3.250000e+02 2.150000e+02 2.727234e+02 0.000000e+00 2.000000e-03
3.350000e+02 2.150000e+02 2.699606e+02 0.000000e+00 2.000000e-03
3.450000e+02 2.150000e+02 2.669565e+02 0.000000e+00 2.000000e-03
3.550000e+02 2.150000e+02 2.637429e+02 0.000000e+00 2.000000e-03
3.650000e+02 2.150000e+02 2.603533e+02 0.000000e+00 2.000000e-03
3.750000e+02 2.150000e+02 2.568218e+02 0.000000e+00 2.000000e-03
1.050000e+02 2.250000e+02 2.501165e+02 0.000000e+00 2.000000e-03
1.150000e+02 2.250000e+02 2.538671e+02 0.000000e+00 2.000000e-03
1.250000e+02 2.250000e+02 2.575429e+02 0.000000e+00 2.000000e-03
1.350000e+02 2.250000e+02 2.611099e+02 0.000000e+00 2.000000e-03
1.450000e+02 2.250000e+02 2.645335e+02 0.000000e+00 2.000000e-03
1.550000e+02 2.250000e+02 2.677793e+02 0.000000e+00 2.000000e-03
1.650000e+02 2.250000e+02 2.708137e+02 0.000000e+00 2.000000e-03
1.750000e+02 2.250000e+02 2.736041e+02 0.000000e+00 2.000000e-03
1.850000e+02 2.250000e+02 2.761201e+02 0.000000e+00 2.000000e-03
1.950000e+02 2.250000e+02 2.783334e+02 0.000000e+00 2.000000e-03
2.050000e+02 2.250000e+02 2.802189e+02 0.000000e+00 2.000000e-03
2.150000e+02 2.250000e+02 2.817546e+02 0.000000e+00 2.000000e-03
2.250000e+02 2.250000e+02 2.829227e+02 0.000000e+00 2.000000e-03
2.350000e+02 2.250000e+02 2.837092e+02 0.000000e+00 2.000000e-03
2.450000e+02 2.250000e+02 2.841048e+02 0.000000e+00 2.000000e-03
2.550000e+02 2.250000e+02 2.841048e+02 0.000000e+00 2.000000e-03
2.650000e+02 2.250000e+02 2.837092e+02 0.000000e+00 2.000000e-03
2.750000e+02 2.250000e+02 2.829227e+02 0.000000e+00 2.000000e-03
2.850000e+02 2.250000e+02 2.817546e+02 0.000000e+00 2.000000e-03
2.950000e+02 2.250000e+02 2.802189e+02 0.000000e+00 2.000000e-03
3.050000e+02 2.250000e+02 2.783334e+02 0.000000e+00 2.000000e-03
3.150000e+02 2.250000e+02 2.761201e+02 0.000000e+00 2.000000e-03
3.250000e+02 2.250000e+02 2.736041e+02 0.000000e+00 2.000000e-03
3.350000e+02 2.250000e+02 2.708137e+02 0.000000e+00 2.000000e-03
3.450000e+02 2.250000e+02 2.677793e+02 0.000000e+00 2.000000e-03
3.550000e+02 2.250000e+02 2.645335e+02 0.000000e+00 2.000000e-03
3.650000e+02 2.250000e+02 2.611099e+02 0.000000e+00 2.000000e-03
3.750000e+02 2.250000e+02 2.575429e+02 0.000000e+00 2.000000e-03
1.050000e+02 2.350000e+02 2.505080e+02 0.000000e+00 2.000000e-03
1.150000e+02 2.350000e+02 2.542811e+02 0.000000e+00 2.000000e-03
1.250000e+02 2.350000e+02 2.579790e+02 0.000000e+00 2.000000e-03
1.350000e+02 2.350000e+02 2.615674e+02 0.000000e+00 2.000000e-03
1.450000e+02 2.350000e+02 2.650117e+02 0.000000e+00 2.000000e-03
1.550000e+02 2.350000e+02 2.682770e+02 0.000000e+00 2.000000e-03
1.650000e+02 2.350000e+02 2.713296e+02 0.000000e+00 2.000000e-03
1.750000e+02 2.350000e+02 2.741368e+02 0.000000e+00 2.000000e-03
1.850000e+02 2.350000e+02 2.766679e+02 0.000000e+00 2.000000e-03
1.950000e+02 2.350000e+02 2.788945e+02 0.000000e+00 2.000000e-03
2.050000e+02 2.350000e+02 2.807913e+02 0.000000e+00 2.000000e-03
2.150000e+02 2.350000e+02 2.823363e+02 0.000000e+00 2.000000e-03
2.250000e+02 2.350000e+02 2.835114e+02 0.000000e+00 2.000000e-03
2.350000e+02 2.350000e+02 2.843026e+02 0.000000e+00 2.000000e-03
2.450000e+02 2.350000e+02 2.847007e+02 0.000000e+00 2.000000e-03
2.550000e+02 2.350000e+02 2.847007e+02 0.000000e+00 2.000000e-03
2.650000e+02 2.350000e+02 2.843026e+02 0.000000e+00 2.000000e-03
2.750000e+02 2.350000e+02 2.835114e+02 0.000000e+00 2.000000e-03
2.850000e+02 2.350000e+02 2.823363e+02 0.000000e+00 2.000000e-03
2.950000e+02 2.350000e+02 2.807913e+02 0.000000e+00 2.000000e-03
3.050000e+02 2.350000e+02 2.788945e+02 0.000000e+00 2.000000e-03
3.150000e+02 2.350000e+02 2.766679e+02 0.000000e+00 2.000000e-03
3.250000e+02 2.350000e+02 2.741368e+02 0.000000e+00 2.000000e-03
3.350000e+02 2.350000e+02 2.713296e+02 0.000000e+00 2.000000e-03
3.450000e+02 2.350000e+02 2.682770e+02 0.000000e+00 2.000000e-03
3.550000e+02 2.350000e+02 2.650117e+02 0.000000e+00 2.000000e-03
3.650000e+02 2.350000e+02 2.615674e+02 0.000000e+00 2.000000e-03
3.750000e+02 2.350000e+02 2.579791e+02 0.000000e+00 2.000000e-03
1.050000e+02 2.450000e+02 2.506390e+02 0.000000e+00 2.000000e-03
1.150000e+02 2.450000e+02 2.544197e+02 0.000000e+00 2.000000e-03
1.250000e+02 2.450000e+02 2.581250e+02 0.000000e+00 2.000000e-03
1.350000e+02 2.450000e+02 2.617206e+02 0.000000e+00 2.000000e-03
1.450000e+02 2.450000e+02 2.651717e+02 0.000000e+00 2.000000e-03
1.550000e+02 2.450000e+02 2.684435e+02 0.000000e+00 2.000000e-03
1.650000e+02 2.450000e+02 2.715022e+02 0.000000e+00 2.000000e-03
1.750000e+02 2.450000e+02 2.743151e+02 0.000000e+00 2.000000e-03
1.850000e+02 2.450000e+02 2.768512e+02 0.000000e+00 2.000000e-03
1.950000e+02 2.450000e+02 2.790823e+02 0.000000e+00 2.000000e-03
2.050000e+02 2.450000e+02 2.809829e+02 0.000000e+00 2.000000e-03
2.150000e+02 2.450000e+02 2.825310e+02 0.000000e+00 2.000000e-03
2.250000e+02 2.450000e+02 2.837084e+02 0.000000e+00 2.000000e-03
2.350000e+02 2.450000e+02 2.845013e+02 0.000000e+00 2.000000e-03
2.450000e+02 2.450000e+02 2.849001e+02 0.000000e+00 2.000000e-03
2.550000e+02 2.450000e+02 2.849001e+02 0.000000e+00 2.000000e-03
2.650000e+02 2.450000e+02 2.845013e+02 0.000000e+00 2.000000e-03
2.750000e+02 2.450000e+02 2.837084e+02 0.000000e+00 2.000000e-03
2.850000e+02 2.450000e+02 2.825310e+02 0.000000e+00 2.000000e-03
2.950000e+02 2.450000e+02 2.809829e+02 0.000000e+00 2.000000e-03
3.050000e+02 2.450000e+02 2.790823e+02 0.000000e+00 2.000000e-03
3.150000e+02 2.450000e+02 2.768512e+02 0.000000e+00 2.000000e-03
3.250000e+02 2.450000e+02 2.743151e+02 0.000000e+00 2.000000e-03
3.350000e+02 2.450000e+02 2.715022e+02 0.000000e+00 2.000000e-03
3.450000e+02 2.450000e+02 2.684435e+02 0.000000e+00 2.000000e-03
3.550000e+02 2.450000e+02 2.651717e+02 0.000000e+00 2.000000e-03
3.650000e+02 2.450000e+02 2.617206e+02 0.000000e+00 2.000000e-03
3.750000e+02 2.450000e+02 2.581250e+02 0.000000e+00 2.000000e-03
1.050000e+02 2.550000e+02 2.505080e+02 0.000000e+00 2.000000e-03
1.150000e+02 2.550000e+02 2.542811e+02 0.000000e+00 2.000000e-03
1.250000e+02 2.550000e+02 2.579790e+02 0.000000e+00 2.000000e-03
1.350000e+02 2.550000e+02 2.615674e+02 0.000000e+00 2.000000e-03
1.450000e+02 2.550000e+02 2.650117e+02 0.000000e+00 2.000000e-03
1.550000e+02 2.550000e+02 2.682770e+02 0.000000e+00 2.000000e-03
1.650000e+02 2.550000e+02 2.713296e+02 0.000000e+00 2.000000e-03
1.750000e+02 2.550000e+02 2.741368e+02 0.000000e+00 2.000000e-03
1.850000e+02 2.550000e+02 2.766679e+02 0.000000e+00 2.000000e-03
1.950000e+02 2.550000e+02 2.788945e+02 0.000000e+00 2.000000e-03
2.050000e+02 2.550000e+02 2.807913e+02 0.000000e+00 2.000000e-03
2.150000e+02 2.550000e+02 2.823363e+02 0.000000e+00 2.000000e-03
2.250000e+02 2.550000e+02 2.835114e+02 0.000000e+00 2.000000e-03
2.350000e+02 2.550000e+02 2.843026e+02 0.000000e+00 2.000000e-03
2.450000e+02 2.550000e+02 2.847007e+02 0.000000e+00 2.000000e-03
2.550000e+02 2.550000e+02 2.847007e+02 0.000000e+00 2.000000e-03
2.650000e+02 2.550000e+02 2.843026e+02 0.000000e+00 2.000000e-03
2.750000e+02 2.550000e+02 2.835114e+02 0.000000e+00 2.000000e-03
2.850000e+02 2.550000e+02 2.823363e+02 0.000000e+00 2.000000e-03
2.950000e+02 2.550000e+02 2.807913e+02 0.000000e+00 2.000000e-03
3.050000e+02 2.550000e+02 2.788945e+02 0.000000e+00 2.000000e-03
3.150000e+02 2.550000e+02 2.766679e+02 0.000000e+00 2.000000e-03
3.250000e+02 2.550000e+02 2.741368e+02 0.000000e+00 2.000000e-03
3.350000e+02 2.550000e+02 2.713296e+02 0.000000e+00 2.000000e-03
3.450000e+02 2.550000e+02 2.682770e+02 0.000000e+00 2.000000e-03
3.550000e+02 2.550000e+02 2.650117e+02 0.000000e+00 2.000000e-03
3.650000e+02 2.550000e+02 2.615674e+02 0.000000e+00 2.000000e-03
3.750000e+02 2.550000e+02 2.579791e+02 0.000000e+00 2.000000e-03
1.050000e+02 2.650000e+02 2.501165e+02 0.000000e+00 2.000000e-03
1.150000e+02 2.650000e+02 2.538671e+02 0.000000e+00 2.000000e-03
1.250000e+02 2.650000e+02 2.575429e+02 0.000000e+00 2.000000e-03
1.350000e+02 2.650000e+02 2.611099e+02 0.000000e+00 2.000000e-03
1.450000e+02 2.650000e+02 2.645335e+02 0.000000e+00 2.000000e-03
1.550000e+02 2.650000e+02 2.677793e+02 0.000000e+00 2.000000e-03
1.650000e+02 2.650000e+02 2.708137e+02 0.000000e+00 2.000000e-03
1.750000e+02 2.650000e+02 2.736041e+02 0.000000e+00 2.000000e-03
1.850000e+02 2.650000e+02 2.761201e+02 0.000000e+00 2.000000e-03
1.950000e+02 2.650000e+02 2.783334e+02 0.000000e+00 2.000000e-03
2.050000e+02 2.650000e+02 2.802189e+02 0.000000e+00 2.000000e-03
2.150000e+02 2.650000e+02 2.817546e+02 0.000000e+00 2.000000e-03
2.250000e+02 2.650000e+02 2.829227e+02 0.000000e+00 2.000000e-03
2.350000e+02 2.650000e+02 2.837092e+02 0.000000e+00 2.000000e-03
2.450000e+02 2.650000e+02 2.841048e+02 0.000000e+00 2.000000e-03
2.550000e+02 2.650000e+02 2.841048e+02 0.000000e+00 2.000000e-03
2.650000e+02 2.650000e+02 2.837092e+02 0.000000e+00 2.000000e-03
2.750000e+02 2.650000e+02 2.829227e+02 0.000000e+00 2.000000e-03
2.850000e+02 2.650000e+02 2.817546e+02 0.000000e+00 2.000000e-03
2.950000e+02 2.650000e+02 2.802189e+02 0.000000e+00 2.000000e-03
3.050000e+02 2.650000e+02 2.783334e+02 0.000000e+00 2.000000e-03
3.150000e+02 2.650000e+02 2.761201e+02 0.000000e+00 2.000000e-03
3.250000e+02 2.650000e+02 2.736041e+02 0.000000e+00 2.000000e-03
3.350000e+02 2.650000e+02 2.708137e+02 0.000000e+00 2.000000e-03
3.450000e+02 2.650000e+02 2.677793e+02 0.000000e+00 2.000000e-03
3.550000e+02 2.650000e+02 2.645335e+02 0.000000e+00 2.000000e-03
3.650000e+02 2.650000e+02 2.611099e+02 0.000000e+00 2.000000e-03
3.750000e+02 2.650000e+02 2.575429e+02 0.000000e+00 2.000000e-03
1.050000e+02 2.750000e+02 2.494692e+02 0.000000e+00 2.000000e-03
1.150000e+02 2.750000e+02 2.531825e+02 0.000000e+00 2.000000e-03
1.250000e+02 2.750000e+02 2.568218e+02 0.000000e+00 2.000000e-03
1.350000e+02 2.750000e+02 2.603533e+02 0.000000e+00 2.000000e-03
1.450000e+02 2.750000e+02 2.637429e+02 0.000000e+00 2.000000e-03
1.550000e+02 2.750000e+02 2.669565e+02 0.000000e+00 2.000000e-03
1.650000e+02 2.750000e+02 2.699606e+02 0.000000e+00 2.000000e-03
1.750000e+02 2.750000e+02 2.727234e+02 0.000000e+00 2.000000e-03
1.850000e+02 2.750000e+02 2.752143e+02 0.000000e+00 2.000000e-03
1.950000e+02 2.750000e+02 2.774056e+02 0.000000e+00 2.000000e-03
2.050000e+02 2.750000e+02 2.792724e+02 0.000000e+00 2.000000e-03
2.150000e+02 2.750000e+02 2.807928e+02 0.000000e+00 2.000000e-03
2.250000e+02 2.750000e+02 2.819493e+02 0.000000e+00 2.000000e-03
2.350000e+02 2.750000e+02 2.827280e+02 0.000000e+00 2.000000e-03
2.450000e+02 2.750000e+02 2.831197e+02 0.000000e+00 2.000000e-03
2.550000e+02 2.750000e+02 2.831197e+02 0.000000e+00 2.000000e-03
2.650000e+02 2.750000e+02 2.827280e+02 0.000000e+00 2.000000e-03
2.750000e+02 2.750000e+02 2.819493e+02 0.000000e+00 2.000000e-03
2.850000e+02 2.750000e+02 2.807928e+02 0.000000e+00 2.000000e-03
2.950000e+02 2.750000e+02 2.792724e+02 0.000000e+00 2.000000e-03
3.050000e+02 2.750000e+02 2.774056e+02 0.000000e+00 2.000000e-03
3.150000e+02 2.750000e+02 2.752143e+02 0.000000e+00 2.000000e-03
3.250000e+02 2.750000e+02 2.727234e+02 0.000000e+00 2.000000e-03
3.350000e+02 2.750000e+02 2.699606e+02 0.000000e+00 2.000000e-03
3.450000e+02 2.750000e+02 2.669565e+02 0.000000e+00 2.000000e-03
3.550000e+02 2.750000e+02 2.637429e+02 0.000000e+00 2.000000e-03
3.650000e+02 2.750000e+02 2.603533e+02 0.000000e+00 2.000000e-03
3.750000e+02 2.750000e+02 2.568218e+02 0.000000e+00 2.000000e-03
1.050000e+02 2.850000e+02 2.485739e+02 0.000000e+00 2.000000e-03
1.150000e+02 2.850000e+02 2.522356e+02 0.000000e+00 2.000000e-03
1.250000e+02 2.850000e+02 2.558243e+02 0.000000e+00 2.000000e-03
1.350000e+02 2.850000e+02 2.593068e+02 0.000000e+00 2.000000e-03
1.450000e+02 2.850000e+02 2.626493e+02 0.000000e+00 2.000000e-03
1.550000e+02 2.850000e+02 2.658182e+02 0.000000e+00 2.000000e-03
1.650000e+02 2.850000e+02 2.687807e+02 0.000000e+00 2.000000e-03
1.750000e+02 2.850000e+02 2.715050e+02 0.000000e+00 2.000000e-03
1.850000e+02 2.850000e+02 2.739614e+02 0.000000e+00 2.000000e-03
1.950000e+02 2.850000e+02 2.761223e+02 0.000000e+00 2.000000e-03
2.050000e+02 2.850000e+02 2.779631e+02 0.000000e+00 2.000000e-03
2.150000e+02 2.850000e+02 2.794624e+02 0.000000e+00 2.000000e-03
2.250000e+02 2.850000e+02 2.806028e+02 0.000000e+00 2.000000e-03
2.350000e+02 2.850000e+02 2.813707e+02 0.000000e+00 2.000000e-03
2.450000e+02 2.850000e+02 2.817570e+02 0.000000e+00 2.000000e-03
2.550000e+02 2.850000e+02 2.817570e+02 0.000000e+00 2.000000e-03
2.650000e+02 2.850000e+02 2.813707e+02 0.000000e+00 2.000000e-03
2.750000e+02 2.850000e+02 2.806028e+02 0.000000e+00 2.000000e-03
2.850000e+02 2.850000e+02 2.794624e+02 0.000000e+00 2.000000e-03
2.950000e+02 2.850000e+02 2.779631e+02 0.000000e+00 2.000000e-03
3.050000e+02 2.850000e+02 2.761223e+02 0.000000e+00 2.000000e-03
3.150000e+02 2.850000e+02 2.739614e+02 0.000000e+00 2.000000e-03
3.250000e+02 2.850000e+02 2.715050e+02 0.000000e+00 2.000000e-03
3.350000e+02 2.850000e+02 2.687807e+02 0.000000e+00 2.000000e-03
3.450000e+02 2.850000e+02 2.658182e+02 0.000000e+00 2.000000e-03
3.550000e+02 2.850000e+02 2.626493e+02 0.000000e+00 2.000000e-03
3.650000e+02 2.850000e+02 2.593068e+02 0.000000e+00 2.000000e-03
3.750000e+02 2.850000e+02 2.558243e+02 0.000000e+00 2.000000e-03
1.050000e+02 2.950000e+02 2.474409e+02 0.000000e+00 2.000000e-03
1.150000e+02 2.950000e+02 2.510374e+02 0.000000e+00 2.000000e-03
1.250000e+02 2.950000e+02 2.545621e+02 0.000000e+00 2.000000e-03
1.350000e+02 2.950000e+02 2.579825e+02 0.000000e+00 2.000000e-03
1.450000e+02 2.950000e+02 2.612655e+02 0.000000e+00 2.000000e-03
1.550000e+02 2.950000e+02 2.643779e+02 0.000000e+00 2.000000e-03
1.650000e+02 2.950000e+02 2.672876e+02 0.000000e+00 2.000000e-03
1.750000e+02 2.950000e+02 2.699634e+02 0.000000e+00 2.000000e-03
1.850000e+02 2.950000e+02 2.723760e+02 0.000000e+00 2.000000e-03
1.950000e+02 2.950000e+02 2.744984e+02 0.000000e+00 2.000000e-03
2.050000e+02 2.950000e+02 2.763063e+02 0.000000e+00 2.000000e-03
2.150000e+02 2.950000e+02 2.777790e+02 0.000000e+00 2.000000e-03
2.250000e+02 2.950000e+02 2.788990e+02 0.000000e+00 2.000000e-03
2.350000e+02 2.950000e+02 2.796533e+02 0.000000e+00 2.000000e-03
2.450000e+02 2.950000e+02 2.800326e+02 0.000000e+00 2.000000e-03
2.550000e+02 2.950000e+02 2.800326e+02 0.000000e+00 2.000000e-03
2.650000e+02 2.950000e+02 2.796533e+02 0.000000e+00 2.000000e-03
2.750000e+02 2.950000e+02 2.788990e+02 0.000000e+00 2.000000e-03
2.850000e+02 2.950000e+02 2.777790e+02 0.000000e+00 2.000000e-03
2.950000e+02 2.950000e+02 2.763063e+02 0.000000e+00 2.000000e-03
3.050000e+02 2.950000e+02 2.744984e+02 0.000000e+00 2.000000e-03
3.150000e+02 2.950000e+02 2.723760e+02 0.000000e+00 2.000000e-03
3.250000e+02 2.950000e+02 2.699634e+02 0.000000e+00 2.000000e-03
3.350000e+02 2.950000e+02 2.672876e+02 0.000000e+00 2.000000e-03
3.450000e+02 2.950000e+02 2.643779e+02 0.000000e+00 2.000000e-03
3.550000e+02 2.950000e+02 2.612655e+02 0.000000e+00 2.000000e-03
3.650000e+02 2.950000e+02 2.579825e+02 0.000000e+00 2.000000e-03
3.750000e+02 2.950000e+02 2.545621e+02 0.000000e+00 2.000000e-03
1.050000e+02 3.050000e+02 2.460835e+02 0.000000e+00 2.000000e-03
1.150000e+02 3.050000e+02 2.496018e+02 0.000000e+00 2.000000e-03
1.250000e+02 3.050000e+02 2.530500e+02 0.000000e+00 2.000000e-03
1.350000e+02 3.050000e+02 2.563960e+02 0.000000e+00 2.000000e-03
1.450000e+02 3.050000e+02 2.596076e+02 0.000000e+00 2.000000e-03
1.550000e+02 3.050000e+02 2.626524e+02 0.000000e+00 2.000000e-03
1.650000e+02 3.050000e+02 2.654988e+02 0.000000e+00 2.000000e-03
1.750000e+02 3.050000e+02 2.681164e+02 0.000000e+00 2.000000e-03
1.850000e+02 3.050000e+02 2.704766e+02 0.000000e+00 2.000000e-03
1.950000e+02 3.050000e+02 2.725528e+02 0.000000e+00 2.000000e-03
2.050000e+02 3.050000e+02 2.743215e+02 0.000000e+00 2.000000e-03
2.150000e+02 3.050000e+02 2.757621e+02 0.000000e+00 2.000000e-03
2.250000e+02 3.050000e+02 2.768578e+02 0.000000e+00 2.000000e-03
2.350000e+02 3.050000e+02 2.775956e+02 0.000000e+00 2.000000e-03
2.450000e+02 3.050000e+02 2.779668e+02 0.000000e+00 2.000000e-03
2.550000e+02 3.050000e+02 2.779668e+02 0.000000e+00 2.000000e-03
2.650000e+02 3.050000e+02 2.775956e+02 0.000000e+00 2.000000e-03
2.750000e+02 3.050000e+02 2.768578e+02 0.000000e+00 2.000000e-03
2.850000e+02 3.050000e+02 2.757621e+02 0.000000e+00 2.000000e-03
2.950000e+02 3.050000e+02 2.743215e+02 0.000000e+00 2.000000e-03
3.050000e+02 3.050000e+02 2.725528e+02 0.000000e+00 2.000000e-03
3.150000e+02 3.050000e+02 2.704766e+02 0.000000e+00 2.000000e-03
3.250000e+02 3.050000e+02 2.681164e+02 0.000000e+00 2.000000e-03
3.350000e+02 3.050000e+02 2.654988e+02 0.000000e+00 2.000000e-03
3.450000e+02 3.050000e+02 2.626524e+02 0.000000e+00 2.000000e-03
3.550000e+02 3.050000e+02 2.596076e+02 0.000000e+00 2.000000e-03
3.650000e+02 3.050000e+02 2.563960e+02 0.000000e+00 2.000000e-03
3.750000e+02 3.050000e+02 2.530500e+02 0.000000e+00 2.000000e-03
1.050000e+02 3.150000e+02 2.445174e+02 0.000000e+00 2.000000e-03
1.150000e+02 3.150000e+02 2.479454e+02 0.000000e+00 2.000000e-03
1.250000e+02 3.150000e+02 2.513052e+02 0.000000e+00 2.000000e-03
1.350000e+02 3.150000e+02 2.545655e+02 0.000000e+00 2.000000e-03
1.450000e+02 3.150000e+02 2.576947e+02 0.000000e+00 2.000000e-03
1.550000e+02 3.150000e+02 2.606614e+02 0.000000e+00 2.000000e-03
1.650000e+02 3.150000e+02 2.634348e+02 0.000000e+00 2.000000e-03
1.750000e+02 3.150000e+02 2.659854e+02 0.000000e+00 2.000000e-03
1.850000e+02 3.150000e+02 2.682850e+02 0.000000e+00 2.000000e-03
1.950000e+02 3.150000e+02 2.703080e+02 0.000000e+00 2.000000e-03
2.050000e+02 3.150000e+02 2.720313e+02 0.000000e+00 2.000000e-03
2.150000e+02 3.150000e+02 2.734350e+02 0.000000e+00 2.000000e-03
2.250000e+02 3.150000e+02 2.745027e+02 0.000000e+00 2.000000e-03
2.350000e+02 3.150000e+02 2.752216e+02 0.000000e+00 2.000000e-03
2.450000e+02 3.150000e+02 2.755831e+02 0.000000e+00 2.000000e-03
2.550000e+02 3.150000e+02 2.755831e+02 0.000000e+00 2.000000e-03
2.650000e+02 3.150000e+02 2.752216e+02 0.000000e+00 2.000000e-03
2.750000e+02 3.150000e+02 2.745027e+02 0.000000e+00 2.000000e-03
2.850000e+02 3.150000e+02 2.734350e+02 0.000000e+00 2.000000e-03
2.950000e+02 3.150000e+02 2.720313e+02 0.000000e+00 2.000000e-03
3.050000e+02 3.150000e+02 2.703080e+02 0.000000e+00 2.000000e-03
3.150000e+02 3.150000e+02 2.682850e+02 0.000000e+00 2.000000e-03
3.250000e+02 3.150000e+02 2.659854e+02 0.000000e+00 2.000000e-03
3.350000e+02 3.150000e+02 2.634348e+02 0.000000e+00 2.000000e-03
3.450000e+02 3.150000e+02 2.606614e+02 0.000000e+00 2.000000e-03
3.550000e+02 3.150000e+02 2.576947e+02 0.000000e+00 2.000000e-03
3.650000e+02 3.150000e+02 2.545655e+02 0.000000e+00 2.000000e-03
3.750000e+02 3.150000e+02 2.513052e+02 0.000000e+00 2.000000e-03
1.050000e+02 3.250000e+02 2.427601e+02 0.000000e+00 2.000000e-03
1.150000e+02 3.250000e+02 2.460869e+02 0.000000e+00 2.000000e-03
1.250000e+02 3.250000e+02 2.493475e+02 0.000000e+00 2.000000e-03
1.350000e+02 3.250000e+02 2.525115e+02 0.000000e+00 2.000000e-03
1.450000e+02 3.250000e+02 2.555483e+02 0.000000e+00 2.000000e-03
1.550000e+02 3.250000e+02 2.584275e+02 0.000000e+00 2.000000e-03
1.650000e+02 3.250000e+02 2.611190e+02 0.000000e+00 2.000000e-03
1.750000e+02 3.250000e+02 2.635942e+02 0.000000e+00 2.000000e-03
1.850000e+02 3.250000e+02 2.658260e+02 0.000000e+00 2.000000e-03
1.950000e+02 3.250000e+02 2.677892e+02 0.000000e+00 2.000000e-03
2.050000e+02 3.250000e+02 2.694617e+02 0.000000e+00 2.000000e-03
2.150000e+02 3.250000e+02 2.708240e+02 0.000000e+00 2.000000e-03
2.250000e+02 3.250000e+02 2.718600e+02 0.000000e+00 2.000000e-03
2.350000e+02 3.250000e+02 2.725577e+02 0.000000e+00 2.000000e-03
2.450000e+02 3.250000e+02 2.729086e+02 0.000000e+00 2.000000e-03
2.550000e+02 3.250000e+02 2.729086e+02 0.000000e+00 2.000000e-03
2.650000e+02 3.250000e+02 2.725577e+02 0.000000e+00 2.000000e-03
2.750000e+02 3.250000e+02 2.718600e+02 0.000000e+00 2.000000e-03
2.850000e+02 3.250000e+02 2.708240e+02 0.000000e+00 2.000000e-03
2.950000e+02 3.250000e+02 2.694617e+02 0.000000e+00 2.000000e-03
3.050000e+02 3.250000e+02 2.677892e+02 0.000000e+00 2.000000e-03
3.150000e+02 3.250000e+02 2.658260e+02 0.000000e+00 2.000000e-03
3.250000e+02 3.250000e+02 2.635942e+02 0.000000e+00 2.000000e-03
3.350000e+02 3.250000e+02 2.611190e+02 0.000000e+00 2.000000e-03
3.450000e+02 3.250000e+02 2.584275e+02 0.000000e+00 2.000000e-03
3.550000e+02 3.250000e+02 2.555483e+02 0.000000e+00 2.000000e-03
3.650000e+02 3.250000e+02 2.525115e+02 0.000000e+00 2.000000e-03
3.750000e+02 3.250000e+02 2.493475e+02 0.000000e+00 2.000000e-03
1.050000e+02 3.350000e+02 2.408312e+02 0.000000e+00 2.000000e-03
1.150000e+02 3.350000e+02 2.440469e+02 0.000000e+00 2.000000e-03
1.250000e+02 3.350000e+02 2.471986e+02 0.000000e+00 2.000000e-03
1.350000e+02 3.350000e+02 2.502569e+02 0.000000e+00 2.000000e-03
1.450000e+02 3.350000e+02 2.531923e+02 0.000000e+00 2.000000e-03
1.550000e+02 3.350000e+02 2.559753e+02 0.000000e+00 2.000000e-03
1.650000e+02 3.350000e+02 2.585770e+02 0.000000e+00 2.000000e-03
1.750000e+02 3.350000e+02 2.609695e+02 0.000000e+00 2.000000e-03
1.850000e+02 3.350000e+02 2.631267e+02 0.000000e+00 2.000000e-03
1.950000e+02 3.350000e+02 2.650244e+02 0.000000e+00 2.000000e-03
2.050000e+02 3.350000e+02 2.666411e+02 0.000000e+00 2.000000e-03
2.150000e+02 3.350000e+02 2.679578e+02 0.000000e+00 2.000000e-03
2.250000e+02 3.350000e+02 2.689593e+02 0.000000e+00 2.000000e-03
2.350000e+02 3.350000e+02 2.696337e+02 0.000000e+00 2.000000e-03
2.450000e+02 3.350000e+02 2.699729e+02 0.000000e+00 2.000000e-03
2.550000e+02 3.350000e+02 2.699729e+02 0.000000e+00 2.000000e-03
2.650000e+02 3.350000e+02 2.696337e+02 0.000000e+00 2.000000e-03
2.750000e+02 3.350000e+02 2.689593e+02 0.000000e+00 2.000000e-03
2.850000e+02 3.350000e+02 2.679578e+02 0.000000e+00 2.000000e-03
2.950000e+02 3.350000e+02 2.666411e+02 0.000000e+00 2.000000e-03
3.050000e+02 3.350000e+02 2.650244e+02 0.000000e+00 2.000000e-03
3.150000e+02 3.350000e+02 2.631267e+02 0.000000e+00 2.000000e-03
3.250000e+02 3.350000e+02 2.609695e+02 0.000000e+00 2.000000e-03
3.350000e+02 3.350000e+02 2.585770e+02 0.000000e+00 2.000000e-03
3.450000e+02 3.350000e+02 2.559753e+02 0.000000e+00 2.000000e-03
3.550000e+02 3.350000e+02 2.531923e+02 0.000000e+00 2.000000e-03
3.650000e+02 3.350000e+02 2.502569e+02 0.000000e+00 2.000000e-03
3.750000e+02 3.350000e+02 2.471986e+02 0.000000e+00 2.000000e-03
1.050000e+02 3.450000e+02 2.387514e+02 0.000000e+00 2.000000e-03
1.150000e+02 3.450000e+02 2.418474e+02 0.000000e+00 2.000000e-03
1.250000e+02 3.450000e+02 2.448817e+02 0.000000e+00 2.000000e-03
1.350000e+02 3.450000e+02 2.478261e+02 0.000000e+00 2.000000e-03
1.450000e+02 3.450000e+02 2.506521e+02 0.000000e+00 2.000000e-03
1.550000e+02 3.450000e+02 2.533315e+02 0.000000e+00 2.000000e-03
1.650000e+02 3.450000e+02 2.558362e+02 0.000000e+00 2.000000e-03
1.750000e+02 3.450000e+02 2.581396e+02 0.000000e+00 2.000000e-03
1.850000e+02 3.450000e+02 2.602165e+02 0.000000e+00 2.000000e-03
1.950000e+02 3.450000e+02 2.620435e+02 0.000000e+00 2.000000e-03
2.050000e+02 3.450000e+02 2.635999e+02 0.000000e+00 2.000000e-03
2.150000e+02 3.450000e+02 2.648676e+02 0.000000e+00 2.000000e-03
2.250000e+02 3.450000e+02 2.658318e+02 0.000000e+00 2.000000e-03
2.350000e+02 3.450000e+02 2.664810e+02 0.000000e+00 2.000000e-03
2.450000e+02 3.450000e+02 2.668076e+02 0.000000e+00 2.000000e-03
2.550000e+02 3.450000e+02 2.668076e+02 0.000000e+00 2.000000e-03
2.650000e+02 3.450000e+02 2.664810e+02 0.000000e+00 2.000000e-03
2.750000e+02 3.450000e+02 2.658318e+02 0.000000e+00 2.000000e-03
2.850000e+02 3.450000e+02 2.648676e+02 0.000000e+00 2.000000e-03
2.950000e+02 3.450000e+02 2.635999e+02 0.000000e+00 2.000000e-03
3.050000e+02 3.450000e+02 2.620435e+02 0.000000e+00 2.000000e-03
3.150000e+02 3.450000e+02 2.602165e+02 0.000000e+00 2.000000e-03
3.250000e+02 3.450000e+02 2.581396e+02 0.000000e+00 2.000000e-03
3.350000e+02 3.450000e+02 2.558362e+02 0.000000e+00 2.000000e-03
3.450000e+02 3.450000e+02 2.533315e+02 0.000000e+00 2.000000e-03
3.550000e+02 3.450000e+02 2.506521e+02 0.000000e+00 2.000000e-03
3.650000e+02 3.450000e+02 2.478261e+02 0.000000e+00 2.000000e-03
3.750000e+02 3.450000e+02 2.448817e+02 0.000000e+00 2.000000e-03
1.050000e+02 3.550000e+02 2.365428e+02 0.000000e+00 2.000000e-03
1.150000e+02 3.550000e+02 2.395115e+02 0.000000e+00 2.000000e-03
1.250000e+02 3.550000e+02 2.424211e+02 0.000000e+00 2.000000e-03
1.350000e+02 3.550000e+02 2.452446e+02 0.000000e+00 2.000000e-03
1.450000e+02 3.550000e+02 2.479545e+02 0.000000e+00 2.000000e-03
1.550000e+02 3.550000e+02 2.505237e+02 0.000000e+00 2.000000e-03
1.650000e+02 3.550000e+02 2.529255e+02 0.000000e+00 2.000000e-03
1.750000e+02 3.550000e+02 2.551343e+02 0.000000e+00 2.000000e-03
1.850000e+02 3.550000e+02 2.571258e+02 0.000000e+00 2.000000e-03
1.950000e+02 3.550000e+02 2.588778e+02 0.000000e+00 2.000000e-03
2.050000e+02 3.550000e+02 2.603702e+02 0.000000e+00 2.000000e-03
2.150000e+02 3.550000e+02 2.615858e+02 0.000000e+00 2.000000e-03
2.250000e+02 3.550000e+02 2.625104e+02 0.000000e+00 2.000000e-03
2.350000e+02 3.550000e+02 2.631330e+02 0.000000e+00 2.000000e-03
2.450000e+02 3.550000e+02 2.634461e+02 0.000000e+00 2.000000e-03
2.550000e+02 3.550000e+02 2.634461e+02 0.000000e+00 2.000000e-03
2.650000e+02 3.550000e+02 2.631330e+02 0.000000e+00 2.000000e-03
2.750000e+02 3.550000e+02 2.625104e+02 0.000000e+00 2.000000e-03
2.850000e+02 3.550000e+02 2.615858e+02 0.000000e+00 2.000000e-03
2.950000e+02 3.550000e+02 2.603702e+02 0.000000e+00 2.000000e-03
3.050000e+02 3.550000e+02 2.588778e+02 0.000000e+00 2.000000e-03
3.150000e+02 3.550000e+02 2.571258e+02 0.000000e+00 2.000000e-03
3.250000e+02 3.550000e+02 2.551343e+02 0.000000e+00 2.000000e-03
3.350000e+02 3.550000e+02 2.529255e+02 0.000000e+00 2.000000e-03
3.450000e+02 3.550000e+02 2.505237e+02 0.000000e+00 2.000000e-03
3.550000e+02 3.550000e+02 2.479545e+02 0.000000e+00 2.000000e-03
3.650000e+02 3.550000e+02 2.452446e+02 0.000000e+00 2.000000e-03
3.750000e+02 3.550000e+02 2.424211e+02 0.000000e+00 2.000000e-03
1.050000e+02 3.650000e+02 2.342278e+02 0.000000e+00 2.000000e-03
1.150000e+02 3.650000e+02 2.370632e+02 0.000000e+00 2.000000e-03
1.250000e+02 3.650000e+02 2.398421e+02 0.000000e+00 2.000000e-03
1.350000e+02 3.650000e+02 2.425387e+02 0.000000e+00 2.000000e-03
1.450000e+02 3.650000e+02 2.451270e+02 0.000000e+00 2.000000e-03
1.550000e+02 3.650000e+02 2.475808e+02 0.000000e+00 2.000000e-03
1.650000e+02 3.650000e+02 2.498748e+02 0.000000e+00 2.000000e-03
1.750000e+02 3.650000e+02 2.519843e+02 0.000000e+00 2.000000e-03
1.850000e+02 3.650000e+02 2.538864e+02 0.000000e+00 2.000000e-03
1.950000e+02 3.650000e+02 2.555596e+02 0.000000e+00 2.000000e-03
2.050000e+02 3.650000e+02 2.569850e+02 0.000000e+00 2.000000e-03
2.150000e+02 3.650000e+02 2.581461e+02 0.000000e+00 2.000000e-03
2.250000e+02 3.650000e+02 2.590291e+02 0.000000e+00 2.000000e-03
2.350000e+02 3.650000e+02 2.596237e+02 0.000000e+00 2.000000e-03
2.450000e+02 3.650000e+02 2.599228e+02 0.000000e+00 2.000000e-03
2.550000e+02 3.650000e+02 2.599228e+02 0.000000e+00 2.000000e-03
2.650000e+02 3.650000e+02 2.596237e+02 0.000000e+00 2.000000e-03
2.750000e+02 3.650000e+02 2.590291e+02 0.000000e+00 2.000000e-03
2.850000e+02 3.650000e+02 2.581461e+02 0.000000e+00 2.000000e-03
2.950000e+02 3.650000e+02 2.569850e+02 0.000000e+00 2.000000e-03
3.050000e+02 3.650000e+02 2.555596e+02 0.000000e+00 2.000000e-03
3.150000e+02 3.650000e+02 2.538864e+02 0.000000e+00 2.000000e-03
3.250000e+02 3.650000e+02 2.519843e+02 0.000000e+00 2.000000e-03
3.350000e+02 3.650000e+02 2.498748e+02 0.000000e+00 2.000000e-03
3.450000e+02 3.650000e+02 2.475808e+02 0.000000e+00 2.000000e-03
3.550000e+02 3.650000e+02 2.451270e+02 0.000000e+00 2.000000e-03
3.650000e+02 3.650000e+02 2.425387e+02 0.000000e+00 2.000000e-03
3.750000e+02 3.650000e+02 2.398421e+02 0.000000e+00 2.000000e-03
simpegPF/Dev/MAG/Intgrl_MAG_Inv_Driver.py
+16 -33
View File
@@ -14,7 +14,6 @@ wd = survey.std
ndata = survey.srcField.rxList[0].locs.shape[0]
beta_in = 1e+5
eps_p = 1e-4
eps_q = 1e-4
@@ -24,8 +23,8 @@ nC = len(actv)
# Create active map to go from reduce set to full
actvMap = Maps.InjectActiveCells(mesh, actv, -100)
# Creat reduced identity map
idenMap = Maps.IdentityMap(nP = nC)
# Create reduced identity map
idenMap = Maps.IdentityMap(nP=nC)
# Get magnetization vector for MOF
M_xyz = driver.magnetizationModel
@@ -34,12 +33,6 @@ M_xyz = driver.magnetizationModel
midx = int(mesh.nCx/2)
midy = int(mesh.nCy/2)
# 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)
#==============================================================================
#%% Plot obs data
PF.Magnetics.plot_obs_2D(rxLoc,d, 'Observed Data')
@@ -47,9 +40,7 @@ PF.Magnetics.plot_obs_2D(rxLoc,d, 'Observed Data')
prob = PF.Magnetics.Problem3D_Integral(mesh, mapping=idenMap, actInd=actv)
prob.solverOpts['accuracyTol'] = 1e-4
survey.pair(prob)
#survey.makeSyntheticData(data, std=0.01)
#survey.dobs=d
#survey.mtrue = model
# Write out the predicted
pred = prob.fields(driver.m0)
PF.Magnetics.writeUBCobs('Pred.dat', survey, pred)
@@ -68,24 +59,19 @@ plt.gca().set_aspect('equal', adjustable='box')
reg = Regularization.Simple(mesh, indActive=actv, mapping=idenMap)
reg.mref = driver.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()
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)
invProb = InvProblem.BaseInvProblem(dmis, reg, opt)
beta = Directives.BetaSchedule(coolingFactor=2, coolingRate=1)
#betaest = Directives.BetaEstimate_ByEig()
beta_init = Directives.BetaEstimate_ByEig()
target = Directives.TargetMisfit()
update_Jacobi = Directives.Update_lin_PreCond(onlyOnStart=True)
inv = Inversion.BaseInversion(invProb, directiveList=[beta,target])
inv = Inversion.BaseInversion(invProb, directiveList=[beta,target,beta_init,update_Jacobi])
# Run inversion
mrec = inv.run(driver.m0)
@@ -94,7 +80,6 @@ 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)
@@ -134,6 +119,7 @@ plt.gca().set_aspect('equal', adjustable='box')
phim = invProb.phi_m_last
phid = invProb.phi_d
# Set parameters for sparsity
reg = Regularization.Sparse(mesh, indActive = actv, mapping = idenMap)
reg.recModel = mrec
reg.mref = driver.mref
@@ -143,26 +129,22 @@ reg.eps_q = eps_q
reg.norms = driver.lpnorms
diagA = np.sum(prob.G**2.,axis=0) + beta_in*(reg.W.T*reg.W).diagonal()
PC = Utils.sdiag(diagA**-1.)
#reg.alpha_s = 1.
dmis = DataMisfit.l2_DataMisfit(survey)
dmis.Wd = wd
opt = Optimization.ProjectedGNCG(maxIter=20 ,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)
# Create inversion directives
beta = Directives.BetaSchedule(coolingFactor=1, coolingRate=1)
update_beta = Directives.Scale_Beta(tol = 0.05)
#betaest = Directives.BetaEstimate_ByEig()
target = Directives.TargetMisfit()
IRLS =Directives.Update_IRLS( phi_m_last = phim, phi_d_last = phid )
update_Jacobi = Directives.Update_lin_PreCond(onlyOnStart=False)
save_log = Directives.SaveOutputEveryIteration()
save_log.fileName = 'LogName_blabla'
inv = Inversion.BaseInversion(invProb, directiveList=[beta,IRLS,update_beta])
inv = Inversion.BaseInversion(invProb, directiveList=[beta,IRLS,update_beta,update_Jacobi,save_log])
m0 = mrec
@@ -171,6 +153,7 @@ mrec = inv.run(m0)
m_out = actvMap*mrec
# Write final model out.
Mesh.TensorMesh.writeModelUBC(mesh,'SimPEG_inv_l0l2.sus',m_out)
pred = prob.fields(mrec)
simpegPF/notebooks/SimPEG Tutorial - MAG Linear Problem.ipynb
+32 -29
View File
File diff suppressed because one or more lines are too long