The document describes the use of fr2in.for
It was written in 2011-12 by Alex Brown with the guidance of Ian Thompson.

It is used for (p,t), (t,p), (p,3He), (3He,p) and (3He,n) reactions
on targets with J=0 (even=even). Any final J is allowed.

If you want to make a *.eps figure you also need
plotsig.for and cps.for.

All of these programs can be compiled with gfortran

for.zip contains the source codes for fr2in, plotsig and cps
ca40.zip contains the input and output files for the example below

-------------------------------------------------------
Sample input file j0.dai for 40Ca 0+ to 38Ca 0+

bbcm00_co00.tna
re = pt
el = 35.

The first line must be the name of the *.tna file from oxbash, nushell or nushellx.
The following lines must all be in the format xx = yy
re is the type of reaction. The options are pt, tp, ph and hp
el is the laboratory projectile energy in MeV

With no other input the program uses the following (default) 
internal optical potentials (for details see below): 

proton:   CH 
deuteron: DG 
triton:   XL 
helium:   BG 

To run fr2in type

fr2in < j0.dai > j0.dao

fr2in makes the fresco input files j0t1.in (direct only) and 
j0t2.in (direct plus sequential) for fresco.

The windows command file fr2.bat given below runs fr2in and fresco and 
then makes a plot.

Other input options are documented in the *.dao file made with fr2in.
These include the following given by example:
nf = 2  for the second 0+ state (the default nf = 1)
ef = 2. excitation energy (MeV) of the final state (the default is 0.)
sp = 0.   to change effective single-particle separation energy, the default is (bef-bei)/2    
rm = 30.  to change rmatch in fresco (in fm) 

pp = CH, BG, ME or PE optical potentials for protons 
dp = DG, PP or LH     optical potentials for deuterons 
tp = XL, BG or PA     optical potentials for tritons 
hp = BG or PA         optical potentials for helium

proton potentials
 [BG] Bechetti-Greenlees (A>40 20<E<50 MeV):    Phys Rev 182, 1190 (1969)
 [CH] Chapel-Hill 89 Global set (A>40 E>10 MeV):Phys Rep 201, 57 (1991)
 [ME] Menet (30<E<60):                          ADNTD 17, 6 (1976)
 [PE] Perey (E<20MeV):                          ADNTD 17, 6 (1976)

deuteron potentials
 [LH] Lohr-Haeberli (A>40 8<E<13 MeV):          ADNTD 17, 6 (1976)
 [PP] Perey-Perey (12<E<25 MeV) no spin-orbit:  ADNTD 17, 6 (1976)
 [DG] Daehnick Global (A>27 12<E<90 MeV)        PRC 21, 2253 (1980)

triton potentials
 [BG] Bechetti-Greenlees (not well determined): ADNTD 17, 6 (1976)
 [XL] X. Li et al Global potential < 40 MeV:    NPA 789, 103 (2007)
 [PA] D.Y. Pang et al. GDP08                    PRC 79, 024615 (2009)

helium potentials
 [BG] Bechetti-Greenlees (not well determined): ADNTD 17, 6 (1976)
 [PA] D.Y. Pang et al. GDP08                    PRC 79, 024615 (2009)

All of these internal options are based on a subroutine provided
by Jeff Tostevin.
When fr2in is run it makes pp.dao, dp.dao, tp.dao and hp.dao files.
These contain the optical potential parameters from the internal options.

These can be modified as in the following example for protons

rename pp.dao to pp*.pot where * is the name of your modification
modify pp*.pot
put the following line in the *.dai file

pp = pp*

then run fr2in over again.   

If one has many *.dai files and would like to change the 
potentials for all of them. Make a pot.dai file. If the pot.dai
file exists it will override anything in the other *.dai files.
pot.dai should contain a list of potential options in the
formats given above.     


----------------------------------------
the bbcm00_co00.tna file if given below this line.
The numbers in this file can be modified, but the format must remain the same.
----------------------------------------
   
pn format, made from nushellx *.trz file
! model-space name = sd       
! interaction name = usdb     
!
!     ----------------------------------
!     |  ai <-> af  =    38  <->   40  |
!     |  zi <-> zf  =    20  <->   20  |
!     | 2ji <-> 2jf =     0  <->    0  |
!     |  pi <->  pf =     +  <->    +  |
!     ----------------------------------
!
! one-body transition densities a(dj) = 
! = <f|| [a+(k1)a(k2)]^(dj) ||i>/sqrt(2dj+1)
! for input to some dwba codes: z = a(dj)/sqrt(2ji+1)
! with edmonds (de-shalit talmi) reduced matrix element convention
!
! for n,l,2j =    1  2  3  label  k =   5
! for n,l,2j =    1  2  5  label  k =   4
! for n,l,2j =    2  0  1  label  k =   6
!
! ji,   jf,   ti,   tf,  ti,  fz,
! dj,   ni,   nf,       ef,       ei,      exi,      exf,
!
  0.0,  0.0,  1.0,  0.0, -1.0,  0.0,
  0.0,   1.,   1.,  -279.961,  -250.864,     0.000,     0.000,
   4,   4,   0.00000,   0.16403,  ! k1,k2,a(pp),a(nn)
   5,   5,   0.00000,   0.97237,
   6,   6,   0.00000,   0.16610,
   0,
  0.0,   2.,   1.,  -279.961,  -244.836,     6.028,     0.000,
   4,   4,   0.00000,   0.07794,
   5,   5,   0.00000,  -0.18065,
   6,   6,   0.00000,   0.98051,
   0,
   
   
--------------------------------------
fr2.bat (given below) is a windows command file
that runs fr2in and fresco and then makes a figure
--------------------------------------

rem - run fr2in to make the fresco inputs 
rem - *t1.in for direct
rem - *t2.in for direct plus sequential

fr2in < %1.dai > %1.dao

rem ----- run fresco for *t1.in
copy babt1.in %1t1.in
del babt1.in
fresco < %1t1.in > %1t1.out

rem ----- copy the fresco output into the *.elas, *.sig files
if exist fort.200 goto aa 
copy fort.201 %1t1.elas
copy fort.202 %1t1.sig
del fort.201
del fort.202
goto bb
:aa
copy fort.200 %1t1.elas
copy fort.201+fort.202 %1t1.sig
del fort.200
del fort.201
del fort.202
:bb
copy fort.13 %1t1.tot

rem ----- run fresco for *t2.in

copy babt2.in %1t2.in
del babt2.in
fresco < %1t2.in > %1t2.out
rem ----- copy the fresco output into the *.elas, *.sig files
if exist fort.200 goto aa 
copy fort.201 %1t2.elas
copy fort.202 %1t2.sig
del fort.201
del fort.202
goto bb
:aa
copy fort.200 %1t2.elas
copy fort.201+fort.202 %1t2.sig
del fort.200
del fort.201
del fort.202
:bb
copy fort.13 %1t2.tot

rem ----- convert the *.sig to *.plt 
plotsig %1

rem ----- make *.eps with log y axis - you must have cps.exe to use this
cps %1 > temp.out

rem ----- make *.eps with decade y axis - you must have cps.exe to use this
cps %1x > temp.out
del temp.out

rem ----- delete the uneeded file
del fort.13
del fort.16
del fort.3
del fort.35
del fort.38
del fort.39
del fort.4
del fort.40
del fort.45
del fort.46
del fort.48
del fort.56
del fort.7
del fort.75