Benesh : Procedure options ( Main ) ;

/* Procedure for calculating */
/* the total nuclear cross sections */
/* according to the formula of Benesh, Cook and Vary */

Declare
(I_A1,I_A2) Bin Fixed(15),
(R_A1,R_A2) Dec Float(15),
bmin Dec Float(15),
pi Dec Float(15) init(3.141592654E0),
R_SigTot Dec Float(15),
C_Dummy Char(8) Var;


Put Edit('Calculation of the total nuclear interaction') (Skip,A);
Put Edit('cross section') (Skip,A);
Put Edit('according to Benesh, Cook and Vary') (Skip,A);
Put Edit('(Phys. Rev. C 40 (1989) 1198)') (Skip,A);

Put Edit('Enter A of projectile') (Skip,Skip,A);
Get List(I_A1);
Put Edit('Enter A of target') (A);
Get List(I_A2);

R_A1 = Float(I_A1,15);
R_A2 = Float(I_A2,15);

bmin = 1.34E0 * (R_A1** (1.E0 / 3.E0)
+ R_A2** (1.E0 / 3.E0)
- .75E0 * (R_A1** (-1.E0 / 3.E0)
+ R_A2** (-1.E0 / 3.E0)));
R_Sigtot = pi * bmin**2 * 10.E0; /* in units of mb */


Put Edit('The cross section (in mb) is ') (Skip,A);
Put Edit(R_SigTot) (F(10));

Put Edit('Check validity of energy range in Benesh et al.!')
(SKIP,A);

Get Edit(C_Dummy) (A(8));

End Benesh ;