imprim.m

来自「德国人开发的地震处理分析软件」· M 代码 · 共 69 行

function imprim

global asciiaxes;
global archsalant;
global exten;
global archX;
global asciicol;
global psdX;
global corX;
global psdXY;
global corXY
global tfXY;
global psdgrafY;

mensaje1=  'NOMBRE DE ARCHIVO INCORRECTO';
mensaje2 = 'ARCHIVO IMPRESO             ';

axes(asciiaxes);hold off;
fill([0 0 1 1 0],[0 1 1 0 0],'k');
axis('off');

sal=[archsalant exten];

fid=fopen(sal,'wt');
if fid < 0;
   text('Position',[0.00 0.00],'Units','normalized','Color','r','String',mensaje1);
   hold off;
   return;
else;
   if exten=='.psd';
      if psdgrafY==[];
         fprintf(fid,'%12.4e %12.4e %12.4e\n',psdX);
         fclose(fid);
         sal=[archsalant,'.cor'];
         fid=fopen(sal,'wt');
         fprintf(fid,'%12.4e\n',corX);
         fclose(fid);
      else;
         fprintf(fid,'%12.4e %12.4e %12.4e %12.4e %12.4e %12.4e\n',psdXY);
         fclose(fid);
         sal=[archsalant,'.tf'];
         fid=fopen(sal,'wt');
         fprintf(fid,'%12.4e %12.4e %12.4e\n',tfXY);
         fclose(fid);
         sal=[archsalant,'.cor'];
         fid=fopen(sal,'wt');
         fprintf(fid,'%12.4e %12.4e %12.4e\n',corXY);
         fclose(fid);
         clear global psdX corX psdXY corXY tfXY; 
      end;
   else 
      if asciicol==1;
         fprintf(fid,'%12.4e\n',archX);
      else
         if asciicol==5;
            fprintf(fid,'%12.4e %12.4e %12.4e %12.4e %12.4e\n',archX);
         else
            if asciicol==3;
               fprintf(fid,'%12.4e %12.4e %12.4e\n',archX);
            end;
         end;
      end;
      clear global archX;
      fclose(fid);
   end;
   text('Position',[0.00 0.00],'Units','normalized','Color','r','String',mensaje2);
   hold off;
end;