![]() % Depths of the lower limits of the model blocksįprintf( fidw, ' %. % Number of nodes between adjacent electrodesįprintf( fidw, ' %g, ', R( 1 : end - 1)) % Position of the first electrode (m), number of blocks, number of modelįprintf( fidw, ' %g, %g, %g\n ', x1, ncol, nr) % "2" indicates user-defined block depths X1 = x( 1) % Position of first electrodeįidw = fopen(, 'wt ') įprintf( fidw,) Nr = length( R) % Number of model resistivities Nlin = size( MM, 1) % Number of model lines Ncol = size( MM, 2) % Number of model columns Nspac = x( 3)- x( 1) % Unit electrode spacing Nelec = ( size( x, 2) 1)/ 2 % Number of electrodes MMnew = MMnew - 1 % Labels start from "0" R = unique( RR) % Resistivities of non-empty bins In ERT-1 we removed data with standard deviation above 5, and in ERT-2 data with a standard deviation above 8. Before processing, we integrated topography data and manually removed the bad data points in the dipoledipole data sets. nfac minM)) % Compute bin resistivities (not needed) The resulting data was processed using Res2DInv inversion software. % replace block resistivities by bin numbers (1 through nr0) % Sort resistivity values into nr0 logarithmically equidistant bins and % Read model data from inversion file (.INV) number of differentĪry = 3 % Array type (1=Wenner alpha, 2=pole-pole, 3=dipole-dipole, 4=Wenner beta, 5=Wenner Gamma, 6=inline pole-dipole, 7=Wenner-Schlumberger, 8=equatorial dipole-dipole ![]() % Model discretization > Use extended model) Inversion must have been carried out with extended model (Inversion > % model blocks between adjacent electrodes) Inversion must have been carried out with refined model (i.e., two Reduces number of different model resistivities to a max.
0 Comments
Leave a Reply. |