vec/del *
fn=beam_energy_63158_63194.out
npa=1.
npc=5
fn=beam_energy_62914_63092.out
npa=1.
npc=4
fn=beam_energy_62873_62912.out
npa=1.
npc=5
fn=beam_energy_63287_63322.out
npa=4.
npc=1.
**
fn=beam_energy_63160_63280.out 
* use for 63160 to 63218
npa=1.
npc=5.
* use for 63219 to 63280
*npa=5.
*npc=5.
**
fn=beam_energy_63287_63322.out
npa=4.
npc=1.
**
fn=beam_energy_63324_63416.out
npa=4.
npc=5.
*
fn=beam_energy_63094_63156.out
npa=1.
npc=2
**
fn=beam_energy_63417_63633.out
* before run 63465
npa=2.
npc=3.
* after run 63465
npa=4.
npc=3.
*
vec/read nrun,ehalla,erra,ehallc,errc [fn]
sigma elinaca=ehalla/(2*[npa]+0.1125)
sigma elinacc=ehallc/(2*[npc]+0.1125)
sigma perd=(elinacc-elinaca)/elinacc*10000
vec/write nrun,elinaca,elinacc,perd ' ' 4f15.5