run2.m

This software gave volume MoM solution by the CG-FFT method

%run2

frq=str2num(get(FRQ,'string'))*1000000;

b=[5e-2 100];
b1=[nx-1 ny-1 nz-1];
b2=[str2num(get(ANGLET,'string')) str2num(get(ANGLEP,'string'))];
b3=[str2num(get(THPH,'string'))+1 str2num(get(STEPS,'string'))];
b4=str2num(get(POL1,'string'));
b5=str2num(get(NOSTEPS,'string'));
b6=[dx*nx*U dy*ny*U dz*nz*U];
b71='S'; 
b72='PEPE';
b73='0';
b74='0';
b75='0';

b8='B';
if (b4==1), b8='B'; end;
if (b4==2), b8='C'; end;
b91='1';
b92='1.';
cd ..
cd fFiles

!del input.DAT
fid=fopen('input.DAT','w');
fwrite(fid,[num2str(frq),10]);
fwrite(fid,[num2str(b1),10]);
fwrite(fid,[num2str(b),10]);
fwrite(fid,[num2str(b2),10]);
fwrite(fid,[num2str(b3),10]);
fwrite(fid,[num2str(b4),10]);
fwrite(fid,[num2str(b5),10]);
fwrite(fid,[num2str(b6),10]);
fwrite(fid,[b71,10]);
fwrite(fid,[b72,10]);
fwrite(fid,[b73,10]);
fwrite(fid,[b74,10]);
fwrite(fid,[b75,10]);
fwrite(fid,[b8,10]);
fwrite(fid,[b91,10]);
fwrite(fid,[b92,10]);
fclose(fid);

%a=ones(nx*ny*nz,1);

if (exist('a1')==0), 
   cd ..
   cd mfiles
   errordlg('You must first create a grid before calculating RCS')
   break;
end;
save input1.DAT a1 -ascii
save input2.DAT a2 -ascii
save input3.DAT a3 -ascii

s=size(gridd);
if (s(1)==2), gridd=gridd'; end;

save inputm.DAT gridd -ascii
s=['        parameter (npx=',num2str(nx),',npy=',num2str(ny),',npz=',num2str(nz),')'];
fid=fopen('PAR.INC','w');
fwrite(fid,[s,10]);
fclose(fid);

errortrap on
c1=clock;
! fl32 mom3d2.f /G5 /Ox /Op
errortrap off
c2=clock
if (etime(c2,c1)<10), 
   errordlg('You do not have the FORTRAN compiler fl32.exe in your path')
   cd ..
   cd mFiles
   break;
end;
! mom3d2.exe
close (2);

load input.RES -ascii
cd ..
cd mFiles

figure(2)
plot(input(:,b3(1)),input(:,7),'r');

ylabel('RCS (dBsm)');
if (b3(1)==1),
  xlabel('Theta, (DEG)');
  title(['RCS (dBsm) at Phi=',num2str(b2(2)),' in different Theta']);
end;
if (b3(1)==2),
  xlabel('Phi, (DEG)');
  title(['RCS (dBsm) at Theta=',num2str(b2(1)),' in different Phi']);
end;
grid;