diff --git a/simpegDCIP/Dev/DC3D_fwr_MtIsa.py b/simpegDCIP/Dev/DC3D_fwr_MtIsa.py index 9dd5905d..854a32ce 100644 --- a/simpegDCIP/Dev/DC3D_fwr_MtIsa.py +++ b/simpegDCIP/Dev/DC3D_fwr_MtIsa.py @@ -38,7 +38,7 @@ import scipy.interpolate as interpolation import os #home_dir = 'C:\\Users\\dominiquef.MIRAGEOSCIENCE\\ownCloud\\Research\\Modelling\\Synthetic\\Two_Sphere' -home_dir ='C:\Users\dominiquef.MIRAGEOSCIENCE\Documents\GIT\SimPEG\simpegdc\simpegDCIP\Dev' +home_dir ='C:\Users\dominiquef.MIRAGEOSCIENCE\ownCloud\Research\MtIsa\Modeling' dsep = '\\' #from scipy.linalg import solve_banded @@ -48,11 +48,13 @@ mesh = Utils.meshutils.readUBCTensorMesh(home_dir + '\MtIsa_20m.msh') #mesh = Utils.meshutils.readUBCTensorMesh(home_dir + '\Mesh_50m.msh') # Load model -model = Utils.meshutils.readUBCTensorModel(home_dir + '\MtIsa_3D.con',mesh) +model = Utils.meshutils.readUBCTensorModel(home_dir + '\MtIsa_20m.con',mesh) #model = Utils.meshutils.readUBCTensorModel(home_dir + '\Synthetic.con',mesh) #model = Utils.meshutils.readUBCTensorModel(home_dir + '\Lalor_model_50m.con',mesh) #model = Utils.meshutils.readUBCTensorModel(home_dir + '\TwoSpheres.con',mesh) +#model[model>1] = 0.08 + #model = model**0 * 1e-2 # Specify survey type stype = 'pdp' @@ -60,7 +62,7 @@ stype = 'pdp' # Survey parameters a = 100 b = 100 -n = 20 +n = 15 # Forward solver slvr = 'BiCGStab' #'LU' @@ -130,7 +132,7 @@ plt.sca(ax_prim) # Takes two points from ginput and create survey #if re.match(stype,'gradient'): -gin = [(400.,12200.), (1400.,12200.)] +gin = [(400.,12200.), (1800.,12200.)] #else: #gin = plt.ginput(2, timeout = 0) @@ -278,10 +280,10 @@ if not re.match(stype,'gradient'): #============================================================================== plt.figure() - axs = plt.subplot() + axs = plt.subplot(2,1,1) plt.xlim([0,nc*dx]) - plt.ylim([mesh2d.vectorNy[-1]-dl_len,mesh2d.vectorNy[-1]]) + plt.ylim([mesh2d.vectorNy[-1]-dl_len/2,mesh2d.vectorNy[-1]]) plt.gca().set_aspect('equal', adjustable='box') plt.pcolormesh(mesh2d.vectorNx,mesh2d.vectorNy,np.log10(m2D),alpha=0.5, cmap='gray')#axes = [mesh2d.vectorNx[0],mesh2d.vectorNx[-1],mesh2d.vectorNy[0],mesh2d.vectorNy[-1]]) @@ -345,11 +347,11 @@ if not re.match(stype,'gradient'): #%% #Load model minv = DC.readUBC_DC2DModel(inv_dir + dsep + 'dcinv2d.con') - plt.figure() - axs = plt.subplot() + #plt.figure() + axs = plt.subplot(2,1,2) plt.xlim([0,nc*dx]) - plt.ylim([mesh2d.vectorNy[-1]-dl_len,mesh2d.vectorNy[-1]]) + plt.ylim([mesh2d.vectorNy[-1]-dl_len/2,mesh2d.vectorNy[-1]]) plt.gca().set_aspect('equal', adjustable='box') minv = np.reshape(minv,(mesh2d.nCy,mesh2d.nCx)) @@ -387,4 +389,15 @@ elif re.match(stype,'gradient'): var = 'Gradient Array - a-spacing: ' + str(a) + ' m' plt.title(var) plt.colorbar() - plt.contour(grid_x,grid_z,grid_rho, colors='k') \ No newline at end of file + plt.contour(grid_x,grid_z,grid_rho, colors='k') + +#%% Load tight model and plot +mesh = Utils.meshutils.readUBCTensorMesh(home_dir + '\MtIsa_5m.msh') + +# Load model +model = Utils.meshutils.readUBCTensorModel(home_dir + '\MtIsa_5m.con',mesh) +model = model.reshape((mesh.nCz,mesh.nCx)) +plt.figure() +plt.imshow(np.log10(model),extent = (125,375,0,75),origin='lower') +plt.colorbar(fraction=0.015) + \ No newline at end of file diff --git a/simpegDCIP/Dev/Inv2D/dcinv2d.log b/simpegDCIP/Dev/Inv2D/dcinv2d.log index 27da5405..5b30b7f7 100644 --- a/simpegDCIP/Dev/Inv2D/dcinv2d.log +++ b/simpegDCIP/Dev/Inv2D/dcinv2d.log @@ -14,7 +14,7 @@ Distributed by: Mira Geoscience Ltd. - DCINV2D started on: 1/12/2016 16:42:18 + DCINV2D started on: 1/12/2016 21:07:02 Reading input file: dcinv2d.inp @@ -93,7 +93,7 @@ model norm change = 0.00000E+00 norm comp Wx = 3.70226E-01 norm comp Wz = 3.77603E-01 -iter cpu time: 0:00:01.20 +iter cpu time: 0:00:01.81 Iteration 2 @@ -114,7 +114,7 @@ model norm change = 2.77567E+00 norm comp Wx = 1.61088E+00 norm comp Wz = 1.40814E+00 -iter cpu time: 0:00:00.78 +iter cpu time: 0:00:00.83 Iteration 3 @@ -135,7 +135,7 @@ model norm change = 1.12370E+00 norm comp Wx = 4.11507E+00 norm comp Wz = 3.02271E+00 -iter cpu time: 0:00:01.00 +iter cpu time: 0:00:00.64 Iteration 4 @@ -156,7 +156,7 @@ model norm change = 6.62554E-01 norm comp Wx = 8.32898E+00 norm comp Wz = 5.18098E+00 -iter cpu time: 0:00:00.86 +iter cpu time: 0:00:01.47 Iteration 5 @@ -177,7 +177,7 @@ model norm change = 4.40344E-01 norm comp Wx = 1.43993E+01 norm comp Wz = 7.86671E+00 -iter cpu time: 0:00:00.88 +iter cpu time: 0:00:01.41 Iteration 6 @@ -199,7 +199,7 @@ model norm change = 2.84090E-01 norm comp Wx = 2.11440E+01 norm comp Wz = 1.08731E+01 -iter cpu time: 0:00:00.81 +iter cpu time: 0:00:01.02 Iteration 7 @@ -221,7 +221,7 @@ model norm change = 1.89686E-01 norm comp Wx = 2.69488E+01 norm comp Wz = 1.38454E+01 -iter cpu time: 0:00:00.89 +iter cpu time: 0:00:01.00 Iteration 8 @@ -243,7 +243,7 @@ model norm change = 1.53314E-01 norm comp Wx = 3.20341E+01 norm comp Wz = 1.68852E+01 -iter cpu time: 0:00:00.91 +iter cpu time: 0:00:01.08 Iteration 9 @@ -264,7 +264,7 @@ model norm change = 1.48742E-01 norm comp Wx = 3.78831E+01 norm comp Wz = 2.05954E+01 -iter cpu time: 0:00:00.78 +iter cpu time: 0:00:00.62 Iteration 10 @@ -285,7 +285,7 @@ model norm change = 1.81682E-01 norm comp Wx = 4.67077E+01 norm comp Wz = 2.57538E+01 -iter cpu time: 0:00:00.66 +iter cpu time: 0:00:01.14 Iteration 11 @@ -306,7 +306,7 @@ model norm change = 2.42021E-01 norm comp Wx = 6.50107E+01 norm comp Wz = 3.14941E+01 -iter cpu time: 0:00:00.81 +iter cpu time: 0:00:01.18 Iteration 12 @@ -327,7 +327,7 @@ model norm change = 1.55454E-01 norm comp Wx = 7.90861E+01 norm comp Wz = 3.81319E+01 -iter cpu time: 0:00:00.72 +iter cpu time: 0:00:01.10 Iteration 13 @@ -349,7 +349,7 @@ model norm change = -6.98309E-03 norm comp Wx = 7.63269E+01 norm comp Wz = 3.53536E+01 -iter cpu time: 0:00:01.02 +iter cpu time: 0:00:00.88 Iteration 14 @@ -371,7 +371,7 @@ model norm change = -3.03571E-02 norm comp Wx = 7.21243E+01 norm comp Wz = 3.19461E+01 -iter cpu time: 0:00:01.21 +iter cpu time: 0:00:00.89 Target misfit achieved. Minimizing model norm. @@ -396,7 +396,7 @@ model norm change = -6.88525E-03 norm comp Wx = 6.90749E+01 norm comp Wz = 3.41159E+01 -iter cpu time: 0:00:00.76 +iter cpu time: 0:00:00.95 Exit at convergence. @@ -404,7 +404,7 @@ iter cpu time: 0:00:00.76 Iterations performed: 15 -total cpu time: 0:00:13.54 +total cpu time: 0:00:16.25 - DCINV2D ended on: 1/12/2016 16:42:31 + DCINV2D ended on: 1/12/2016 21:07:18