Shell Programming and Scripting

      character*70 infile,outfile
      integer run_mode
      common/mode/run_mode,infile

      write(*,*) 'input event file: (e.g. events.lhe)'
      read(*,'(a)')  infile
c
      write(*,*) 'name for output file: (e.g. dec-events.lhe)'
      read(*,'(a)')  outfile

#!/usr/bin/env bash

# @(#) do-all-files	Processing with static data, but different pairs of files.

pe() { for _i;do printf "%s" "$_i";done; printf "\n"; }
pl() { pe;pe "-----" ;pe "$*"; }
db() { ( printf " db, ";for _i;do printf "%s" "$_i";done;printf "\n" ) >&2 ; }
db() { : ; }
C=$HOME/bin/context && [ -f $C ] && $C gfortran

FILE=${1-"file-pairs"}

pl " File of input, output filename pairs, $FILE:"
head $FILE

i=1
pl " Processing pairs of files:"
while read input output
do
echo 
echo "-----"
echo " Iteration $i"
echo
./a.out <<EOF
3
$input
$output
rho
5
EOF
i=$( expr $i + 1 )
done < <( cat $FILE )

exit 0

% ./do-all-files 

Environment: LC_ALL = C, LANG = C
(Versions displayed with local utility "version")
OS, ker|rel, machine: Linux, 2.6.26-2-amd64, x86_64
Distribution        : Debian GNU/Linux 5.0.8 (lenny) 
bash GNU bash 3.2.39
gfortran GNU Fortran (Debian 4.3.2-1.1) 4.3.2

-----
 File of input, output filename pairs, file-pairs:
i1 o1
input-2 output-2
in_3.txt out_3.txt

-----
 Processing pairs of files:

-----
 Iteration 1


   Input run mode:
  Run mode is            3

 input event file: (e.g. events.lhe)
  input is i1                            

 name for output file: (e.g. dec-events.lhe)
  output is o1                            

  Input particle to be decayed (e.g. t~):
  particle is rho                           

 1   ta- -> vt   e- ve~           
 2   ta- -> vt  mu- vm~           
 3   ta- -> vt   l- vl~ (e+mu)    
 4   ta- -> vt   pi-              
 5   ta- -> vt  rho(770)-         
 your choice is:
  mode is:           5

-----
 Iteration 2


   Input run mode:
  Run mode is            3

 input event file: (e.g. events.lhe)
  input is input-2                       

 name for output file: (e.g. dec-events.lhe)
  output is output-2                      

  Input particle to be decayed (e.g. t~):
  particle is rho                           

 1   ta- -> vt   e- ve~           
 2   ta- -> vt  mu- vm~           
 3   ta- -> vt   l- vl~ (e+mu)    
 4   ta- -> vt   pi-              
 5   ta- -> vt  rho(770)-         
 your choice is:
  mode is:           5

-----
 Iteration 3


   Input run mode:
  Run mode is            3

 input event file: (e.g. events.lhe)
  input is in_3.txt                      

 name for output file: (e.g. dec-events.lhe)
  output is out_3.txt                     

  Input particle to be decayed (e.g. t~):
  particle is rho                           

 1   ta- -> vt   e- ve~           
 2   ta- -> vt  mu- vm~           
 3   ta- -> vt   l- vl~ (e+mu)    
 4   ta- -> vt   pi-              
 5   ta- -> vt  rho(770)-         
 your choice is:
  mode is:           5