Run TRIM in bath mode, but get Run-time error '62' \n Input past end of file

Question

Edoardo Renaldin 2021-11-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1595774-run-trim-in-bath-mode-but-get-run-time-error-62-n-input-past-end-of-file

编辑： dpb 2021-11-27

Hello, everyone! I am writing a simple code to run TRIM in batch mode to calculate the average proton energy exiting from a target material. I developed the following script where I change the inlet proton energy in the file TRIM.IN, but once I run TRIM.exe it give me error '62', do you know how to solve this problem?

function [] = SRIM_automatic(energy_next,energy_prev)
% Enter units: energy [MeV], step [MeV]
energy_next = energy_next * 10^3; %[keV] energy to insert
energy_prev = energy_prev * 10^3; %[keV] energy to substitute
%system('start /D C:\Users\renald_e\Desktop\SRIM TRIM.exe');
fid = fopen('H:\folder...\TRIM.IN','r');
i = 1;
while ~feof(fid)
    temp_line = fgetl(fid);
    trim_file{i} = temp_line;
    i = i + 1;
end
fclose(fid);
% Edit the initial proton energy: third element on third line
str_loc = strfind(trim_file{3},num2str(energy_prev));
str_length = length(num2str(energy_next));
% Substitute the proton energy
switch str_length
    case 4
        trim_file{3}(str_loc:str_loc+str_length-1) = sprintf('%c%s',' ',num2str(energy_next));
    case 5
        trim_file{3}(str_loc:str_loc+str_length-1) = sprintf('%s',num2str(energy_next));
end
fid = fopen('H:\folder...\TRIM.IN','w');
for i=1:numel(trim_file)
    if i == numel(trim_file)
        fprintf(fid,'%s',trim_file{i});
        break
    else
        fprintf(fid,'%s\n',trim_file{i});
    end
end
fclose(fid);
end

5 个评论
显示 3更早的评论隐藏 3更早的评论

dpb 2021-11-25

编辑：dpb 2021-11-25

在 MATLAB Online 中打开

Debug it at the command line first so you at least can tell which end it is that has the problem with the full error message.

Either way, it's far easier to just read the full file if you're going to bring it all into memory.

If all you're going to do is to modify the third record, as this does, there's really no reason to bring the file into memory, just work on each line -- if don't modify it, then just echo it back out...of course, this needs a second temporary output file to save the modifications into and then copy it back over the input when done.

function [] = SRIM_automatic(energy_next,energy_prev)
  % Enter units: energy [MeV], step [MeV]
  energy_next = energy_next * 10^3; %[keV] energy to insert
  energy_prev = energy_prev * 10^3; %[keV] energy to substitute
  fid = fopen('H:\folder...\TRIM.IN','r');
  fod = fopen('H:\folder...\TRIM.OUT','w');
  i=0;
  while ~feof(fid)
    i=i+1;                  % increment first, possibly where your issue is...
    l=fgets(fid);           % fgets() is mandatory here to return \n as well
    if i==3                 % this is the magic line to fix up, make the substitution
      l=strrep(l,num2str(energy_prev,energy_next));
    end
    fwrite(fod,l)
end
fclose('all');
movefile('H:\folder...\TRIM.OUT','H:\folder...\TRIM.IN')

This has the one possible advantage that you do have a temporary file that could be verified before destroying the input file if you gave yourself the opportunity to look at the .OUT file content before executing movefile

Edoardo Renaldin 2021-11-26

@dpb Thank you very much, it worked!

Sorry if I did not answered to the question about the "fixed field widths", it seems that the spacing does not affect the proper functioning of the program.

Do you have any suggestion about a website or book where I can improve my knowledge about MATLAB? Like this stuff of algorythm optimization, avoiding to store the whole file...?

Thank you again.

dpb 2021-11-26

编辑：dpb 2021-11-27

Glad to help...coaching/teaching is something I try to do as well as just answer a specific Q? You note I never did address your initial problem of why your code ran in to the EOF() issue but steered you in another direction entirely. :)

As for books, I'm sorry, no...I don't know of any that are either still in print or of a level for novice; I'm an old geezer and all my reference material other than the current documentation is at least 40 years old or, if more recent, very specialized stuff I bought while still in the active consulting gig -- but it's approaching 20 year now since I gave that up, too.

I know there is a link at the MathWorks site on MATLAB books, probably start there...

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

dpb 2021-11-25

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1595774-run-trim-in-bath-mode-but-get-run-time-error-62-n-input-past-end-of-file#answer_840294

编辑：dpb 2021-11-25

在 MATLAB Online 中打开

The above looks like could be much simpler if it is just doing a string substitution on a given line...try something like

function [] = SRIM_automatic(energy_next,energy_prev)
  % Enter units: energy [MeV], step [MeV]
  energy_next = energy_next * 10^3; %[keV] energy to insert
  energy_prev = energy_prev * 10^3; %[keV] energy to substitute
  %system('start /D C:\Users\renald_e\Desktop\SRIM TRIM.exe');
  trim_file=textread('H:\folder...\TRIM.IN','%s','whitespace','','delimiter',newline);
  % Edit the initial proton energy: third element on third line
  trim_file(3)=strrep(num2str(energy_prev),num2str(energy_next));
  % Write the new file
  writecell('H:\folder...\TRIM.IN',trim_file,'FileType','text','QuoteStrings',0);
end

NB: Air code, not tested -- as does yours, this code overwrites the original file, be SURE to have backups before testing.

NB2: Above replaces the prev string with next string; it does not worry about the spacing. If it is critical, then

newstr=num2str(energy_next);                % get the new conversion as char
lnew=length(newstr);                        % and it's length
newstr=strcat({blanks(4-lnew)}, newstr);    % pad new string to four if needed

3 个评论
显示 1更早的评论隐藏 1更早的评论

dpb 2021-11-25

Here's a convenient place for venerable textread.

You've not answered the Q? though, of whether the fixed width fields are significant or not -- if they are, your code isn't compensating for it as is...you test for the length of the new number, but not for the length of the text it is replacing.

I modified the above a little...can deal with the length Q? if know what is needed.

dpb 2021-11-25

编辑：dpb 2021-11-25

在 MATLAB Online 中打开

Illustration the above (modified) code logic works as intended -- I downloaded your file, changed the 13000 to 12980 and read the modified file back in to check...all in three lines of code...

txt=textread('trim.txt','%s','whitespace','','delimiter',newline);   % read as cellstr array
txt(3)=strrep(txt(3),'13000','12980');                               % make a substitution 
writecell(txt,'trim.txt','FileType','text','QuoteStrings',0)         % write the modified file

Let's see what we got...type the file back out at the command line...

>> type trim.txt
==> SRIM-2013.00 This file controls TRIM Calculations.
Ion: Z1 ,  M1,  Energy (keV), Angle,Number,Bragg Corr,AutoSave Number.
     1   1.008       12980       0     400      .94      400
Cascades(1=No;2=Full;3=Sputt;4-5=Ions;6-7=Neutrons), Random Number Seed, Reminders
                      1                                   0       0
Diskfiles (0=no,1=yes): Ranges, Backscatt, Transmit, Sputtered, Collisions(1=Ion;2=Ion+Recoils), Special EXYZ.txt file
                          0       0           1       0               0                               0
Target material : Number of Elements & Layers
"H (72000) into Target                   "       2               1
PlotType (0-5); Plot Depths: Xmin, Xmax(Ang.) [=0 0 for Viewing Full Target]
       5                         0         3000000
Target Elements:    Z   Mass(amu)
Atom 1 = Er =       68  167.26
Atom 2 = O =         8  15.999
Layer   Layer Name /               Width Density    Er(68)    O(8)
Numb.   Description                (Ang) (g/cm3)    Stoich  Stoich
 1      "Target"           3000000.11920929  4.72      .4      .6
0  Target layer phases (0=Solid, 1=Gas)
0 
Target Compound Corrections (Bragg)
 1  
Individual target atom displacement energies (eV)
      25      28
Individual target atom lattice binding energies (eV)
       3       3
Individual target atom surface binding energies (eV)
    3.05       2
Stopping Power Version (1=2011, 0=2011)
 0 
>> 