fillsig.m

生成三维模型异常体填充文件

      stg=1:3;            % List of stages to execute
    wordy=2;              % Messages: 0-silence, 1-basic messages, 2-detailed messages
     mfit=1e-4;           % accuracy level for iterative full IE solution
srcpar{1}=[1 1];          % source parameters
srcpar{2}=[1 2];          % source parameters
     sig0=[0.01 0.03];         % vector of background sigma
      hh0=[500];             % vector of background layer thicknesses
      an0=[1];            % vector of background layer anisotropies
        x=-1550:100:1550;   % Cell center x-coordinates of the anomalous area (m)
        y=-950:100:950;   % Cell center y-coordinates of the anomalous area (m)
        z=[250:100:1050 ]; % Cell center z-levels
       dz=[100 100 100 100 100 100 100 100 100];  % vertical cell sizes 
  solflag= 4;             % Solver flag.


Nxyz = length(x)*length(y)*length(z);
sb = bgcond(hh0,sig0,x,y,z);
ts = sb;			   % total conductivity
ds = zeros(Nxyz,1);		   % anomalous conductivity


[ds,ts] = rectbody(x,y,z,sb,ts,1  ,-1600, 1600,-1000,-900,200,1100);
[ds,ts] = rectbody(x,y,z,sb,ts,1  ,-1600, 1600,900,1000,200,1100);
[ds,ts] = rectbody(x,y,z,sb,ts,1  ,-1600, -1500,-1000,1000,200,1100);
[ds,ts] = rectbody(x,y,z,sb,ts,1  ,1500, 1600,-1000,1000,200,1100);
[ds,ts] = rectbody(x,y,z,sb,ts,1  ,-1600, 1600,-1000,1000,1000,1100);


ds = [ds zeros(Nxyz,1)];
save sigbody.dat ts -ascii