write_parmela.m

使用matlab编辑文本并输出

clear all;
n=5000; % particle number.
x = normrnd(0,0.0004,[1 n]); % particle position distribution sigam=0.002 cm
x_p=normrnd(0,0.0005,[1 n]); % particle position prime distribution,sigma=0.002 rad
y=normrnd(0,0.0004,[1 n]);
y_p=normrnd(0,0.0005,[1 n]);

%rn = random('Normal',0,1,1,n);
%plot(x,x_p,'.')
fid = fopen('C:\LANL\mywork\Pepper_Pot.acc','w');
fprintf(fid,'run 1 1 2856. -.000001 150\n');
fprintf(fid,'title \n');
fprintf(fid,'Pepper Pot  \n');
fprintf(fid,'\n');
fprintf(fid,';start particle parameters input\n');
for i=1:n
fprintf(fid,'input 4, %f %f %f %f %f %f\n',x(i),x_p(i),y(i),y_p(i),0,0);
end
fprintf(fid,';end particle parameters input\n');
fprintf(fid,'\n');

fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');

fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');

fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');

fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');

fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');

fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'\n');

%fprintf(fid,'rfqout 0 5 1 60 10 5\n');
fprintf(fid,'\n');

fprintf(fid,'quad  12.0 5. 1 2000\n');

fprintf(fid,'\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');

fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');

fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');

fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');

fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');

fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');

fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');

fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');

fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');

fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');

fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');

fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');

fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');

fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');

fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');

fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');

fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');

fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');
fprintf(fid,'drift 1 5 1\n');

fprintf(fid,'drift 10 5 1\n');
fprintf(fid,'drift 10 5 1\n');
fprintf(fid,'drift 10 5 1\n');
fprintf(fid,'drift 10 5 1\n');
fprintf(fid,'drift 10 5 1\n');

fprintf(fid,'output 5\n');
fprintf(fid,'\n');
fprintf(fid,'start 40. 5 9200000 1 20\n');
fprintf(fid,'\n');
fprintf(fid,'end\n');
fclose(fid);
dos('RUN_INP');% before runing, add c:\LANL\mywork\ to path