borra.m

德国人开发的地震处理分析软件

function borra

global A;
global Apsd;
global Atf;
global Acor;
global Acorxy;
global Afft;
global Afil;
global Ares;
global Aine;
global Aria;
global Aamp;
global archivo;
global deltat;
global psddeltaf;
global tfdeltaf;
global fftdeltaf;
global fildeltaf;
global resdeltat;
global inedeltat;
global ampdeltaf;
global numero;
global numeropsd;
global numerotf;
global numerocor;
global numerocorxy;
global numerofft;
global numerofil;
global numerores;
global numeroine;
global numeroaria;
global numeroamp;
global nlbgraf;
global MaxA;
global ngraf;
global npsdgraf;
global ntfgraf;
global ncorgraf;
global ncorxygraf;
global nfftgraf;
global nfilgraf;
global nresgraf;
global ninegraf;
global nariagraf;
global nampgraf;
global tmax;
global linutil;
global letras;
global grafArange;

ngraf=[];

for i=1:length(deltat);
    ngraf=[ngraf,i];
end;

dif      =  numero(nlbgraf+1)-numero(nlbgraf);
final    =  length(deltat);
A        =  [A(1:numero(nlbgraf)-1);A(numero(nlbgraf+1):length(A))]; 
numero   =  [numero(1:nlbgraf);numero(nlbgraf+2:final+1)-dif];
linutil  =  [linutil(1:nlbgraf-1);linutil(nlbgraf+1:final)];
letras   =  [letras(1:nlbgraf-1);letras(nlbgraf+1:final)];
deltat   =  [deltat(1:nlbgraf-1);deltat(nlbgraf+1:final)];
archivo  =  [archivo(1:nlbgraf-1,:);archivo(nlbgraf+1:final,:)];
MaxA     =  [MaxA(1:nlbgraf-1);MaxA(nlbgraf+1:final)];

Dngraf = ngraf;
j = 0;
for I=Dngraf;
    j=j+1;
    if I==nlbgraf;
       ngraf=[ngraf(1:j-1), ngraf(j+1:length(Dngraf))];
       break;
    end;
end;

j = 0;
for I=ngraf;
    j=j+1;
    if I>nlbgraf;
       ngraf(j)=ngraf(j)-1;
    end;
end;

j=0;
for i=nfftgraf;
    j=j+1;
    if i==nlbgraf;
       dif      =  numerofft(j+1)-numerofft(j);
       final    =  length(fftdeltaf);
       Afft     =  [Afft(1:numerofft(j)-1);Afft(numerofft(j+1):length(Afft))]; 
       numerofft=  [numerofft(1:j);numerofft(j+2:final+1)-dif];
       fftdeltaf=  [fftdeltaf(1:j-1);fftdeltaf(j+1:final)];
       break;
    end;
end;

Dngraf = nfftgraf;
j = 0;
for I=Dngraf;
    j=j+1;
    if I==nlbgraf;
       nfftgraf=[nfftgraf(1:j-1), nfftgraf(j+1:length(Dngraf))];
       break;
    end;
end;

j = 0;
for I=nfftgraf;
    j=j+1;
    if I>nlbgraf;
       nfftgraf(j)=nfftgraf(j)-1;
    end;
end;

j=0;
for i=nfilgraf;
    j=j+1;
    if i==nlbgraf;
       dif      =  numerofil(j+1)-numerofil(j);
       final    =  length(fildeltaf);
       Afil     =  [Afil(1:numerofil(j)-1);Afil(numerofil(j+1):length(Afil))]; 
       numerofil=  [numerofil(1:j);numerofil(j+2:final+1)-dif];
       fildeltaf=  [fildeltaf(1:j-1);fildeltaf(j+1:final)];
       break;
    end;
end;

Dngraf = nfilgraf;
j = 0;
for I=Dngraf;
    j=j+1;
    if I==nlbgraf;
       nfilgraf=[nfilgraf(1:j-1), nfilgraf(j+1:length(Dngraf))];
       break;
    end;
end;

j = 0;
for I=nfilgraf;
    j=j+1;
    if I>nlbgraf;
       nfilgraf(j)=nfilgraf(j)-1;
    end;
end;

j=0;
for i=nampgraf;
    j=j+1;
    if i==nlbgraf;
       dif      =  numeroamp(j+1)-numeroamp(j);
       final    =  length(ampdeltaf);
       Aamp     =  [Aamp(1:numeroamp(j)-1);Aamp(numeroamp(j+1):length(Aamp))]; 
       numeroamp=  [numeroamp(1:j);numeroamp(j+2:final+1)-dif];
       ampdeltaf=  [ampdeltaf(1:j-1);ampdeltaf(j+1:final)];
       break;
    end;
end;

Dngraf = nampgraf;
j = 0;
for I=Dngraf;
    j=j+1;
    if I==nlbgraf;
       nampgraf=[nampgraf(1:j-1), nampgraf(j+1:length(Dngraf))];
       break;
    end;
end;

j = 0;
for I=nampgraf;
    j=j+1;
    if I>nlbgraf;
       nampgraf(j)=nampgraf(j)-1;
    end;
end;

j=0;
for i=nresgraf;
    j=j+1;
    if i==nlbgraf;
       dif      =  numerores(j+1)-numerores(j);
       final    =  length(resdeltat);
       Ares     =  [Ares(1:numerores(j)-1,:);Ares(numerores(j+1):length(Ares),:)]; 
       numerores=  [numerores(1:j);numerores(j+2:final+1)-dif];
       resdeltat=  [resdeltat(1:j-1);resdeltat(j+1:final)];
       break;
    end;
end;

Dngraf = nresgraf;
j = 0;
for I=Dngraf;
    j=j+1;
    if I==nlbgraf;
       nresgraf=[nresgraf(1:j-1), nresgraf(j+1:length(Dngraf))];
       break;
    end;
end;

j = 0;
for I=nresgraf;
    j=j+1;
    if I>nlbgraf;
       nresgraf(j)=nresgraf(j)-1;
    end;
end;

j=0;
for i=ninegraf;
    j=j+1;
    if i==nlbgraf;
       dif      =  numeroine(j+1)-numeroine(j);
       final    =  length(inedeltat);
       Aine     =  [Aine(1:numeroine(j)-1,:);Aine(numeroine(j+1):length(Aine),:)]; 
       numeroine=  [numeroine(1:j);numeroine(j+2:final+1)-dif];
       inedeltat=  [inedeltat(1:j-1);inedeltat(j+1:final)];
       break;
    end;
end;

Dngraf = ninegraf;
j = 0;
for I=Dngraf;
    j=j+1;
    if I==nlbgraf;
       ninegraf=[ninegraf(1:j-1), ninegraf(j+1:length(Dngraf))];
       break;
    end;
end;

j = 0;
for I=ninegraf;
    j=j+1;
    if I>nlbgraf;
       ninegraf(j)=ninegraf(j)-1;
    end;
end;

j=0;
for i=nariagraf;
    j=j+1;
    if i==nlbgraf;
       dif      =  numeroaria(j+1)-numeroaria(j);
       final    =  length(numeroaria)-1;
       Aria     =  [Aria(1:numeroaria(j)-1,:);Aria(numeroaria(j+1):length(Aria),:)]; 
       numeroaria=  [numeroaria(1:j);numeroaria(j+2:final+1)-dif];
       break;
    end;
end;

Dngraf = nariagraf;
j = 0;
for I=Dngraf;
    j=j+1;
    if I==nlbgraf;
       nariagraf=[nariagraf(1:j-1), nariagraf(j+1:length(Dngraf))];
       break;
    end;
end;

j = 0;
for I=nariagraf;
    j=j+1;
    if I>nlbgraf;
       nariagraf(j)=nariagraf(j)-1;
    end;
end;

j=0;
for i=npsdgraf;
    j=j+1;
    if i==nlbgraf;
       dif      =  numeropsd(j+1)-numeropsd(j);
       final    =  length(psddeltaf);
       Apsd     =  [Apsd(1:numeropsd(j)-1,:);Apsd(numeropsd(j+1):length(Apsd),:)]; 
       numeropsd=  [numeropsd(1:j);numeropsd(j+2:final+1)-dif];
       psddeltaf=  [psddeltaf(1:j-1);psddeltaf(j+1:final)];
       break;
    end;
end;

Dngraf = npsdgraf;
j = 0;
for I=Dngraf;
    j=j+1;
    if I==nlbgraf;
       npsdgraf=[npsdgraf(1:j-1), npsdgraf(j+1:length(Dngraf))];
       break;
    end;
end;

j = 0;
for I=npsdgraf;
    j=j+1;
    if I>nlbgraf;
       npsdgraf(j)=npsdgraf(j)-1;
    end;
end;

j=0;
for i=ncorgraf;
    j=j+1;
    if i==nlbgraf;
       dif      =  numerocor(j+1)-numerocor(j);
       final    =  length(ncorgraf);
       Acor     =  [Acor(1:numerocor(j)-1);Acor(numerocor(j+1):length(Acor))]; 
       numerocor=  [numerocor(1:j);numerocor(j+2:final+1)-dif];
       break;
    end;
end;

Dngraf = ncorgraf;
j = 0;
for I=Dngraf;
    j=j+1;
    if I==nlbgraf;
       ncorgraf=[ncorgraf(1:j-1), ncorgraf(j+1:length(Dngraf))];
       break;
    end;
end;

j = 0;
for I=ncorgraf;
    j=j+1;
    if I>nlbgraf;
       ncorgraf(j)=ncorgraf(j)-1;
    end;
end;

if ncorxygraf~=[];

   j=0;
   for i=ncorxygraf(1,:);
       j=j+1;
       k=ncorxygraf(2,j);
       if i==nlbgraf | k==nlbgraf;
          dif        =  numerocorxy(j+1)-numerocorxy(j);
          final      =  length(ncorxygraf(1,:));
          Acorxy     =  [Acorxy(1:numerocorxy(j)-1);Acorxy(numerocorxy(j+1):length(Acorxy))]; 
          numerocorxy=  [numerocorxy(1:j);numerocorxy(j+2:final+1)-dif];
          j=j-1;
       end;
   end;

   Dngraf = ncorxygraf;
   j = 0;
   for I=Dngraf(1,:);
       j=j+1;
       k=ncorxygraf(2,j);
       if I==nlbgraf | k==nlbgraf;
          ncorxygraf=[ncorxygraf(:,1:j-1), ncorxygraf(:,j+1:length(Dngraf(1,:)))];
          j=j-1;
       end;
   end;

   if ncorxygraf~=[];
    
   j = 0;
   for I=ncorxygraf(1,:);
       j=j+1;
       k=ncorxygraf(2,j);
       if I>nlbgraf;
          ncorxygraf(1,j)=ncorxygraf(1,j)-1;
       end;
       if k>nlbgraf;
          ncorxygraf(2,j)=ncorxygraf(2,j)-1;
       end;     
   end;
   end;
end;

if ntfgraf~=[];

   j=0;
   for i=ntfgraf(1,:);
       j=j+1;
       k=ntfgraf(2,j);
       if i==nlbgraf | k==nlbgraf;
          dif     =  numerotf(j+1)-numerotf(j);
          final   =  length(ntfgraf(1,:));
          Atf     =  [Atf(1:numerotf(j)-1,:);Atf(numerotf(j+1):length(Atf),:)]; 
          numerotf=  [numerotf(1:j);numerotf(j+2:final+1)-dif];
          j=j-1;
       end;
   end;

   Dngraf = ntfgraf;
   j = 0;
   for I=Dngraf(1,:);
       j=j+1;
       k=ntfgraf(2,j);
       if I==nlbgraf | k==nlbgraf;
          ntfgraf=[ntfgraf(:,1:j-1), ntfgraf(:,j+1:length(Dngraf(1,:)))];
          j=j-1;
       end; 
   end;

   if ntfgraf~=[];

   j = 0;
   for I=ntfgraf(1,:);
       j=j+1;
       k=ntfgraf(2,j);
       if I>nlbgraf;
          ntfgraf(1,j)=ntfgraf(1,j)-1;
       end;
       if k>nlbgraf;
          ntfgraf(2,j)=ntfgraf(2,j)-1;
       end;     
   end;
   end;
end;

if length(deltat)~=0
   tmax=max(deltat(ngraf).*(numero(ngraf+1)-numero(ngraf)));
   Amax=max(MaxA(ngraf));
   grafArange=[-Amax Amax];
   grafica;
else;
   close;
end;